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