#include #include #include #include #include char *name = "NMSE example 3.9"; double f_if(float x) { float r27666 = 1.0f; float r27667 = x; float r27668 = r27666 / r27667; float r27669 = tan(r27667); float r27670 = r27666 / r27669; float r27671 = r27668 - r27670; return r27671; } double f_id(double x) { double r27672 = 1.0; double r27673 = x; double r27674 = r27672 / r27673; double r27675 = tan(r27673); double r27676 = r27672 / r27675; double r27677 = r27674 - r27676; return r27677; } double f_of(float x) { float r27678 = x; float r27679 = -6.667433711053936e-08f; bool r27680 = r27678 <= r27679; float r27681 = 1.0f; float r27682 = r27681 / r27678; float r27683 = tan(r27678); float r27684 = r27681 / r27683; float r27685 = r27682 - r27684; float r27686 = 1183236732512159.2f; bool r27687 = r27678 <= r27686; float r27688 = 0.022222222222222223f; float r27689 = 3.0f; float r27690 = pow(r27678, r27689); float r27691 = r27688 * r27690; float r27692 = 0.0021164021164021165f; float r27693 = 5.0f; float r27694 = pow(r27678, r27693); float r27695 = r27692 * r27694; float r27696 = 0.3333333333333333f; float r27697 = r27696 * r27678; float r27698 = r27695 + r27697; float r27699 = r27691 + r27698; float r27700 = r27687 ? r27699 : r27685; float r27701 = r27680 ? r27685 : r27700; return r27701; } double f_od(double x) { double r27702 = x; double r27703 = -6.667433711053936e-08; bool r27704 = r27702 <= r27703; double r27705 = 1.0; double r27706 = r27705 / r27702; double r27707 = tan(r27702); double r27708 = r27705 / r27707; double r27709 = r27706 - r27708; double r27710 = 1183236732512159.2; bool r27711 = r27702 <= r27710; double r27712 = 0.022222222222222223; double r27713 = 3.0; double r27714 = pow(r27702, r27713); double r27715 = r27712 * r27714; double r27716 = 0.0021164021164021165; double r27717 = 5.0; double r27718 = pow(r27702, r27717); double r27719 = r27716 * r27718; double r27720 = 0.3333333333333333; double r27721 = r27720 * r27702; double r27722 = r27719 + r27721; double r27723 = r27715 + r27722; double r27724 = r27711 ? r27723 : r27709; double r27725 = r27704 ? r27709 : r27724; return r27725; } 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 r27726, r27727, r27728, r27729, r27730, r27731; void setup_mpfr_f_im() { mpfr_set_default_prec(2448); mpfr_init_set_str(r27726, "1", 10, MPFR_RNDN); mpfr_init(r27727); mpfr_init(r27728); mpfr_init(r27729); mpfr_init(r27730); mpfr_init(r27731); } double f_im(double x) { ; mpfr_set_d(r27727, x, MPFR_RNDN); mpfr_div(r27728, r27726, r27727, MPFR_RNDN); mpfr_tan(r27729, r27727, MPFR_RNDN); mpfr_div(r27730, r27726, r27729, MPFR_RNDN); mpfr_sub(r27731, r27728, r27730, MPFR_RNDN); return mpfr_get_d(r27731, MPFR_RNDN); } static mpfr_t r27732, r27733, r27734, r27735, r27736, r27737, r27738, r27739, r27740, r27741, r27742, r27743, r27744, r27745, r27746, r27747, r27748, r27749, r27750, r27751, r27752, r27753, r27754, r27755; void setup_mpfr_f_fm() { mpfr_set_default_prec(2448); mpfr_init(r27732); mpfr_init_set_str(r27733, "-6.667433711053936e-08", 10, MPFR_RNDN); mpfr_init(r27734); mpfr_init_set_str(r27735, "1", 10, MPFR_RNDN); mpfr_init(r27736); mpfr_init(r27737); mpfr_init(r27738); mpfr_init(r27739); mpfr_init_set_str(r27740, "1183236732512159.2", 10, MPFR_RNDN); mpfr_init(r27741); mpfr_init_set_str(r27742, "1/45", 10, MPFR_RNDN); mpfr_init_set_str(r27743, "3", 10, MPFR_RNDN); mpfr_init(r27744); mpfr_init(r27745); mpfr_init_set_str(r27746, "2/945", 10, MPFR_RNDN); mpfr_init_set_str(r27747, "5", 10, MPFR_RNDN); mpfr_init(r27748); mpfr_init(r27749); mpfr_init_set_str(r27750, "1/3", 10, MPFR_RNDN); mpfr_init(r27751); mpfr_init(r27752); mpfr_init(r27753); mpfr_init(r27754); mpfr_init(r27755); } double f_fm(double x) { mpfr_set_d(r27732, x, MPFR_RNDN); ; mpfr_set_si(r27734, mpfr_cmp(r27732, r27733) <= 0, MPFR_RNDN); ; mpfr_div(r27736, r27735, r27732, MPFR_RNDN); mpfr_tan(r27737, r27732, MPFR_RNDN); mpfr_div(r27738, r27735, r27737, MPFR_RNDN); mpfr_sub(r27739, r27736, r27738, MPFR_RNDN); ; mpfr_set_si(r27741, mpfr_cmp(r27732, r27740) <= 0, MPFR_RNDN); ; ; mpfr_pow(r27744, r27732, r27743, MPFR_RNDN); mpfr_mul(r27745, r27742, r27744, MPFR_RNDN); ; ; mpfr_pow(r27748, r27732, r27747, MPFR_RNDN); mpfr_mul(r27749, r27746, r27748, MPFR_RNDN); ; mpfr_mul(r27751, r27750, r27732, MPFR_RNDN); mpfr_add(r27752, r27749, r27751, MPFR_RNDN); mpfr_add(r27753, r27745, r27752, MPFR_RNDN); if (mpfr_get_si(r27741, MPFR_RNDN)) { mpfr_set(r27754, r27753, MPFR_RNDN); } else { mpfr_set(r27754, r27739, MPFR_RNDN); }; if (mpfr_get_si(r27734, MPFR_RNDN)) { mpfr_set(r27755, r27739, MPFR_RNDN); } else { mpfr_set(r27755, r27754, MPFR_RNDN); }; return mpfr_get_d(r27755, MPFR_RNDN); } static mpfr_t r27756, r27757, r27758, r27759, r27760, r27761, r27762, r27763, r27764, r27765, r27766, r27767, r27768, r27769, r27770, r27771, r27772, r27773, r27774, r27775, r27776, r27777, r27778, r27779; void setup_mpfr_f_dm() { mpfr_set_default_prec(2448); mpfr_init(r27756); mpfr_init_set_str(r27757, "-6.667433711053936e-08", 10, MPFR_RNDN); mpfr_init(r27758); mpfr_init_set_str(r27759, "1", 10, MPFR_RNDN); mpfr_init(r27760); mpfr_init(r27761); mpfr_init(r27762); mpfr_init(r27763); mpfr_init_set_str(r27764, "1183236732512159.2", 10, MPFR_RNDN); mpfr_init(r27765); mpfr_init_set_str(r27766, "1/45", 10, MPFR_RNDN); mpfr_init_set_str(r27767, "3", 10, MPFR_RNDN); mpfr_init(r27768); mpfr_init(r27769); mpfr_init_set_str(r27770, "2/945", 10, MPFR_RNDN); mpfr_init_set_str(r27771, "5", 10, MPFR_RNDN); mpfr_init(r27772); mpfr_init(r27773); mpfr_init_set_str(r27774, "1/3", 10, MPFR_RNDN); mpfr_init(r27775); mpfr_init(r27776); mpfr_init(r27777); mpfr_init(r27778); mpfr_init(r27779); } double f_dm(double x) { mpfr_set_d(r27756, x, MPFR_RNDN); ; mpfr_set_si(r27758, mpfr_cmp(r27756, r27757) <= 0, MPFR_RNDN); ; mpfr_div(r27760, r27759, r27756, MPFR_RNDN); mpfr_tan(r27761, r27756, MPFR_RNDN); mpfr_div(r27762, r27759, r27761, MPFR_RNDN); mpfr_sub(r27763, r27760, r27762, MPFR_RNDN); ; mpfr_set_si(r27765, mpfr_cmp(r27756, r27764) <= 0, MPFR_RNDN); ; ; mpfr_pow(r27768, r27756, r27767, MPFR_RNDN); mpfr_mul(r27769, r27766, r27768, MPFR_RNDN); ; ; mpfr_pow(r27772, r27756, r27771, MPFR_RNDN); mpfr_mul(r27773, r27770, r27772, MPFR_RNDN); ; mpfr_mul(r27775, r27774, r27756, MPFR_RNDN); mpfr_add(r27776, r27773, r27775, MPFR_RNDN); mpfr_add(r27777, r27769, r27776, MPFR_RNDN); if (mpfr_get_si(r27765, MPFR_RNDN)) { mpfr_set(r27778, r27777, MPFR_RNDN); } else { mpfr_set(r27778, r27763, MPFR_RNDN); }; if (mpfr_get_si(r27758, MPFR_RNDN)) { mpfr_set(r27779, r27763, MPFR_RNDN); } else { mpfr_set(r27779, r27778, MPFR_RNDN); }; return mpfr_get_d(r27779, MPFR_RNDN); }