#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r22723 = i; float r22724 = alpha; float r22725 = beta; float r22726 = r22724 + r22725; float r22727 = r22726 + r22723; float r22728 = r22723 * r22727; float r22729 = r22725 * r22724; float r22730 = r22729 + r22728; float r22731 = r22728 * r22730; float r22732 = 2; float r22733 = r22732 * r22723; float r22734 = r22726 + r22733; float r22735 = r22734 * r22734; float r22736 = r22731 / r22735; float r22737 = 1.0; float r22738 = r22735 - r22737; float r22739 = r22736 / r22738; return r22739; } double f_id(double alpha, double beta, double i) { double r22740 = i; double r22741 = alpha; double r22742 = beta; double r22743 = r22741 + r22742; double r22744 = r22743 + r22740; double r22745 = r22740 * r22744; double r22746 = r22742 * r22741; double r22747 = r22746 + r22745; double r22748 = r22745 * r22747; double r22749 = 2; double r22750 = r22749 * r22740; double r22751 = r22743 + r22750; double r22752 = r22751 * r22751; double r22753 = r22748 / r22752; double r22754 = 1.0; double r22755 = r22752 - r22754; double r22756 = r22753 / r22755; return r22756; } double f_of(float alpha, float beta, float i) { float r22757 = alpha; float r22758 = 1.0859168924204577e+148; bool r22759 = r22757 <= r22758; float r22760 = i; float r22761 = beta; float r22762 = r22757 + r22761; float r22763 = r22762 + r22760; float r22764 = r22760 * r22763; float r22765 = 2; float r22766 = r22765 * r22760; float r22767 = r22762 + r22766; float r22768 = r22764 / r22767; float r22769 = r22761 * r22757; float r22770 = r22769 + r22764; float r22771 = r22770 / r22767; float r22772 = r22768 * r22771; float r22773 = r22767 * r22767; float r22774 = 1.0; float r22775 = r22773 - r22774; float r22776 = r22772 / r22775; float r22777 = 0; float r22778 = r22759 ? r22776 : r22777; return r22778; } double f_od(double alpha, double beta, double i) { double r22779 = alpha; double r22780 = 1.0859168924204577e+148; bool r22781 = r22779 <= r22780; double r22782 = i; double r22783 = beta; double r22784 = r22779 + r22783; double r22785 = r22784 + r22782; double r22786 = r22782 * r22785; double r22787 = 2; double r22788 = r22787 * r22782; double r22789 = r22784 + r22788; double r22790 = r22786 / r22789; double r22791 = r22783 * r22779; double r22792 = r22791 + r22786; double r22793 = r22792 / r22789; double r22794 = r22790 * r22793; double r22795 = r22789 * r22789; double r22796 = 1.0; double r22797 = r22795 - r22796; double r22798 = r22794 / r22797; double r22799 = 0; double r22800 = r22781 ? r22798 : r22799; return r22800; } 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 r22801, r22802, r22803, r22804, r22805, r22806, r22807, r22808, r22809, r22810, r22811, r22812, r22813, r22814, r22815, r22816, r22817; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22801); mpfr_init(r22802); mpfr_init(r22803); mpfr_init(r22804); mpfr_init(r22805); mpfr_init(r22806); mpfr_init(r22807); mpfr_init(r22808); mpfr_init(r22809); mpfr_init_set_str(r22810, "2", 10, MPFR_RNDN); mpfr_init(r22811); mpfr_init(r22812); mpfr_init(r22813); mpfr_init(r22814); mpfr_init_set_str(r22815, "1.0", 10, MPFR_RNDN); mpfr_init(r22816); mpfr_init(r22817); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r22801, i, MPFR_RNDN); mpfr_set_d(r22802, alpha, MPFR_RNDN); mpfr_set_d(r22803, beta, MPFR_RNDN); mpfr_add(r22804, r22802, r22803, MPFR_RNDN); mpfr_add(r22805, r22804, r22801, MPFR_RNDN); mpfr_mul(r22806, r22801, r22805, MPFR_RNDN); mpfr_mul(r22807, r22803, r22802, MPFR_RNDN); mpfr_add(r22808, r22807, r22806, MPFR_RNDN); mpfr_mul(r22809, r22806, r22808, MPFR_RNDN); ; mpfr_mul(r22811, r22810, r22801, MPFR_RNDN); mpfr_add(r22812, r22804, r22811, MPFR_RNDN); mpfr_mul(r22813, r22812, r22812, MPFR_RNDN); mpfr_div(r22814, r22809, r22813, MPFR_RNDN); ; mpfr_sub(r22816, r22813, r22815, MPFR_RNDN); mpfr_div(r22817, r22814, r22816, MPFR_RNDN); return mpfr_get_d(r22817, MPFR_RNDN); } static mpfr_t r22818, r22819, r22820, r22821, r22822, r22823, r22824, r22825, r22826, r22827, r22828, r22829, r22830, r22831, r22832, r22833, r22834, r22835, r22836, r22837, r22838, r22839; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22818); mpfr_init_set_str(r22819, "1.0859168924204577e+148", 10, MPFR_RNDN); mpfr_init(r22820); mpfr_init(r22821); mpfr_init(r22822); mpfr_init(r22823); mpfr_init(r22824); mpfr_init(r22825); mpfr_init_set_str(r22826, "2", 10, MPFR_RNDN); mpfr_init(r22827); mpfr_init(r22828); mpfr_init(r22829); mpfr_init(r22830); mpfr_init(r22831); mpfr_init(r22832); mpfr_init(r22833); mpfr_init(r22834); mpfr_init_set_str(r22835, "1.0", 10, MPFR_RNDN); mpfr_init(r22836); mpfr_init(r22837); mpfr_init_set_str(r22838, "0", 10, MPFR_RNDN); mpfr_init(r22839); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r22818, alpha, MPFR_RNDN); ; mpfr_set_si(r22820, mpfr_cmp(r22818, r22819) <= 0, MPFR_RNDN); mpfr_set_d(r22821, i, MPFR_RNDN); mpfr_set_d(r22822, beta, MPFR_RNDN); mpfr_add(r22823, r22818, r22822, MPFR_RNDN); mpfr_add(r22824, r22823, r22821, MPFR_RNDN); mpfr_mul(r22825, r22821, r22824, MPFR_RNDN); ; mpfr_mul(r22827, r22826, r22821, MPFR_RNDN); mpfr_add(r22828, r22823, r22827, MPFR_RNDN); mpfr_div(r22829, r22825, r22828, MPFR_RNDN); mpfr_mul(r22830, r22822, r22818, MPFR_RNDN); mpfr_add(r22831, r22830, r22825, MPFR_RNDN); mpfr_div(r22832, r22831, r22828, MPFR_RNDN); mpfr_mul(r22833, r22829, r22832, MPFR_RNDN); mpfr_mul(r22834, r22828, r22828, MPFR_RNDN); ; mpfr_sub(r22836, r22834, r22835, MPFR_RNDN); mpfr_div(r22837, r22833, r22836, MPFR_RNDN); ; if (mpfr_get_si(r22820, MPFR_RNDN)) { mpfr_set(r22839, r22837, MPFR_RNDN); } else { mpfr_set(r22839, r22838, MPFR_RNDN); }; return mpfr_get_d(r22839, MPFR_RNDN); } static mpfr_t r22840, r22841, r22842, r22843, r22844, r22845, r22846, r22847, r22848, r22849, r22850, r22851, r22852, r22853, r22854, r22855, r22856, r22857, r22858, r22859, r22860, r22861; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22840); mpfr_init_set_str(r22841, "1.0859168924204577e+148", 10, MPFR_RNDN); mpfr_init(r22842); mpfr_init(r22843); mpfr_init(r22844); mpfr_init(r22845); mpfr_init(r22846); mpfr_init(r22847); mpfr_init_set_str(r22848, "2", 10, MPFR_RNDN); mpfr_init(r22849); mpfr_init(r22850); mpfr_init(r22851); mpfr_init(r22852); mpfr_init(r22853); mpfr_init(r22854); mpfr_init(r22855); mpfr_init(r22856); mpfr_init_set_str(r22857, "1.0", 10, MPFR_RNDN); mpfr_init(r22858); mpfr_init(r22859); mpfr_init_set_str(r22860, "0", 10, MPFR_RNDN); mpfr_init(r22861); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r22840, alpha, MPFR_RNDN); ; mpfr_set_si(r22842, mpfr_cmp(r22840, r22841) <= 0, MPFR_RNDN); mpfr_set_d(r22843, i, MPFR_RNDN); mpfr_set_d(r22844, beta, MPFR_RNDN); mpfr_add(r22845, r22840, r22844, MPFR_RNDN); mpfr_add(r22846, r22845, r22843, MPFR_RNDN); mpfr_mul(r22847, r22843, r22846, MPFR_RNDN); ; mpfr_mul(r22849, r22848, r22843, MPFR_RNDN); mpfr_add(r22850, r22845, r22849, MPFR_RNDN); mpfr_div(r22851, r22847, r22850, MPFR_RNDN); mpfr_mul(r22852, r22844, r22840, MPFR_RNDN); mpfr_add(r22853, r22852, r22847, MPFR_RNDN); mpfr_div(r22854, r22853, r22850, MPFR_RNDN); mpfr_mul(r22855, r22851, r22854, MPFR_RNDN); mpfr_mul(r22856, r22850, r22850, MPFR_RNDN); ; mpfr_sub(r22858, r22856, r22857, MPFR_RNDN); mpfr_div(r22859, r22855, r22858, MPFR_RNDN); ; if (mpfr_get_si(r22842, MPFR_RNDN)) { mpfr_set(r22861, r22859, MPFR_RNDN); } else { mpfr_set(r22861, r22860, MPFR_RNDN); }; return mpfr_get_d(r22861, MPFR_RNDN); }