#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r21097 = atan2(1.0, 0.0); float r21098 = 2; float r21099 = r21097 / r21098; float r21100 = 1; float r21101 = b; float r21102 = r21101 * r21101; float r21103 = a; float r21104 = r21103 * r21103; float r21105 = r21102 - r21104; float r21106 = r21100 / r21105; float r21107 = r21099 * r21106; float r21108 = r21100 / r21103; float r21109 = r21100 / r21101; float r21110 = r21108 - r21109; float r21111 = r21107 * r21110; return r21111; } double f_id(double a, double b) { double r21112 = atan2(1.0, 0.0); double r21113 = 2; double r21114 = r21112 / r21113; double r21115 = 1; double r21116 = b; double r21117 = r21116 * r21116; double r21118 = a; double r21119 = r21118 * r21118; double r21120 = r21117 - r21119; double r21121 = r21115 / r21120; double r21122 = r21114 * r21121; double r21123 = r21115 / r21118; double r21124 = r21115 / r21116; double r21125 = r21123 - r21124; double r21126 = r21122 * r21125; return r21126; } double f_of(float a, float b) { float r21127 = 1; float r21128 = a; float r21129 = r21127 / r21128; float r21130 = b; float r21131 = r21127 / r21130; float r21132 = r21129 - r21131; float r21133 = r21130 + r21128; float r21134 = r21132 / r21133; float r21135 = r21130 - r21128; float r21136 = r21134 / r21135; float r21137 = atan2(1.0, 0.0); float r21138 = 2; float r21139 = r21137 / r21138; float r21140 = r21136 * r21139; return r21140; } double f_od(double a, double b) { double r21141 = 1; double r21142 = a; double r21143 = r21141 / r21142; double r21144 = b; double r21145 = r21141 / r21144; double r21146 = r21143 - r21145; double r21147 = r21144 + r21142; double r21148 = r21146 / r21147; double r21149 = r21144 - r21142; double r21150 = r21148 / r21149; double r21151 = atan2(1.0, 0.0); double r21152 = 2; double r21153 = r21151 / r21152; double r21154 = r21150 * r21153; return r21154; } 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 r21155, r21156, r21157, r21158, r21159, r21160, r21161, r21162, r21163, r21164, r21165, r21166, r21167, r21168, r21169; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21155); mpfr_init_set_str(r21156, "2", 10, MPFR_RNDN); mpfr_init(r21157); mpfr_init_set_str(r21158, "1", 10, MPFR_RNDN); mpfr_init(r21159); mpfr_init(r21160); mpfr_init(r21161); mpfr_init(r21162); mpfr_init(r21163); mpfr_init(r21164); mpfr_init(r21165); mpfr_init(r21166); mpfr_init(r21167); mpfr_init(r21168); mpfr_init(r21169); } double f_im(double a, double b) { mpfr_const_pi(r21155, MPFR_RNDN); ; mpfr_div(r21157, r21155, r21156, MPFR_RNDN); ; mpfr_set_d(r21159, b, MPFR_RNDN); mpfr_mul(r21160, r21159, r21159, MPFR_RNDN); mpfr_set_d(r21161, a, MPFR_RNDN); mpfr_mul(r21162, r21161, r21161, MPFR_RNDN); mpfr_sub(r21163, r21160, r21162, MPFR_RNDN); mpfr_div(r21164, r21158, r21163, MPFR_RNDN); mpfr_mul(r21165, r21157, r21164, MPFR_RNDN); mpfr_div(r21166, r21158, r21161, MPFR_RNDN); mpfr_div(r21167, r21158, r21159, MPFR_RNDN); mpfr_sub(r21168, r21166, r21167, MPFR_RNDN); mpfr_mul(r21169, r21165, r21168, MPFR_RNDN); return mpfr_get_d(r21169, MPFR_RNDN); } static mpfr_t r21170, r21171, r21172, r21173, r21174, r21175, r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21170, "1", 10, MPFR_RNDN); mpfr_init(r21171); mpfr_init(r21172); mpfr_init(r21173); mpfr_init(r21174); mpfr_init(r21175); mpfr_init(r21176); mpfr_init(r21177); mpfr_init(r21178); mpfr_init(r21179); mpfr_init(r21180); mpfr_init_set_str(r21181, "2", 10, MPFR_RNDN); mpfr_init(r21182); mpfr_init(r21183); } double f_fm(double a, double b) { ; mpfr_set_d(r21171, a, MPFR_RNDN); mpfr_div(r21172, r21170, r21171, MPFR_RNDN); mpfr_set_d(r21173, b, MPFR_RNDN); mpfr_div(r21174, r21170, r21173, MPFR_RNDN); mpfr_sub(r21175, r21172, r21174, MPFR_RNDN); mpfr_add(r21176, r21173, r21171, MPFR_RNDN); mpfr_div(r21177, r21175, r21176, MPFR_RNDN); mpfr_sub(r21178, r21173, r21171, MPFR_RNDN); mpfr_div(r21179, r21177, r21178, MPFR_RNDN); mpfr_const_pi(r21180, MPFR_RNDN); ; mpfr_div(r21182, r21180, r21181, MPFR_RNDN); mpfr_mul(r21183, r21179, r21182, MPFR_RNDN); return mpfr_get_d(r21183, MPFR_RNDN); } static mpfr_t r21184, r21185, r21186, r21187, r21188, r21189, r21190, r21191, r21192, r21193, r21194, r21195, r21196, r21197; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21184, "1", 10, MPFR_RNDN); mpfr_init(r21185); mpfr_init(r21186); mpfr_init(r21187); mpfr_init(r21188); mpfr_init(r21189); mpfr_init(r21190); mpfr_init(r21191); mpfr_init(r21192); mpfr_init(r21193); mpfr_init(r21194); mpfr_init_set_str(r21195, "2", 10, MPFR_RNDN); mpfr_init(r21196); mpfr_init(r21197); } double f_dm(double a, double b) { ; mpfr_set_d(r21185, a, MPFR_RNDN); mpfr_div(r21186, r21184, r21185, MPFR_RNDN); mpfr_set_d(r21187, b, MPFR_RNDN); mpfr_div(r21188, r21184, r21187, MPFR_RNDN); mpfr_sub(r21189, r21186, r21188, MPFR_RNDN); mpfr_add(r21190, r21187, r21185, MPFR_RNDN); mpfr_div(r21191, r21189, r21190, MPFR_RNDN); mpfr_sub(r21192, r21187, r21185, MPFR_RNDN); mpfr_div(r21193, r21191, r21192, MPFR_RNDN); mpfr_const_pi(r21194, MPFR_RNDN); ; mpfr_div(r21196, r21194, r21195, MPFR_RNDN); mpfr_mul(r21197, r21193, r21196, MPFR_RNDN); return mpfr_get_d(r21197, MPFR_RNDN); }