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