#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r26748 = alpha; float r26749 = beta; float r26750 = r26748 + r26749; float r26751 = r26749 * r26748; float r26752 = r26750 + r26751; float r26753 = 1.0; float r26754 = r26752 + r26753; float r26755 = 2; float r26756 = 1; float r26757 = r26755 * r26756; float r26758 = r26750 + r26757; float r26759 = r26754 / r26758; float r26760 = r26759 / r26758; float r26761 = r26758 + r26753; float r26762 = r26760 / r26761; return r26762; } double f_id(double alpha, double beta) { double r26763 = alpha; double r26764 = beta; double r26765 = r26763 + r26764; double r26766 = r26764 * r26763; double r26767 = r26765 + r26766; double r26768 = 1.0; double r26769 = r26767 + r26768; double r26770 = 2; double r26771 = 1; double r26772 = r26770 * r26771; double r26773 = r26765 + r26772; double r26774 = r26769 / r26773; double r26775 = r26774 / r26773; double r26776 = r26773 + r26768; double r26777 = r26775 / r26776; return r26777; } double f_of(float alpha, float beta) { float r26778 = beta; float r26779 = 4.193960641989379e+198; bool r26780 = r26778 <= r26779; float r26781 = 1; float r26782 = sqrt(r26781); float r26783 = alpha; float r26784 = r26778 + r26783; float r26785 = 1.0; float r26786 = 2; float r26787 = r26785 + r26786; float r26788 = r26784 + r26787; float r26789 = r26782 / r26788; float r26790 = r26778 * r26783; float r26791 = r26790 + r26778; float r26792 = r26783 + r26785; float r26793 = r26791 + r26792; float r26794 = sqrt(r26793); float r26795 = r26786 + r26783; float r26796 = r26795 + r26778; float r26797 = r26794 / r26796; float r26798 = r26789 * r26797; float r26799 = r26798 * r26797; float r26800 = 0; float r26801 = r26780 ? r26799 : r26800; return r26801; } double f_od(double alpha, double beta) { double r26802 = beta; double r26803 = 4.193960641989379e+198; bool r26804 = r26802 <= r26803; double r26805 = 1; double r26806 = sqrt(r26805); double r26807 = alpha; double r26808 = r26802 + r26807; double r26809 = 1.0; double r26810 = 2; double r26811 = r26809 + r26810; double r26812 = r26808 + r26811; double r26813 = r26806 / r26812; double r26814 = r26802 * r26807; double r26815 = r26814 + r26802; double r26816 = r26807 + r26809; double r26817 = r26815 + r26816; double r26818 = sqrt(r26817); double r26819 = r26810 + r26807; double r26820 = r26819 + r26802; double r26821 = r26818 / r26820; double r26822 = r26813 * r26821; double r26823 = r26822 * r26821; double r26824 = 0; double r26825 = r26804 ? r26823 : r26824; return r26825; } 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 r26826, r26827, r26828, r26829, r26830, r26831, r26832, r26833, r26834, r26835, r26836, r26837, r26838, r26839, r26840; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r26826); mpfr_init(r26827); mpfr_init(r26828); mpfr_init(r26829); mpfr_init(r26830); mpfr_init_set_str(r26831, "1.0", 10, MPFR_RNDN); mpfr_init(r26832); mpfr_init_set_str(r26833, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26834, "1", 10, MPFR_RNDN); mpfr_init(r26835); mpfr_init(r26836); mpfr_init(r26837); mpfr_init(r26838); mpfr_init(r26839); mpfr_init(r26840); } double f_im(double alpha, double beta) { mpfr_set_d(r26826, alpha, MPFR_RNDN); mpfr_set_d(r26827, beta, MPFR_RNDN); mpfr_add(r26828, r26826, r26827, MPFR_RNDN); mpfr_mul(r26829, r26827, r26826, MPFR_RNDN); mpfr_add(r26830, r26828, r26829, MPFR_RNDN); ; mpfr_add(r26832, r26830, r26831, MPFR_RNDN); ; ; mpfr_mul(r26835, r26833, r26834, MPFR_RNDN); mpfr_add(r26836, r26828, r26835, MPFR_RNDN); mpfr_div(r26837, r26832, r26836, MPFR_RNDN); mpfr_div(r26838, r26837, r26836, MPFR_RNDN); mpfr_add(r26839, r26836, r26831, MPFR_RNDN); mpfr_div(r26840, r26838, r26839, MPFR_RNDN); return mpfr_get_d(r26840, MPFR_RNDN); } static mpfr_t r26841, r26842, r26843, r26844, r26845, r26846, r26847, r26848, r26849, r26850, r26851, r26852, r26853, r26854, r26855, r26856, r26857, r26858, r26859, r26860, r26861, r26862, r26863, r26864; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r26841); mpfr_init_set_str(r26842, "4.193960641989379e+198", 10, MPFR_RNDN); mpfr_init(r26843); mpfr_init_set_str(r26844, "1", 10, MPFR_RNDN); mpfr_init(r26845); mpfr_init(r26846); mpfr_init(r26847); mpfr_init_set_str(r26848, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r26849, "2", 10, MPFR_RNDN); mpfr_init(r26850); mpfr_init(r26851); mpfr_init(r26852); mpfr_init(r26853); mpfr_init(r26854); mpfr_init(r26855); mpfr_init(r26856); mpfr_init(r26857); mpfr_init(r26858); mpfr_init(r26859); mpfr_init(r26860); mpfr_init(r26861); mpfr_init(r26862); mpfr_init_set_str(r26863, "0", 10, MPFR_RNDN); mpfr_init(r26864); } double f_fm(double alpha, double beta) { mpfr_set_d(r26841, beta, MPFR_RNDN); ; mpfr_set_si(r26843, mpfr_cmp(r26841, r26842) <= 0, MPFR_RNDN); ; mpfr_sqrt(r26845, r26844, MPFR_RNDN); mpfr_set_d(r26846, alpha, MPFR_RNDN); mpfr_add(r26847, r26841, r26846, MPFR_RNDN); ; ; mpfr_add(r26850, r26848, r26849, MPFR_RNDN); mpfr_add(r26851, r26847, r26850, MPFR_RNDN); mpfr_div(r26852, r26845, r26851, MPFR_RNDN); mpfr_mul(r26853, r26841, r26846, MPFR_RNDN); mpfr_add(r26854, r26853, r26841, MPFR_RNDN); mpfr_add(r26855, r26846, r26848, MPFR_RNDN); mpfr_add(r26856, r26854, r26855, MPFR_RNDN); mpfr_sqrt(r26857, r26856, MPFR_RNDN); mpfr_add(r26858, r26849, r26846, MPFR_RNDN); mpfr_add(r26859, r26858, r26841, MPFR_RNDN); mpfr_div(r26860, r26857, r26859, MPFR_RNDN); mpfr_mul(r26861, r26852, r26860, MPFR_RNDN); mpfr_mul(r26862, r26861, r26860, MPFR_RNDN); ; if (mpfr_get_si(r26843, MPFR_RNDN)) { mpfr_set(r26864, r26862, MPFR_RNDN); } else { mpfr_set(r26864, r26863, MPFR_RNDN); }; return mpfr_get_d(r26864, MPFR_RNDN); } static mpfr_t r26865, r26866, r26867, r26868, r26869, r26870, r26871, r26872, r26873, r26874, r26875, r26876, r26877, r26878, r26879, r26880, r26881, r26882, r26883, r26884, r26885, r26886, r26887, r26888; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r26865); mpfr_init_set_str(r26866, "4.193960641989379e+198", 10, MPFR_RNDN); mpfr_init(r26867); mpfr_init_set_str(r26868, "1", 10, MPFR_RNDN); mpfr_init(r26869); mpfr_init(r26870); mpfr_init(r26871); mpfr_init_set_str(r26872, "1.0", 10, MPFR_RNDN); mpfr_init_set_str(r26873, "2", 10, MPFR_RNDN); mpfr_init(r26874); mpfr_init(r26875); mpfr_init(r26876); mpfr_init(r26877); mpfr_init(r26878); mpfr_init(r26879); mpfr_init(r26880); mpfr_init(r26881); mpfr_init(r26882); mpfr_init(r26883); mpfr_init(r26884); mpfr_init(r26885); mpfr_init(r26886); mpfr_init_set_str(r26887, "0", 10, MPFR_RNDN); mpfr_init(r26888); } double f_dm(double alpha, double beta) { mpfr_set_d(r26865, beta, MPFR_RNDN); ; mpfr_set_si(r26867, mpfr_cmp(r26865, r26866) <= 0, MPFR_RNDN); ; mpfr_sqrt(r26869, r26868, MPFR_RNDN); mpfr_set_d(r26870, alpha, MPFR_RNDN); mpfr_add(r26871, r26865, r26870, MPFR_RNDN); ; ; mpfr_add(r26874, r26872, r26873, MPFR_RNDN); mpfr_add(r26875, r26871, r26874, MPFR_RNDN); mpfr_div(r26876, r26869, r26875, MPFR_RNDN); mpfr_mul(r26877, r26865, r26870, MPFR_RNDN); mpfr_add(r26878, r26877, r26865, MPFR_RNDN); mpfr_add(r26879, r26870, r26872, MPFR_RNDN); mpfr_add(r26880, r26878, r26879, MPFR_RNDN); mpfr_sqrt(r26881, r26880, MPFR_RNDN); mpfr_add(r26882, r26873, r26870, MPFR_RNDN); mpfr_add(r26883, r26882, r26865, MPFR_RNDN); mpfr_div(r26884, r26881, r26883, MPFR_RNDN); mpfr_mul(r26885, r26876, r26884, MPFR_RNDN); mpfr_mul(r26886, r26885, r26884, MPFR_RNDN); ; if (mpfr_get_si(r26867, MPFR_RNDN)) { mpfr_set(r26888, r26886, MPFR_RNDN); } else { mpfr_set(r26888, r26887, MPFR_RNDN); }; return mpfr_get_d(r26888, MPFR_RNDN); }