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