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