#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 r22502 = b; float r22503 = -r22502; float r22504 = r22502 * r22502; float r22505 = 3; float r22506 = a; float r22507 = r22505 * r22506; float r22508 = c; float r22509 = r22507 * r22508; float r22510 = r22504 - r22509; float r22511 = sqrt(r22510); float r22512 = r22503 + r22511; float r22513 = r22512 / r22507; return r22513; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r22514 = b; double r22515 = -r22514; double r22516 = r22514 * r22514; double r22517 = 3; double r22518 = a; double r22519 = r22517 * r22518; double r22520 = c; double r22521 = r22519 * r22520; double r22522 = r22516 - r22521; double r22523 = sqrt(r22522); double r22524 = r22515 + r22523; double r22525 = r22524 / r22519; return r22525; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r22526 = c; float r22527 = a; float r22528 = 3; float r22529 = r22527 * r22528; float r22530 = r22526 * r22529; float r22531 = b; float r22532 = -r22531; float r22533 = r22531 * r22531; float r22534 = r22528 * r22527; float r22535 = r22534 * r22526; float r22536 = r22533 - r22535; float r22537 = sqrt(r22536); float r22538 = r22532 - r22537; float r22539 = r22530 / r22538; float r22540 = 1; float r22541 = r22540 / r22534; float r22542 = r22539 * r22541; return r22542; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r22543 = c; double r22544 = a; double r22545 = 3; double r22546 = r22544 * r22545; double r22547 = r22543 * r22546; double r22548 = b; double r22549 = -r22548; double r22550 = r22548 * r22548; double r22551 = r22545 * r22544; double r22552 = r22551 * r22543; double r22553 = r22550 - r22552; double r22554 = sqrt(r22553); double r22555 = r22549 - r22554; double r22556 = r22547 / r22555; double r22557 = 1; double r22558 = r22557 / r22551; double r22559 = r22556 * r22558; return r22559; } 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 r22560, r22561, r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r22560); mpfr_init(r22561); mpfr_init(r22562); mpfr_init_set_str(r22563, "3", 10, MPFR_RNDN); mpfr_init(r22564); mpfr_init(r22565); mpfr_init(r22566); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init(r22570); mpfr_init(r22571); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r22560, b, MPFR_RNDN); mpfr_neg(r22561, r22560, MPFR_RNDN); mpfr_mul(r22562, r22560, r22560, MPFR_RNDN); ; mpfr_set_d(r22564, a, MPFR_RNDN); mpfr_mul(r22565, r22563, r22564, MPFR_RNDN); mpfr_set_d(r22566, c, MPFR_RNDN); mpfr_mul(r22567, r22565, r22566, MPFR_RNDN); mpfr_sub(r22568, r22562, r22567, MPFR_RNDN); mpfr_sqrt(r22569, r22568, MPFR_RNDN); mpfr_add(r22570, r22561, r22569, MPFR_RNDN); mpfr_div(r22571, r22570, r22565, MPFR_RNDN); return mpfr_get_d(r22571, MPFR_RNDN); } static mpfr_t r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r22572); mpfr_init(r22573); mpfr_init_set_str(r22574, "3", 10, MPFR_RNDN); mpfr_init(r22575); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); mpfr_init(r22580); mpfr_init(r22581); mpfr_init(r22582); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init_set_str(r22586, "1", 10, MPFR_RNDN); mpfr_init(r22587); mpfr_init(r22588); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r22572, c, MPFR_RNDN); mpfr_set_d(r22573, a, MPFR_RNDN); ; mpfr_mul(r22575, r22573, r22574, MPFR_RNDN); mpfr_mul(r22576, r22572, r22575, MPFR_RNDN); mpfr_set_d(r22577, b, MPFR_RNDN); mpfr_neg(r22578, r22577, MPFR_RNDN); mpfr_mul(r22579, r22577, r22577, MPFR_RNDN); mpfr_mul(r22580, r22574, r22573, MPFR_RNDN); mpfr_mul(r22581, r22580, r22572, MPFR_RNDN); mpfr_sub(r22582, r22579, r22581, MPFR_RNDN); mpfr_sqrt(r22583, r22582, MPFR_RNDN); mpfr_sub(r22584, r22578, r22583, MPFR_RNDN); mpfr_div(r22585, r22576, r22584, MPFR_RNDN); ; mpfr_div(r22587, r22586, r22580, MPFR_RNDN); mpfr_mul(r22588, r22585, r22587, MPFR_RNDN); return mpfr_get_d(r22588, MPFR_RNDN); } static mpfr_t r22589, r22590, r22591, r22592, r22593, r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r22589); mpfr_init(r22590); mpfr_init_set_str(r22591, "3", 10, MPFR_RNDN); mpfr_init(r22592); mpfr_init(r22593); mpfr_init(r22594); mpfr_init(r22595); mpfr_init(r22596); mpfr_init(r22597); mpfr_init(r22598); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init_set_str(r22603, "1", 10, MPFR_RNDN); mpfr_init(r22604); mpfr_init(r22605); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r22589, c, MPFR_RNDN); mpfr_set_d(r22590, a, MPFR_RNDN); ; mpfr_mul(r22592, r22590, r22591, MPFR_RNDN); mpfr_mul(r22593, r22589, r22592, MPFR_RNDN); mpfr_set_d(r22594, b, MPFR_RNDN); mpfr_neg(r22595, r22594, MPFR_RNDN); mpfr_mul(r22596, r22594, r22594, MPFR_RNDN); mpfr_mul(r22597, r22591, r22590, MPFR_RNDN); mpfr_mul(r22598, r22597, r22589, MPFR_RNDN); mpfr_sub(r22599, r22596, r22598, MPFR_RNDN); mpfr_sqrt(r22600, r22599, MPFR_RNDN); mpfr_sub(r22601, r22595, r22600, MPFR_RNDN); mpfr_div(r22602, r22593, r22601, MPFR_RNDN); ; mpfr_div(r22604, r22603, r22597, MPFR_RNDN); mpfr_mul(r22605, r22602, r22604, MPFR_RNDN); return mpfr_get_d(r22605, MPFR_RNDN); }