#include #include #include #include #include char *name = "NMSE problem 3.4.2"; double f_if(float a, float b, float eps) { float r15510 = eps; float r15511 = a; float r15512 = b; float r15513 = r15511 + r15512; float r15514 = r15513 * r15510; float r15515 = exp(r15514); float r15516 = 1.0f; float r15517 = r15515 - r15516; float r15518 = r15510 * r15517; float r15519 = r15511 * r15510; float r15520 = exp(r15519); float r15521 = r15520 - r15516; float r15522 = r15512 * r15510; float r15523 = exp(r15522); float r15524 = r15523 - r15516; float r15525 = r15521 * r15524; float r15526 = r15518 / r15525; return r15526; } double f_id(double a, double b, double eps) { double r15527 = eps; double r15528 = a; double r15529 = b; double r15530 = r15528 + r15529; double r15531 = r15530 * r15527; double r15532 = exp(r15531); double r15533 = 1.0; double r15534 = r15532 - r15533; double r15535 = r15527 * r15534; double r15536 = r15528 * r15527; double r15537 = exp(r15536); double r15538 = r15537 - r15533; double r15539 = r15529 * r15527; double r15540 = exp(r15539); double r15541 = r15540 - r15533; double r15542 = r15538 * r15541; double r15543 = r15535 / r15542; return r15543; } double f_of(float a, float b, float eps) { float r15544 = eps; float r15545 = b; float r15546 = a; float r15547 = r15545 + r15546; float r15548 = r15544 * r15547; float r15549 = expm1(r15548); float r15550 = r15546 * r15544; float r15551 = expm1(r15550); float r15552 = r15549 / r15551; float r15553 = 1.0f; float r15554 = r15545 * r15544; float r15555 = expm1(r15554); float r15556 = r15555 / r15544; float r15557 = r15553 / r15556; float r15558 = r15552 * r15557; return r15558; } double f_od(double a, double b, double eps) { double r15559 = eps; double r15560 = b; double r15561 = a; double r15562 = r15560 + r15561; double r15563 = r15559 * r15562; double r15564 = expm1(r15563); double r15565 = r15561 * r15559; double r15566 = expm1(r15565); double r15567 = r15564 / r15566; double r15568 = 1.0; double r15569 = r15560 * r15559; double r15570 = expm1(r15569); double r15571 = r15570 / r15559; double r15572 = r15568 / r15571; double r15573 = r15567 * r15572; return r15573; } 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 r15574, r15575, r15576, r15577, r15578, r15579, r15580, r15581, r15582, r15583, r15584, r15585, r15586, r15587, r15588, r15589, r15590; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r15574); mpfr_init(r15575); mpfr_init(r15576); mpfr_init(r15577); mpfr_init(r15578); mpfr_init(r15579); mpfr_init_set_str(r15580, "1", 10, MPFR_RNDN); mpfr_init(r15581); mpfr_init(r15582); mpfr_init(r15583); mpfr_init(r15584); mpfr_init(r15585); mpfr_init(r15586); mpfr_init(r15587); mpfr_init(r15588); mpfr_init(r15589); mpfr_init(r15590); } double f_im(double a, double b, double eps) { mpfr_set_d(r15574, eps, MPFR_RNDN); mpfr_set_d(r15575, a, MPFR_RNDN); mpfr_set_d(r15576, b, MPFR_RNDN); mpfr_add(r15577, r15575, r15576, MPFR_RNDN); mpfr_mul(r15578, r15577, r15574, MPFR_RNDN); mpfr_exp(r15579, r15578, MPFR_RNDN); ; mpfr_sub(r15581, r15579, r15580, MPFR_RNDN); mpfr_mul(r15582, r15574, r15581, MPFR_RNDN); mpfr_mul(r15583, r15575, r15574, MPFR_RNDN); mpfr_exp(r15584, r15583, MPFR_RNDN); mpfr_sub(r15585, r15584, r15580, MPFR_RNDN); mpfr_mul(r15586, r15576, r15574, MPFR_RNDN); mpfr_exp(r15587, r15586, MPFR_RNDN); mpfr_sub(r15588, r15587, r15580, MPFR_RNDN); mpfr_mul(r15589, r15585, r15588, MPFR_RNDN); mpfr_div(r15590, r15582, r15589, MPFR_RNDN); return mpfr_get_d(r15590, MPFR_RNDN); } static mpfr_t r15591, r15592, r15593, r15594, r15595, r15596, r15597, r15598, r15599, r15600, r15601, r15602, r15603, r15604, r15605; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r15591); mpfr_init(r15592); mpfr_init(r15593); mpfr_init(r15594); mpfr_init(r15595); mpfr_init(r15596); mpfr_init(r15597); mpfr_init(r15598); mpfr_init(r15599); mpfr_init_set_str(r15600, "1", 10, MPFR_RNDN); mpfr_init(r15601); mpfr_init(r15602); mpfr_init(r15603); mpfr_init(r15604); mpfr_init(r15605); } double f_fm(double a, double b, double eps) { mpfr_set_d(r15591, eps, MPFR_RNDN); mpfr_set_d(r15592, b, MPFR_RNDN); mpfr_set_d(r15593, a, MPFR_RNDN); mpfr_add(r15594, r15592, r15593, MPFR_RNDN); mpfr_mul(r15595, r15591, r15594, MPFR_RNDN); mpfr_expm1(r15596, r15595, MPFR_RNDN); mpfr_mul(r15597, r15593, r15591, MPFR_RNDN); mpfr_expm1(r15598, r15597, MPFR_RNDN); mpfr_div(r15599, r15596, r15598, MPFR_RNDN); ; mpfr_mul(r15601, r15592, r15591, MPFR_RNDN); mpfr_expm1(r15602, r15601, MPFR_RNDN); mpfr_div(r15603, r15602, r15591, MPFR_RNDN); mpfr_div(r15604, r15600, r15603, MPFR_RNDN); mpfr_mul(r15605, r15599, r15604, MPFR_RNDN); return mpfr_get_d(r15605, MPFR_RNDN); } static mpfr_t r15606, r15607, r15608, r15609, r15610, r15611, r15612, r15613, r15614, r15615, r15616, r15617, r15618, r15619, r15620; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r15606); mpfr_init(r15607); mpfr_init(r15608); mpfr_init(r15609); mpfr_init(r15610); mpfr_init(r15611); mpfr_init(r15612); mpfr_init(r15613); mpfr_init(r15614); mpfr_init_set_str(r15615, "1", 10, MPFR_RNDN); mpfr_init(r15616); mpfr_init(r15617); mpfr_init(r15618); mpfr_init(r15619); mpfr_init(r15620); } double f_dm(double a, double b, double eps) { mpfr_set_d(r15606, eps, MPFR_RNDN); mpfr_set_d(r15607, b, MPFR_RNDN); mpfr_set_d(r15608, a, MPFR_RNDN); mpfr_add(r15609, r15607, r15608, MPFR_RNDN); mpfr_mul(r15610, r15606, r15609, MPFR_RNDN); mpfr_expm1(r15611, r15610, MPFR_RNDN); mpfr_mul(r15612, r15608, r15606, MPFR_RNDN); mpfr_expm1(r15613, r15612, MPFR_RNDN); mpfr_div(r15614, r15611, r15613, MPFR_RNDN); ; mpfr_mul(r15616, r15607, r15606, MPFR_RNDN); mpfr_expm1(r15617, r15616, MPFR_RNDN); mpfr_div(r15618, r15617, r15606, MPFR_RNDN); mpfr_div(r15619, r15615, r15618, MPFR_RNDN); mpfr_mul(r15620, r15614, r15619, MPFR_RNDN); return mpfr_get_d(r15620, MPFR_RNDN); }