#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 r23525 = b; float r23526 = -r23525; float r23527 = r23525 * r23525; float r23528 = 3; float r23529 = a; float r23530 = r23528 * r23529; float r23531 = c; float r23532 = r23530 * r23531; float r23533 = r23527 - r23532; float r23534 = sqrt(r23533); float r23535 = r23526 + r23534; float r23536 = r23535 / r23530; return r23536; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23537 = b; double r23538 = -r23537; double r23539 = r23537 * r23537; double r23540 = 3; double r23541 = a; double r23542 = r23540 * r23541; double r23543 = c; double r23544 = r23542 * r23543; double r23545 = r23539 - r23544; double r23546 = sqrt(r23545); double r23547 = r23538 + r23546; double r23548 = r23547 / r23542; return r23548; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23549 = c; float r23550 = a; float r23551 = r23549 * r23550; float r23552 = 1; float r23553 = r23551 / r23552; float r23554 = 3; float r23555 = -r23554; float r23556 = r23554 * r23550; float r23557 = -r23549; float r23558 = b; float r23559 = r23558 * r23558; float r23560 = fma(r23556, r23557, r23559); float r23561 = sqrt(r23560); float r23562 = r23561 + r23558; float r23563 = r23555 / r23562; float r23564 = r23553 * r23563; float r23565 = r23564 / r23556; return r23565; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23566 = c; double r23567 = a; double r23568 = r23566 * r23567; double r23569 = 1; double r23570 = r23568 / r23569; double r23571 = 3; double r23572 = -r23571; double r23573 = r23571 * r23567; double r23574 = -r23566; double r23575 = b; double r23576 = r23575 * r23575; double r23577 = fma(r23573, r23574, r23576); double r23578 = sqrt(r23577); double r23579 = r23578 + r23575; double r23580 = r23572 / r23579; double r23581 = r23570 * r23580; double r23582 = r23581 / r23573; return r23582; } 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 r23583, r23584, r23585, r23586, r23587, r23588, r23589, r23590, r23591, r23592, r23593, r23594; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r23583); mpfr_init(r23584); mpfr_init(r23585); mpfr_init_set_str(r23586, "3", 10, MPFR_RNDN); mpfr_init(r23587); mpfr_init(r23588); mpfr_init(r23589); mpfr_init(r23590); mpfr_init(r23591); mpfr_init(r23592); mpfr_init(r23593); mpfr_init(r23594); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23583, b, MPFR_RNDN); mpfr_neg(r23584, r23583, MPFR_RNDN); mpfr_mul(r23585, r23583, r23583, MPFR_RNDN); ; mpfr_set_d(r23587, a, MPFR_RNDN); mpfr_mul(r23588, r23586, r23587, MPFR_RNDN); mpfr_set_d(r23589, c, MPFR_RNDN); mpfr_mul(r23590, r23588, r23589, MPFR_RNDN); mpfr_sub(r23591, r23585, r23590, MPFR_RNDN); mpfr_sqrt(r23592, r23591, MPFR_RNDN); mpfr_add(r23593, r23584, r23592, MPFR_RNDN); mpfr_div(r23594, r23593, r23588, MPFR_RNDN); return mpfr_get_d(r23594, MPFR_RNDN); } static mpfr_t r23595, r23596, r23597, r23598, r23599, r23600, r23601, r23602, r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r23595); mpfr_init(r23596); mpfr_init(r23597); mpfr_init_set_str(r23598, "1", 10, MPFR_RNDN); mpfr_init(r23599); mpfr_init_set_str(r23600, "3", 10, MPFR_RNDN); mpfr_init(r23601); mpfr_init(r23602); mpfr_init(r23603); mpfr_init(r23604); mpfr_init(r23605); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init(r23609); mpfr_init(r23610); mpfr_init(r23611); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23595, c, MPFR_RNDN); mpfr_set_d(r23596, a, MPFR_RNDN); mpfr_mul(r23597, r23595, r23596, MPFR_RNDN); ; mpfr_div(r23599, r23597, r23598, MPFR_RNDN); ; mpfr_neg(r23601, r23600, MPFR_RNDN); mpfr_mul(r23602, r23600, r23596, MPFR_RNDN); mpfr_neg(r23603, r23595, MPFR_RNDN); mpfr_set_d(r23604, b, MPFR_RNDN); mpfr_mul(r23605, r23604, r23604, MPFR_RNDN); mpfr_fma(r23606, r23602, r23603, r23605, MPFR_RNDN); mpfr_sqrt(r23607, r23606, MPFR_RNDN); mpfr_add(r23608, r23607, r23604, MPFR_RNDN); mpfr_div(r23609, r23601, r23608, MPFR_RNDN); mpfr_mul(r23610, r23599, r23609, MPFR_RNDN); mpfr_div(r23611, r23610, r23602, MPFR_RNDN); return mpfr_get_d(r23611, MPFR_RNDN); } static mpfr_t r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626, r23627, r23628; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r23612); mpfr_init(r23613); mpfr_init(r23614); mpfr_init_set_str(r23615, "1", 10, MPFR_RNDN); mpfr_init(r23616); mpfr_init_set_str(r23617, "3", 10, MPFR_RNDN); mpfr_init(r23618); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); mpfr_init(r23624); mpfr_init(r23625); mpfr_init(r23626); mpfr_init(r23627); mpfr_init(r23628); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23612, c, MPFR_RNDN); mpfr_set_d(r23613, a, MPFR_RNDN); mpfr_mul(r23614, r23612, r23613, MPFR_RNDN); ; mpfr_div(r23616, r23614, r23615, MPFR_RNDN); ; mpfr_neg(r23618, r23617, MPFR_RNDN); mpfr_mul(r23619, r23617, r23613, MPFR_RNDN); mpfr_neg(r23620, r23612, MPFR_RNDN); mpfr_set_d(r23621, b, MPFR_RNDN); mpfr_mul(r23622, r23621, r23621, MPFR_RNDN); mpfr_fma(r23623, r23619, r23620, r23622, MPFR_RNDN); mpfr_sqrt(r23624, r23623, MPFR_RNDN); mpfr_add(r23625, r23624, r23621, MPFR_RNDN); mpfr_div(r23626, r23618, r23625, MPFR_RNDN); mpfr_mul(r23627, r23616, r23626, MPFR_RNDN); mpfr_div(r23628, r23627, r23619, MPFR_RNDN); return mpfr_get_d(r23628, MPFR_RNDN); }