#include #include #include #include #include char *name = "Cubic critical, narrow range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r26016 = b; float r26017 = -r26016; float r26018 = r26016 * r26016; float r26019 = 3; float r26020 = a; float r26021 = r26019 * r26020; float r26022 = c; float r26023 = r26021 * r26022; float r26024 = r26018 - r26023; float r26025 = sqrt(r26024); float r26026 = r26017 + r26025; float r26027 = r26026 / r26021; return r26027; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r26028 = b; double r26029 = -r26028; double r26030 = r26028 * r26028; double r26031 = 3; double r26032 = a; double r26033 = r26031 * r26032; double r26034 = c; double r26035 = r26033 * r26034; double r26036 = r26030 - r26035; double r26037 = sqrt(r26036); double r26038 = r26029 + r26037; double r26039 = r26038 / r26033; return r26039; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r26040 = a; float r26041 = 3; float r26042 = r26040 / r26041; float r26043 = c; float r26044 = r26042 * r26043; float r26045 = -r26041; float r26046 = r26041 * r26040; float r26047 = -r26043; float r26048 = b; float r26049 = r26048 * r26048; float r26050 = fma(r26046, r26047, r26049); float r26051 = sqrt(r26050); float r26052 = r26051 + r26048; float r26053 = r26045 / r26052; float r26054 = r26053 / r26040; float r26055 = r26044 * r26054; return r26055; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r26056 = a; double r26057 = 3; double r26058 = r26056 / r26057; double r26059 = c; double r26060 = r26058 * r26059; double r26061 = -r26057; double r26062 = r26057 * r26056; double r26063 = -r26059; double r26064 = b; double r26065 = r26064 * r26064; double r26066 = fma(r26062, r26063, r26065); double r26067 = sqrt(r26066); double r26068 = r26067 + r26064; double r26069 = r26061 / r26068; double r26070 = r26069 / r26056; double r26071 = r26060 * r26070; return r26071; } 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 r26072, r26073, r26074, r26075, r26076, r26077, r26078, r26079, r26080, r26081, r26082, r26083; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r26072); mpfr_init(r26073); mpfr_init(r26074); mpfr_init_set_str(r26075, "3", 10, MPFR_RNDN); mpfr_init(r26076); mpfr_init(r26077); mpfr_init(r26078); mpfr_init(r26079); mpfr_init(r26080); mpfr_init(r26081); mpfr_init(r26082); mpfr_init(r26083); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26072, b, MPFR_RNDN); mpfr_neg(r26073, r26072, MPFR_RNDN); mpfr_mul(r26074, r26072, r26072, MPFR_RNDN); ; mpfr_set_d(r26076, a, MPFR_RNDN); mpfr_mul(r26077, r26075, r26076, MPFR_RNDN); mpfr_set_d(r26078, c, MPFR_RNDN); mpfr_mul(r26079, r26077, r26078, MPFR_RNDN); mpfr_sub(r26080, r26074, r26079, MPFR_RNDN); mpfr_sqrt(r26081, r26080, MPFR_RNDN); mpfr_add(r26082, r26073, r26081, MPFR_RNDN); mpfr_div(r26083, r26082, r26077, MPFR_RNDN); return mpfr_get_d(r26083, MPFR_RNDN); } static mpfr_t r26084, r26085, r26086, r26087, r26088, r26089, r26090, r26091, r26092, r26093, r26094, r26095, r26096, r26097, r26098, r26099; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r26084); 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); mpfr_init(r26095); mpfr_init(r26096); mpfr_init(r26097); mpfr_init(r26098); mpfr_init(r26099); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26084, a, MPFR_RNDN); ; mpfr_div(r26086, r26084, r26085, MPFR_RNDN); mpfr_set_d(r26087, c, MPFR_RNDN); mpfr_mul(r26088, r26086, r26087, MPFR_RNDN); mpfr_neg(r26089, r26085, MPFR_RNDN); mpfr_mul(r26090, r26085, r26084, MPFR_RNDN); mpfr_neg(r26091, r26087, MPFR_RNDN); mpfr_set_d(r26092, b, MPFR_RNDN); mpfr_mul(r26093, r26092, r26092, MPFR_RNDN); mpfr_fma(r26094, r26090, r26091, r26093, MPFR_RNDN); mpfr_sqrt(r26095, r26094, MPFR_RNDN); mpfr_add(r26096, r26095, r26092, MPFR_RNDN); mpfr_div(r26097, r26089, r26096, MPFR_RNDN); mpfr_div(r26098, r26097, r26084, MPFR_RNDN); mpfr_mul(r26099, r26088, r26098, MPFR_RNDN); return mpfr_get_d(r26099, MPFR_RNDN); } static mpfr_t r26100, r26101, r26102, r26103, r26104, r26105, r26106, r26107, r26108, r26109, r26110, r26111, r26112, r26113, r26114, r26115; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r26100); mpfr_init_set_str(r26101, "3", 10, MPFR_RNDN); mpfr_init(r26102); mpfr_init(r26103); mpfr_init(r26104); mpfr_init(r26105); mpfr_init(r26106); mpfr_init(r26107); mpfr_init(r26108); mpfr_init(r26109); mpfr_init(r26110); mpfr_init(r26111); mpfr_init(r26112); mpfr_init(r26113); mpfr_init(r26114); mpfr_init(r26115); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r26100, a, MPFR_RNDN); ; mpfr_div(r26102, r26100, r26101, MPFR_RNDN); mpfr_set_d(r26103, c, MPFR_RNDN); mpfr_mul(r26104, r26102, r26103, MPFR_RNDN); mpfr_neg(r26105, r26101, MPFR_RNDN); mpfr_mul(r26106, r26101, r26100, MPFR_RNDN); mpfr_neg(r26107, r26103, MPFR_RNDN); mpfr_set_d(r26108, b, MPFR_RNDN); mpfr_mul(r26109, r26108, r26108, MPFR_RNDN); mpfr_fma(r26110, r26106, r26107, r26109, MPFR_RNDN); mpfr_sqrt(r26111, r26110, MPFR_RNDN); mpfr_add(r26112, r26111, r26108, MPFR_RNDN); mpfr_div(r26113, r26105, r26112, MPFR_RNDN); mpfr_div(r26114, r26113, r26100, MPFR_RNDN); mpfr_mul(r26115, r26104, r26114, MPFR_RNDN); return mpfr_get_d(r26115, MPFR_RNDN); }