#include #include #include #include #include char *name = "NMSE problem 3.3.3"; double f_if(float x) { float r26616 = 1.0f; float r26617 = x; float r26618 = r26617 + r26616; float r26619 = r26616 / r26618; float r26620 = 2.0f; float r26621 = r26620 / r26617; float r26622 = r26619 - r26621; float r26623 = r26617 - r26616; float r26624 = r26616 / r26623; float r26625 = r26622 + r26624; return r26625; } double f_id(double x) { double r26626 = 1.0; double r26627 = x; double r26628 = r26627 + r26626; double r26629 = r26626 / r26628; double r26630 = 2.0; double r26631 = r26630 / r26627; double r26632 = r26629 - r26631; double r26633 = r26627 - r26626; double r26634 = r26626 / r26633; double r26635 = r26632 + r26634; return r26635; } double f_of(float x) { float r26636 = x; float r26637 = -6.667433711053936e-08f; bool r26638 = r26636 <= r26637; float r26639 = 2.0f; float r26640 = 7.0f; float r26641 = pow(r26636, r26640); float r26642 = r26639 / r26641; float r26643 = 5.0f; float r26644 = pow(r26636, r26643); float r26645 = r26639 / r26644; float r26646 = r26642 + r26645; float r26647 = 1.0f; float r26648 = r26647 / r26636; float r26649 = r26636 * r26636; float r26650 = r26639 / r26649; float r26651 = r26648 * r26650; float r26652 = r26646 + r26651; float r26653 = 1183236732512159.2f; bool r26654 = r26636 <= r26653; float r26655 = r26636 + r26647; float r26656 = r26647 / r26655; float r26657 = r26639 / r26636; float r26658 = r26656 - r26657; float r26659 = r26636 - r26647; float r26660 = r26647 / r26659; float r26661 = r26658 + r26660; float r26662 = r26654 ? r26661 : r26652; float r26663 = r26638 ? r26652 : r26662; return r26663; } double f_od(double x) { double r26664 = x; double r26665 = -6.667433711053936e-08; bool r26666 = r26664 <= r26665; double r26667 = 2.0; double r26668 = 7.0; double r26669 = pow(r26664, r26668); double r26670 = r26667 / r26669; double r26671 = 5.0; double r26672 = pow(r26664, r26671); double r26673 = r26667 / r26672; double r26674 = r26670 + r26673; double r26675 = 1.0; double r26676 = r26675 / r26664; double r26677 = r26664 * r26664; double r26678 = r26667 / r26677; double r26679 = r26676 * r26678; double r26680 = r26674 + r26679; double r26681 = 1183236732512159.2; bool r26682 = r26664 <= r26681; double r26683 = r26664 + r26675; double r26684 = r26675 / r26683; double r26685 = r26667 / r26664; double r26686 = r26684 - r26685; double r26687 = r26664 - r26675; double r26688 = r26675 / r26687; double r26689 = r26686 + r26688; double r26690 = r26682 ? r26689 : r26680; double r26691 = r26666 ? r26680 : r26690; return r26691; } 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 r26692, r26693, r26694, r26695, r26696, r26697, r26698, r26699, r26700, r26701; void setup_mpfr_f_im() { mpfr_set_default_prec(1168); mpfr_init_set_str(r26692, "1", 10, MPFR_RNDN); mpfr_init(r26693); mpfr_init(r26694); mpfr_init(r26695); mpfr_init_set_str(r26696, "2", 10, MPFR_RNDN); mpfr_init(r26697); mpfr_init(r26698); mpfr_init(r26699); mpfr_init(r26700); mpfr_init(r26701); } double f_im(double x) { ; mpfr_set_d(r26693, x, MPFR_RNDN); mpfr_add(r26694, r26693, r26692, MPFR_RNDN); mpfr_div(r26695, r26692, r26694, MPFR_RNDN); ; mpfr_div(r26697, r26696, r26693, MPFR_RNDN); mpfr_sub(r26698, r26695, r26697, MPFR_RNDN); mpfr_sub(r26699, r26693, r26692, MPFR_RNDN); mpfr_div(r26700, r26692, r26699, MPFR_RNDN); mpfr_add(r26701, r26698, r26700, MPFR_RNDN); return mpfr_get_d(r26701, MPFR_RNDN); } static mpfr_t r26702, r26703, r26704, r26705, r26706, r26707, r26708, r26709, r26710, r26711, r26712, r26713, r26714, r26715, r26716, r26717, r26718, r26719, r26720, r26721, r26722, r26723, r26724, r26725, r26726, r26727, r26728, r26729; void setup_mpfr_f_fm() { mpfr_set_default_prec(1168); mpfr_init(r26702); mpfr_init_set_str(r26703, "-6.667433711053936e-08", 10, MPFR_RNDN); mpfr_init(r26704); mpfr_init_set_str(r26705, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26706, "7", 10, MPFR_RNDN); mpfr_init(r26707); mpfr_init(r26708); mpfr_init_set_str(r26709, "5", 10, MPFR_RNDN); mpfr_init(r26710); mpfr_init(r26711); mpfr_init(r26712); mpfr_init_set_str(r26713, "1", 10, MPFR_RNDN); mpfr_init(r26714); mpfr_init(r26715); mpfr_init(r26716); mpfr_init(r26717); mpfr_init(r26718); mpfr_init_set_str(r26719, "1183236732512159.2", 10, MPFR_RNDN); mpfr_init(r26720); mpfr_init(r26721); mpfr_init(r26722); mpfr_init(r26723); mpfr_init(r26724); mpfr_init(r26725); mpfr_init(r26726); mpfr_init(r26727); mpfr_init(r26728); mpfr_init(r26729); } double f_fm(double x) { mpfr_set_d(r26702, x, MPFR_RNDN); ; mpfr_set_si(r26704, mpfr_cmp(r26702, r26703) <= 0, MPFR_RNDN); ; ; mpfr_pow(r26707, r26702, r26706, MPFR_RNDN); mpfr_div(r26708, r26705, r26707, MPFR_RNDN); ; mpfr_pow(r26710, r26702, r26709, MPFR_RNDN); mpfr_div(r26711, r26705, r26710, MPFR_RNDN); mpfr_add(r26712, r26708, r26711, MPFR_RNDN); ; mpfr_div(r26714, r26713, r26702, MPFR_RNDN); mpfr_sqr(r26715, r26702, MPFR_RNDN); mpfr_div(r26716, r26705, r26715, MPFR_RNDN); mpfr_mul(r26717, r26714, r26716, MPFR_RNDN); mpfr_add(r26718, r26712, r26717, MPFR_RNDN); ; mpfr_set_si(r26720, mpfr_cmp(r26702, r26719) <= 0, MPFR_RNDN); mpfr_add(r26721, r26702, r26713, MPFR_RNDN); mpfr_div(r26722, r26713, r26721, MPFR_RNDN); mpfr_div(r26723, r26705, r26702, MPFR_RNDN); mpfr_sub(r26724, r26722, r26723, MPFR_RNDN); mpfr_sub(r26725, r26702, r26713, MPFR_RNDN); mpfr_div(r26726, r26713, r26725, MPFR_RNDN); mpfr_add(r26727, r26724, r26726, MPFR_RNDN); if (mpfr_get_si(r26720, MPFR_RNDN)) { mpfr_set(r26728, r26727, MPFR_RNDN); } else { mpfr_set(r26728, r26718, MPFR_RNDN); }; if (mpfr_get_si(r26704, MPFR_RNDN)) { mpfr_set(r26729, r26718, MPFR_RNDN); } else { mpfr_set(r26729, r26728, MPFR_RNDN); }; return mpfr_get_d(r26729, MPFR_RNDN); } static mpfr_t r26730, r26731, r26732, r26733, r26734, r26735, r26736, r26737, r26738, r26739, r26740, r26741, r26742, r26743, r26744, r26745, r26746, r26747, r26748, r26749, r26750, r26751, r26752, r26753, r26754, r26755, r26756, r26757; void setup_mpfr_f_dm() { mpfr_set_default_prec(1168); mpfr_init(r26730); mpfr_init_set_str(r26731, "-6.667433711053936e-08", 10, MPFR_RNDN); mpfr_init(r26732); mpfr_init_set_str(r26733, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26734, "7", 10, MPFR_RNDN); mpfr_init(r26735); mpfr_init(r26736); mpfr_init_set_str(r26737, "5", 10, MPFR_RNDN); mpfr_init(r26738); mpfr_init(r26739); mpfr_init(r26740); mpfr_init_set_str(r26741, "1", 10, MPFR_RNDN); mpfr_init(r26742); mpfr_init(r26743); mpfr_init(r26744); mpfr_init(r26745); mpfr_init(r26746); mpfr_init_set_str(r26747, "1183236732512159.2", 10, MPFR_RNDN); mpfr_init(r26748); mpfr_init(r26749); mpfr_init(r26750); mpfr_init(r26751); mpfr_init(r26752); mpfr_init(r26753); mpfr_init(r26754); mpfr_init(r26755); mpfr_init(r26756); mpfr_init(r26757); } double f_dm(double x) { mpfr_set_d(r26730, x, MPFR_RNDN); ; mpfr_set_si(r26732, mpfr_cmp(r26730, r26731) <= 0, MPFR_RNDN); ; ; mpfr_pow(r26735, r26730, r26734, MPFR_RNDN); mpfr_div(r26736, r26733, r26735, MPFR_RNDN); ; mpfr_pow(r26738, r26730, r26737, MPFR_RNDN); mpfr_div(r26739, r26733, r26738, MPFR_RNDN); mpfr_add(r26740, r26736, r26739, MPFR_RNDN); ; mpfr_div(r26742, r26741, r26730, MPFR_RNDN); mpfr_sqr(r26743, r26730, MPFR_RNDN); mpfr_div(r26744, r26733, r26743, MPFR_RNDN); mpfr_mul(r26745, r26742, r26744, MPFR_RNDN); mpfr_add(r26746, r26740, r26745, MPFR_RNDN); ; mpfr_set_si(r26748, mpfr_cmp(r26730, r26747) <= 0, MPFR_RNDN); mpfr_add(r26749, r26730, r26741, MPFR_RNDN); mpfr_div(r26750, r26741, r26749, MPFR_RNDN); mpfr_div(r26751, r26733, r26730, MPFR_RNDN); mpfr_sub(r26752, r26750, r26751, MPFR_RNDN); mpfr_sub(r26753, r26730, r26741, MPFR_RNDN); mpfr_div(r26754, r26741, r26753, MPFR_RNDN); mpfr_add(r26755, r26752, r26754, MPFR_RNDN); if (mpfr_get_si(r26748, MPFR_RNDN)) { mpfr_set(r26756, r26755, MPFR_RNDN); } else { mpfr_set(r26756, r26746, MPFR_RNDN); }; if (mpfr_get_si(r26732, MPFR_RNDN)) { mpfr_set(r26757, r26746, MPFR_RNDN); } else { mpfr_set(r26757, r26756, MPFR_RNDN); }; return mpfr_get_d(r26757, MPFR_RNDN); }