#include #include #include #include #include char *name = "Bouland and Aaronson, Equation (24)"; double f_if(float a, float b) { float r25954 = a; float r25955 = r25954 * r25954; float r25956 = b; float r25957 = r25956 * r25956; float r25958 = r25955 + r25957; float r25959 = 2; float r25960 = pow(r25958, r25959); float r25961 = 4; float r25962 = 1; float r25963 = r25962 - r25954; float r25964 = r25955 * r25963; float r25965 = 3; float r25966 = r25965 + r25954; float r25967 = r25957 * r25966; float r25968 = r25964 + r25967; float r25969 = r25961 * r25968; float r25970 = r25960 + r25969; float r25971 = r25970 - r25962; return r25971; } double f_id(double a, double b) { double r25972 = a; double r25973 = r25972 * r25972; double r25974 = b; double r25975 = r25974 * r25974; double r25976 = r25973 + r25975; double r25977 = 2; double r25978 = pow(r25976, r25977); double r25979 = 4; double r25980 = 1; double r25981 = r25980 - r25972; double r25982 = r25973 * r25981; double r25983 = 3; double r25984 = r25983 + r25972; double r25985 = r25975 * r25984; double r25986 = r25982 + r25985; double r25987 = r25979 * r25986; double r25988 = r25978 + r25987; double r25989 = r25988 - r25980; return r25989; } double f_of(float a, float b) { float r25990 = 4; float r25991 = a; float r25992 = 1; float r25993 = r25992 - r25991; float r25994 = r25991 * r25993; float r25995 = 3; float r25996 = r25991 + r25995; float r25997 = b; float r25998 = r25997 * r25997; float r25999 = r25996 * r25998; float r26000 = fma(r25991, r25994, r25999); float r26001 = r25991 * r25997; float r26002 = r26001 * r26001; float r26003 = 2; float r26004 = pow(r25997, r25990); float r26005 = pow(r25991, r25990); float r26006 = r26004 + r26005; float r26007 = fma(r26002, r26003, r26006); float r26008 = fma(r25990, r26000, r26007); float r26009 = r26008 - r25992; return r26009; } double f_od(double a, double b) { double r26010 = 4; double r26011 = a; double r26012 = 1; double r26013 = r26012 - r26011; double r26014 = r26011 * r26013; double r26015 = 3; double r26016 = r26011 + r26015; double r26017 = b; double r26018 = r26017 * r26017; double r26019 = r26016 * r26018; double r26020 = fma(r26011, r26014, r26019); double r26021 = r26011 * r26017; double r26022 = r26021 * r26021; double r26023 = 2; double r26024 = pow(r26017, r26010); double r26025 = pow(r26011, r26010); double r26026 = r26024 + r26025; double r26027 = fma(r26022, r26023, r26026); double r26028 = fma(r26010, r26020, r26027); double r26029 = r26028 - r26012; return r26029; } 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 r26030, r26031, r26032, r26033, r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r26030); mpfr_init(r26031); mpfr_init(r26032); mpfr_init(r26033); mpfr_init(r26034); mpfr_init_set_str(r26035, "2", 10, MPFR_RNDN); mpfr_init(r26036); mpfr_init_set_str(r26037, "4", 10, MPFR_RNDN); mpfr_init_set_str(r26038, "1", 10, MPFR_RNDN); mpfr_init(r26039); mpfr_init(r26040); mpfr_init_set_str(r26041, "3", 10, MPFR_RNDN); mpfr_init(r26042); mpfr_init(r26043); mpfr_init(r26044); mpfr_init(r26045); mpfr_init(r26046); mpfr_init(r26047); } double f_im(double a, double b) { mpfr_set_d(r26030, a, MPFR_RNDN); mpfr_mul(r26031, r26030, r26030, MPFR_RNDN); mpfr_set_d(r26032, b, MPFR_RNDN); mpfr_mul(r26033, r26032, r26032, MPFR_RNDN); mpfr_add(r26034, r26031, r26033, MPFR_RNDN); ; mpfr_pow(r26036, r26034, r26035, MPFR_RNDN); ; ; mpfr_sub(r26039, r26038, r26030, MPFR_RNDN); mpfr_mul(r26040, r26031, r26039, MPFR_RNDN); ; mpfr_add(r26042, r26041, r26030, MPFR_RNDN); mpfr_mul(r26043, r26033, r26042, MPFR_RNDN); mpfr_add(r26044, r26040, r26043, MPFR_RNDN); mpfr_mul(r26045, r26037, r26044, MPFR_RNDN); mpfr_add(r26046, r26036, r26045, MPFR_RNDN); mpfr_sub(r26047, r26046, r26038, MPFR_RNDN); return mpfr_get_d(r26047, MPFR_RNDN); } static mpfr_t r26048, r26049, r26050, r26051, r26052, r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064, r26065, r26066, r26067; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26048, "4", 10, MPFR_RNDN); mpfr_init(r26049); mpfr_init_set_str(r26050, "1", 10, MPFR_RNDN); mpfr_init(r26051); mpfr_init(r26052); mpfr_init_set_str(r26053, "3", 10, MPFR_RNDN); mpfr_init(r26054); mpfr_init(r26055); mpfr_init(r26056); mpfr_init(r26057); mpfr_init(r26058); mpfr_init(r26059); mpfr_init(r26060); mpfr_init_set_str(r26061, "2", 10, MPFR_RNDN); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); mpfr_init(r26065); mpfr_init(r26066); mpfr_init(r26067); } double f_fm(double a, double b) { ; mpfr_set_d(r26049, a, MPFR_RNDN); ; mpfr_sub(r26051, r26050, r26049, MPFR_RNDN); mpfr_mul(r26052, r26049, r26051, MPFR_RNDN); ; mpfr_add(r26054, r26049, r26053, MPFR_RNDN); mpfr_set_d(r26055, b, MPFR_RNDN); mpfr_mul(r26056, r26055, r26055, MPFR_RNDN); mpfr_mul(r26057, r26054, r26056, MPFR_RNDN); mpfr_fma(r26058, r26049, r26052, r26057, MPFR_RNDN); mpfr_mul(r26059, r26049, r26055, MPFR_RNDN); mpfr_mul(r26060, r26059, r26059, MPFR_RNDN); ; mpfr_pow(r26062, r26055, r26048, MPFR_RNDN); mpfr_pow(r26063, r26049, r26048, MPFR_RNDN); mpfr_add(r26064, r26062, r26063, MPFR_RNDN); mpfr_fma(r26065, r26060, r26061, r26064, MPFR_RNDN); mpfr_fma(r26066, r26048, r26058, r26065, MPFR_RNDN); mpfr_sub(r26067, r26066, r26050, MPFR_RNDN); return mpfr_get_d(r26067, MPFR_RNDN); } static mpfr_t r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init_set_str(r26068, "4", 10, MPFR_RNDN); mpfr_init(r26069); mpfr_init_set_str(r26070, "1", 10, MPFR_RNDN); mpfr_init(r26071); mpfr_init(r26072); mpfr_init_set_str(r26073, "3", 10, MPFR_RNDN); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init_set_str(r26081, "2", 10, MPFR_RNDN); mpfr_init(r26082); mpfr_init(r26083); mpfr_init(r26084); mpfr_init(r26085); mpfr_init(r26086); mpfr_init(r26087); } double f_dm(double a, double b) { ; mpfr_set_d(r26069, a, MPFR_RNDN); ; mpfr_sub(r26071, r26070, r26069, MPFR_RNDN); mpfr_mul(r26072, r26069, r26071, MPFR_RNDN); ; mpfr_add(r26074, r26069, r26073, MPFR_RNDN); mpfr_set_d(r26075, b, MPFR_RNDN); mpfr_mul(r26076, r26075, r26075, MPFR_RNDN); mpfr_mul(r26077, r26074, r26076, MPFR_RNDN); mpfr_fma(r26078, r26069, r26072, r26077, MPFR_RNDN); mpfr_mul(r26079, r26069, r26075, MPFR_RNDN); mpfr_mul(r26080, r26079, r26079, MPFR_RNDN); ; mpfr_pow(r26082, r26075, r26068, MPFR_RNDN); mpfr_pow(r26083, r26069, r26068, MPFR_RNDN); mpfr_add(r26084, r26082, r26083, MPFR_RNDN); mpfr_fma(r26085, r26080, r26081, r26084, MPFR_RNDN); mpfr_fma(r26086, r26068, r26078, r26085, MPFR_RNDN); mpfr_sub(r26087, r26086, r26070, MPFR_RNDN); return mpfr_get_d(r26087, MPFR_RNDN); }