#include #include #include #include #include char *name = "NMSE example 3.6"; double f_if(float x) { float r17130 = 1.0f; float r17131 = x; float r17132 = sqrt(r17131); float r17133 = r17130 / r17132; float r17134 = r17131 + r17130; float r17135 = sqrt(r17134); float r17136 = r17130 / r17135; float r17137 = r17133 - r17136; return r17137; } double f_id(double x) { double r17138 = 1.0; double r17139 = x; double r17140 = sqrt(r17139); double r17141 = r17138 / r17140; double r17142 = r17139 + r17138; double r17143 = sqrt(r17142); double r17144 = r17138 / r17143; double r17145 = r17141 - r17144; return r17145; } double f_of(float x) { float r17146 = x; float r17147 = 139.7125244140625f; bool r17148 = r17146 <= r17147; float r17149 = 1.0f; float r17150 = sqrt(r17146); float r17151 = r17149 / r17150; float r17152 = r17146 + r17149; float r17153 = sqrt(r17152); float r17154 = log(r17153); float r17155 = -r17154; float r17156 = exp(r17155); float r17157 = r17151 - r17156; float r17158 = 4.0f; float r17159 = pow(r17146, r17158); float r17160 = r17149 / r17159; float r17161 = r17149 / r17146; float r17162 = r17161 / r17146; float r17163 = r17160 + r17162; float r17164 = r17146 * (r17146 * r17146); float r17165 = r17149 / r17164; float r17166 = r17163 - r17165; float r17167 = r17149 + r17146; float r17168 = sqrt(r17167); float r17169 = r17149 / r17168; float r17170 = r17151 + r17169; float r17171 = r17166 / r17170; float r17172 = r17148 ? r17157 : r17171; return r17172; } double f_od(double x) { double r17173 = x; double r17174 = 139.7125244140625; bool r17175 = r17173 <= r17174; double r17176 = 1.0; double r17177 = sqrt(r17173); double r17178 = r17176 / r17177; double r17179 = r17173 + r17176; double r17180 = sqrt(r17179); double r17181 = log(r17180); double r17182 = -r17181; double r17183 = exp(r17182); double r17184 = r17178 - r17183; double r17185 = 4.0; double r17186 = pow(r17173, r17185); double r17187 = r17176 / r17186; double r17188 = r17176 / r17173; double r17189 = r17188 / r17173; double r17190 = r17187 + r17189; double r17191 = r17173 * (r17173 * r17173); double r17192 = r17176 / r17191; double r17193 = r17190 - r17192; double r17194 = r17176 + r17173; double r17195 = sqrt(r17194); double r17196 = r17176 / r17195; double r17197 = r17178 + r17196; double r17198 = r17193 / r17197; double r17199 = r17175 ? r17184 : r17198; return r17199; } 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 r17200, r17201, r17202, r17203, r17204, r17205, r17206, r17207; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r17200, "1", 10, MPFR_RNDN); mpfr_init(r17201); mpfr_init(r17202); mpfr_init(r17203); mpfr_init(r17204); mpfr_init(r17205); mpfr_init(r17206); mpfr_init(r17207); } double f_im(double x) { ; mpfr_set_d(r17201, x, MPFR_RNDN); mpfr_sqrt(r17202, r17201, MPFR_RNDN); mpfr_div(r17203, r17200, r17202, MPFR_RNDN); mpfr_add(r17204, r17201, r17200, MPFR_RNDN); mpfr_sqrt(r17205, r17204, MPFR_RNDN); mpfr_div(r17206, r17200, r17205, MPFR_RNDN); mpfr_sub(r17207, r17203, r17206, MPFR_RNDN); return mpfr_get_d(r17207, MPFR_RNDN); } static mpfr_t r17208, r17209, r17210, r17211, r17212, r17213, r17214, r17215, r17216, r17217, r17218, r17219, r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227, r17228, r17229, r17230, r17231, r17232, r17233, r17234; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17208); mpfr_init_set_str(r17209, "139.71252f0", 10, MPFR_RNDN); mpfr_init(r17210); mpfr_init_set_str(r17211, "1", 10, MPFR_RNDN); mpfr_init(r17212); mpfr_init(r17213); mpfr_init(r17214); mpfr_init(r17215); mpfr_init(r17216); mpfr_init(r17217); mpfr_init(r17218); mpfr_init(r17219); mpfr_init_set_str(r17220, "4", 10, MPFR_RNDN); mpfr_init(r17221); mpfr_init(r17222); mpfr_init(r17223); mpfr_init(r17224); mpfr_init(r17225); mpfr_init(r17226); mpfr_init(r17227); mpfr_init(r17228); mpfr_init(r17229); mpfr_init(r17230); mpfr_init(r17231); mpfr_init(r17232); mpfr_init(r17233); mpfr_init(r17234); } double f_fm(double x) { mpfr_set_d(r17208, x, MPFR_RNDN); ; mpfr_set_si(r17210, mpfr_cmp(r17208, r17209) <= 0, MPFR_RNDN); ; mpfr_sqrt(r17212, r17208, MPFR_RNDN); mpfr_div(r17213, r17211, r17212, MPFR_RNDN); mpfr_add(r17214, r17208, r17211, MPFR_RNDN); mpfr_sqrt(r17215, r17214, MPFR_RNDN); mpfr_log(r17216, r17215, MPFR_RNDN); mpfr_neg(r17217, r17216, MPFR_RNDN); mpfr_exp(r17218, r17217, MPFR_RNDN); mpfr_sub(r17219, r17213, r17218, MPFR_RNDN); ; mpfr_pow(r17221, r17208, r17220, MPFR_RNDN); mpfr_div(r17222, r17211, r17221, MPFR_RNDN); mpfr_div(r17223, r17211, r17208, MPFR_RNDN); mpfr_div(r17224, r17223, r17208, MPFR_RNDN); mpfr_add(r17225, r17222, r17224, MPFR_RNDN); mpfr_mul(r17226, r17208, r17208, MPFR_RNDN); mpfr_mul(r17226, r17226, r17208, MPFR_RNDN); mpfr_div(r17227, r17211, r17226, MPFR_RNDN); mpfr_sub(r17228, r17225, r17227, MPFR_RNDN); mpfr_add(r17229, r17211, r17208, MPFR_RNDN); mpfr_sqrt(r17230, r17229, MPFR_RNDN); mpfr_div(r17231, r17211, r17230, MPFR_RNDN); mpfr_add(r17232, r17213, r17231, MPFR_RNDN); mpfr_div(r17233, r17228, r17232, MPFR_RNDN); if (mpfr_get_si(r17210, MPFR_RNDN)) { mpfr_set(r17234, r17219, MPFR_RNDN); } else { mpfr_set(r17234, r17233, MPFR_RNDN); }; return mpfr_get_d(r17234, MPFR_RNDN); } static mpfr_t r17235, r17236, r17237, r17238, r17239, r17240, r17241, r17242, r17243, r17244, r17245, r17246, r17247, r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255, r17256, r17257, r17258, r17259, r17260, r17261; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17235); mpfr_init_set_str(r17236, "139.71252f0", 10, MPFR_RNDN); mpfr_init(r17237); mpfr_init_set_str(r17238, "1", 10, MPFR_RNDN); mpfr_init(r17239); mpfr_init(r17240); mpfr_init(r17241); mpfr_init(r17242); mpfr_init(r17243); mpfr_init(r17244); mpfr_init(r17245); mpfr_init(r17246); mpfr_init_set_str(r17247, "4", 10, MPFR_RNDN); mpfr_init(r17248); mpfr_init(r17249); mpfr_init(r17250); mpfr_init(r17251); mpfr_init(r17252); mpfr_init(r17253); mpfr_init(r17254); mpfr_init(r17255); mpfr_init(r17256); mpfr_init(r17257); mpfr_init(r17258); mpfr_init(r17259); mpfr_init(r17260); mpfr_init(r17261); } double f_dm(double x) { mpfr_set_d(r17235, x, MPFR_RNDN); ; mpfr_set_si(r17237, mpfr_cmp(r17235, r17236) <= 0, MPFR_RNDN); ; mpfr_sqrt(r17239, r17235, MPFR_RNDN); mpfr_div(r17240, r17238, r17239, MPFR_RNDN); mpfr_add(r17241, r17235, r17238, MPFR_RNDN); mpfr_sqrt(r17242, r17241, MPFR_RNDN); mpfr_log(r17243, r17242, MPFR_RNDN); mpfr_neg(r17244, r17243, MPFR_RNDN); mpfr_exp(r17245, r17244, MPFR_RNDN); mpfr_sub(r17246, r17240, r17245, MPFR_RNDN); ; mpfr_pow(r17248, r17235, r17247, MPFR_RNDN); mpfr_div(r17249, r17238, r17248, MPFR_RNDN); mpfr_div(r17250, r17238, r17235, MPFR_RNDN); mpfr_div(r17251, r17250, r17235, MPFR_RNDN); mpfr_add(r17252, r17249, r17251, MPFR_RNDN); mpfr_mul(r17253, r17235, r17235, MPFR_RNDN); mpfr_mul(r17253, r17253, r17235, MPFR_RNDN); mpfr_div(r17254, r17238, r17253, MPFR_RNDN); mpfr_sub(r17255, r17252, r17254, MPFR_RNDN); mpfr_add(r17256, r17238, r17235, MPFR_RNDN); mpfr_sqrt(r17257, r17256, MPFR_RNDN); mpfr_div(r17258, r17238, r17257, MPFR_RNDN); mpfr_add(r17259, r17240, r17258, MPFR_RNDN); mpfr_div(r17260, r17255, r17259, MPFR_RNDN); if (mpfr_get_si(r17237, MPFR_RNDN)) { mpfr_set(r17261, r17246, MPFR_RNDN); } else { mpfr_set(r17261, r17260, MPFR_RNDN); }; return mpfr_get_d(r17261, MPFR_RNDN); }