#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 r25999 = b; float r26000 = -r25999; float r26001 = r25999 * r25999; float r26002 = 3; float r26003 = a; float r26004 = r26002 * r26003; float r26005 = c; float r26006 = r26004 * r26005; float r26007 = r26001 - r26006; float r26008 = sqrt(r26007); float r26009 = r26000 + r26008; float r26010 = r26009 / r26004; return r26010; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26011 = b; double r26012 = -r26011; double r26013 = r26011 * r26011; double r26014 = 3; double r26015 = a; double r26016 = r26014 * r26015; double r26017 = c; double r26018 = r26016 * r26017; double r26019 = r26013 - r26018; double r26020 = sqrt(r26019); double r26021 = r26012 + r26020; double r26022 = r26021 / r26016; return r26022; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26023 = c; float r26024 = -r26023; float r26025 = a; float r26026 = r26024 * r26025; float r26027 = 1; float r26028 = 3; float r26029 = r26028 * r26025; float r26030 = b; float r26031 = r26030 * r26030; float r26032 = fma(r26029, r26024, r26031); float r26033 = sqrt(r26032); float r26034 = r26033 + r26030; float r26035 = r26027 / r26034; float r26036 = r26035 / r26025; float r26037 = r26026 * r26036; return r26037; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26038 = c; double r26039 = -r26038; double r26040 = a; double r26041 = r26039 * r26040; double r26042 = 1; double r26043 = 3; double r26044 = r26043 * r26040; double r26045 = b; double r26046 = r26045 * r26045; double r26047 = fma(r26044, r26039, r26046); double r26048 = sqrt(r26047); double r26049 = r26048 + r26045; double r26050 = r26042 / r26049; double r26051 = r26050 / r26040; double r26052 = r26041 * r26051; return r26052; } 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 r26053, r26054, r26055, r26056, r26057, r26058, r26059, r26060, r26061, r26062, r26063, r26064; void setup_mpfr_f_im() { mpfr_set_default_prec(848); mpfr_init(r26053); mpfr_init(r26054); mpfr_init(r26055); mpfr_init_set_str(r26056, "3", 10, MPFR_RNDN); mpfr_init(r26057); mpfr_init(r26058); mpfr_init(r26059); mpfr_init(r26060); mpfr_init(r26061); mpfr_init(r26062); mpfr_init(r26063); mpfr_init(r26064); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26053, b, MPFR_RNDN); mpfr_neg(r26054, r26053, MPFR_RNDN); mpfr_mul(r26055, r26053, r26053, MPFR_RNDN); ; mpfr_set_d(r26057, a, MPFR_RNDN); mpfr_mul(r26058, r26056, r26057, MPFR_RNDN); mpfr_set_d(r26059, c, MPFR_RNDN); mpfr_mul(r26060, r26058, r26059, MPFR_RNDN); mpfr_sub(r26061, r26055, r26060, MPFR_RNDN); mpfr_sqrt(r26062, r26061, MPFR_RNDN); mpfr_add(r26063, r26054, r26062, MPFR_RNDN); mpfr_div(r26064, r26063, r26058, MPFR_RNDN); return mpfr_get_d(r26064, MPFR_RNDN); } static mpfr_t r26065, r26066, r26067, r26068, r26069, r26070, r26071, r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079; void setup_mpfr_f_fm() { mpfr_set_default_prec(848); mpfr_init(r26065); mpfr_init(r26066); mpfr_init(r26067); mpfr_init(r26068); mpfr_init_set_str(r26069, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26070, "3", 10, MPFR_RNDN); mpfr_init(r26071); mpfr_init(r26072); mpfr_init(r26073); mpfr_init(r26074); mpfr_init(r26075); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26065, c, MPFR_RNDN); mpfr_neg(r26066, r26065, MPFR_RNDN); mpfr_set_d(r26067, a, MPFR_RNDN); mpfr_mul(r26068, r26066, r26067, MPFR_RNDN); ; ; mpfr_mul(r26071, r26070, r26067, MPFR_RNDN); mpfr_set_d(r26072, b, MPFR_RNDN); mpfr_mul(r26073, r26072, r26072, MPFR_RNDN); mpfr_fma(r26074, r26071, r26066, r26073, MPFR_RNDN); mpfr_sqrt(r26075, r26074, MPFR_RNDN); mpfr_add(r26076, r26075, r26072, MPFR_RNDN); mpfr_div(r26077, r26069, r26076, MPFR_RNDN); mpfr_div(r26078, r26077, r26067, MPFR_RNDN); mpfr_mul(r26079, r26068, r26078, MPFR_RNDN); return mpfr_get_d(r26079, MPFR_RNDN); } static mpfr_t r26080, r26081, r26082, r26083, r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094; void setup_mpfr_f_dm() { mpfr_set_default_prec(848); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); mpfr_init_set_str(r26084, "1", 10, MPFR_RNDN); mpfr_init_set_str(r26085, "3", 10, MPFR_RNDN); mpfr_init(r26086); mpfr_init(r26087); mpfr_init(r26088); mpfr_init(r26089); mpfr_init(r26090); mpfr_init(r26091); mpfr_init(r26092); mpfr_init(r26093); mpfr_init(r26094); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26080, c, MPFR_RNDN); mpfr_neg(r26081, r26080, MPFR_RNDN); mpfr_set_d(r26082, a, MPFR_RNDN); mpfr_mul(r26083, r26081, r26082, MPFR_RNDN); ; ; mpfr_mul(r26086, r26085, r26082, MPFR_RNDN); mpfr_set_d(r26087, b, MPFR_RNDN); mpfr_mul(r26088, r26087, r26087, MPFR_RNDN); mpfr_fma(r26089, r26086, r26081, r26088, MPFR_RNDN); mpfr_sqrt(r26090, r26089, MPFR_RNDN); mpfr_add(r26091, r26090, r26087, MPFR_RNDN); mpfr_div(r26092, r26084, r26091, MPFR_RNDN); mpfr_div(r26093, r26092, r26082, MPFR_RNDN); mpfr_mul(r26094, r26083, r26093, MPFR_RNDN); return mpfr_get_d(r26094, MPFR_RNDN); }