#include #include #include #include #include char *name = "NMSE problem 3.3.6"; double f_if(float N) { float r15233 = N; float r15234 = 1.0f; float r15235 = r15233 + r15234; float r15236 = log(r15235); float r15237 = log(r15233); float r15238 = r15236 - r15237; return r15238; } double f_id(double N) { double r15239 = N; double r15240 = 1.0; double r15241 = r15239 + r15240; double r15242 = log(r15241); double r15243 = log(r15239); double r15244 = r15242 - r15243; return r15244; } double f_of(float N) { float r15245 = N; float r15246 = 344.2880859375f; bool r15247 = r15245 <= r15246; float r15248 = 1.0f; float r15249 = r15245 + r15248; float r15250 = log(r15249); float r15251 = cbrt(r15250); float r15252 = r15251 * (r15251 * r15251); float r15253 = log(r15245); float r15254 = r15252 - r15253; float r15255 = r15248 / r15245; float r15256 = 0.3333333432674408f; float r15257 = r15245 * (r15245 * r15245); float r15258 = r15256 / r15257; float r15259 = 0.5f; float r15260 = r15245 * r15245; float r15261 = r15259 / r15260; float r15262 = r15258 - r15261; float r15263 = r15255 + r15262; float r15264 = r15247 ? r15254 : r15263; return r15264; } double f_od(double N) { double r15265 = N; double r15266 = 344.2880859375; bool r15267 = r15265 <= r15266; double r15268 = 1.0; double r15269 = r15265 + r15268; double r15270 = log(r15269); double r15271 = cbrt(r15270); double r15272 = r15271 * (r15271 * r15271); double r15273 = log(r15265); double r15274 = r15272 - r15273; double r15275 = r15268 / r15265; double r15276 = 0.3333333432674408; double r15277 = r15265 * (r15265 * r15265); double r15278 = r15276 / r15277; double r15279 = 0.5; double r15280 = r15265 * r15265; double r15281 = r15279 / r15280; double r15282 = r15278 - r15281; double r15283 = r15275 + r15282; double r15284 = r15267 ? r15274 : r15283; return r15284; } 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 r15285, r15286, r15287, r15288, r15289, r15290; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15285); mpfr_init_set_str(r15286, "1", 10, MPFR_RNDN); mpfr_init(r15287); mpfr_init(r15288); mpfr_init(r15289); mpfr_init(r15290); } double f_im(double N) { mpfr_set_d(r15285, N, MPFR_RNDN); ; mpfr_add(r15287, r15285, r15286, MPFR_RNDN); mpfr_log(r15288, r15287, MPFR_RNDN); mpfr_log(r15289, r15285, MPFR_RNDN); mpfr_sub(r15290, r15288, r15289, MPFR_RNDN); return mpfr_get_d(r15290, MPFR_RNDN); } static mpfr_t r15291, r15292, r15293, r15294, r15295, r15296, r15297, r15298, r15299, r15300, r15301, r15302, r15303, r15304, r15305, r15306, r15307, r15308, r15309, r15310; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15291); mpfr_init_set_str(r15292, "344.2881f0", 10, MPFR_RNDN); mpfr_init(r15293); mpfr_init_set_str(r15294, "1", 10, MPFR_RNDN); mpfr_init(r15295); mpfr_init(r15296); mpfr_init(r15297); mpfr_init(r15298); mpfr_init(r15299); mpfr_init(r15300); mpfr_init(r15301); mpfr_init_set_str(r15302, "1/3", 10, MPFR_RNDN); mpfr_init(r15303); mpfr_init(r15304); mpfr_init_set_str(r15305, "1/2", 10, MPFR_RNDN); mpfr_init(r15306); mpfr_init(r15307); mpfr_init(r15308); mpfr_init(r15309); mpfr_init(r15310); } double f_fm(double N) { mpfr_set_d(r15291, N, MPFR_RNDN); ; mpfr_set_si(r15293, mpfr_cmp(r15291, r15292) <= 0, MPFR_RNDN); ; mpfr_add(r15295, r15291, r15294, MPFR_RNDN); mpfr_log(r15296, r15295, MPFR_RNDN); mpfr_cbrt(r15297, r15296, MPFR_RNDN); mpfr_mul(r15298, r15297, r15297, MPFR_RNDN); mpfr_mul(r15298, r15298, r15297, MPFR_RNDN); mpfr_log(r15299, r15291, MPFR_RNDN); mpfr_sub(r15300, r15298, r15299, MPFR_RNDN); mpfr_div(r15301, r15294, r15291, MPFR_RNDN); ; mpfr_mul(r15303, r15291, r15291, MPFR_RNDN); mpfr_mul(r15303, r15303, r15291, MPFR_RNDN); mpfr_div(r15304, r15302, r15303, MPFR_RNDN); ; mpfr_mul(r15306, r15291, r15291, MPFR_RNDN); mpfr_div(r15307, r15305, r15306, MPFR_RNDN); mpfr_sub(r15308, r15304, r15307, MPFR_RNDN); mpfr_add(r15309, r15301, r15308, MPFR_RNDN); if (mpfr_get_si(r15293, MPFR_RNDN)) { mpfr_set(r15310, r15300, MPFR_RNDN); } else { mpfr_set(r15310, r15309, MPFR_RNDN); }; return mpfr_get_d(r15310, MPFR_RNDN); } static mpfr_t r15311, r15312, r15313, r15314, r15315, r15316, r15317, r15318, r15319, r15320, r15321, r15322, r15323, r15324, r15325, r15326, r15327, r15328, r15329, r15330; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15311); mpfr_init_set_str(r15312, "344.2881f0", 10, MPFR_RNDN); mpfr_init(r15313); mpfr_init_set_str(r15314, "1", 10, MPFR_RNDN); mpfr_init(r15315); mpfr_init(r15316); mpfr_init(r15317); mpfr_init(r15318); mpfr_init(r15319); mpfr_init(r15320); mpfr_init(r15321); mpfr_init_set_str(r15322, "1/3", 10, MPFR_RNDN); mpfr_init(r15323); mpfr_init(r15324); mpfr_init_set_str(r15325, "1/2", 10, MPFR_RNDN); mpfr_init(r15326); mpfr_init(r15327); mpfr_init(r15328); mpfr_init(r15329); mpfr_init(r15330); } double f_dm(double N) { mpfr_set_d(r15311, N, MPFR_RNDN); ; mpfr_set_si(r15313, mpfr_cmp(r15311, r15312) <= 0, MPFR_RNDN); ; mpfr_add(r15315, r15311, r15314, MPFR_RNDN); mpfr_log(r15316, r15315, MPFR_RNDN); mpfr_cbrt(r15317, r15316, MPFR_RNDN); mpfr_mul(r15318, r15317, r15317, MPFR_RNDN); mpfr_mul(r15318, r15318, r15317, MPFR_RNDN); mpfr_log(r15319, r15311, MPFR_RNDN); mpfr_sub(r15320, r15318, r15319, MPFR_RNDN); mpfr_div(r15321, r15314, r15311, MPFR_RNDN); ; mpfr_mul(r15323, r15311, r15311, MPFR_RNDN); mpfr_mul(r15323, r15323, r15311, MPFR_RNDN); mpfr_div(r15324, r15322, r15323, MPFR_RNDN); ; mpfr_mul(r15326, r15311, r15311, MPFR_RNDN); mpfr_div(r15327, r15325, r15326, MPFR_RNDN); mpfr_sub(r15328, r15324, r15327, MPFR_RNDN); mpfr_add(r15329, r15321, r15328, MPFR_RNDN); if (mpfr_get_si(r15313, MPFR_RNDN)) { mpfr_set(r15330, r15320, MPFR_RNDN); } else { mpfr_set(r15330, r15329, MPFR_RNDN); }; return mpfr_get_d(r15330, MPFR_RNDN); }