#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 r23697 = b; float r23698 = -r23697; float r23699 = r23697 * r23697; float r23700 = 3; float r23701 = a; float r23702 = r23700 * r23701; float r23703 = c; float r23704 = r23702 * r23703; float r23705 = r23699 - r23704; float r23706 = sqrt(r23705); float r23707 = r23698 + r23706; float r23708 = r23707 / r23702; return r23708; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23709 = b; double r23710 = -r23709; double r23711 = r23709 * r23709; double r23712 = 3; double r23713 = a; double r23714 = r23712 * r23713; double r23715 = c; double r23716 = r23714 * r23715; double r23717 = r23711 - r23716; double r23718 = sqrt(r23717); double r23719 = r23710 + r23718; double r23720 = r23719 / r23714; return r23720; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23721 = c; float r23722 = -r23721; float r23723 = a; float r23724 = r23722 * r23723; float r23725 = 1; float r23726 = 3; float r23727 = r23726 * r23723; float r23728 = b; float r23729 = r23728 * r23728; float r23730 = fma(r23727, r23722, r23729); float r23731 = sqrt(r23730); float r23732 = r23731 + r23728; float r23733 = r23725 / r23732; float r23734 = r23733 / r23723; float r23735 = r23724 * r23734; return r23735; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23736 = c; double r23737 = -r23736; double r23738 = a; double r23739 = r23737 * r23738; double r23740 = 1; double r23741 = 3; double r23742 = r23741 * r23738; double r23743 = b; double r23744 = r23743 * r23743; double r23745 = fma(r23742, r23737, r23744); double r23746 = sqrt(r23745); double r23747 = r23746 + r23743; double r23748 = r23740 / r23747; double r23749 = r23748 / r23738; double r23750 = r23739 * r23749; return r23750; } 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 r23751, r23752, r23753, r23754, r23755, r23756, r23757, r23758, r23759, r23760, r23761, r23762; void setup_mpfr_f_im() { mpfr_set_default_prec(656); mpfr_init(r23751); mpfr_init(r23752); mpfr_init(r23753); mpfr_init_set_str(r23754, "3", 10, MPFR_RNDN); mpfr_init(r23755); mpfr_init(r23756); mpfr_init(r23757); mpfr_init(r23758); mpfr_init(r23759); mpfr_init(r23760); mpfr_init(r23761); mpfr_init(r23762); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23751, b, MPFR_RNDN); mpfr_neg(r23752, r23751, MPFR_RNDN); mpfr_mul(r23753, r23751, r23751, MPFR_RNDN); ; mpfr_set_d(r23755, a, MPFR_RNDN); mpfr_mul(r23756, r23754, r23755, MPFR_RNDN); mpfr_set_d(r23757, c, MPFR_RNDN); mpfr_mul(r23758, r23756, r23757, MPFR_RNDN); mpfr_sub(r23759, r23753, r23758, MPFR_RNDN); mpfr_sqrt(r23760, r23759, MPFR_RNDN); mpfr_add(r23761, r23752, r23760, MPFR_RNDN); mpfr_div(r23762, r23761, r23756, MPFR_RNDN); return mpfr_get_d(r23762, MPFR_RNDN); } static mpfr_t r23763, r23764, r23765, r23766, r23767, r23768, r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777; void setup_mpfr_f_fm() { mpfr_set_default_prec(656); mpfr_init(r23763); mpfr_init(r23764); mpfr_init(r23765); mpfr_init(r23766); mpfr_init_set_str(r23767, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23768, "3", 10, MPFR_RNDN); mpfr_init(r23769); mpfr_init(r23770); mpfr_init(r23771); mpfr_init(r23772); mpfr_init(r23773); mpfr_init(r23774); mpfr_init(r23775); mpfr_init(r23776); mpfr_init(r23777); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23763, c, MPFR_RNDN); mpfr_neg(r23764, r23763, MPFR_RNDN); mpfr_set_d(r23765, a, MPFR_RNDN); mpfr_mul(r23766, r23764, r23765, MPFR_RNDN); ; ; mpfr_mul(r23769, r23768, r23765, MPFR_RNDN); mpfr_set_d(r23770, b, MPFR_RNDN); mpfr_mul(r23771, r23770, r23770, MPFR_RNDN); mpfr_fma(r23772, r23769, r23764, r23771, MPFR_RNDN); mpfr_sqrt(r23773, r23772, MPFR_RNDN); mpfr_add(r23774, r23773, r23770, MPFR_RNDN); mpfr_div(r23775, r23767, r23774, MPFR_RNDN); mpfr_div(r23776, r23775, r23765, MPFR_RNDN); mpfr_mul(r23777, r23766, r23776, MPFR_RNDN); return mpfr_get_d(r23777, MPFR_RNDN); } static mpfr_t r23778, r23779, r23780, r23781, r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789, r23790, r23791, r23792; void setup_mpfr_f_dm() { mpfr_set_default_prec(656); mpfr_init(r23778); mpfr_init(r23779); mpfr_init(r23780); mpfr_init(r23781); mpfr_init_set_str(r23782, "1", 10, MPFR_RNDN); mpfr_init_set_str(r23783, "3", 10, MPFR_RNDN); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); mpfr_init(r23790); mpfr_init(r23791); mpfr_init(r23792); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23778, c, MPFR_RNDN); mpfr_neg(r23779, r23778, MPFR_RNDN); mpfr_set_d(r23780, a, MPFR_RNDN); mpfr_mul(r23781, r23779, r23780, MPFR_RNDN); ; ; mpfr_mul(r23784, r23783, r23780, MPFR_RNDN); mpfr_set_d(r23785, b, MPFR_RNDN); mpfr_mul(r23786, r23785, r23785, MPFR_RNDN); mpfr_fma(r23787, r23784, r23779, r23786, MPFR_RNDN); mpfr_sqrt(r23788, r23787, MPFR_RNDN); mpfr_add(r23789, r23788, r23785, MPFR_RNDN); mpfr_div(r23790, r23782, r23789, MPFR_RNDN); mpfr_div(r23791, r23790, r23780, MPFR_RNDN); mpfr_mul(r23792, r23781, r23791, MPFR_RNDN); return mpfr_get_d(r23792, MPFR_RNDN); }