#include #include #include #include #include char *name = "NMSE problem 3.4.5"; double f_if(float x) { float r14767 = x; float r14768 = sin(r14767); float r14769 = r14767 - r14768; float r14770 = tan(r14767); float r14771 = r14767 - r14770; float r14772 = r14769 / r14771; return r14772; } double f_id(double x) { double r14773 = x; double r14774 = sin(r14773); double r14775 = r14773 - r14774; double r14776 = tan(r14773); double r14777 = r14773 - r14776; double r14778 = r14775 / r14777; return r14778; } double f_of(float x) { float r14779 = x; float r14780 = r14779 * r14779; float r14781 = 0.225f; float r14782 = r14780 * r14781; float r14783 = r14782 * (r14782 * r14782); float r14784 = 0.009642857142857142f; float r14785 = 4.0f; float r14786 = pow(r14779, r14785); float r14787 = 0.5f; float r14788 = fma(r14784, r14786, r14787); float r14789 = r14788 * (r14788 * r14788); float r14790 = r14783 - r14789; float r14791 = r14781 * r14780; float r14792 = r14791 * r14791; float r14793 = r14784 * r14786; float r14794 = r14793 + r14787; float r14795 = r14794 * r14794; float r14796 = r14791 * r14794; float r14797 = r14795 + r14796; float r14798 = r14792 + r14797; float r14799 = r14790 / r14798; return r14799; } double f_od(double x) { double r14800 = x; double r14801 = r14800 * r14800; double r14802 = 0.225; double r14803 = r14801 * r14802; double r14804 = r14803 * (r14803 * r14803); double r14805 = 0.009642857142857142; double r14806 = 4.0; double r14807 = pow(r14800, r14806); double r14808 = 0.5; double r14809 = fma(r14805, r14807, r14808); double r14810 = r14809 * (r14809 * r14809); double r14811 = r14804 - r14810; double r14812 = r14802 * r14801; double r14813 = r14812 * r14812; double r14814 = r14805 * r14807; double r14815 = r14814 + r14808; double r14816 = r14815 * r14815; double r14817 = r14812 * r14815; double r14818 = r14816 + r14817; double r14819 = r14813 + r14818; double r14820 = r14811 / r14819; return r14820; } 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 r14821, r14822, r14823, r14824, r14825, r14826; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r14821); mpfr_init(r14822); mpfr_init(r14823); mpfr_init(r14824); mpfr_init(r14825); mpfr_init(r14826); } double f_im(double x) { mpfr_set_d(r14821, x, MPFR_RNDN); mpfr_sin(r14822, r14821, MPFR_RNDN); mpfr_sub(r14823, r14821, r14822, MPFR_RNDN); mpfr_tan(r14824, r14821, MPFR_RNDN); mpfr_sub(r14825, r14821, r14824, MPFR_RNDN); mpfr_div(r14826, r14823, r14825, MPFR_RNDN); return mpfr_get_d(r14826, MPFR_RNDN); } static mpfr_t r14827, r14828, r14829, r14830, r14831, r14832, r14833, r14834, r14835, r14836, r14837, r14838, r14839, r14840, r14841, r14842, r14843, r14844, r14845, r14846, r14847; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r14827); mpfr_init(r14828); mpfr_init_set_str(r14829, "9/40", 10, MPFR_RNDN); mpfr_init(r14830); mpfr_init(r14831); mpfr_init_set_str(r14832, "27/2800", 10, MPFR_RNDN); mpfr_init_set_str(r14833, "4", 10, MPFR_RNDN); mpfr_init(r14834); mpfr_init_set_str(r14835, "1/2", 10, MPFR_RNDN); mpfr_init(r14836); mpfr_init(r14837); mpfr_init(r14838); mpfr_init(r14839); mpfr_init(r14840); mpfr_init(r14841); mpfr_init(r14842); mpfr_init(r14843); mpfr_init(r14844); mpfr_init(r14845); mpfr_init(r14846); mpfr_init(r14847); } double f_fm(double x) { mpfr_set_d(r14827, x, MPFR_RNDN); mpfr_sqr(r14828, r14827, MPFR_RNDN); ; mpfr_mul(r14830, r14828, r14829, MPFR_RNDN); mpfr_mul(r14831, r14830, r14830, MPFR_RNDN); mpfr_mul(r14831, r14831, r14830, MPFR_RNDN); ; ; mpfr_pow(r14834, r14827, r14833, MPFR_RNDN); ; mpfr_fma(r14836, r14832, r14834, r14835, MPFR_RNDN); mpfr_mul(r14837, r14836, r14836, MPFR_RNDN); mpfr_mul(r14837, r14837, r14836, MPFR_RNDN); mpfr_sub(r14838, r14831, r14837, MPFR_RNDN); mpfr_mul(r14839, r14829, r14828, MPFR_RNDN); mpfr_sqr(r14840, r14839, MPFR_RNDN); mpfr_mul(r14841, r14832, r14834, MPFR_RNDN); mpfr_add(r14842, r14841, r14835, MPFR_RNDN); mpfr_sqr(r14843, r14842, MPFR_RNDN); mpfr_mul(r14844, r14839, r14842, MPFR_RNDN); mpfr_add(r14845, r14843, r14844, MPFR_RNDN); mpfr_add(r14846, r14840, r14845, MPFR_RNDN); mpfr_div(r14847, r14838, r14846, MPFR_RNDN); return mpfr_get_d(r14847, MPFR_RNDN); } static mpfr_t r14848, r14849, r14850, r14851, r14852, r14853, r14854, r14855, r14856, r14857, r14858, r14859, r14860, r14861, r14862, r14863, r14864, r14865, r14866, r14867, r14868; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r14848); mpfr_init(r14849); mpfr_init_set_str(r14850, "9/40", 10, MPFR_RNDN); mpfr_init(r14851); mpfr_init(r14852); mpfr_init_set_str(r14853, "27/2800", 10, MPFR_RNDN); mpfr_init_set_str(r14854, "4", 10, MPFR_RNDN); mpfr_init(r14855); mpfr_init_set_str(r14856, "1/2", 10, MPFR_RNDN); mpfr_init(r14857); mpfr_init(r14858); mpfr_init(r14859); mpfr_init(r14860); mpfr_init(r14861); mpfr_init(r14862); mpfr_init(r14863); mpfr_init(r14864); mpfr_init(r14865); mpfr_init(r14866); mpfr_init(r14867); mpfr_init(r14868); } double f_dm(double x) { mpfr_set_d(r14848, x, MPFR_RNDN); mpfr_sqr(r14849, r14848, MPFR_RNDN); ; mpfr_mul(r14851, r14849, r14850, MPFR_RNDN); mpfr_mul(r14852, r14851, r14851, MPFR_RNDN); mpfr_mul(r14852, r14852, r14851, MPFR_RNDN); ; ; mpfr_pow(r14855, r14848, r14854, MPFR_RNDN); ; mpfr_fma(r14857, r14853, r14855, r14856, MPFR_RNDN); mpfr_mul(r14858, r14857, r14857, MPFR_RNDN); mpfr_mul(r14858, r14858, r14857, MPFR_RNDN); mpfr_sub(r14859, r14852, r14858, MPFR_RNDN); mpfr_mul(r14860, r14850, r14849, MPFR_RNDN); mpfr_sqr(r14861, r14860, MPFR_RNDN); mpfr_mul(r14862, r14853, r14855, MPFR_RNDN); mpfr_add(r14863, r14862, r14856, MPFR_RNDN); mpfr_sqr(r14864, r14863, MPFR_RNDN); mpfr_mul(r14865, r14860, r14863, MPFR_RNDN); mpfr_add(r14866, r14864, r14865, MPFR_RNDN); mpfr_add(r14867, r14861, r14866, MPFR_RNDN); mpfr_div(r14868, r14859, r14867, MPFR_RNDN); return mpfr_get_d(r14868, MPFR_RNDN); }