#include #include #include #include #include char *name = "Henrywood and Agarwal, Equation (13)"; double f_if(float c0, float w, float h, float D, float d, float M) { float r23679 = c0; float r23680 = 2; float r23681 = w; float r23682 = r23680 * r23681; float r23683 = r23679 / r23682; float r23684 = d; float r23685 = r23684 * r23684; float r23686 = r23679 * r23685; float r23687 = h; float r23688 = r23681 * r23687; float r23689 = D; float r23690 = r23689 * r23689; float r23691 = r23688 * r23690; float r23692 = r23686 / r23691; float r23693 = r23692 * r23692; float r23694 = M; float r23695 = r23694 * r23694; float r23696 = r23693 - r23695; float r23697 = sqrt(r23696); float r23698 = r23692 + r23697; float r23699 = r23683 * r23698; return r23699; } double f_id(double c0, double w, double h, double D, double d, double M) { double r23700 = c0; double r23701 = 2; double r23702 = w; double r23703 = r23701 * r23702; double r23704 = r23700 / r23703; double r23705 = d; double r23706 = r23705 * r23705; double r23707 = r23700 * r23706; double r23708 = h; double r23709 = r23702 * r23708; double r23710 = D; double r23711 = r23710 * r23710; double r23712 = r23709 * r23711; double r23713 = r23707 / r23712; double r23714 = r23713 * r23713; double r23715 = M; double r23716 = r23715 * r23715; double r23717 = r23714 - r23716; double r23718 = sqrt(r23717); double r23719 = r23713 + r23718; double r23720 = r23704 * r23719; return r23720; } double f_of(float c0, float w, float h, float D, float d, float __attribute__((unused)) M) { float r23721 = c0; float r23722 = 2; float r23723 = w; float r23724 = r23722 * r23723; float r23725 = r23721 / r23724; float r23726 = d; float r23727 = pow(r23726, r23722); float r23728 = r23727 * r23721; float r23729 = D; float r23730 = pow(r23729, r23722); float r23731 = h; float r23732 = r23731 * r23723; float r23733 = r23730 * r23732; float r23734 = r23728 / r23733; float r23735 = r23725 * r23734; float r23736 = -inf.0; bool r23737 = r23735 <= r23736; float r23738 = 0; float r23739 = 143688.22359878916; bool r23740 = r23735 <= r23739; float r23741 = r23722 * r23734; float r23742 = r23725 * r23741; float r23743 = r23740 ? r23742 : r23738; float r23744 = r23737 ? r23738 : r23743; return r23744; } double f_od(double c0, double w, double h, double D, double d, double __attribute__((unused)) M) { double r23745 = c0; double r23746 = 2; double r23747 = w; double r23748 = r23746 * r23747; double r23749 = r23745 / r23748; double r23750 = d; double r23751 = pow(r23750, r23746); double r23752 = r23751 * r23745; double r23753 = D; double r23754 = pow(r23753, r23746); double r23755 = h; double r23756 = r23755 * r23747; double r23757 = r23754 * r23756; double r23758 = r23752 / r23757; double r23759 = r23749 * r23758; double r23760 = -inf.0; bool r23761 = r23759 <= r23760; double r23762 = 0; double r23763 = 143688.22359878916; bool r23764 = r23759 <= r23763; double r23765 = r23746 * r23758; double r23766 = r23749 * r23765; double r23767 = r23764 ? r23766 : r23762; double r23768 = r23761 ? r23762 : r23767; return r23768; } 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 r23769, r23770, r23771, r23772, r23773, r23774, r23775, r23776, r23777, r23778, r23779, r23780, r23781, r23782, r23783, r23784, r23785, r23786, r23787, r23788, r23789; void setup_mpfr_f_im() { mpfr_set_default_prec(7312); mpfr_init(r23769); mpfr_init_set_str(r23770, "2", 10, MPFR_RNDN); 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); mpfr_init(r23781); mpfr_init(r23782); mpfr_init(r23783); mpfr_init(r23784); mpfr_init(r23785); mpfr_init(r23786); mpfr_init(r23787); mpfr_init(r23788); mpfr_init(r23789); } double f_im(double c0, double w, double h, double D, double d, double M) { mpfr_set_d(r23769, c0, MPFR_RNDN); ; mpfr_set_d(r23771, w, MPFR_RNDN); mpfr_mul(r23772, r23770, r23771, MPFR_RNDN); mpfr_div(r23773, r23769, r23772, MPFR_RNDN); mpfr_set_d(r23774, d, MPFR_RNDN); mpfr_mul(r23775, r23774, r23774, MPFR_RNDN); mpfr_mul(r23776, r23769, r23775, MPFR_RNDN); mpfr_set_d(r23777, h, MPFR_RNDN); mpfr_mul(r23778, r23771, r23777, MPFR_RNDN); mpfr_set_d(r23779, D, MPFR_RNDN); mpfr_mul(r23780, r23779, r23779, MPFR_RNDN); mpfr_mul(r23781, r23778, r23780, MPFR_RNDN); mpfr_div(r23782, r23776, r23781, MPFR_RNDN); mpfr_mul(r23783, r23782, r23782, MPFR_RNDN); mpfr_set_d(r23784, M, MPFR_RNDN); mpfr_mul(r23785, r23784, r23784, MPFR_RNDN); mpfr_sub(r23786, r23783, r23785, MPFR_RNDN); mpfr_sqrt(r23787, r23786, MPFR_RNDN); mpfr_add(r23788, r23782, r23787, MPFR_RNDN); mpfr_mul(r23789, r23773, r23788, MPFR_RNDN); return mpfr_get_d(r23789, MPFR_RNDN); } static mpfr_t r23790, r23791, r23792, r23793, r23794, r23795, r23796, r23797, r23798, r23799, r23800, r23801, r23802, r23803, r23804, r23805, r23806, r23807, r23808, r23809, r23810, r23811, r23812, r23813; void setup_mpfr_f_fm() { mpfr_set_default_prec(7312); mpfr_init(r23790); mpfr_init_set_str(r23791, "2", 10, MPFR_RNDN); mpfr_init(r23792); mpfr_init(r23793); mpfr_init(r23794); mpfr_init(r23795); mpfr_init(r23796); mpfr_init(r23797); mpfr_init(r23798); mpfr_init(r23799); mpfr_init(r23800); mpfr_init(r23801); mpfr_init(r23802); mpfr_init(r23803); mpfr_init(r23804); mpfr_init_set_str(r23805, "-inf.0", 10, MPFR_RNDN); mpfr_init(r23806); mpfr_init_set_str(r23807, "0", 10, MPFR_RNDN); mpfr_init_set_str(r23808, "143688.22359878916", 10, MPFR_RNDN); mpfr_init(r23809); mpfr_init(r23810); mpfr_init(r23811); mpfr_init(r23812); mpfr_init(r23813); } double f_fm(double c0, double w, double h, double D, double d, double __attribute__((unused)) M) { mpfr_set_d(r23790, c0, MPFR_RNDN); ; mpfr_set_d(r23792, w, MPFR_RNDN); mpfr_mul(r23793, r23791, r23792, MPFR_RNDN); mpfr_div(r23794, r23790, r23793, MPFR_RNDN); mpfr_set_d(r23795, d, MPFR_RNDN); mpfr_pow(r23796, r23795, r23791, MPFR_RNDN); mpfr_mul(r23797, r23796, r23790, MPFR_RNDN); mpfr_set_d(r23798, D, MPFR_RNDN); mpfr_pow(r23799, r23798, r23791, MPFR_RNDN); mpfr_set_d(r23800, h, MPFR_RNDN); mpfr_mul(r23801, r23800, r23792, MPFR_RNDN); mpfr_mul(r23802, r23799, r23801, MPFR_RNDN); mpfr_div(r23803, r23797, r23802, MPFR_RNDN); mpfr_mul(r23804, r23794, r23803, MPFR_RNDN); ; mpfr_set_si(r23806, mpfr_cmp(r23804, r23805) <= 0, MPFR_RNDN); ; ; mpfr_set_si(r23809, mpfr_cmp(r23804, r23808) <= 0, MPFR_RNDN); mpfr_mul(r23810, r23791, r23803, MPFR_RNDN); mpfr_mul(r23811, r23794, r23810, MPFR_RNDN); if (mpfr_get_si(r23809, MPFR_RNDN)) { mpfr_set(r23812, r23811, MPFR_RNDN); } else { mpfr_set(r23812, r23807, MPFR_RNDN); }; if (mpfr_get_si(r23806, MPFR_RNDN)) { mpfr_set(r23813, r23807, MPFR_RNDN); } else { mpfr_set(r23813, r23812, MPFR_RNDN); }; return mpfr_get_d(r23813, MPFR_RNDN); } static mpfr_t r23814, r23815, r23816, r23817, r23818, r23819, r23820, r23821, r23822, r23823, r23824, r23825, r23826, r23827, r23828, r23829, r23830, r23831, r23832, r23833, r23834, r23835, r23836, r23837; void setup_mpfr_f_dm() { mpfr_set_default_prec(7312); mpfr_init(r23814); mpfr_init_set_str(r23815, "2", 10, MPFR_RNDN); mpfr_init(r23816); 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); mpfr_init(r23828); mpfr_init_set_str(r23829, "-inf.0", 10, MPFR_RNDN); mpfr_init(r23830); mpfr_init_set_str(r23831, "0", 10, MPFR_RNDN); mpfr_init_set_str(r23832, "143688.22359878916", 10, MPFR_RNDN); mpfr_init(r23833); mpfr_init(r23834); mpfr_init(r23835); mpfr_init(r23836); mpfr_init(r23837); } double f_dm(double c0, double w, double h, double D, double d, double __attribute__((unused)) M) { mpfr_set_d(r23814, c0, MPFR_RNDN); ; mpfr_set_d(r23816, w, MPFR_RNDN); mpfr_mul(r23817, r23815, r23816, MPFR_RNDN); mpfr_div(r23818, r23814, r23817, MPFR_RNDN); mpfr_set_d(r23819, d, MPFR_RNDN); mpfr_pow(r23820, r23819, r23815, MPFR_RNDN); mpfr_mul(r23821, r23820, r23814, MPFR_RNDN); mpfr_set_d(r23822, D, MPFR_RNDN); mpfr_pow(r23823, r23822, r23815, MPFR_RNDN); mpfr_set_d(r23824, h, MPFR_RNDN); mpfr_mul(r23825, r23824, r23816, MPFR_RNDN); mpfr_mul(r23826, r23823, r23825, MPFR_RNDN); mpfr_div(r23827, r23821, r23826, MPFR_RNDN); mpfr_mul(r23828, r23818, r23827, MPFR_RNDN); ; mpfr_set_si(r23830, mpfr_cmp(r23828, r23829) <= 0, MPFR_RNDN); ; ; mpfr_set_si(r23833, mpfr_cmp(r23828, r23832) <= 0, MPFR_RNDN); mpfr_mul(r23834, r23815, r23827, MPFR_RNDN); mpfr_mul(r23835, r23818, r23834, MPFR_RNDN); if (mpfr_get_si(r23833, MPFR_RNDN)) { mpfr_set(r23836, r23835, MPFR_RNDN); } else { mpfr_set(r23836, r23831, MPFR_RNDN); }; if (mpfr_get_si(r23830, MPFR_RNDN)) { mpfr_set(r23837, r23831, MPFR_RNDN); } else { mpfr_set(r23837, r23836, MPFR_RNDN); }; return mpfr_get_d(r23837, MPFR_RNDN); }