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