#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (25)"; double f_if(float a, float b) { float r17998 = a; float r17999 = r17998 * r17998; float r18000 = b; float r18001 = r18000 * r18000; float r18002 = r17999 + r18001; float r18003 = r18002 * r18002; float r18004 = 4.0f; float r18005 = 1.0f; float r18006 = r18005 + r17998; float r18007 = r17999 * r18006; float r18008 = 3.0f; float r18009 = r18008 * r17998; float r18010 = r18005 - r18009; float r18011 = r18001 * r18010; float r18012 = r18007 + r18011; float r18013 = r18004 * r18012; float r18014 = r18003 + r18013; float r18015 = r18014 - r18005; return r18015; } double f_id(double a, double b) { double r18016 = a; double r18017 = r18016 * r18016; double r18018 = b; double r18019 = r18018 * r18018; double r18020 = r18017 + r18019; double r18021 = r18020 * r18020; double r18022 = 4.0; double r18023 = 1.0; double r18024 = r18023 + r18016; double r18025 = r18017 * r18024; double r18026 = 3.0; double r18027 = r18026 * r18016; double r18028 = r18023 - r18027; double r18029 = r18019 * r18028; double r18030 = r18025 + r18029; double r18031 = r18022 * r18030; double r18032 = r18021 + r18031; double r18033 = r18032 - r18023; return r18033; } double f_of(float a, float b) { float r18034 = b; float r18035 = 4.0f; float r18036 = pow(r18034, r18035); float r18037 = 2.0f; float r18038 = r18034 * r18034; float r18039 = a; float r18040 = r18039 * r18039; float r18041 = r18038 * r18040; float r18042 = r18037 * r18041; float r18043 = pow(r18039, r18035); float r18044 = r18042 + r18043; float r18045 = r18036 + r18044; float r18046 = 1.0f; float r18047 = r18046 + r18039; float r18048 = r18040 * r18047; float r18049 = 3.0f; float r18050 = r18049 * r18039; float r18051 = r18046 - r18050; float r18052 = r18038 * r18051; float r18053 = r18048 + r18052; float r18054 = r18035 * r18053; float r18055 = r18045 + r18054; float r18056 = r18055 - r18046; return r18056; } double f_od(double a, double b) { double r18057 = b; double r18058 = 4.0; double r18059 = pow(r18057, r18058); double r18060 = 2.0; double r18061 = r18057 * r18057; double r18062 = a; double r18063 = r18062 * r18062; double r18064 = r18061 * r18063; double r18065 = r18060 * r18064; double r18066 = pow(r18062, r18058); double r18067 = r18065 + r18066; double r18068 = r18059 + r18067; double r18069 = 1.0; double r18070 = r18069 + r18062; double r18071 = r18063 * r18070; double r18072 = 3.0; double r18073 = r18072 * r18062; double r18074 = r18069 - r18073; double r18075 = r18061 * r18074; double r18076 = r18071 + r18075; double r18077 = r18058 * r18076; double r18078 = r18068 + r18077; double r18079 = r18078 - r18069; return r18079; } 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 r18080, r18081, r18082, r18083, r18084, r18085, r18086, r18087, r18088, r18089, r18090, r18091, r18092, r18093, r18094, r18095, r18096, r18097; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r18080); mpfr_init(r18081); mpfr_init(r18082); mpfr_init(r18083); mpfr_init(r18084); mpfr_init(r18085); mpfr_init_set_str(r18086, "4", 10, MPFR_RNDN); mpfr_init_set_str(r18087, "1", 10, MPFR_RNDN); mpfr_init(r18088); mpfr_init(r18089); mpfr_init_set_str(r18090, "3", 10, MPFR_RNDN); mpfr_init(r18091); mpfr_init(r18092); mpfr_init(r18093); mpfr_init(r18094); mpfr_init(r18095); mpfr_init(r18096); mpfr_init(r18097); } double f_im(double a, double b) { mpfr_set_d(r18080, a, MPFR_RNDN); mpfr_sqr(r18081, r18080, MPFR_RNDN); mpfr_set_d(r18082, b, MPFR_RNDN); mpfr_sqr(r18083, r18082, MPFR_RNDN); mpfr_add(r18084, r18081, r18083, MPFR_RNDN); mpfr_sqr(r18085, r18084, MPFR_RNDN); ; ; mpfr_add(r18088, r18087, r18080, MPFR_RNDN); mpfr_mul(r18089, r18081, r18088, MPFR_RNDN); ; mpfr_mul(r18091, r18090, r18080, MPFR_RNDN); mpfr_sub(r18092, r18087, r18091, MPFR_RNDN); mpfr_mul(r18093, r18083, r18092, MPFR_RNDN); mpfr_add(r18094, r18089, r18093, MPFR_RNDN); mpfr_mul(r18095, r18086, r18094, MPFR_RNDN); mpfr_add(r18096, r18085, r18095, MPFR_RNDN); mpfr_sub(r18097, r18096, r18087, MPFR_RNDN); return mpfr_get_d(r18097, MPFR_RNDN); } static mpfr_t r18098, r18099, r18100, r18101, r18102, r18103, r18104, r18105, r18106, r18107, r18108, r18109, r18110, r18111, r18112, r18113, r18114, r18115, r18116, r18117, r18118, r18119, r18120; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r18098); mpfr_init_set_str(r18099, "4", 10, MPFR_RNDN); mpfr_init(r18100); mpfr_init_set_str(r18101, "2", 10, MPFR_RNDN); mpfr_init(r18102); mpfr_init(r18103); mpfr_init(r18104); mpfr_init(r18105); mpfr_init(r18106); mpfr_init(r18107); mpfr_init(r18108); mpfr_init(r18109); mpfr_init_set_str(r18110, "1", 10, MPFR_RNDN); mpfr_init(r18111); mpfr_init(r18112); mpfr_init_set_str(r18113, "3", 10, MPFR_RNDN); mpfr_init(r18114); mpfr_init(r18115); mpfr_init(r18116); mpfr_init(r18117); mpfr_init(r18118); mpfr_init(r18119); mpfr_init(r18120); } double f_fm(double a, double b) { mpfr_set_d(r18098, b, MPFR_RNDN); ; mpfr_pow(r18100, r18098, r18099, MPFR_RNDN); ; mpfr_sqr(r18102, r18098, MPFR_RNDN); mpfr_set_d(r18103, a, MPFR_RNDN); mpfr_sqr(r18104, r18103, MPFR_RNDN); mpfr_mul(r18105, r18102, r18104, MPFR_RNDN); mpfr_mul(r18106, r18101, r18105, MPFR_RNDN); mpfr_pow(r18107, r18103, r18099, MPFR_RNDN); mpfr_add(r18108, r18106, r18107, MPFR_RNDN); mpfr_add(r18109, r18100, r18108, MPFR_RNDN); ; mpfr_add(r18111, r18110, r18103, MPFR_RNDN); mpfr_mul(r18112, r18104, r18111, MPFR_RNDN); ; mpfr_mul(r18114, r18113, r18103, MPFR_RNDN); mpfr_sub(r18115, r18110, r18114, MPFR_RNDN); mpfr_mul(r18116, r18102, r18115, MPFR_RNDN); mpfr_add(r18117, r18112, r18116, MPFR_RNDN); mpfr_mul(r18118, r18099, r18117, MPFR_RNDN); mpfr_add(r18119, r18109, r18118, MPFR_RNDN); mpfr_sub(r18120, r18119, r18110, MPFR_RNDN); return mpfr_get_d(r18120, MPFR_RNDN); } static mpfr_t r18121, r18122, r18123, r18124, r18125, r18126, r18127, r18128, r18129, r18130, r18131, r18132, r18133, r18134, r18135, r18136, r18137, r18138, r18139, r18140, r18141, r18142, r18143; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r18121); mpfr_init_set_str(r18122, "4", 10, MPFR_RNDN); mpfr_init(r18123); mpfr_init_set_str(r18124, "2", 10, MPFR_RNDN); mpfr_init(r18125); mpfr_init(r18126); mpfr_init(r18127); mpfr_init(r18128); mpfr_init(r18129); mpfr_init(r18130); mpfr_init(r18131); mpfr_init(r18132); mpfr_init_set_str(r18133, "1", 10, MPFR_RNDN); mpfr_init(r18134); mpfr_init(r18135); mpfr_init_set_str(r18136, "3", 10, MPFR_RNDN); mpfr_init(r18137); mpfr_init(r18138); mpfr_init(r18139); mpfr_init(r18140); mpfr_init(r18141); mpfr_init(r18142); mpfr_init(r18143); } double f_dm(double a, double b) { mpfr_set_d(r18121, b, MPFR_RNDN); ; mpfr_pow(r18123, r18121, r18122, MPFR_RNDN); ; mpfr_sqr(r18125, r18121, MPFR_RNDN); mpfr_set_d(r18126, a, MPFR_RNDN); mpfr_sqr(r18127, r18126, MPFR_RNDN); mpfr_mul(r18128, r18125, r18127, MPFR_RNDN); mpfr_mul(r18129, r18124, r18128, MPFR_RNDN); mpfr_pow(r18130, r18126, r18122, MPFR_RNDN); mpfr_add(r18131, r18129, r18130, MPFR_RNDN); mpfr_add(r18132, r18123, r18131, MPFR_RNDN); ; mpfr_add(r18134, r18133, r18126, MPFR_RNDN); mpfr_mul(r18135, r18127, r18134, MPFR_RNDN); ; mpfr_mul(r18137, r18136, r18126, MPFR_RNDN); mpfr_sub(r18138, r18133, r18137, MPFR_RNDN); mpfr_mul(r18139, r18125, r18138, MPFR_RNDN); mpfr_add(r18140, r18135, r18139, MPFR_RNDN); mpfr_mul(r18141, r18122, r18140, MPFR_RNDN); mpfr_add(r18142, r18132, r18141, MPFR_RNDN); mpfr_sub(r18143, r18142, r18133, MPFR_RNDN); return mpfr_get_d(r18143, MPFR_RNDN); }