#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r26780 = alpha; float r26781 = beta; float r26782 = r26780 + r26781; float r26783 = r26781 * r26780; float r26784 = r26782 + r26783; float r26785 = 1.0; float r26786 = r26784 + r26785; float r26787 = 2; float r26788 = 1; float r26789 = r26787 * r26788; float r26790 = r26782 + r26789; float r26791 = r26786 / r26790; float r26792 = r26791 / r26790; float r26793 = r26790 + r26785; float r26794 = r26792 / r26793; return r26794; } double f_id(double alpha, double beta) { double r26795 = alpha; double r26796 = beta; double r26797 = r26795 + r26796; double r26798 = r26796 * r26795; double r26799 = r26797 + r26798; double r26800 = 1.0; double r26801 = r26799 + r26800; double r26802 = 2; double r26803 = 1; double r26804 = r26802 * r26803; double r26805 = r26797 + r26804; double r26806 = r26801 / r26805; double r26807 = r26806 / r26805; double r26808 = r26805 + r26800; double r26809 = r26807 / r26808; return r26809; } double f_of(float alpha, float beta) { float r26810 = alpha; float r26811 = 1.5394648535800867e+190; bool r26812 = r26810 <= r26811; float r26813 = beta; float r26814 = r26813 * r26810; float r26815 = r26814 + r26813; float r26816 = 1.0; float r26817 = r26810 + r26816; float r26818 = r26815 + r26817; float r26819 = sqrt(r26818); float r26820 = sqrt(r26819); float r26821 = r26813 + r26810; float r26822 = 2; float r26823 = r26816 + r26822; float r26824 = r26821 + r26823; float r26825 = r26820 / r26824; float r26826 = r26822 + r26810; float r26827 = r26826 + r26813; float r26828 = r26820 / r26827; float r26829 = r26825 * r26828; float r26830 = r26819 / r26827; float r26831 = r26829 * r26830; float r26832 = 0; float r26833 = r26812 ? r26831 : r26832; return r26833; } double f_od(double alpha, double beta) { double r26834 = alpha; double r26835 = 1.5394648535800867e+190; bool r26836 = r26834 <= r26835; double r26837 = beta; double r26838 = r26837 * r26834; double r26839 = r26838 + r26837; double r26840 = 1.0; double r26841 = r26834 + r26840; double r26842 = r26839 + r26841; double r26843 = sqrt(r26842); double r26844 = sqrt(r26843); double r26845 = r26837 + r26834; double r26846 = 2; double r26847 = r26840 + r26846; double r26848 = r26845 + r26847; double r26849 = r26844 / r26848; double r26850 = r26846 + r26834; double r26851 = r26850 + r26837; double r26852 = r26844 / r26851; double r26853 = r26849 * r26852; double r26854 = r26843 / r26851; double r26855 = r26853 * r26854; double r26856 = 0; double r26857 = r26836 ? r26855 : r26856; return r26857; } 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 r26858, r26859, r26860, r26861, r26862, r26863, r26864, r26865, r26866, r26867, r26868, r26869, r26870, r26871, r26872; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r26858); mpfr_init(r26859); mpfr_init(r26860); mpfr_init(r26861); mpfr_init(r26862); mpfr_init_set_str(r26863, "1.0", 10, MPFR_RNDN); mpfr_init(r26864); mpfr_init_set_str(r26865, "2", 10, MPFR_RNDN); mpfr_init_set_str(r26866, "1", 10, MPFR_RNDN); mpfr_init(r26867); mpfr_init(r26868); mpfr_init(r26869); mpfr_init(r26870); mpfr_init(r26871); mpfr_init(r26872); } double f_im(double alpha, double beta) { mpfr_set_d(r26858, alpha, MPFR_RNDN); mpfr_set_d(r26859, beta, MPFR_RNDN); mpfr_add(r26860, r26858, r26859, MPFR_RNDN); mpfr_mul(r26861, r26859, r26858, MPFR_RNDN); mpfr_add(r26862, r26860, r26861, MPFR_RNDN); ; mpfr_add(r26864, r26862, r26863, MPFR_RNDN); ; ; mpfr_mul(r26867, r26865, r26866, MPFR_RNDN); mpfr_add(r26868, r26860, r26867, MPFR_RNDN); mpfr_div(r26869, r26864, r26868, MPFR_RNDN); mpfr_div(r26870, r26869, r26868, MPFR_RNDN); mpfr_add(r26871, r26868, r26863, MPFR_RNDN); mpfr_div(r26872, r26870, r26871, MPFR_RNDN); return mpfr_get_d(r26872, MPFR_RNDN); } static mpfr_t r26873, r26874, r26875, r26876, r26877, r26878, r26879, r26880, r26881, r26882, r26883, r26884, r26885, r26886, r26887, r26888, r26889, r26890, r26891, r26892, r26893, r26894, r26895, r26896; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r26873); mpfr_init_set_str(r26874, "1.5394648535800867e+190", 10, MPFR_RNDN); mpfr_init(r26875); mpfr_init(r26876); mpfr_init(r26877); mpfr_init(r26878); mpfr_init_set_str(r26879, "1.0", 10, MPFR_RNDN); mpfr_init(r26880); mpfr_init(r26881); mpfr_init(r26882); mpfr_init(r26883); mpfr_init(r26884); mpfr_init_set_str(r26885, "2", 10, MPFR_RNDN); mpfr_init(r26886); mpfr_init(r26887); mpfr_init(r26888); mpfr_init(r26889); mpfr_init(r26890); mpfr_init(r26891); mpfr_init(r26892); mpfr_init(r26893); mpfr_init(r26894); mpfr_init_set_str(r26895, "0", 10, MPFR_RNDN); mpfr_init(r26896); } double f_fm(double alpha, double beta) { mpfr_set_d(r26873, alpha, MPFR_RNDN); ; mpfr_set_si(r26875, mpfr_cmp(r26873, r26874) <= 0, MPFR_RNDN); mpfr_set_d(r26876, beta, MPFR_RNDN); mpfr_mul(r26877, r26876, r26873, MPFR_RNDN); mpfr_add(r26878, r26877, r26876, MPFR_RNDN); ; mpfr_add(r26880, r26873, r26879, MPFR_RNDN); mpfr_add(r26881, r26878, r26880, MPFR_RNDN); mpfr_sqrt(r26882, r26881, MPFR_RNDN); mpfr_sqrt(r26883, r26882, MPFR_RNDN); mpfr_add(r26884, r26876, r26873, MPFR_RNDN); ; mpfr_add(r26886, r26879, r26885, MPFR_RNDN); mpfr_add(r26887, r26884, r26886, MPFR_RNDN); mpfr_div(r26888, r26883, r26887, MPFR_RNDN); mpfr_add(r26889, r26885, r26873, MPFR_RNDN); mpfr_add(r26890, r26889, r26876, MPFR_RNDN); mpfr_div(r26891, r26883, r26890, MPFR_RNDN); mpfr_mul(r26892, r26888, r26891, MPFR_RNDN); mpfr_div(r26893, r26882, r26890, MPFR_RNDN); mpfr_mul(r26894, r26892, r26893, MPFR_RNDN); ; if (mpfr_get_si(r26875, MPFR_RNDN)) { mpfr_set(r26896, r26894, MPFR_RNDN); } else { mpfr_set(r26896, r26895, MPFR_RNDN); }; return mpfr_get_d(r26896, MPFR_RNDN); } static mpfr_t r26897, r26898, r26899, r26900, r26901, r26902, r26903, r26904, r26905, r26906, r26907, r26908, r26909, r26910, r26911, r26912, r26913, r26914, r26915, r26916, r26917, r26918, r26919, r26920; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r26897); mpfr_init_set_str(r26898, "1.5394648535800867e+190", 10, MPFR_RNDN); mpfr_init(r26899); mpfr_init(r26900); mpfr_init(r26901); mpfr_init(r26902); mpfr_init_set_str(r26903, "1.0", 10, MPFR_RNDN); mpfr_init(r26904); mpfr_init(r26905); mpfr_init(r26906); mpfr_init(r26907); mpfr_init(r26908); mpfr_init_set_str(r26909, "2", 10, MPFR_RNDN); mpfr_init(r26910); mpfr_init(r26911); mpfr_init(r26912); mpfr_init(r26913); mpfr_init(r26914); mpfr_init(r26915); mpfr_init(r26916); mpfr_init(r26917); mpfr_init(r26918); mpfr_init_set_str(r26919, "0", 10, MPFR_RNDN); mpfr_init(r26920); } double f_dm(double alpha, double beta) { mpfr_set_d(r26897, alpha, MPFR_RNDN); ; mpfr_set_si(r26899, mpfr_cmp(r26897, r26898) <= 0, MPFR_RNDN); mpfr_set_d(r26900, beta, MPFR_RNDN); mpfr_mul(r26901, r26900, r26897, MPFR_RNDN); mpfr_add(r26902, r26901, r26900, MPFR_RNDN); ; mpfr_add(r26904, r26897, r26903, MPFR_RNDN); mpfr_add(r26905, r26902, r26904, MPFR_RNDN); mpfr_sqrt(r26906, r26905, MPFR_RNDN); mpfr_sqrt(r26907, r26906, MPFR_RNDN); mpfr_add(r26908, r26900, r26897, MPFR_RNDN); ; mpfr_add(r26910, r26903, r26909, MPFR_RNDN); mpfr_add(r26911, r26908, r26910, MPFR_RNDN); mpfr_div(r26912, r26907, r26911, MPFR_RNDN); mpfr_add(r26913, r26909, r26897, MPFR_RNDN); mpfr_add(r26914, r26913, r26900, MPFR_RNDN); mpfr_div(r26915, r26907, r26914, MPFR_RNDN); mpfr_mul(r26916, r26912, r26915, MPFR_RNDN); mpfr_div(r26917, r26906, r26914, MPFR_RNDN); mpfr_mul(r26918, r26916, r26917, MPFR_RNDN); ; if (mpfr_get_si(r26899, MPFR_RNDN)) { mpfr_set(r26920, r26918, MPFR_RNDN); } else { mpfr_set(r26920, r26919, MPFR_RNDN); }; return mpfr_get_d(r26920, MPFR_RNDN); }