#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 r23567 = b; float r23568 = -r23567; float r23569 = r23567 * r23567; float r23570 = 3; float r23571 = a; float r23572 = r23570 * r23571; float r23573 = c; float r23574 = r23572 * r23573; float r23575 = r23569 - r23574; float r23576 = sqrt(r23575); float r23577 = r23568 + r23576; float r23578 = r23577 / r23572; return r23578; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23579 = b; double r23580 = -r23579; double r23581 = r23579 * r23579; double r23582 = 3; double r23583 = a; double r23584 = r23582 * r23583; double r23585 = c; double r23586 = r23584 * r23585; double r23587 = r23581 - r23586; double r23588 = sqrt(r23587); double r23589 = r23580 + r23588; double r23590 = r23589 / r23584; return r23590; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23591 = c; float r23592 = a; float r23593 = r23591 * r23592; float r23594 = 3; float r23595 = -r23594; float r23596 = r23593 * r23595; float r23597 = r23594 * r23592; float r23598 = -r23591; float r23599 = b; float r23600 = r23599 * r23599; float r23601 = fma(r23597, r23598, r23600); float r23602 = sqrt(r23601); float r23603 = r23602 + r23599; float r23604 = r23596 / r23603; float r23605 = r23604 / r23597; return r23605; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23606 = c; double r23607 = a; double r23608 = r23606 * r23607; double r23609 = 3; double r23610 = -r23609; double r23611 = r23608 * r23610; double r23612 = r23609 * r23607; double r23613 = -r23606; double r23614 = b; double r23615 = r23614 * r23614; double r23616 = fma(r23612, r23613, r23615); double r23617 = sqrt(r23616); double r23618 = r23617 + r23614; double r23619 = r23611 / r23618; double r23620 = r23619 / r23612; return r23620; } 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 r23621, r23622, r23623, r23624, r23625, r23626, r23627, r23628, r23629, r23630, r23631, r23632; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); mpfr_init_set_str(r23624, "3", 10, MPFR_RNDN); mpfr_init(r23625); mpfr_init(r23626); mpfr_init(r23627); mpfr_init(r23628); mpfr_init(r23629); mpfr_init(r23630); mpfr_init(r23631); mpfr_init(r23632); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23621, b, MPFR_RNDN); mpfr_neg(r23622, r23621, MPFR_RNDN); mpfr_mul(r23623, r23621, r23621, MPFR_RNDN); ; mpfr_set_d(r23625, a, MPFR_RNDN); mpfr_mul(r23626, r23624, r23625, MPFR_RNDN); mpfr_set_d(r23627, c, MPFR_RNDN); mpfr_mul(r23628, r23626, r23627, MPFR_RNDN); mpfr_sub(r23629, r23623, r23628, MPFR_RNDN); mpfr_sqrt(r23630, r23629, MPFR_RNDN); mpfr_add(r23631, r23622, r23630, MPFR_RNDN); mpfr_div(r23632, r23631, r23626, MPFR_RNDN); return mpfr_get_d(r23632, MPFR_RNDN); } static mpfr_t r23633, r23634, r23635, r23636, r23637, r23638, r23639, r23640, r23641, r23642, r23643, r23644, r23645, r23646, r23647; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r23633); mpfr_init(r23634); mpfr_init(r23635); mpfr_init_set_str(r23636, "3", 10, MPFR_RNDN); mpfr_init(r23637); mpfr_init(r23638); mpfr_init(r23639); mpfr_init(r23640); mpfr_init(r23641); mpfr_init(r23642); mpfr_init(r23643); mpfr_init(r23644); mpfr_init(r23645); mpfr_init(r23646); mpfr_init(r23647); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23633, c, MPFR_RNDN); mpfr_set_d(r23634, a, MPFR_RNDN); mpfr_mul(r23635, r23633, r23634, MPFR_RNDN); ; mpfr_neg(r23637, r23636, MPFR_RNDN); mpfr_mul(r23638, r23635, r23637, MPFR_RNDN); mpfr_mul(r23639, r23636, r23634, MPFR_RNDN); mpfr_neg(r23640, r23633, MPFR_RNDN); mpfr_set_d(r23641, b, MPFR_RNDN); mpfr_mul(r23642, r23641, r23641, MPFR_RNDN); mpfr_fma(r23643, r23639, r23640, r23642, MPFR_RNDN); mpfr_sqrt(r23644, r23643, MPFR_RNDN); mpfr_add(r23645, r23644, r23641, MPFR_RNDN); mpfr_div(r23646, r23638, r23645, MPFR_RNDN); mpfr_div(r23647, r23646, r23639, MPFR_RNDN); return mpfr_get_d(r23647, MPFR_RNDN); } static mpfr_t r23648, r23649, r23650, r23651, r23652, r23653, r23654, r23655, r23656, r23657, r23658, r23659, r23660, r23661, r23662; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r23648); mpfr_init(r23649); mpfr_init(r23650); mpfr_init_set_str(r23651, "3", 10, MPFR_RNDN); mpfr_init(r23652); mpfr_init(r23653); mpfr_init(r23654); mpfr_init(r23655); mpfr_init(r23656); mpfr_init(r23657); mpfr_init(r23658); mpfr_init(r23659); mpfr_init(r23660); mpfr_init(r23661); mpfr_init(r23662); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23648, c, MPFR_RNDN); mpfr_set_d(r23649, a, MPFR_RNDN); mpfr_mul(r23650, r23648, r23649, MPFR_RNDN); ; mpfr_neg(r23652, r23651, MPFR_RNDN); mpfr_mul(r23653, r23650, r23652, MPFR_RNDN); mpfr_mul(r23654, r23651, r23649, MPFR_RNDN); mpfr_neg(r23655, r23648, MPFR_RNDN); mpfr_set_d(r23656, b, MPFR_RNDN); mpfr_mul(r23657, r23656, r23656, MPFR_RNDN); mpfr_fma(r23658, r23654, r23655, r23657, MPFR_RNDN); mpfr_sqrt(r23659, r23658, MPFR_RNDN); mpfr_add(r23660, r23659, r23656, MPFR_RNDN); mpfr_div(r23661, r23653, r23660, MPFR_RNDN); mpfr_div(r23662, r23661, r23654, MPFR_RNDN); return mpfr_get_d(r23662, MPFR_RNDN); }