#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r21144 = atan2(1.0, 0.0); float r21145 = 2; float r21146 = r21144 / r21145; float r21147 = 1; float r21148 = b; float r21149 = r21148 * r21148; float r21150 = a; float r21151 = r21150 * r21150; float r21152 = r21149 - r21151; float r21153 = r21147 / r21152; float r21154 = r21146 * r21153; float r21155 = r21147 / r21150; float r21156 = r21147 / r21148; float r21157 = r21155 - r21156; float r21158 = r21154 * r21157; return r21158; } double f_id(double a, double b) { double r21159 = atan2(1.0, 0.0); double r21160 = 2; double r21161 = r21159 / r21160; double r21162 = 1; double r21163 = b; double r21164 = r21163 * r21163; double r21165 = a; double r21166 = r21165 * r21165; double r21167 = r21164 - r21166; double r21168 = r21162 / r21167; double r21169 = r21161 * r21168; double r21170 = r21162 / r21165; double r21171 = r21162 / r21163; double r21172 = r21170 - r21171; double r21173 = r21169 * r21172; return r21173; } double f_of(float a, float b) { float r21174 = atan2(1.0, 0.0); float r21175 = b; float r21176 = a; float r21177 = r21175 + r21176; float r21178 = r21174 / r21177; float r21179 = 1; float r21180 = r21179 / r21176; float r21181 = r21179 / r21175; float r21182 = r21180 - r21181; float r21183 = 2; float r21184 = r21175 - r21176; float r21185 = r21183 * r21184; float r21186 = r21182 / r21185; float r21187 = r21178 * r21186; float r21188 = pow(r21187, r21179); return r21188; } double f_od(double a, double b) { double r21189 = atan2(1.0, 0.0); double r21190 = b; double r21191 = a; double r21192 = r21190 + r21191; double r21193 = r21189 / r21192; double r21194 = 1; double r21195 = r21194 / r21191; double r21196 = r21194 / r21190; double r21197 = r21195 - r21196; double r21198 = 2; double r21199 = r21190 - r21191; double r21200 = r21198 * r21199; double r21201 = r21197 / r21200; double r21202 = r21193 * r21201; double r21203 = pow(r21202, r21194); return r21203; } void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) { mpfr_fmod(r, n, d, rmd); if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd); } static mpfr_t r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216, r21217, r21218; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21204); mpfr_init_set_str(r21205, "2", 10, MPFR_RNDN); mpfr_init(r21206); mpfr_init_set_str(r21207, "1", 10, MPFR_RNDN); mpfr_init(r21208); mpfr_init(r21209); mpfr_init(r21210); mpfr_init(r21211); mpfr_init(r21212); mpfr_init(r21213); mpfr_init(r21214); mpfr_init(r21215); mpfr_init(r21216); mpfr_init(r21217); mpfr_init(r21218); } double f_im(double a, double b) { mpfr_const_pi(r21204, MPFR_RNDN); ; mpfr_div(r21206, r21204, r21205, MPFR_RNDN); ; mpfr_set_d(r21208, b, MPFR_RNDN); mpfr_mul(r21209, r21208, r21208, MPFR_RNDN); mpfr_set_d(r21210, a, MPFR_RNDN); mpfr_mul(r21211, r21210, r21210, MPFR_RNDN); mpfr_sub(r21212, r21209, r21211, MPFR_RNDN); mpfr_div(r21213, r21207, r21212, MPFR_RNDN); mpfr_mul(r21214, r21206, r21213, MPFR_RNDN); mpfr_div(r21215, r21207, r21210, MPFR_RNDN); mpfr_div(r21216, r21207, r21208, MPFR_RNDN); mpfr_sub(r21217, r21215, r21216, MPFR_RNDN); mpfr_mul(r21218, r21214, r21217, MPFR_RNDN); return mpfr_get_d(r21218, MPFR_RNDN); } static mpfr_t r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227, r21228, r21229, r21230, r21231, r21232, r21233; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r21219); mpfr_init(r21220); mpfr_init(r21221); mpfr_init(r21222); mpfr_init(r21223); mpfr_init_set_str(r21224, "1", 10, MPFR_RNDN); mpfr_init(r21225); mpfr_init(r21226); mpfr_init(r21227); mpfr_init_set_str(r21228, "2", 10, MPFR_RNDN); mpfr_init(r21229); mpfr_init(r21230); mpfr_init(r21231); mpfr_init(r21232); mpfr_init(r21233); } double f_fm(double a, double b) { mpfr_const_pi(r21219, MPFR_RNDN); mpfr_set_d(r21220, b, MPFR_RNDN); mpfr_set_d(r21221, a, MPFR_RNDN); mpfr_add(r21222, r21220, r21221, MPFR_RNDN); mpfr_div(r21223, r21219, r21222, MPFR_RNDN); ; mpfr_div(r21225, r21224, r21221, MPFR_RNDN); mpfr_div(r21226, r21224, r21220, MPFR_RNDN); mpfr_sub(r21227, r21225, r21226, MPFR_RNDN); ; mpfr_sub(r21229, r21220, r21221, MPFR_RNDN); mpfr_mul(r21230, r21228, r21229, MPFR_RNDN); mpfr_div(r21231, r21227, r21230, MPFR_RNDN); mpfr_mul(r21232, r21223, r21231, MPFR_RNDN); mpfr_pow(r21233, r21232, r21224, MPFR_RNDN); return mpfr_get_d(r21233, MPFR_RNDN); } static mpfr_t r21234, r21235, r21236, r21237, r21238, r21239, r21240, r21241, r21242, r21243, r21244, r21245, r21246, r21247, r21248; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r21234); mpfr_init(r21235); mpfr_init(r21236); mpfr_init(r21237); mpfr_init(r21238); mpfr_init_set_str(r21239, "1", 10, MPFR_RNDN); mpfr_init(r21240); mpfr_init(r21241); mpfr_init(r21242); mpfr_init_set_str(r21243, "2", 10, MPFR_RNDN); mpfr_init(r21244); mpfr_init(r21245); mpfr_init(r21246); mpfr_init(r21247); mpfr_init(r21248); } double f_dm(double a, double b) { mpfr_const_pi(r21234, MPFR_RNDN); mpfr_set_d(r21235, b, MPFR_RNDN); mpfr_set_d(r21236, a, MPFR_RNDN); mpfr_add(r21237, r21235, r21236, MPFR_RNDN); mpfr_div(r21238, r21234, r21237, MPFR_RNDN); ; mpfr_div(r21240, r21239, r21236, MPFR_RNDN); mpfr_div(r21241, r21239, r21235, MPFR_RNDN); mpfr_sub(r21242, r21240, r21241, MPFR_RNDN); ; mpfr_sub(r21244, r21235, r21236, MPFR_RNDN); mpfr_mul(r21245, r21243, r21244, MPFR_RNDN); mpfr_div(r21246, r21242, r21245, MPFR_RNDN); mpfr_mul(r21247, r21238, r21246, MPFR_RNDN); mpfr_pow(r21248, r21247, r21239, MPFR_RNDN); return mpfr_get_d(r21248, MPFR_RNDN); }