#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r25725 = i; float r25726 = alpha; float r25727 = beta; float r25728 = r25726 + r25727; float r25729 = r25728 + r25725; float r25730 = r25725 * r25729; float r25731 = r25727 * r25726; float r25732 = r25731 + r25730; float r25733 = r25730 * r25732; float r25734 = 2; float r25735 = r25734 * r25725; float r25736 = r25728 + r25735; float r25737 = r25736 * r25736; float r25738 = r25733 / r25737; float r25739 = 1.0; float r25740 = r25737 - r25739; float r25741 = r25738 / r25740; return r25741; } double f_id(double alpha, double beta, double i) { double r25742 = i; double r25743 = alpha; double r25744 = beta; double r25745 = r25743 + r25744; double r25746 = r25745 + r25742; double r25747 = r25742 * r25746; double r25748 = r25744 * r25743; double r25749 = r25748 + r25747; double r25750 = r25747 * r25749; double r25751 = 2; double r25752 = r25751 * r25742; double r25753 = r25745 + r25752; double r25754 = r25753 * r25753; double r25755 = r25750 / r25754; double r25756 = 1.0; double r25757 = r25754 - r25756; double r25758 = r25755 / r25757; return r25758; } double f_of(float alpha, float beta, float i) { float r25759 = alpha; float r25760 = i; float r25761 = r25759 + r25760; float r25762 = beta; float r25763 = r25761 + r25762; float r25764 = r25762 * r25759; float r25765 = fma(r25763, r25760, r25764); float r25766 = 2; float r25767 = fma(r25760, r25766, r25759); float r25768 = r25767 + r25762; float r25769 = 1.0; float r25770 = -r25769; float r25771 = fma(r25768, r25768, r25770); float r25772 = r25765 / r25771; float r25773 = r25763 / r25768; float r25774 = r25760 / r25768; float r25775 = r25773 * r25774; float r25776 = r25772 * r25775; float r25777 = cbrt(r25776); float r25778 = r25777 * r25777; float r25779 = r25778 * r25777; float r25780 = 0.0633766128278581; bool r25781 = r25779 <= r25780; float r25782 = expm1(r25773); float r25783 = log1p(r25782); float r25784 = r25783 * r25774; float r25785 = r25772 * r25784; float r25786 = 0; float r25787 = r25781 ? r25785 : r25786; return r25787; } double f_od(double alpha, double beta, double i) { double r25788 = alpha; double r25789 = i; double r25790 = r25788 + r25789; double r25791 = beta; double r25792 = r25790 + r25791; double r25793 = r25791 * r25788; double r25794 = fma(r25792, r25789, r25793); double r25795 = 2; double r25796 = fma(r25789, r25795, r25788); double r25797 = r25796 + r25791; double r25798 = 1.0; double r25799 = -r25798; double r25800 = fma(r25797, r25797, r25799); double r25801 = r25794 / r25800; double r25802 = r25792 / r25797; double r25803 = r25789 / r25797; double r25804 = r25802 * r25803; double r25805 = r25801 * r25804; double r25806 = cbrt(r25805); double r25807 = r25806 * r25806; double r25808 = r25807 * r25806; double r25809 = 0.0633766128278581; bool r25810 = r25808 <= r25809; double r25811 = expm1(r25802); double r25812 = log1p(r25811); double r25813 = r25812 * r25803; double r25814 = r25801 * r25813; double r25815 = 0; double r25816 = r25810 ? r25814 : r25815; return r25816; } 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 r25817, r25818, r25819, r25820, r25821, r25822, r25823, r25824, r25825, r25826, r25827, r25828, r25829, r25830, r25831, r25832, r25833; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25817); mpfr_init(r25818); mpfr_init(r25819); mpfr_init(r25820); mpfr_init(r25821); mpfr_init(r25822); mpfr_init(r25823); mpfr_init(r25824); mpfr_init(r25825); mpfr_init_set_str(r25826, "2", 10, MPFR_RNDN); mpfr_init(r25827); mpfr_init(r25828); mpfr_init(r25829); mpfr_init(r25830); mpfr_init_set_str(r25831, "1.0", 10, MPFR_RNDN); mpfr_init(r25832); mpfr_init(r25833); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r25817, i, MPFR_RNDN); mpfr_set_d(r25818, alpha, MPFR_RNDN); mpfr_set_d(r25819, beta, MPFR_RNDN); mpfr_add(r25820, r25818, r25819, MPFR_RNDN); mpfr_add(r25821, r25820, r25817, MPFR_RNDN); mpfr_mul(r25822, r25817, r25821, MPFR_RNDN); mpfr_mul(r25823, r25819, r25818, MPFR_RNDN); mpfr_add(r25824, r25823, r25822, MPFR_RNDN); mpfr_mul(r25825, r25822, r25824, MPFR_RNDN); ; mpfr_mul(r25827, r25826, r25817, MPFR_RNDN); mpfr_add(r25828, r25820, r25827, MPFR_RNDN); mpfr_mul(r25829, r25828, r25828, MPFR_RNDN); mpfr_div(r25830, r25825, r25829, MPFR_RNDN); ; mpfr_sub(r25832, r25829, r25831, MPFR_RNDN); mpfr_div(r25833, r25830, r25832, MPFR_RNDN); return mpfr_get_d(r25833, MPFR_RNDN); } static mpfr_t r25834, r25835, r25836, r25837, r25838, r25839, r25840, r25841, r25842, r25843, r25844, r25845, r25846, r25847, r25848, r25849, r25850, r25851, r25852, r25853, r25854, r25855, r25856, r25857, r25858, r25859, r25860, r25861, r25862; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25834); mpfr_init(r25835); mpfr_init(r25836); mpfr_init(r25837); mpfr_init(r25838); mpfr_init(r25839); mpfr_init(r25840); mpfr_init_set_str(r25841, "2", 10, MPFR_RNDN); mpfr_init(r25842); mpfr_init(r25843); mpfr_init_set_str(r25844, "1.0", 10, MPFR_RNDN); mpfr_init(r25845); mpfr_init(r25846); mpfr_init(r25847); mpfr_init(r25848); mpfr_init(r25849); mpfr_init(r25850); mpfr_init(r25851); mpfr_init(r25852); mpfr_init(r25853); mpfr_init(r25854); mpfr_init_set_str(r25855, "0.0633766128278581", 10, MPFR_RNDN); mpfr_init(r25856); mpfr_init(r25857); mpfr_init(r25858); mpfr_init(r25859); mpfr_init(r25860); mpfr_init_set_str(r25861, "0", 10, MPFR_RNDN); mpfr_init(r25862); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r25834, alpha, MPFR_RNDN); mpfr_set_d(r25835, i, MPFR_RNDN); mpfr_add(r25836, r25834, r25835, MPFR_RNDN); mpfr_set_d(r25837, beta, MPFR_RNDN); mpfr_add(r25838, r25836, r25837, MPFR_RNDN); mpfr_mul(r25839, r25837, r25834, MPFR_RNDN); mpfr_fma(r25840, r25838, r25835, r25839, MPFR_RNDN); ; mpfr_fma(r25842, r25835, r25841, r25834, MPFR_RNDN); mpfr_add(r25843, r25842, r25837, MPFR_RNDN); ; mpfr_neg(r25845, r25844, MPFR_RNDN); mpfr_fma(r25846, r25843, r25843, r25845, MPFR_RNDN); mpfr_div(r25847, r25840, r25846, MPFR_RNDN); mpfr_div(r25848, r25838, r25843, MPFR_RNDN); mpfr_div(r25849, r25835, r25843, MPFR_RNDN); mpfr_mul(r25850, r25848, r25849, MPFR_RNDN); mpfr_mul(r25851, r25847, r25850, MPFR_RNDN); mpfr_cbrt(r25852, r25851, MPFR_RNDN); mpfr_mul(r25853, r25852, r25852, MPFR_RNDN); mpfr_mul(r25854, r25853, r25852, MPFR_RNDN); ; mpfr_set_si(r25856, mpfr_cmp(r25854, r25855) <= 0, MPFR_RNDN); mpfr_expm1(r25857, r25848, MPFR_RNDN); mpfr_log1p(r25858, r25857, MPFR_RNDN); mpfr_mul(r25859, r25858, r25849, MPFR_RNDN); mpfr_mul(r25860, r25847, r25859, MPFR_RNDN); ; if (mpfr_get_si(r25856, MPFR_RNDN)) { mpfr_set(r25862, r25860, MPFR_RNDN); } else { mpfr_set(r25862, r25861, MPFR_RNDN); }; return mpfr_get_d(r25862, MPFR_RNDN); } static mpfr_t r25863, r25864, r25865, r25866, r25867, r25868, r25869, r25870, r25871, r25872, r25873, r25874, r25875, r25876, r25877, r25878, r25879, r25880, r25881, r25882, r25883, r25884, r25885, r25886, r25887, r25888, r25889, r25890, r25891; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25863); mpfr_init(r25864); mpfr_init(r25865); mpfr_init(r25866); mpfr_init(r25867); mpfr_init(r25868); mpfr_init(r25869); mpfr_init_set_str(r25870, "2", 10, MPFR_RNDN); mpfr_init(r25871); mpfr_init(r25872); mpfr_init_set_str(r25873, "1.0", 10, MPFR_RNDN); mpfr_init(r25874); mpfr_init(r25875); mpfr_init(r25876); mpfr_init(r25877); mpfr_init(r25878); mpfr_init(r25879); mpfr_init(r25880); mpfr_init(r25881); mpfr_init(r25882); mpfr_init(r25883); mpfr_init_set_str(r25884, "0.0633766128278581", 10, MPFR_RNDN); mpfr_init(r25885); mpfr_init(r25886); mpfr_init(r25887); mpfr_init(r25888); mpfr_init(r25889); mpfr_init_set_str(r25890, "0", 10, MPFR_RNDN); mpfr_init(r25891); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r25863, alpha, MPFR_RNDN); mpfr_set_d(r25864, i, MPFR_RNDN); mpfr_add(r25865, r25863, r25864, MPFR_RNDN); mpfr_set_d(r25866, beta, MPFR_RNDN); mpfr_add(r25867, r25865, r25866, MPFR_RNDN); mpfr_mul(r25868, r25866, r25863, MPFR_RNDN); mpfr_fma(r25869, r25867, r25864, r25868, MPFR_RNDN); ; mpfr_fma(r25871, r25864, r25870, r25863, MPFR_RNDN); mpfr_add(r25872, r25871, r25866, MPFR_RNDN); ; mpfr_neg(r25874, r25873, MPFR_RNDN); mpfr_fma(r25875, r25872, r25872, r25874, MPFR_RNDN); mpfr_div(r25876, r25869, r25875, MPFR_RNDN); mpfr_div(r25877, r25867, r25872, MPFR_RNDN); mpfr_div(r25878, r25864, r25872, MPFR_RNDN); mpfr_mul(r25879, r25877, r25878, MPFR_RNDN); mpfr_mul(r25880, r25876, r25879, MPFR_RNDN); mpfr_cbrt(r25881, r25880, MPFR_RNDN); mpfr_mul(r25882, r25881, r25881, MPFR_RNDN); mpfr_mul(r25883, r25882, r25881, MPFR_RNDN); ; mpfr_set_si(r25885, mpfr_cmp(r25883, r25884) <= 0, MPFR_RNDN); mpfr_expm1(r25886, r25877, MPFR_RNDN); mpfr_log1p(r25887, r25886, MPFR_RNDN); mpfr_mul(r25888, r25887, r25878, MPFR_RNDN); mpfr_mul(r25889, r25876, r25888, MPFR_RNDN); ; if (mpfr_get_si(r25885, MPFR_RNDN)) { mpfr_set(r25891, r25889, MPFR_RNDN); } else { mpfr_set(r25891, r25890, MPFR_RNDN); }; return mpfr_get_d(r25891, MPFR_RNDN); }