#include #include #include #include #include char *name = "NMSE section 3.11"; double f_if(float x) { float r26183 = x; float r26184 = exp(r26183); float r26185 = 1; float r26186 = r26184 - r26185; float r26187 = r26184 / r26186; return r26187; } double f_id(double x) { double r26188 = x; double r26189 = exp(r26188); double r26190 = 1; double r26191 = r26189 - r26190; double r26192 = r26189 / r26191; return r26192; } double f_of(float x) { float r26193 = x; float r26194 = exp(r26193); float r26195 = 1; float r26196 = r26194 - r26195; float r26197 = r26194 / r26196; float r26198 = -0.06048001426015734; bool r26199 = r26197 <= r26198; float r26200 = 1/2; float r26201 = r26195 / r26193; float r26202 = 1/12; float r26203 = r26202 * r26193; float r26204 = r26201 + r26203; float r26205 = r26200 + r26204; float r26206 = 4644336287.380058; bool r26207 = r26197 <= r26206; float r26208 = -r26193; float r26209 = exp(r26208); float r26210 = r26195 - r26209; float r26211 = r26195 / r26210; float r26212 = r26207 ? r26211 : r26205; float r26213 = r26199 ? r26205 : r26212; return r26213; } double f_od(double x) { double r26214 = x; double r26215 = exp(r26214); double r26216 = 1; double r26217 = r26215 - r26216; double r26218 = r26215 / r26217; double r26219 = -0.06048001426015734; bool r26220 = r26218 <= r26219; double r26221 = 1/2; double r26222 = r26216 / r26214; double r26223 = 1/12; double r26224 = r26223 * r26214; double r26225 = r26222 + r26224; double r26226 = r26221 + r26225; double r26227 = 4644336287.380058; bool r26228 = r26218 <= r26227; double r26229 = -r26214; double r26230 = exp(r26229); double r26231 = r26216 - r26230; double r26232 = r26216 / r26231; double r26233 = r26228 ? r26232 : r26226; double r26234 = r26220 ? r26226 : r26233; return r26234; } 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 r26235, r26236, r26237, r26238, r26239; void setup_mpfr_f_im() { mpfr_set_default_prec(1424); mpfr_init(r26235); mpfr_init(r26236); mpfr_init_set_str(r26237, "1", 10, MPFR_RNDN); mpfr_init(r26238); mpfr_init(r26239); } double f_im(double x) { mpfr_set_d(r26235, x, MPFR_RNDN); mpfr_exp(r26236, r26235, MPFR_RNDN); ; mpfr_sub(r26238, r26236, r26237, MPFR_RNDN); mpfr_div(r26239, r26236, r26238, MPFR_RNDN); return mpfr_get_d(r26239, MPFR_RNDN); } static mpfr_t r26240, r26241, r26242, r26243, r26244, r26245, r26246, r26247, r26248, r26249, r26250, r26251, r26252, r26253, r26254, r26255, r26256, r26257, r26258, r26259, r26260; void setup_mpfr_f_fm() { mpfr_set_default_prec(1424); mpfr_init(r26240); mpfr_init(r26241); mpfr_init_set_str(r26242, "1", 10, MPFR_RNDN); mpfr_init(r26243); mpfr_init(r26244); mpfr_init_set_str(r26245, "-0.06048001426015734", 10, MPFR_RNDN); mpfr_init(r26246); mpfr_init_set_str(r26247, "1/2", 10, MPFR_RNDN); mpfr_init(r26248); mpfr_init_set_str(r26249, "1/12", 10, MPFR_RNDN); mpfr_init(r26250); mpfr_init(r26251); mpfr_init(r26252); mpfr_init_set_str(r26253, "4644336287.380058", 10, MPFR_RNDN); mpfr_init(r26254); mpfr_init(r26255); mpfr_init(r26256); mpfr_init(r26257); mpfr_init(r26258); mpfr_init(r26259); mpfr_init(r26260); } double f_fm(double x) { mpfr_set_d(r26240, x, MPFR_RNDN); mpfr_exp(r26241, r26240, MPFR_RNDN); ; mpfr_sub(r26243, r26241, r26242, MPFR_RNDN); mpfr_div(r26244, r26241, r26243, MPFR_RNDN); ; mpfr_set_si(r26246, mpfr_cmp(r26244, r26245) <= 0, MPFR_RNDN); ; mpfr_div(r26248, r26242, r26240, MPFR_RNDN); ; mpfr_mul(r26250, r26249, r26240, MPFR_RNDN); mpfr_add(r26251, r26248, r26250, MPFR_RNDN); mpfr_add(r26252, r26247, r26251, MPFR_RNDN); ; mpfr_set_si(r26254, mpfr_cmp(r26244, r26253) <= 0, MPFR_RNDN); mpfr_neg(r26255, r26240, MPFR_RNDN); mpfr_exp(r26256, r26255, MPFR_RNDN); mpfr_sub(r26257, r26242, r26256, MPFR_RNDN); mpfr_div(r26258, r26242, r26257, MPFR_RNDN); if (mpfr_get_si(r26254, MPFR_RNDN)) { mpfr_set(r26259, r26258, MPFR_RNDN); } else { mpfr_set(r26259, r26252, MPFR_RNDN); }; if (mpfr_get_si(r26246, MPFR_RNDN)) { mpfr_set(r26260, r26252, MPFR_RNDN); } else { mpfr_set(r26260, r26259, MPFR_RNDN); }; return mpfr_get_d(r26260, MPFR_RNDN); } static mpfr_t r26261, r26262, r26263, r26264, r26265, r26266, r26267, r26268, r26269, r26270, r26271, r26272, r26273, r26274, r26275, r26276, r26277, r26278, r26279, r26280, r26281; void setup_mpfr_f_dm() { mpfr_set_default_prec(1424); mpfr_init(r26261); mpfr_init(r26262); mpfr_init_set_str(r26263, "1", 10, MPFR_RNDN); mpfr_init(r26264); mpfr_init(r26265); mpfr_init_set_str(r26266, "-0.06048001426015734", 10, MPFR_RNDN); mpfr_init(r26267); mpfr_init_set_str(r26268, "1/2", 10, MPFR_RNDN); mpfr_init(r26269); mpfr_init_set_str(r26270, "1/12", 10, MPFR_RNDN); mpfr_init(r26271); mpfr_init(r26272); mpfr_init(r26273); mpfr_init_set_str(r26274, "4644336287.380058", 10, MPFR_RNDN); mpfr_init(r26275); mpfr_init(r26276); mpfr_init(r26277); mpfr_init(r26278); mpfr_init(r26279); mpfr_init(r26280); mpfr_init(r26281); } double f_dm(double x) { mpfr_set_d(r26261, x, MPFR_RNDN); mpfr_exp(r26262, r26261, MPFR_RNDN); ; mpfr_sub(r26264, r26262, r26263, MPFR_RNDN); mpfr_div(r26265, r26262, r26264, MPFR_RNDN); ; mpfr_set_si(r26267, mpfr_cmp(r26265, r26266) <= 0, MPFR_RNDN); ; mpfr_div(r26269, r26263, r26261, MPFR_RNDN); ; mpfr_mul(r26271, r26270, r26261, MPFR_RNDN); mpfr_add(r26272, r26269, r26271, MPFR_RNDN); mpfr_add(r26273, r26268, r26272, MPFR_RNDN); ; mpfr_set_si(r26275, mpfr_cmp(r26265, r26274) <= 0, MPFR_RNDN); mpfr_neg(r26276, r26261, MPFR_RNDN); mpfr_exp(r26277, r26276, MPFR_RNDN); mpfr_sub(r26278, r26263, r26277, MPFR_RNDN); mpfr_div(r26279, r26263, r26278, MPFR_RNDN); if (mpfr_get_si(r26275, MPFR_RNDN)) { mpfr_set(r26280, r26279, MPFR_RNDN); } else { mpfr_set(r26280, r26273, MPFR_RNDN); }; if (mpfr_get_si(r26267, MPFR_RNDN)) { mpfr_set(r26281, r26273, MPFR_RNDN); } else { mpfr_set(r26281, r26280, MPFR_RNDN); }; return mpfr_get_d(r26281, MPFR_RNDN); }