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