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