#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 r23732 = b; float r23733 = -r23732; float r23734 = r23732 * r23732; float r23735 = 3; float r23736 = a; float r23737 = r23735 * r23736; float r23738 = c; float r23739 = r23737 * r23738; float r23740 = r23734 - r23739; float r23741 = sqrt(r23740); float r23742 = r23733 + r23741; float r23743 = r23742 / r23737; return r23743; } double f_id(double a, double b, double c, double __attribute__((unused)) d) { double r23744 = b; double r23745 = -r23744; double r23746 = r23744 * r23744; double r23747 = 3; double r23748 = a; double r23749 = r23747 * r23748; double r23750 = c; double r23751 = r23749 * r23750; double r23752 = r23746 - r23751; double r23753 = sqrt(r23752); double r23754 = r23745 + r23753; double r23755 = r23754 / r23749; return r23755; } double f_of(float a, float b, float c, float __attribute__((unused)) d) { float r23756 = b; float r23757 = -r23756; float r23758 = r23756 * r23756; float r23759 = 3; float r23760 = a; float r23761 = r23759 * r23760; float r23762 = c; float r23763 = r23761 * r23762; float r23764 = cbrt(r23763); float r23765 = r23764 * r23764; float r23766 = r23765 * r23764; float r23767 = r23758 - r23766; float r23768 = sqrt(r23767); float r23769 = r23757 + r23768; float r23770 = r23769 / r23761; return r23770; } double f_od(double a, double b, double c, double __attribute__((unused)) d) { double r23771 = b; double r23772 = -r23771; double r23773 = r23771 * r23771; double r23774 = 3; double r23775 = a; double r23776 = r23774 * r23775; double r23777 = c; double r23778 = r23776 * r23777; double r23779 = cbrt(r23778); double r23780 = r23779 * r23779; double r23781 = r23780 * r23779; double r23782 = r23773 - r23781; double r23783 = sqrt(r23782); double r23784 = r23772 + r23783; double r23785 = r23784 / r23776; return r23785; } 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 r23786, r23787, r23788, r23789, r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init_set_str(r23789, "3", 10, MPFR_RNDN); mpfr_init(r23790); mpfr_init(r23791); mpfr_init(r23792); mpfr_init(r23793); mpfr_init(r23794); mpfr_init(r23795); mpfr_init(r23796); mpfr_init(r23797); } double f_im(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23786, b, MPFR_RNDN); mpfr_neg(r23787, r23786, MPFR_RNDN); mpfr_mul(r23788, r23786, r23786, MPFR_RNDN); ; mpfr_set_d(r23790, a, MPFR_RNDN); mpfr_mul(r23791, r23789, r23790, MPFR_RNDN); mpfr_set_d(r23792, c, MPFR_RNDN); mpfr_mul(r23793, r23791, r23792, MPFR_RNDN); mpfr_sub(r23794, r23788, r23793, MPFR_RNDN); mpfr_sqrt(r23795, r23794, MPFR_RNDN); mpfr_add(r23796, r23787, r23795, MPFR_RNDN); mpfr_div(r23797, r23796, r23791, MPFR_RNDN); return mpfr_get_d(r23797, MPFR_RNDN); } static mpfr_t r23798, r23799, r23800, r23801, r23802, r23803, r23804, r23805, r23806, r23807, r23808, r23809, r23810, r23811, r23812; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r23798); mpfr_init(r23799); mpfr_init(r23800); mpfr_init_set_str(r23801, "3", 10, MPFR_RNDN); mpfr_init(r23802); mpfr_init(r23803); mpfr_init(r23804); mpfr_init(r23805); mpfr_init(r23806); mpfr_init(r23807); mpfr_init(r23808); mpfr_init(r23809); mpfr_init(r23810); mpfr_init(r23811); mpfr_init(r23812); } double f_fm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23798, b, MPFR_RNDN); mpfr_neg(r23799, r23798, MPFR_RNDN); mpfr_mul(r23800, r23798, r23798, MPFR_RNDN); ; mpfr_set_d(r23802, a, MPFR_RNDN); mpfr_mul(r23803, r23801, r23802, MPFR_RNDN); mpfr_set_d(r23804, c, MPFR_RNDN); mpfr_mul(r23805, r23803, r23804, MPFR_RNDN); mpfr_cbrt(r23806, r23805, MPFR_RNDN); mpfr_mul(r23807, r23806, r23806, MPFR_RNDN); mpfr_mul(r23808, r23807, r23806, MPFR_RNDN); mpfr_sub(r23809, r23800, r23808, MPFR_RNDN); mpfr_sqrt(r23810, r23809, MPFR_RNDN); mpfr_add(r23811, r23799, r23810, MPFR_RNDN); mpfr_div(r23812, r23811, r23803, MPFR_RNDN); return mpfr_get_d(r23812, MPFR_RNDN); } static mpfr_t r23813, r23814, r23815, r23816, r23817, r23818, r23819, r23820, r23821, r23822, r23823, r23824, r23825, r23826, r23827; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r23813); mpfr_init(r23814); mpfr_init(r23815); mpfr_init_set_str(r23816, "3", 10, MPFR_RNDN); mpfr_init(r23817); mpfr_init(r23818); mpfr_init(r23819); mpfr_init(r23820); mpfr_init(r23821); mpfr_init(r23822); mpfr_init(r23823); mpfr_init(r23824); mpfr_init(r23825); mpfr_init(r23826); mpfr_init(r23827); } double f_dm(double a, double b, double c, double __attribute__((unused)) d) { mpfr_set_d(r23813, b, MPFR_RNDN); mpfr_neg(r23814, r23813, MPFR_RNDN); mpfr_mul(r23815, r23813, r23813, MPFR_RNDN); ; mpfr_set_d(r23817, a, MPFR_RNDN); mpfr_mul(r23818, r23816, r23817, MPFR_RNDN); mpfr_set_d(r23819, c, MPFR_RNDN); mpfr_mul(r23820, r23818, r23819, MPFR_RNDN); mpfr_cbrt(r23821, r23820, MPFR_RNDN); mpfr_mul(r23822, r23821, r23821, MPFR_RNDN); mpfr_mul(r23823, r23822, r23821, MPFR_RNDN); mpfr_sub(r23824, r23815, r23823, MPFR_RNDN); mpfr_sqrt(r23825, r23824, MPFR_RNDN); mpfr_add(r23826, r23814, r23825, MPFR_RNDN); mpfr_div(r23827, r23826, r23818, MPFR_RNDN); return mpfr_get_d(r23827, MPFR_RNDN); }