#include #include #include #include #include char *name = "Cubic critical, wide range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r25950 = b; float r25951 = -r25950; float r25952 = r25950 * r25950; float r25953 = 3; float r25954 = a; float r25955 = r25953 * r25954; float r25956 = c; float r25957 = r25955 * r25956; float r25958 = r25952 - r25957; float r25959 = sqrt(r25958); float r25960 = r25951 + r25959; float r25961 = r25960 / r25955; return r25961; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r25962 = b; double r25963 = -r25962; double r25964 = r25962 * r25962; double r25965 = 3; double r25966 = a; double r25967 = r25965 * r25966; double r25968 = c; double r25969 = r25967 * r25968; double r25970 = r25964 - r25969; double r25971 = sqrt(r25970); double r25972 = r25963 + r25971; double r25973 = r25972 / r25967; return r25973; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r25974 = a; float r25975 = -r25974; float r25976 = c; float r25977 = r25975 * r25976; float r25978 = 1; float r25979 = b; float r25980 = r25979 * r25979; float r25981 = r25974 * r25976; float r25982 = 3; float r25983 = r25981 * r25982; float r25984 = r25980 - r25983; float r25985 = sqrt(r25984); float r25986 = r25985 + r25979; float r25987 = r25978 / r25986; float r25988 = r25987 / r25974; float r25989 = r25977 * r25988; return r25989; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r25990 = a; double r25991 = -r25990; double r25992 = c; double r25993 = r25991 * r25992; double r25994 = 1; double r25995 = b; double r25996 = r25995 * r25995; double r25997 = r25990 * r25992; double r25998 = 3; double r25999 = r25997 * r25998; double r26000 = r25996 - r25999; double r26001 = sqrt(r26000); double r26002 = r26001 + r25995; double r26003 = r25994 / r26002; double r26004 = r26003 / r25990; double r26005 = r25993 * r26004; return r26005; } 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 r26006, r26007, r26008, r26009, r26010, r26011, r26012, r26013, r26014, r26015, r26016, r26017; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r26006); mpfr_init(r26007); mpfr_init(r26008); mpfr_init_set_str(r26009, "3", 10, MPFR_RNDN); mpfr_init(r26010); mpfr_init(r26011); mpfr_init(r26012); mpfr_init(r26013); mpfr_init(r26014); mpfr_init(r26015); mpfr_init(r26016); mpfr_init(r26017); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26006, b, MPFR_RNDN); mpfr_neg(r26007, r26006, MPFR_RNDN); mpfr_mul(r26008, r26006, r26006, MPFR_RNDN); ; mpfr_set_d(r26010, a, MPFR_RNDN); mpfr_mul(r26011, r26009, r26010, MPFR_RNDN); mpfr_set_d(r26012, c, MPFR_RNDN); mpfr_mul(r26013, r26011, r26012, MPFR_RNDN); mpfr_sub(r26014, r26008, r26013, MPFR_RNDN); mpfr_sqrt(r26015, r26014, MPFR_RNDN); mpfr_add(r26016, r26007, r26015, MPFR_RNDN); mpfr_div(r26017, r26016, r26011, MPFR_RNDN); return mpfr_get_d(r26017, MPFR_RNDN); } static mpfr_t r26018, r26019, r26020, r26021, r26022, r26023, r26024, r26025, r26026, r26027, r26028, r26029, r26030, r26031, r26032, r26033; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r26018); mpfr_init(r26019); mpfr_init(r26020); mpfr_init(r26021); mpfr_init_set_str(r26022, "1", 10, MPFR_RNDN); mpfr_init(r26023); mpfr_init(r26024); mpfr_init(r26025); mpfr_init_set_str(r26026, "3", 10, MPFR_RNDN); mpfr_init(r26027); mpfr_init(r26028); mpfr_init(r26029); mpfr_init(r26030); mpfr_init(r26031); mpfr_init(r26032); mpfr_init(r26033); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26018, a, MPFR_RNDN); mpfr_neg(r26019, r26018, MPFR_RNDN); mpfr_set_d(r26020, c, MPFR_RNDN); mpfr_mul(r26021, r26019, r26020, MPFR_RNDN); ; mpfr_set_d(r26023, b, MPFR_RNDN); mpfr_mul(r26024, r26023, r26023, MPFR_RNDN); mpfr_mul(r26025, r26018, r26020, MPFR_RNDN); ; mpfr_mul(r26027, r26025, r26026, MPFR_RNDN); mpfr_sub(r26028, r26024, r26027, MPFR_RNDN); mpfr_sqrt(r26029, r26028, MPFR_RNDN); mpfr_add(r26030, r26029, r26023, MPFR_RNDN); mpfr_div(r26031, r26022, r26030, MPFR_RNDN); mpfr_div(r26032, r26031, r26018, MPFR_RNDN); mpfr_mul(r26033, r26021, r26032, MPFR_RNDN); return mpfr_get_d(r26033, MPFR_RNDN); } static mpfr_t r26034, r26035, r26036, r26037, r26038, r26039, r26040, r26041, r26042, r26043, r26044, r26045, r26046, r26047, r26048, r26049; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r26034); mpfr_init(r26035); mpfr_init(r26036); mpfr_init(r26037); mpfr_init_set_str(r26038, "1", 10, MPFR_RNDN); mpfr_init(r26039); mpfr_init(r26040); mpfr_init(r26041); mpfr_init_set_str(r26042, "3", 10, MPFR_RNDN); mpfr_init(r26043); mpfr_init(r26044); mpfr_init(r26045); mpfr_init(r26046); mpfr_init(r26047); mpfr_init(r26048); mpfr_init(r26049); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26034, a, MPFR_RNDN); mpfr_neg(r26035, r26034, MPFR_RNDN); mpfr_set_d(r26036, c, MPFR_RNDN); mpfr_mul(r26037, r26035, r26036, MPFR_RNDN); ; mpfr_set_d(r26039, b, MPFR_RNDN); mpfr_mul(r26040, r26039, r26039, MPFR_RNDN); mpfr_mul(r26041, r26034, r26036, MPFR_RNDN); ; mpfr_mul(r26043, r26041, r26042, MPFR_RNDN); mpfr_sub(r26044, r26040, r26043, MPFR_RNDN); mpfr_sqrt(r26045, r26044, MPFR_RNDN); mpfr_add(r26046, r26045, r26039, MPFR_RNDN); mpfr_div(r26047, r26038, r26046, MPFR_RNDN); mpfr_div(r26048, r26047, r26034, MPFR_RNDN); mpfr_mul(r26049, r26037, r26048, MPFR_RNDN); return mpfr_get_d(r26049, MPFR_RNDN); }