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