#include #include #include #include #include char *name = "NMSE section 3.11"; double f_if(float x) { float r16097 = x; float r16098 = exp(r16097); float r16099 = 1.0f; float r16100 = r16098 - r16099; float r16101 = r16098 / r16100; return r16101; } double f_id(double x) { double r16102 = x; double r16103 = exp(r16102); double r16104 = 1.0; double r16105 = r16103 - r16104; double r16106 = r16103 / r16105; return r16106; } double f_of(float x) { float r16107 = x; float r16108 = 5.398963451385498f; bool r16109 = r16107 <= r16108; float r16110 = exp(r16107); float r16111 = r16107 * r16107; float r16112 = 0.5f; float r16113 = 0.1666666716337204f; float r16114 = r16107 * r16113; float r16115 = r16112 + r16114; float r16116 = r16111 * r16115; float r16117 = r16107 + r16116; float r16118 = r16110 / r16117; float r16119 = cbrt(r16118); float r16120 = r16119 * (r16119 * r16119); float r16121 = 1.0f; float r16122 = -r16107; float r16123 = exp(r16122); float r16124 = r16121 - r16123; float r16125 = r16121 / r16124; float r16126 = r16109 ? r16120 : r16125; return r16126; } double f_od(double x) { double r16127 = x; double r16128 = 5.398963451385498; bool r16129 = r16127 <= r16128; double r16130 = exp(r16127); double r16131 = r16127 * r16127; double r16132 = 0.5; double r16133 = 0.1666666716337204; double r16134 = r16127 * r16133; double r16135 = r16132 + r16134; double r16136 = r16131 * r16135; double r16137 = r16127 + r16136; double r16138 = r16130 / r16137; double r16139 = cbrt(r16138); double r16140 = r16139 * (r16139 * r16139); double r16141 = 1.0; double r16142 = -r16127; double r16143 = exp(r16142); double r16144 = r16141 - r16143; double r16145 = r16141 / r16144; double r16146 = r16129 ? r16140 : r16145; return r16146; } 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 r16147, r16148, r16149, r16150, r16151; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16147); mpfr_init(r16148); mpfr_init_set_str(r16149, "1", 10, MPFR_RNDN); mpfr_init(r16150); mpfr_init(r16151); } double f_im(double x) { mpfr_set_d(r16147, x, MPFR_RNDN); mpfr_exp(r16148, r16147, MPFR_RNDN); ; mpfr_sub(r16150, r16148, r16149, MPFR_RNDN); mpfr_div(r16151, r16148, r16150, MPFR_RNDN); return mpfr_get_d(r16151, MPFR_RNDN); } static mpfr_t r16152, r16153, r16154, r16155, r16156, r16157, r16158, r16159, r16160, r16161, r16162, r16163, r16164, r16165, r16166, r16167, r16168, r16169, r16170, r16171; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16152); mpfr_init_set_str(r16153, "5.3989635f0", 10, MPFR_RNDN); mpfr_init(r16154); mpfr_init(r16155); mpfr_init(r16156); mpfr_init_set_str(r16157, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r16158, "1/6", 10, MPFR_RNDN); mpfr_init(r16159); mpfr_init(r16160); mpfr_init(r16161); mpfr_init(r16162); mpfr_init(r16163); mpfr_init(r16164); mpfr_init(r16165); mpfr_init_set_str(r16166, "1", 10, MPFR_RNDN); mpfr_init(r16167); mpfr_init(r16168); mpfr_init(r16169); mpfr_init(r16170); mpfr_init(r16171); } double f_fm(double x) { mpfr_set_d(r16152, x, MPFR_RNDN); ; mpfr_set_si(r16154, mpfr_cmp(r16152, r16153) <= 0, MPFR_RNDN); mpfr_exp(r16155, r16152, MPFR_RNDN); mpfr_sqr(r16156, r16152, MPFR_RNDN); ; ; mpfr_mul(r16159, r16152, r16158, MPFR_RNDN); mpfr_add(r16160, r16157, r16159, MPFR_RNDN); mpfr_mul(r16161, r16156, r16160, MPFR_RNDN); mpfr_add(r16162, r16152, r16161, MPFR_RNDN); mpfr_div(r16163, r16155, r16162, MPFR_RNDN); mpfr_cbrt(r16164, r16163, MPFR_RNDN); mpfr_mul(r16165, r16164, r16164, MPFR_RNDN); mpfr_mul(r16165, r16165, r16164, MPFR_RNDN); ; mpfr_neg(r16167, r16152, MPFR_RNDN); mpfr_exp(r16168, r16167, MPFR_RNDN); mpfr_sub(r16169, r16166, r16168, MPFR_RNDN); mpfr_div(r16170, r16166, r16169, MPFR_RNDN); if (mpfr_get_si(r16154, MPFR_RNDN)) { mpfr_set(r16171, r16165, MPFR_RNDN); } else { mpfr_set(r16171, r16170, MPFR_RNDN); }; return mpfr_get_d(r16171, MPFR_RNDN); } static mpfr_t r16172, r16173, r16174, r16175, r16176, r16177, r16178, r16179, r16180, r16181, r16182, r16183, r16184, r16185, r16186, r16187, r16188, r16189, r16190, r16191; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16172); mpfr_init_set_str(r16173, "5.3989635f0", 10, MPFR_RNDN); mpfr_init(r16174); mpfr_init(r16175); mpfr_init(r16176); mpfr_init_set_str(r16177, "1/2", 10, MPFR_RNDN); mpfr_init_set_str(r16178, "1/6", 10, MPFR_RNDN); mpfr_init(r16179); mpfr_init(r16180); mpfr_init(r16181); mpfr_init(r16182); mpfr_init(r16183); mpfr_init(r16184); mpfr_init(r16185); mpfr_init_set_str(r16186, "1", 10, MPFR_RNDN); mpfr_init(r16187); mpfr_init(r16188); mpfr_init(r16189); mpfr_init(r16190); mpfr_init(r16191); } double f_dm(double x) { mpfr_set_d(r16172, x, MPFR_RNDN); ; mpfr_set_si(r16174, mpfr_cmp(r16172, r16173) <= 0, MPFR_RNDN); mpfr_exp(r16175, r16172, MPFR_RNDN); mpfr_sqr(r16176, r16172, MPFR_RNDN); ; ; mpfr_mul(r16179, r16172, r16178, MPFR_RNDN); mpfr_add(r16180, r16177, r16179, MPFR_RNDN); mpfr_mul(r16181, r16176, r16180, MPFR_RNDN); mpfr_add(r16182, r16172, r16181, MPFR_RNDN); mpfr_div(r16183, r16175, r16182, MPFR_RNDN); mpfr_cbrt(r16184, r16183, MPFR_RNDN); mpfr_mul(r16185, r16184, r16184, MPFR_RNDN); mpfr_mul(r16185, r16185, r16184, MPFR_RNDN); ; mpfr_neg(r16187, r16172, MPFR_RNDN); mpfr_exp(r16188, r16187, MPFR_RNDN); mpfr_sub(r16189, r16186, r16188, MPFR_RNDN); mpfr_div(r16190, r16186, r16189, MPFR_RNDN); if (mpfr_get_si(r16174, MPFR_RNDN)) { mpfr_set(r16191, r16185, MPFR_RNDN); } else { mpfr_set(r16191, r16190, MPFR_RNDN); }; return mpfr_get_d(r16191, MPFR_RNDN); }