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