#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 r22478 = b; float r22479 = -r22478; float r22480 = r22478 * r22478; float r22481 = 3; float r22482 = a; float r22483 = r22481 * r22482; float r22484 = c; float r22485 = r22483 * r22484; float r22486 = r22480 - r22485; float r22487 = sqrt(r22486); float r22488 = r22479 + r22487; float r22489 = r22488 / r22483; return r22489; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r22490 = b; double r22491 = -r22490; double r22492 = r22490 * r22490; double r22493 = 3; double r22494 = a; double r22495 = r22493 * r22494; double r22496 = c; double r22497 = r22495 * r22496; double r22498 = r22492 - r22497; double r22499 = sqrt(r22498); double r22500 = r22491 + r22499; double r22501 = r22500 / r22495; return r22501; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r22502 = c; float r22503 = 1; float r22504 = r22502 / r22503; float r22505 = a; float r22506 = 3; float r22507 = r22505 * r22506; float r22508 = b; float r22509 = -r22508; float r22510 = r22508 * r22508; float r22511 = r22506 * r22505; float r22512 = r22511 * r22502; float r22513 = r22510 - r22512; float r22514 = sqrt(r22513); float r22515 = r22509 - r22514; float r22516 = r22507 / r22515; float r22517 = r22504 * r22516; float r22518 = r22517 / r22511; return r22518; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r22519 = c; double r22520 = 1; double r22521 = r22519 / r22520; double r22522 = a; double r22523 = 3; double r22524 = r22522 * r22523; double r22525 = b; double r22526 = -r22525; double r22527 = r22525 * r22525; double r22528 = r22523 * r22522; double r22529 = r22528 * r22519; double r22530 = r22527 - r22529; double r22531 = sqrt(r22530); double r22532 = r22526 - r22531; double r22533 = r22524 / r22532; double r22534 = r22521 * r22533; double r22535 = r22534 / r22528; return r22535; } 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 r22536, r22537, r22538, r22539, r22540, r22541, r22542, r22543, r22544, r22545, r22546, r22547; void setup_mpfr_f_im() { mpfr_set_default_prec(912); mpfr_init(r22536); mpfr_init(r22537); mpfr_init(r22538); mpfr_init_set_str(r22539, "3", 10, MPFR_RNDN); mpfr_init(r22540); mpfr_init(r22541); mpfr_init(r22542); mpfr_init(r22543); mpfr_init(r22544); mpfr_init(r22545); mpfr_init(r22546); mpfr_init(r22547); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r22536, b, MPFR_RNDN); mpfr_neg(r22537, r22536, MPFR_RNDN); mpfr_mul(r22538, r22536, r22536, MPFR_RNDN); ; mpfr_set_d(r22540, a, MPFR_RNDN); mpfr_mul(r22541, r22539, r22540, MPFR_RNDN); mpfr_set_d(r22542, c, MPFR_RNDN); mpfr_mul(r22543, r22541, r22542, MPFR_RNDN); mpfr_sub(r22544, r22538, r22543, MPFR_RNDN); mpfr_sqrt(r22545, r22544, MPFR_RNDN); mpfr_add(r22546, r22537, r22545, MPFR_RNDN); mpfr_div(r22547, r22546, r22541, MPFR_RNDN); return mpfr_get_d(r22547, MPFR_RNDN); } static mpfr_t r22548, r22549, r22550, r22551, r22552, r22553, r22554, r22555, r22556, r22557, r22558, r22559, r22560, r22561, r22562, r22563, r22564; void setup_mpfr_f_fm() { mpfr_set_default_prec(912); mpfr_init(r22548); mpfr_init_set_str(r22549, "1", 10, MPFR_RNDN); mpfr_init(r22550); mpfr_init(r22551); mpfr_init_set_str(r22552, "3", 10, MPFR_RNDN); mpfr_init(r22553); mpfr_init(r22554); mpfr_init(r22555); mpfr_init(r22556); mpfr_init(r22557); mpfr_init(r22558); mpfr_init(r22559); mpfr_init(r22560); mpfr_init(r22561); mpfr_init(r22562); mpfr_init(r22563); mpfr_init(r22564); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r22548, c, MPFR_RNDN); ; mpfr_div(r22550, r22548, r22549, MPFR_RNDN); mpfr_set_d(r22551, a, MPFR_RNDN); ; mpfr_mul(r22553, r22551, r22552, MPFR_RNDN); mpfr_set_d(r22554, b, MPFR_RNDN); mpfr_neg(r22555, r22554, MPFR_RNDN); mpfr_mul(r22556, r22554, r22554, MPFR_RNDN); mpfr_mul(r22557, r22552, r22551, MPFR_RNDN); mpfr_mul(r22558, r22557, r22548, MPFR_RNDN); mpfr_sub(r22559, r22556, r22558, MPFR_RNDN); mpfr_sqrt(r22560, r22559, MPFR_RNDN); mpfr_sub(r22561, r22555, r22560, MPFR_RNDN); mpfr_div(r22562, r22553, r22561, MPFR_RNDN); mpfr_mul(r22563, r22550, r22562, MPFR_RNDN); mpfr_div(r22564, r22563, r22557, MPFR_RNDN); return mpfr_get_d(r22564, MPFR_RNDN); } static mpfr_t r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577, r22578, r22579, r22580, r22581; void setup_mpfr_f_dm() { mpfr_set_default_prec(912); mpfr_init(r22565); mpfr_init_set_str(r22566, "1", 10, MPFR_RNDN); mpfr_init(r22567); mpfr_init(r22568); mpfr_init_set_str(r22569, "3", 10, MPFR_RNDN); mpfr_init(r22570); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); mpfr_init(r22576); mpfr_init(r22577); mpfr_init(r22578); mpfr_init(r22579); mpfr_init(r22580); mpfr_init(r22581); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r22565, c, MPFR_RNDN); ; mpfr_div(r22567, r22565, r22566, MPFR_RNDN); mpfr_set_d(r22568, a, MPFR_RNDN); ; mpfr_mul(r22570, r22568, r22569, MPFR_RNDN); mpfr_set_d(r22571, b, MPFR_RNDN); mpfr_neg(r22572, r22571, MPFR_RNDN); mpfr_mul(r22573, r22571, r22571, MPFR_RNDN); mpfr_mul(r22574, r22569, r22568, MPFR_RNDN); mpfr_mul(r22575, r22574, r22565, MPFR_RNDN); mpfr_sub(r22576, r22573, r22575, MPFR_RNDN); mpfr_sqrt(r22577, r22576, MPFR_RNDN); mpfr_sub(r22578, r22572, r22577, MPFR_RNDN); mpfr_div(r22579, r22570, r22578, MPFR_RNDN); mpfr_mul(r22580, r22567, r22579, MPFR_RNDN); mpfr_div(r22581, r22580, r22574, MPFR_RNDN); return mpfr_get_d(r22581, MPFR_RNDN); }