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