#include #include #include #include #include char *name = "NMSE problem 3.3.4"; double f_if(float x) { float r15173 = x; float r15174 = 1.0f; float r15175 = r15173 + r15174; float r15176 = 3.0f; float r15177 = r15174 / r15176; float r15178 = pow(r15175, r15177); float r15179 = pow(r15173, r15177); float r15180 = r15178 - r15179; return r15180; } double f_id(double x) { double r15181 = x; double r15182 = 1.0; double r15183 = r15181 + r15182; double r15184 = 3.0; double r15185 = r15182 / r15184; double r15186 = pow(r15183, r15185); double r15187 = pow(r15181, r15185); double r15188 = r15186 - r15187; return r15188; } double f_of(float x) { float r15189 = x; float r15190 = 1.0f; float r15191 = r15189 + r15190; float r15192 = 3.0f; float r15193 = r15190 / r15192; float r15194 = pow(r15191, r15193); float r15195 = r15194 * (r15194 * r15194); float r15196 = pow(r15189, r15193); float r15197 = r15196 * (r15196 * r15196); float r15198 = sqrt(r15197); float r15199 = r15198 * r15198; float r15200 = r15195 - r15199; float r15201 = r15194 * r15194; float r15202 = r15196 * r15196; float r15203 = r15194 * r15196; float r15204 = r15202 + r15203; float r15205 = r15201 + r15204; float r15206 = r15200 / r15205; return r15206; } double f_od(double x) { double r15207 = x; double r15208 = 1.0; double r15209 = r15207 + r15208; double r15210 = 3.0; double r15211 = r15208 / r15210; double r15212 = pow(r15209, r15211); double r15213 = r15212 * (r15212 * r15212); double r15214 = pow(r15207, r15211); double r15215 = r15214 * (r15214 * r15214); double r15216 = sqrt(r15215); double r15217 = r15216 * r15216; double r15218 = r15213 - r15217; double r15219 = r15212 * r15212; double r15220 = r15214 * r15214; double r15221 = r15212 * r15214; double r15222 = r15220 + r15221; double r15223 = r15219 + r15222; double r15224 = r15218 / r15223; return r15224; } 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 r15225, r15226, r15227, r15228, r15229, r15230, r15231, r15232; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15225); mpfr_init_set_str(r15226, "1", 10, MPFR_RNDN); mpfr_init(r15227); mpfr_init_set_str(r15228, "3", 10, MPFR_RNDN); mpfr_init(r15229); mpfr_init(r15230); mpfr_init(r15231); mpfr_init(r15232); } double f_im(double x) { mpfr_set_d(r15225, x, MPFR_RNDN); ; mpfr_add(r15227, r15225, r15226, MPFR_RNDN); ; mpfr_div(r15229, r15226, r15228, MPFR_RNDN); mpfr_pow(r15230, r15227, r15229, MPFR_RNDN); mpfr_pow(r15231, r15225, r15229, MPFR_RNDN); mpfr_sub(r15232, r15230, r15231, MPFR_RNDN); return mpfr_get_d(r15232, MPFR_RNDN); } static mpfr_t r15233, r15234, r15235, r15236, r15237, r15238, r15239, r15240, r15241, r15242, r15243, r15244, r15245, r15246, r15247, r15248, r15249, r15250; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15233); mpfr_init_set_str(r15234, "1", 10, MPFR_RNDN); mpfr_init(r15235); mpfr_init_set_str(r15236, "3", 10, MPFR_RNDN); mpfr_init(r15237); mpfr_init(r15238); mpfr_init(r15239); mpfr_init(r15240); mpfr_init(r15241); mpfr_init(r15242); mpfr_init(r15243); mpfr_init(r15244); mpfr_init(r15245); mpfr_init(r15246); mpfr_init(r15247); mpfr_init(r15248); mpfr_init(r15249); mpfr_init(r15250); } double f_fm(double x) { mpfr_set_d(r15233, x, MPFR_RNDN); ; mpfr_add(r15235, r15233, r15234, MPFR_RNDN); ; mpfr_div(r15237, r15234, r15236, MPFR_RNDN); mpfr_pow(r15238, r15235, r15237, MPFR_RNDN); mpfr_mul(r15239, r15238, r15238, MPFR_RNDN); mpfr_mul(r15239, r15239, r15238, MPFR_RNDN); mpfr_pow(r15240, r15233, r15237, MPFR_RNDN); mpfr_mul(r15241, r15240, r15240, MPFR_RNDN); mpfr_mul(r15241, r15241, r15240, MPFR_RNDN); mpfr_sqrt(r15242, r15241, MPFR_RNDN); mpfr_sqr(r15243, r15242, MPFR_RNDN); mpfr_sub(r15244, r15239, r15243, MPFR_RNDN); mpfr_sqr(r15245, r15238, MPFR_RNDN); mpfr_sqr(r15246, r15240, MPFR_RNDN); mpfr_mul(r15247, r15238, r15240, MPFR_RNDN); mpfr_add(r15248, r15246, r15247, MPFR_RNDN); mpfr_add(r15249, r15245, r15248, MPFR_RNDN); mpfr_div(r15250, r15244, r15249, MPFR_RNDN); return mpfr_get_d(r15250, MPFR_RNDN); } static mpfr_t r15251, r15252, r15253, r15254, r15255, r15256, r15257, r15258, r15259, r15260, r15261, r15262, r15263, r15264, r15265, r15266, r15267, r15268; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15251); mpfr_init_set_str(r15252, "1", 10, MPFR_RNDN); mpfr_init(r15253); mpfr_init_set_str(r15254, "3", 10, MPFR_RNDN); mpfr_init(r15255); mpfr_init(r15256); mpfr_init(r15257); mpfr_init(r15258); mpfr_init(r15259); mpfr_init(r15260); mpfr_init(r15261); mpfr_init(r15262); mpfr_init(r15263); mpfr_init(r15264); mpfr_init(r15265); mpfr_init(r15266); mpfr_init(r15267); mpfr_init(r15268); } double f_dm(double x) { mpfr_set_d(r15251, x, MPFR_RNDN); ; mpfr_add(r15253, r15251, r15252, MPFR_RNDN); ; mpfr_div(r15255, r15252, r15254, MPFR_RNDN); mpfr_pow(r15256, r15253, r15255, MPFR_RNDN); mpfr_mul(r15257, r15256, r15256, MPFR_RNDN); mpfr_mul(r15257, r15257, r15256, MPFR_RNDN); mpfr_pow(r15258, r15251, r15255, MPFR_RNDN); mpfr_mul(r15259, r15258, r15258, MPFR_RNDN); mpfr_mul(r15259, r15259, r15258, MPFR_RNDN); mpfr_sqrt(r15260, r15259, MPFR_RNDN); mpfr_sqr(r15261, r15260, MPFR_RNDN); mpfr_sub(r15262, r15257, r15261, MPFR_RNDN); mpfr_sqr(r15263, r15256, MPFR_RNDN); mpfr_sqr(r15264, r15258, MPFR_RNDN); mpfr_mul(r15265, r15256, r15258, MPFR_RNDN); mpfr_add(r15266, r15264, r15265, MPFR_RNDN); mpfr_add(r15267, r15263, r15266, MPFR_RNDN); mpfr_div(r15268, r15262, r15267, MPFR_RNDN); return mpfr_get_d(r15268, MPFR_RNDN); }