#include #include #include #include #include char *name = "NMSE problem 3.3.5"; double f_if(float x, float eps) { float r15123 = x; float r15124 = eps; float r15125 = r15123 + r15124; float r15126 = cos(r15125); float r15127 = cos(r15123); float r15128 = r15126 - r15127; return r15128; } double f_id(double x, double eps) { double r15129 = x; double r15130 = eps; double r15131 = r15129 + r15130; double r15132 = cos(r15131); double r15133 = cos(r15129); double r15134 = r15132 - r15133; return r15134; } double f_of(float x, float eps) { float r15135 = eps; float r15136 = -2.4270278086576645e-11f; bool r15137 = r15135 <= r15136; float r15138 = x; float r15139 = cos(r15138); float r15140 = cos(r15135); float r15141 = r15139 * r15140; float r15142 = sin(r15138); float r15143 = sin(r15135); float r15144 = r15142 * r15143; float r15145 = r15141 - r15144; float r15146 = r15145 - r15139; float r15147 = 1.3114572333402165e-07f; bool r15148 = r15135 <= r15147; float r15149 = 0.16666666666666666f; float r15150 = r15135 * r15149; float r15151 = r15138 * (r15138 * r15138); float r15152 = r15150 * r15151; float r15153 = 0.5f; float r15154 = fma(r15153, r15135, r15138); float r15155 = r15135 * r15154; float r15156 = r15152 - r15155; float r15157 = r15148 ? r15156 : r15146; float r15158 = r15137 ? r15146 : r15157; return r15158; } double f_od(double x, double eps) { double r15159 = eps; double r15160 = -2.4270278086576645e-11; bool r15161 = r15159 <= r15160; double r15162 = x; double r15163 = cos(r15162); double r15164 = cos(r15159); double r15165 = r15163 * r15164; double r15166 = sin(r15162); double r15167 = sin(r15159); double r15168 = r15166 * r15167; double r15169 = r15165 - r15168; double r15170 = r15169 - r15163; double r15171 = 1.3114572333402165e-07; bool r15172 = r15159 <= r15171; double r15173 = 0.16666666666666666; double r15174 = r15159 * r15173; double r15175 = r15162 * (r15162 * r15162); double r15176 = r15174 * r15175; double r15177 = 0.5; double r15178 = fma(r15177, r15159, r15162); double r15179 = r15159 * r15178; double r15180 = r15176 - r15179; double r15181 = r15172 ? r15180 : r15170; double r15182 = r15161 ? r15170 : r15181; return r15182; } 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 r15183, r15184, r15185, r15186, r15187, r15188; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15183); mpfr_init(r15184); mpfr_init(r15185); mpfr_init(r15186); mpfr_init(r15187); mpfr_init(r15188); } double f_im(double x, double eps) { mpfr_set_d(r15183, x, MPFR_RNDN); mpfr_set_d(r15184, eps, MPFR_RNDN); mpfr_add(r15185, r15183, r15184, MPFR_RNDN); mpfr_cos(r15186, r15185, MPFR_RNDN); mpfr_cos(r15187, r15183, MPFR_RNDN); mpfr_sub(r15188, r15186, r15187, MPFR_RNDN); return mpfr_get_d(r15188, MPFR_RNDN); } static mpfr_t r15189, r15190, r15191, r15192, r15193, r15194, r15195, r15196, r15197, r15198, r15199, r15200, r15201, r15202, r15203, r15204, r15205, r15206, r15207, r15208, r15209, r15210, r15211, r15212; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15189); mpfr_init_set_str(r15190, "-2.4270278086576645e-11", 10, MPFR_RNDN); mpfr_init(r15191); mpfr_init(r15192); mpfr_init(r15193); mpfr_init(r15194); mpfr_init(r15195); mpfr_init(r15196); mpfr_init(r15197); mpfr_init(r15198); mpfr_init(r15199); mpfr_init(r15200); mpfr_init_set_str(r15201, "1.3114572333402165e-07", 10, MPFR_RNDN); mpfr_init(r15202); mpfr_init_set_str(r15203, "1/6", 10, MPFR_RNDN); mpfr_init(r15204); mpfr_init(r15205); mpfr_init(r15206); mpfr_init_set_str(r15207, "1/2", 10, MPFR_RNDN); mpfr_init(r15208); mpfr_init(r15209); mpfr_init(r15210); mpfr_init(r15211); mpfr_init(r15212); } double f_fm(double x, double eps) { mpfr_set_d(r15189, eps, MPFR_RNDN); ; mpfr_set_si(r15191, mpfr_cmp(r15189, r15190) <= 0, MPFR_RNDN); mpfr_set_d(r15192, x, MPFR_RNDN); mpfr_cos(r15193, r15192, MPFR_RNDN); mpfr_cos(r15194, r15189, MPFR_RNDN); mpfr_mul(r15195, r15193, r15194, MPFR_RNDN); mpfr_sin(r15196, r15192, MPFR_RNDN); mpfr_sin(r15197, r15189, MPFR_RNDN); mpfr_mul(r15198, r15196, r15197, MPFR_RNDN); mpfr_sub(r15199, r15195, r15198, MPFR_RNDN); mpfr_sub(r15200, r15199, r15193, MPFR_RNDN); ; mpfr_set_si(r15202, mpfr_cmp(r15189, r15201) <= 0, MPFR_RNDN); ; mpfr_mul(r15204, r15189, r15203, MPFR_RNDN); mpfr_mul(r15205, r15192, r15192, MPFR_RNDN); mpfr_mul(r15205, r15205, r15192, MPFR_RNDN); mpfr_mul(r15206, r15204, r15205, MPFR_RNDN); ; mpfr_fma(r15208, r15207, r15189, r15192, MPFR_RNDN); mpfr_mul(r15209, r15189, r15208, MPFR_RNDN); mpfr_sub(r15210, r15206, r15209, MPFR_RNDN); if (mpfr_get_si(r15202, MPFR_RNDN)) { mpfr_set(r15211, r15210, MPFR_RNDN); } else { mpfr_set(r15211, r15200, MPFR_RNDN); }; if (mpfr_get_si(r15191, MPFR_RNDN)) { mpfr_set(r15212, r15200, MPFR_RNDN); } else { mpfr_set(r15212, r15211, MPFR_RNDN); }; return mpfr_get_d(r15212, MPFR_RNDN); } static mpfr_t r15213, r15214, r15215, r15216, r15217, r15218, r15219, r15220, r15221, r15222, r15223, r15224, r15225, r15226, r15227, r15228, r15229, r15230, r15231, r15232, r15233, r15234, r15235, r15236; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15213); mpfr_init_set_str(r15214, "-2.4270278086576645e-11", 10, MPFR_RNDN); mpfr_init(r15215); mpfr_init(r15216); mpfr_init(r15217); mpfr_init(r15218); mpfr_init(r15219); mpfr_init(r15220); mpfr_init(r15221); mpfr_init(r15222); mpfr_init(r15223); mpfr_init(r15224); mpfr_init_set_str(r15225, "1.3114572333402165e-07", 10, MPFR_RNDN); mpfr_init(r15226); mpfr_init_set_str(r15227, "1/6", 10, MPFR_RNDN); mpfr_init(r15228); mpfr_init(r15229); mpfr_init(r15230); mpfr_init_set_str(r15231, "1/2", 10, MPFR_RNDN); mpfr_init(r15232); mpfr_init(r15233); mpfr_init(r15234); mpfr_init(r15235); mpfr_init(r15236); } double f_dm(double x, double eps) { mpfr_set_d(r15213, eps, MPFR_RNDN); ; mpfr_set_si(r15215, mpfr_cmp(r15213, r15214) <= 0, MPFR_RNDN); mpfr_set_d(r15216, x, MPFR_RNDN); mpfr_cos(r15217, r15216, MPFR_RNDN); mpfr_cos(r15218, r15213, MPFR_RNDN); mpfr_mul(r15219, r15217, r15218, MPFR_RNDN); mpfr_sin(r15220, r15216, MPFR_RNDN); mpfr_sin(r15221, r15213, MPFR_RNDN); mpfr_mul(r15222, r15220, r15221, MPFR_RNDN); mpfr_sub(r15223, r15219, r15222, MPFR_RNDN); mpfr_sub(r15224, r15223, r15217, MPFR_RNDN); ; mpfr_set_si(r15226, mpfr_cmp(r15213, r15225) <= 0, MPFR_RNDN); ; mpfr_mul(r15228, r15213, r15227, MPFR_RNDN); mpfr_mul(r15229, r15216, r15216, MPFR_RNDN); mpfr_mul(r15229, r15229, r15216, MPFR_RNDN); mpfr_mul(r15230, r15228, r15229, MPFR_RNDN); ; mpfr_fma(r15232, r15231, r15213, r15216, MPFR_RNDN); mpfr_mul(r15233, r15213, r15232, MPFR_RNDN); mpfr_sub(r15234, r15230, r15233, MPFR_RNDN); if (mpfr_get_si(r15226, MPFR_RNDN)) { mpfr_set(r15235, r15234, MPFR_RNDN); } else { mpfr_set(r15235, r15224, MPFR_RNDN); }; if (mpfr_get_si(r15215, MPFR_RNDN)) { mpfr_set(r15236, r15224, MPFR_RNDN); } else { mpfr_set(r15236, r15235, MPFR_RNDN); }; return mpfr_get_d(r15236, MPFR_RNDN); }