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