#include #include #include #include #include char *name = "NMSE example 3.6"; double f_if(float x) { float r17148 = 1.0f; float r17149 = x; float r17150 = sqrt(r17149); float r17151 = r17148 / r17150; float r17152 = r17149 + r17148; float r17153 = sqrt(r17152); float r17154 = r17148 / r17153; float r17155 = r17151 - r17154; return r17155; } double f_id(double x) { double r17156 = 1.0; double r17157 = x; double r17158 = sqrt(r17157); double r17159 = r17156 / r17158; double r17160 = r17157 + r17156; double r17161 = sqrt(r17160); double r17162 = r17156 / r17161; double r17163 = r17159 - r17162; return r17163; } double f_of(float x) { float r17164 = x; float r17165 = 1.1034063900193587e+17f; bool r17166 = r17164 <= r17165; float r17167 = 1.0f; float r17168 = fma(r17164, r17164, r17164); float r17169 = r17167 / r17168; float r17170 = sqrt(r17164); float r17171 = r17167 / r17170; float r17172 = r17164 + r17167; float r17173 = sqrt(r17172); float r17174 = r17167 / r17173; float r17175 = r17171 + r17174; float r17176 = r17169 / r17175; float r17177 = 4.0f; float r17178 = pow(r17164, r17177); float r17179 = r17167 / r17178; float r17180 = r17167 / r17164; float r17181 = r17180 / r17164; float r17182 = r17179 + r17181; float r17183 = r17164 * (r17164 * r17164); float r17184 = r17167 / r17183; float r17185 = r17182 - r17184; float r17186 = r17167 + r17164; float r17187 = sqrt(r17186); float r17188 = r17167 / r17187; float r17189 = r17171 + r17188; float r17190 = r17185 / r17189; float r17191 = r17166 ? r17176 : r17190; return r17191; } double f_od(double x) { double r17192 = x; double r17193 = 1.1034063900193587e+17; bool r17194 = r17192 <= r17193; double r17195 = 1.0; double r17196 = fma(r17192, r17192, r17192); double r17197 = r17195 / r17196; double r17198 = sqrt(r17192); double r17199 = r17195 / r17198; double r17200 = r17192 + r17195; double r17201 = sqrt(r17200); double r17202 = r17195 / r17201; double r17203 = r17199 + r17202; double r17204 = r17197 / r17203; double r17205 = 4.0; double r17206 = pow(r17192, r17205); double r17207 = r17195 / r17206; double r17208 = r17195 / r17192; double r17209 = r17208 / r17192; double r17210 = r17207 + r17209; double r17211 = r17192 * (r17192 * r17192); double r17212 = r17195 / r17211; double r17213 = r17210 - r17212; double r17214 = r17195 + r17192; double r17215 = sqrt(r17214); double r17216 = r17195 / r17215; double r17217 = r17199 + r17216; double r17218 = r17213 / r17217; double r17219 = r17194 ? r17204 : r17218; return r17219; } 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 r17220, r17221, r17222, r17223, r17224, r17225, r17226, r17227; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r17220, "1", 10, MPFR_RNDN); mpfr_init(r17221); mpfr_init(r17222); mpfr_init(r17223); mpfr_init(r17224); mpfr_init(r17225); mpfr_init(r17226); mpfr_init(r17227); } double f_im(double x) { ; mpfr_set_d(r17221, x, MPFR_RNDN); mpfr_sqrt(r17222, r17221, MPFR_RNDN); mpfr_div(r17223, r17220, r17222, MPFR_RNDN); mpfr_add(r17224, r17221, r17220, MPFR_RNDN); mpfr_sqrt(r17225, r17224, MPFR_RNDN); mpfr_div(r17226, r17220, r17225, MPFR_RNDN); mpfr_sub(r17227, r17223, r17226, MPFR_RNDN); return mpfr_get_d(r17227, MPFR_RNDN); } static mpfr_t r17228, r17229, r17230, r17231, r17232, r17233, r17234, r17235, r17236, r17237, r17238, r17239, r17240, r17241, r17242, r17243, r17244, r17245, r17246, r17247, r17248, r17249, r17250, r17251, r17252, r17253, r17254, r17255; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r17228); mpfr_init_set_str(r17229, "1.1034064f+17", 10, MPFR_RNDN); mpfr_init(r17230); mpfr_init_set_str(r17231, "1", 10, MPFR_RNDN); mpfr_init(r17232); mpfr_init(r17233); mpfr_init(r17234); mpfr_init(r17235); mpfr_init(r17236); mpfr_init(r17237); mpfr_init(r17238); mpfr_init(r17239); mpfr_init(r17240); mpfr_init_set_str(r17241, "4", 10, MPFR_RNDN); mpfr_init(r17242); mpfr_init(r17243); mpfr_init(r17244); mpfr_init(r17245); mpfr_init(r17246); mpfr_init(r17247); mpfr_init(r17248); mpfr_init(r17249); mpfr_init(r17250); mpfr_init(r17251); mpfr_init(r17252); mpfr_init(r17253); mpfr_init(r17254); mpfr_init(r17255); } double f_fm(double x) { mpfr_set_d(r17228, x, MPFR_RNDN); ; mpfr_set_si(r17230, mpfr_cmp(r17228, r17229) <= 0, MPFR_RNDN); ; mpfr_fma(r17232, r17228, r17228, r17228, MPFR_RNDN); mpfr_div(r17233, r17231, r17232, MPFR_RNDN); mpfr_sqrt(r17234, r17228, MPFR_RNDN); mpfr_div(r17235, r17231, r17234, MPFR_RNDN); mpfr_add(r17236, r17228, r17231, MPFR_RNDN); mpfr_sqrt(r17237, r17236, MPFR_RNDN); mpfr_div(r17238, r17231, r17237, MPFR_RNDN); mpfr_add(r17239, r17235, r17238, MPFR_RNDN); mpfr_div(r17240, r17233, r17239, MPFR_RNDN); ; mpfr_pow(r17242, r17228, r17241, MPFR_RNDN); mpfr_div(r17243, r17231, r17242, MPFR_RNDN); mpfr_div(r17244, r17231, r17228, MPFR_RNDN); mpfr_div(r17245, r17244, r17228, MPFR_RNDN); mpfr_add(r17246, r17243, r17245, MPFR_RNDN); mpfr_mul(r17247, r17228, r17228, MPFR_RNDN); mpfr_mul(r17247, r17247, r17228, MPFR_RNDN); mpfr_div(r17248, r17231, r17247, MPFR_RNDN); mpfr_sub(r17249, r17246, r17248, MPFR_RNDN); mpfr_add(r17250, r17231, r17228, MPFR_RNDN); mpfr_sqrt(r17251, r17250, MPFR_RNDN); mpfr_div(r17252, r17231, r17251, MPFR_RNDN); mpfr_add(r17253, r17235, r17252, MPFR_RNDN); mpfr_div(r17254, r17249, r17253, MPFR_RNDN); if (mpfr_get_si(r17230, MPFR_RNDN)) { mpfr_set(r17255, r17240, MPFR_RNDN); } else { mpfr_set(r17255, r17254, MPFR_RNDN); }; return mpfr_get_d(r17255, MPFR_RNDN); } static mpfr_t r17256, r17257, r17258, r17259, r17260, r17261, r17262, r17263, r17264, r17265, r17266, r17267, r17268, r17269, r17270, r17271, r17272, r17273, r17274, r17275, r17276, r17277, r17278, r17279, r17280, r17281, r17282, r17283; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r17256); mpfr_init_set_str(r17257, "1.1034064f+17", 10, MPFR_RNDN); mpfr_init(r17258); mpfr_init_set_str(r17259, "1", 10, MPFR_RNDN); mpfr_init(r17260); mpfr_init(r17261); mpfr_init(r17262); mpfr_init(r17263); mpfr_init(r17264); mpfr_init(r17265); mpfr_init(r17266); mpfr_init(r17267); mpfr_init(r17268); mpfr_init_set_str(r17269, "4", 10, MPFR_RNDN); mpfr_init(r17270); mpfr_init(r17271); mpfr_init(r17272); mpfr_init(r17273); mpfr_init(r17274); mpfr_init(r17275); mpfr_init(r17276); mpfr_init(r17277); mpfr_init(r17278); mpfr_init(r17279); mpfr_init(r17280); mpfr_init(r17281); mpfr_init(r17282); mpfr_init(r17283); } double f_dm(double x) { mpfr_set_d(r17256, x, MPFR_RNDN); ; mpfr_set_si(r17258, mpfr_cmp(r17256, r17257) <= 0, MPFR_RNDN); ; mpfr_fma(r17260, r17256, r17256, r17256, MPFR_RNDN); mpfr_div(r17261, r17259, r17260, MPFR_RNDN); mpfr_sqrt(r17262, r17256, MPFR_RNDN); mpfr_div(r17263, r17259, r17262, MPFR_RNDN); mpfr_add(r17264, r17256, r17259, MPFR_RNDN); mpfr_sqrt(r17265, r17264, MPFR_RNDN); mpfr_div(r17266, r17259, r17265, MPFR_RNDN); mpfr_add(r17267, r17263, r17266, MPFR_RNDN); mpfr_div(r17268, r17261, r17267, MPFR_RNDN); ; mpfr_pow(r17270, r17256, r17269, MPFR_RNDN); mpfr_div(r17271, r17259, r17270, MPFR_RNDN); mpfr_div(r17272, r17259, r17256, MPFR_RNDN); mpfr_div(r17273, r17272, r17256, MPFR_RNDN); mpfr_add(r17274, r17271, r17273, MPFR_RNDN); mpfr_mul(r17275, r17256, r17256, MPFR_RNDN); mpfr_mul(r17275, r17275, r17256, MPFR_RNDN); mpfr_div(r17276, r17259, r17275, MPFR_RNDN); mpfr_sub(r17277, r17274, r17276, MPFR_RNDN); mpfr_add(r17278, r17259, r17256, MPFR_RNDN); mpfr_sqrt(r17279, r17278, MPFR_RNDN); mpfr_div(r17280, r17259, r17279, MPFR_RNDN); mpfr_add(r17281, r17263, r17280, MPFR_RNDN); mpfr_div(r17282, r17277, r17281, MPFR_RNDN); if (mpfr_get_si(r17258, MPFR_RNDN)) { mpfr_set(r17283, r17268, MPFR_RNDN); } else { mpfr_set(r17283, r17282, MPFR_RNDN); }; return mpfr_get_d(r17283, MPFR_RNDN); }