#include #include #include #include #include char *name = "NMSE Section 6.1 mentioned, B"; double f_if(float a, float b) { float r20792 = atan2(1.0, 0.0); float r20793 = 2; float r20794 = r20792 / r20793; float r20795 = 1; float r20796 = b; float r20797 = r20796 * r20796; float r20798 = a; float r20799 = r20798 * r20798; float r20800 = r20797 - r20799; float r20801 = r20795 / r20800; float r20802 = r20794 * r20801; float r20803 = r20795 / r20798; float r20804 = r20795 / r20796; float r20805 = r20803 - r20804; float r20806 = r20802 * r20805; return r20806; } double f_id(double a, double b) { double r20807 = atan2(1.0, 0.0); double r20808 = 2; double r20809 = r20807 / r20808; double r20810 = 1; double r20811 = b; double r20812 = r20811 * r20811; double r20813 = a; double r20814 = r20813 * r20813; double r20815 = r20812 - r20814; double r20816 = r20810 / r20815; double r20817 = r20809 * r20816; double r20818 = r20810 / r20813; double r20819 = r20810 / r20811; double r20820 = r20818 - r20819; double r20821 = r20817 * r20820; return r20821; } double f_of(float a, float b) { float r20822 = 1; float r20823 = a; float r20824 = r20822 / r20823; float r20825 = b; float r20826 = r20822 / r20825; float r20827 = r20824 - r20826; float r20828 = r20825 + r20823; float r20829 = r20827 / r20828; float r20830 = atan2(1.0, 0.0); float r20831 = r20829 * r20830; float r20832 = r20825 - r20823; float r20833 = 2; float r20834 = r20832 * r20833; float r20835 = r20831 / r20834; return r20835; } double f_od(double a, double b) { double r20836 = 1; double r20837 = a; double r20838 = r20836 / r20837; double r20839 = b; double r20840 = r20836 / r20839; double r20841 = r20838 - r20840; double r20842 = r20839 + r20837; double r20843 = r20841 / r20842; double r20844 = atan2(1.0, 0.0); double r20845 = r20843 * r20844; double r20846 = r20839 - r20837; double r20847 = 2; double r20848 = r20846 * r20847; double r20849 = r20845 / r20848; return r20849; } 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 r20850, r20851, r20852, r20853, r20854, r20855, r20856, r20857, r20858, r20859, r20860, r20861, r20862, r20863, r20864; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r20850); mpfr_init_set_str(r20851, "2", 10, MPFR_RNDN); mpfr_init(r20852); mpfr_init_set_str(r20853, "1", 10, MPFR_RNDN); mpfr_init(r20854); mpfr_init(r20855); mpfr_init(r20856); mpfr_init(r20857); mpfr_init(r20858); mpfr_init(r20859); mpfr_init(r20860); mpfr_init(r20861); mpfr_init(r20862); mpfr_init(r20863); mpfr_init(r20864); } double f_im(double a, double b) { mpfr_const_pi(r20850, MPFR_RNDN); ; mpfr_div(r20852, r20850, r20851, MPFR_RNDN); ; mpfr_set_d(r20854, b, MPFR_RNDN); mpfr_mul(r20855, r20854, r20854, MPFR_RNDN); mpfr_set_d(r20856, a, MPFR_RNDN); mpfr_mul(r20857, r20856, r20856, MPFR_RNDN); mpfr_sub(r20858, r20855, r20857, MPFR_RNDN); mpfr_div(r20859, r20853, r20858, MPFR_RNDN); mpfr_mul(r20860, r20852, r20859, MPFR_RNDN); mpfr_div(r20861, r20853, r20856, MPFR_RNDN); mpfr_div(r20862, r20853, r20854, MPFR_RNDN); mpfr_sub(r20863, r20861, r20862, MPFR_RNDN); mpfr_mul(r20864, r20860, r20863, MPFR_RNDN); return mpfr_get_d(r20864, MPFR_RNDN); } static mpfr_t r20865, r20866, r20867, r20868, r20869, r20870, r20871, r20872, r20873, r20874, r20875, r20876, r20877, r20878; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20865, "1", 10, MPFR_RNDN); mpfr_init(r20866); mpfr_init(r20867); mpfr_init(r20868); mpfr_init(r20869); mpfr_init(r20870); mpfr_init(r20871); mpfr_init(r20872); mpfr_init(r20873); mpfr_init(r20874); mpfr_init(r20875); mpfr_init_set_str(r20876, "2", 10, MPFR_RNDN); mpfr_init(r20877); mpfr_init(r20878); } double f_fm(double a, double b) { ; mpfr_set_d(r20866, a, MPFR_RNDN); mpfr_div(r20867, r20865, r20866, MPFR_RNDN); mpfr_set_d(r20868, b, MPFR_RNDN); mpfr_div(r20869, r20865, r20868, MPFR_RNDN); mpfr_sub(r20870, r20867, r20869, MPFR_RNDN); mpfr_add(r20871, r20868, r20866, MPFR_RNDN); mpfr_div(r20872, r20870, r20871, MPFR_RNDN); mpfr_const_pi(r20873, MPFR_RNDN); mpfr_mul(r20874, r20872, r20873, MPFR_RNDN); mpfr_sub(r20875, r20868, r20866, MPFR_RNDN); ; mpfr_mul(r20877, r20875, r20876, MPFR_RNDN); mpfr_div(r20878, r20874, r20877, MPFR_RNDN); return mpfr_get_d(r20878, MPFR_RNDN); } static mpfr_t r20879, r20880, r20881, r20882, r20883, r20884, r20885, r20886, r20887, r20888, r20889, r20890, r20891, r20892; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20879, "1", 10, MPFR_RNDN); mpfr_init(r20880); mpfr_init(r20881); mpfr_init(r20882); mpfr_init(r20883); mpfr_init(r20884); mpfr_init(r20885); mpfr_init(r20886); mpfr_init(r20887); mpfr_init(r20888); mpfr_init(r20889); mpfr_init_set_str(r20890, "2", 10, MPFR_RNDN); mpfr_init(r20891); mpfr_init(r20892); } double f_dm(double a, double b) { ; mpfr_set_d(r20880, a, MPFR_RNDN); mpfr_div(r20881, r20879, r20880, MPFR_RNDN); mpfr_set_d(r20882, b, MPFR_RNDN); mpfr_div(r20883, r20879, r20882, MPFR_RNDN); mpfr_sub(r20884, r20881, r20883, MPFR_RNDN); mpfr_add(r20885, r20882, r20880, MPFR_RNDN); mpfr_div(r20886, r20884, r20885, MPFR_RNDN); mpfr_const_pi(r20887, MPFR_RNDN); mpfr_mul(r20888, r20886, r20887, MPFR_RNDN); mpfr_sub(r20889, r20882, r20880, MPFR_RNDN); ; mpfr_mul(r20891, r20889, r20890, MPFR_RNDN); mpfr_div(r20892, r20888, r20891, MPFR_RNDN); return mpfr_get_d(r20892, MPFR_RNDN); }