#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r9757 = i; float r9758 = alpha; float r9759 = beta; float r9760 = r9758 + r9759; float r9761 = r9760 + r9757; float r9762 = r9757 * r9761; float r9763 = r9759 * r9758; float r9764 = r9763 + r9762; float r9765 = r9762 * r9764; float r9766 = 2.0f; float r9767 = r9766 * r9757; float r9768 = r9760 + r9767; float r9769 = r9768 * r9768; float r9770 = r9765 / r9769; float r9771 = 1.0f; float r9772 = r9769 - r9771; float r9773 = r9770 / r9772; return r9773; } double f_id(double alpha, double beta, double i) { double r9774 = i; double r9775 = alpha; double r9776 = beta; double r9777 = r9775 + r9776; double r9778 = r9777 + r9774; double r9779 = r9774 * r9778; double r9780 = r9776 * r9775; double r9781 = r9780 + r9779; double r9782 = r9779 * r9781; double r9783 = 2.0; double r9784 = r9783 * r9774; double r9785 = r9777 + r9784; double r9786 = r9785 * r9785; double r9787 = r9782 / r9786; double r9788 = 1.0; double r9789 = r9786 - r9788; double r9790 = r9787 / r9789; return r9790; } double f_of(float alpha, float beta, float i) { float r9791 = i; float r9792 = 2.0f; float r9793 = r9792 * r9791; float r9794 = beta; float r9795 = alpha; float r9796 = r9794 + r9795; float r9797 = r9793 + r9796; float r9798 = r9791 / r9797; float r9799 = r9795 + r9791; float r9800 = sqrt(r9797); float r9801 = r9799 / r9800; float r9802 = r9794 + r9791; float r9803 = r9802 / r9800; float r9804 = r9801 * r9803; float r9805 = r9798 * r9804; float r9806 = r9794 + r9799; float r9807 = r9797 * r9797; float r9808 = 1.0f; float r9809 = r9807 - r9808; float r9810 = r9806 / r9809; float r9811 = r9805 * r9810; return r9811; } double f_od(double alpha, double beta, double i) { double r9812 = i; double r9813 = 2.0; double r9814 = r9813 * r9812; double r9815 = beta; double r9816 = alpha; double r9817 = r9815 + r9816; double r9818 = r9814 + r9817; double r9819 = r9812 / r9818; double r9820 = r9816 + r9812; double r9821 = sqrt(r9818); double r9822 = r9820 / r9821; double r9823 = r9815 + r9812; double r9824 = r9823 / r9821; double r9825 = r9822 * r9824; double r9826 = r9819 * r9825; double r9827 = r9815 + r9820; double r9828 = r9818 * r9818; double r9829 = 1.0; double r9830 = r9828 - r9829; double r9831 = r9827 / r9830; double r9832 = r9826 * r9831; return r9832; } 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 r9833, r9834, r9835, r9836, r9837, r9838, r9839, r9840, r9841, r9842, r9843, r9844, r9845, r9846, r9847, r9848, r9849; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r9833); mpfr_init(r9834); mpfr_init(r9835); mpfr_init(r9836); mpfr_init(r9837); mpfr_init(r9838); mpfr_init(r9839); mpfr_init(r9840); mpfr_init(r9841); mpfr_init_set_str(r9842, "2", 10, MPFR_RNDN); mpfr_init(r9843); mpfr_init(r9844); mpfr_init(r9845); mpfr_init(r9846); mpfr_init_set_str(r9847, "1.0", 10, MPFR_RNDN); mpfr_init(r9848); mpfr_init(r9849); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r9833, i, MPFR_RNDN); mpfr_set_d(r9834, alpha, MPFR_RNDN); mpfr_set_d(r9835, beta, MPFR_RNDN); mpfr_add(r9836, r9834, r9835, MPFR_RNDN); mpfr_add(r9837, r9836, r9833, MPFR_RNDN); mpfr_mul(r9838, r9833, r9837, MPFR_RNDN); mpfr_mul(r9839, r9835, r9834, MPFR_RNDN); mpfr_add(r9840, r9839, r9838, MPFR_RNDN); mpfr_mul(r9841, r9838, r9840, MPFR_RNDN); ; mpfr_mul(r9843, r9842, r9833, MPFR_RNDN); mpfr_add(r9844, r9836, r9843, MPFR_RNDN); mpfr_mul(r9845, r9844, r9844, MPFR_RNDN); mpfr_div(r9846, r9841, r9845, MPFR_RNDN); ; mpfr_sub(r9848, r9845, r9847, MPFR_RNDN); mpfr_div(r9849, r9846, r9848, MPFR_RNDN); return mpfr_get_d(r9849, MPFR_RNDN); } static mpfr_t r9850, r9851, r9852, r9853, r9854, r9855, r9856, r9857, r9858, r9859, r9860, r9861, r9862, r9863, r9864, r9865, r9866, r9867, r9868, r9869, r9870; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r9850); mpfr_init_set_str(r9851, "2", 10, MPFR_RNDN); mpfr_init(r9852); mpfr_init(r9853); mpfr_init(r9854); mpfr_init(r9855); mpfr_init(r9856); mpfr_init(r9857); mpfr_init(r9858); mpfr_init(r9859); mpfr_init(r9860); mpfr_init(r9861); mpfr_init(r9862); mpfr_init(r9863); mpfr_init(r9864); mpfr_init(r9865); mpfr_init(r9866); mpfr_init_set_str(r9867, "1.0", 10, MPFR_RNDN); mpfr_init(r9868); mpfr_init(r9869); mpfr_init(r9870); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r9850, i, MPFR_RNDN); ; mpfr_mul(r9852, r9851, r9850, MPFR_RNDN); mpfr_set_d(r9853, beta, MPFR_RNDN); mpfr_set_d(r9854, alpha, MPFR_RNDN); mpfr_add(r9855, r9853, r9854, MPFR_RNDN); mpfr_add(r9856, r9852, r9855, MPFR_RNDN); mpfr_div(r9857, r9850, r9856, MPFR_RNDN); mpfr_add(r9858, r9854, r9850, MPFR_RNDN); mpfr_sqrt(r9859, r9856, MPFR_RNDN); mpfr_div(r9860, r9858, r9859, MPFR_RNDN); mpfr_add(r9861, r9853, r9850, MPFR_RNDN); mpfr_div(r9862, r9861, r9859, MPFR_RNDN); mpfr_mul(r9863, r9860, r9862, MPFR_RNDN); mpfr_mul(r9864, r9857, r9863, MPFR_RNDN); mpfr_add(r9865, r9853, r9858, MPFR_RNDN); mpfr_mul(r9866, r9856, r9856, MPFR_RNDN); ; mpfr_sub(r9868, r9866, r9867, MPFR_RNDN); mpfr_div(r9869, r9865, r9868, MPFR_RNDN); mpfr_mul(r9870, r9864, r9869, MPFR_RNDN); return mpfr_get_d(r9870, MPFR_RNDN); } static mpfr_t r9871, r9872, r9873, r9874, r9875, r9876, r9877, r9878, r9879, r9880, r9881, r9882, r9883, r9884, r9885, r9886, r9887, r9888, r9889, r9890, r9891; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r9871); mpfr_init_set_str(r9872, "2", 10, MPFR_RNDN); mpfr_init(r9873); mpfr_init(r9874); mpfr_init(r9875); mpfr_init(r9876); mpfr_init(r9877); mpfr_init(r9878); mpfr_init(r9879); mpfr_init(r9880); mpfr_init(r9881); mpfr_init(r9882); mpfr_init(r9883); mpfr_init(r9884); mpfr_init(r9885); mpfr_init(r9886); mpfr_init(r9887); mpfr_init_set_str(r9888, "1.0", 10, MPFR_RNDN); mpfr_init(r9889); mpfr_init(r9890); mpfr_init(r9891); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r9871, i, MPFR_RNDN); ; mpfr_mul(r9873, r9872, r9871, MPFR_RNDN); mpfr_set_d(r9874, beta, MPFR_RNDN); mpfr_set_d(r9875, alpha, MPFR_RNDN); mpfr_add(r9876, r9874, r9875, MPFR_RNDN); mpfr_add(r9877, r9873, r9876, MPFR_RNDN); mpfr_div(r9878, r9871, r9877, MPFR_RNDN); mpfr_add(r9879, r9875, r9871, MPFR_RNDN); mpfr_sqrt(r9880, r9877, MPFR_RNDN); mpfr_div(r9881, r9879, r9880, MPFR_RNDN); mpfr_add(r9882, r9874, r9871, MPFR_RNDN); mpfr_div(r9883, r9882, r9880, MPFR_RNDN); mpfr_mul(r9884, r9881, r9883, MPFR_RNDN); mpfr_mul(r9885, r9878, r9884, MPFR_RNDN); mpfr_add(r9886, r9874, r9879, MPFR_RNDN); mpfr_mul(r9887, r9877, r9877, MPFR_RNDN); ; mpfr_sub(r9889, r9887, r9888, MPFR_RNDN); mpfr_div(r9890, r9886, r9889, MPFR_RNDN); mpfr_mul(r9891, r9885, r9890, MPFR_RNDN); return mpfr_get_d(r9891, MPFR_RNDN); }