#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r22115 = atan2(1.0, 0.0); float r22116 = 2; float r22117 = r22115 / r22116; float r22118 = 1; float r22119 = b; float r22120 = r22119 * r22119; float r22121 = a; float r22122 = r22121 * r22121; float r22123 = r22120 - r22122; float r22124 = r22118 / r22123; float r22125 = r22117 * r22124; float r22126 = r22118 / r22121; float r22127 = r22118 / r22119; float r22128 = r22126 - r22127; float r22129 = r22125 * r22128; return r22129; } double f_id(double a, double b) { double r22130 = atan2(1.0, 0.0); double r22131 = 2; double r22132 = r22130 / r22131; double r22133 = 1; double r22134 = b; double r22135 = r22134 * r22134; double r22136 = a; double r22137 = r22136 * r22136; double r22138 = r22135 - r22137; double r22139 = r22133 / r22138; double r22140 = r22132 * r22139; double r22141 = r22133 / r22136; double r22142 = r22133 / r22134; double r22143 = r22141 - r22142; double r22144 = r22140 * r22143; return r22144; } double f_of(float a, float b) { float r22145 = 1; float r22146 = a; float r22147 = r22145 / r22146; float r22148 = b; float r22149 = r22145 / r22148; float r22150 = r22147 - r22149; float r22151 = r22148 + r22146; float r22152 = r22150 / r22151; float r22153 = r22148 - r22146; float r22154 = r22152 / r22153; float r22155 = atan2(1.0, 0.0); float r22156 = 2; float r22157 = r22155 / r22156; float r22158 = r22154 * r22157; return r22158; } double f_od(double a, double b) { double r22159 = 1; double r22160 = a; double r22161 = r22159 / r22160; double r22162 = b; double r22163 = r22159 / r22162; double r22164 = r22161 - r22163; double r22165 = r22162 + r22160; double r22166 = r22164 / r22165; double r22167 = r22162 - r22160; double r22168 = r22166 / r22167; double r22169 = atan2(1.0, 0.0); double r22170 = 2; double r22171 = r22169 / r22170; double r22172 = r22168 * r22171; return r22172; } 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 r22173, r22174, r22175, r22176, r22177, r22178, r22179, r22180, r22181, r22182, r22183, r22184, r22185, r22186, r22187; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r22173); mpfr_init_set_str(r22174, "2", 10, MPFR_RNDN); mpfr_init(r22175); mpfr_init_set_str(r22176, "1", 10, MPFR_RNDN); mpfr_init(r22177); mpfr_init(r22178); mpfr_init(r22179); mpfr_init(r22180); mpfr_init(r22181); mpfr_init(r22182); mpfr_init(r22183); mpfr_init(r22184); mpfr_init(r22185); mpfr_init(r22186); mpfr_init(r22187); } double f_im(double a, double b) { mpfr_const_pi(r22173, MPFR_RNDN); ; mpfr_div(r22175, r22173, r22174, MPFR_RNDN); ; mpfr_set_d(r22177, b, MPFR_RNDN); mpfr_mul(r22178, r22177, r22177, MPFR_RNDN); mpfr_set_d(r22179, a, MPFR_RNDN); mpfr_mul(r22180, r22179, r22179, MPFR_RNDN); mpfr_sub(r22181, r22178, r22180, MPFR_RNDN); mpfr_div(r22182, r22176, r22181, MPFR_RNDN); mpfr_mul(r22183, r22175, r22182, MPFR_RNDN); mpfr_div(r22184, r22176, r22179, MPFR_RNDN); mpfr_div(r22185, r22176, r22177, MPFR_RNDN); mpfr_sub(r22186, r22184, r22185, MPFR_RNDN); mpfr_mul(r22187, r22183, r22186, MPFR_RNDN); return mpfr_get_d(r22187, MPFR_RNDN); } static mpfr_t r22188, r22189, r22190, r22191, r22192, r22193, r22194, r22195, r22196, r22197, r22198, r22199, r22200, r22201; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22188, "1", 10, MPFR_RNDN); mpfr_init(r22189); mpfr_init(r22190); mpfr_init(r22191); mpfr_init(r22192); mpfr_init(r22193); mpfr_init(r22194); mpfr_init(r22195); mpfr_init(r22196); mpfr_init(r22197); mpfr_init(r22198); mpfr_init_set_str(r22199, "2", 10, MPFR_RNDN); mpfr_init(r22200); mpfr_init(r22201); } double f_fm(double a, double b) { ; mpfr_set_d(r22189, a, MPFR_RNDN); mpfr_div(r22190, r22188, r22189, MPFR_RNDN); mpfr_set_d(r22191, b, MPFR_RNDN); mpfr_div(r22192, r22188, r22191, MPFR_RNDN); mpfr_sub(r22193, r22190, r22192, MPFR_RNDN); mpfr_add(r22194, r22191, r22189, MPFR_RNDN); mpfr_div(r22195, r22193, r22194, MPFR_RNDN); mpfr_sub(r22196, r22191, r22189, MPFR_RNDN); mpfr_div(r22197, r22195, r22196, MPFR_RNDN); mpfr_const_pi(r22198, MPFR_RNDN); ; mpfr_div(r22200, r22198, r22199, MPFR_RNDN); mpfr_mul(r22201, r22197, r22200, MPFR_RNDN); return mpfr_get_d(r22201, MPFR_RNDN); } static mpfr_t r22202, r22203, r22204, r22205, r22206, r22207, r22208, r22209, r22210, r22211, r22212, r22213, r22214, r22215; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22202, "1", 10, MPFR_RNDN); mpfr_init(r22203); mpfr_init(r22204); mpfr_init(r22205); mpfr_init(r22206); mpfr_init(r22207); mpfr_init(r22208); mpfr_init(r22209); mpfr_init(r22210); mpfr_init(r22211); mpfr_init(r22212); mpfr_init_set_str(r22213, "2", 10, MPFR_RNDN); mpfr_init(r22214); mpfr_init(r22215); } double f_dm(double a, double b) { ; mpfr_set_d(r22203, a, MPFR_RNDN); mpfr_div(r22204, r22202, r22203, MPFR_RNDN); mpfr_set_d(r22205, b, MPFR_RNDN); mpfr_div(r22206, r22202, r22205, MPFR_RNDN); mpfr_sub(r22207, r22204, r22206, MPFR_RNDN); mpfr_add(r22208, r22205, r22203, MPFR_RNDN); mpfr_div(r22209, r22207, r22208, MPFR_RNDN); mpfr_sub(r22210, r22205, r22203, MPFR_RNDN); mpfr_div(r22211, r22209, r22210, MPFR_RNDN); mpfr_const_pi(r22212, MPFR_RNDN); ; mpfr_div(r22214, r22212, r22213, MPFR_RNDN); mpfr_mul(r22215, r22211, r22214, MPFR_RNDN); return mpfr_get_d(r22215, MPFR_RNDN); }