#include #include #include #include #include char *name = "Cubic critical, narrow range"; double f_if(float a, float b, float c, float __attribute__((unused)) d) { float r9551 = b; float r9552 = -r9551; float r9553 = r9551 * r9551; float r9554 = 3.0f; float r9555 = a; float r9556 = r9554 * r9555; float r9557 = c; float r9558 = r9556 * r9557; float r9559 = r9553 - r9558; float r9560 = sqrt(r9559); float r9561 = r9552 + r9560; float r9562 = r9561 / r9556; return r9562; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r9563 = b; double r9564 = -r9563; double r9565 = r9563 * r9563; double r9566 = 3.0; double r9567 = a; double r9568 = r9566 * r9567; double r9569 = c; double r9570 = r9568 * r9569; double r9571 = r9565 - r9570; double r9572 = sqrt(r9571); double r9573 = r9564 + r9572; double r9574 = r9573 / r9568; return r9574; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r9575 = c; float r9576 = -r9575; float r9577 = b; float r9578 = r9577 * r9577; float r9579 = 3.0f; float r9580 = pow(r9578, r9579); float r9581 = a; float r9582 = r9575 * r9581; float r9583 = r9582 * r9579; float r9584 = pow(r9583, r9579); float r9585 = r9580 - r9584; float r9586 = r9578 * r9583; float r9587 = r9583 * r9583; float r9588 = r9586 + r9587; float r9589 = r9578 * r9578; float r9590 = r9589 * r9589; float r9591 = r9590 * r9589; float r9592 = cbrt(r9591); float r9593 = r9588 + r9592; float r9594 = r9585 / r9593; float r9595 = sqrt(r9594); float r9596 = r9577 + r9595; float r9597 = r9576 / r9596; return r9597; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r9598 = c; double r9599 = -r9598; double r9600 = b; double r9601 = r9600 * r9600; double r9602 = 3.0; double r9603 = pow(r9601, r9602); double r9604 = a; double r9605 = r9598 * r9604; double r9606 = r9605 * r9602; double r9607 = pow(r9606, r9602); double r9608 = r9603 - r9607; double r9609 = r9601 * r9606; double r9610 = r9606 * r9606; double r9611 = r9609 + r9610; double r9612 = r9601 * r9601; double r9613 = r9612 * r9612; double r9614 = r9613 * r9612; double r9615 = cbrt(r9614); double r9616 = r9611 + r9615; double r9617 = r9608 / r9616; double r9618 = sqrt(r9617); double r9619 = r9600 + r9618; double r9620 = r9599 / r9619; return r9620; } 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 r9621, r9622, r9623, r9624, r9625, r9626, r9627, r9628, r9629, r9630, r9631, r9632; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r9621); mpfr_init(r9622); mpfr_init(r9623); mpfr_init_set_str(r9624, "3", 10, MPFR_RNDN); mpfr_init(r9625); mpfr_init(r9626); mpfr_init(r9627); mpfr_init(r9628); mpfr_init(r9629); mpfr_init(r9630); mpfr_init(r9631); mpfr_init(r9632); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r9621, b, MPFR_RNDN); mpfr_neg(r9622, r9621, MPFR_RNDN); mpfr_mul(r9623, r9621, r9621, MPFR_RNDN); ; mpfr_set_d(r9625, a, MPFR_RNDN); mpfr_mul(r9626, r9624, r9625, MPFR_RNDN); mpfr_set_d(r9627, c, MPFR_RNDN); mpfr_mul(r9628, r9626, r9627, MPFR_RNDN); mpfr_sub(r9629, r9623, r9628, MPFR_RNDN); mpfr_sqrt(r9630, r9629, MPFR_RNDN); mpfr_add(r9631, r9622, r9630, MPFR_RNDN); mpfr_div(r9632, r9631, r9626, MPFR_RNDN); return mpfr_get_d(r9632, MPFR_RNDN); } static mpfr_t r9633, r9634, r9635, r9636, r9637, r9638, r9639, r9640, r9641, r9642, r9643, r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r9633); mpfr_init(r9634); mpfr_init(r9635); mpfr_init(r9636); mpfr_init_set_str(r9637, "3", 10, MPFR_RNDN); mpfr_init(r9638); mpfr_init(r9639); mpfr_init(r9640); mpfr_init(r9641); mpfr_init(r9642); mpfr_init(r9643); mpfr_init(r9644); mpfr_init(r9645); mpfr_init(r9646); mpfr_init(r9647); mpfr_init(r9648); mpfr_init(r9649); mpfr_init(r9650); mpfr_init(r9651); mpfr_init(r9652); mpfr_init(r9653); mpfr_init(r9654); mpfr_init(r9655); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r9633, c, MPFR_RNDN); mpfr_neg(r9634, r9633, MPFR_RNDN); mpfr_set_d(r9635, b, MPFR_RNDN); mpfr_mul(r9636, r9635, r9635, MPFR_RNDN); ; mpfr_pow(r9638, r9636, r9637, MPFR_RNDN); mpfr_set_d(r9639, a, MPFR_RNDN); mpfr_mul(r9640, r9633, r9639, MPFR_RNDN); mpfr_mul(r9641, r9640, r9637, MPFR_RNDN); mpfr_pow(r9642, r9641, r9637, MPFR_RNDN); mpfr_sub(r9643, r9638, r9642, MPFR_RNDN); mpfr_mul(r9644, r9636, r9641, MPFR_RNDN); mpfr_mul(r9645, r9641, r9641, MPFR_RNDN); mpfr_add(r9646, r9644, r9645, MPFR_RNDN); mpfr_mul(r9647, r9636, r9636, MPFR_RNDN); mpfr_mul(r9648, r9647, r9647, MPFR_RNDN); mpfr_mul(r9649, r9648, r9647, MPFR_RNDN); mpfr_cbrt(r9650, r9649, MPFR_RNDN); mpfr_add(r9651, r9646, r9650, MPFR_RNDN); mpfr_div(r9652, r9643, r9651, MPFR_RNDN); mpfr_sqrt(r9653, r9652, MPFR_RNDN); mpfr_add(r9654, r9635, r9653, MPFR_RNDN); mpfr_div(r9655, r9634, r9654, MPFR_RNDN); return mpfr_get_d(r9655, MPFR_RNDN); } static mpfr_t r9656, r9657, r9658, r9659, r9660, r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669, r9670, r9671, r9672, r9673, r9674, r9675, r9676, r9677, r9678; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r9656); mpfr_init(r9657); mpfr_init(r9658); mpfr_init(r9659); mpfr_init_set_str(r9660, "3", 10, MPFR_RNDN); mpfr_init(r9661); mpfr_init(r9662); mpfr_init(r9663); mpfr_init(r9664); mpfr_init(r9665); mpfr_init(r9666); mpfr_init(r9667); mpfr_init(r9668); mpfr_init(r9669); mpfr_init(r9670); mpfr_init(r9671); mpfr_init(r9672); mpfr_init(r9673); mpfr_init(r9674); mpfr_init(r9675); mpfr_init(r9676); mpfr_init(r9677); mpfr_init(r9678); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r9656, c, MPFR_RNDN); mpfr_neg(r9657, r9656, MPFR_RNDN); mpfr_set_d(r9658, b, MPFR_RNDN); mpfr_mul(r9659, r9658, r9658, MPFR_RNDN); ; mpfr_pow(r9661, r9659, r9660, MPFR_RNDN); mpfr_set_d(r9662, a, MPFR_RNDN); mpfr_mul(r9663, r9656, r9662, MPFR_RNDN); mpfr_mul(r9664, r9663, r9660, MPFR_RNDN); mpfr_pow(r9665, r9664, r9660, MPFR_RNDN); mpfr_sub(r9666, r9661, r9665, MPFR_RNDN); mpfr_mul(r9667, r9659, r9664, MPFR_RNDN); mpfr_mul(r9668, r9664, r9664, MPFR_RNDN); mpfr_add(r9669, r9667, r9668, MPFR_RNDN); mpfr_mul(r9670, r9659, r9659, MPFR_RNDN); mpfr_mul(r9671, r9670, r9670, MPFR_RNDN); mpfr_mul(r9672, r9671, r9670, MPFR_RNDN); mpfr_cbrt(r9673, r9672, MPFR_RNDN); mpfr_add(r9674, r9669, r9673, MPFR_RNDN); mpfr_div(r9675, r9666, r9674, MPFR_RNDN); mpfr_sqrt(r9676, r9675, MPFR_RNDN); mpfr_add(r9677, r9658, r9676, MPFR_RNDN); mpfr_div(r9678, r9657, r9677, MPFR_RNDN); return mpfr_get_d(r9678, MPFR_RNDN); }