#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r25789 = i; float r25790 = alpha; float r25791 = beta; float r25792 = r25790 + r25791; float r25793 = r25792 + r25789; float r25794 = r25789 * r25793; float r25795 = r25791 * r25790; float r25796 = r25795 + r25794; float r25797 = r25794 * r25796; float r25798 = 2; float r25799 = r25798 * r25789; float r25800 = r25792 + r25799; float r25801 = r25800 * r25800; float r25802 = r25797 / r25801; float r25803 = 1.0; float r25804 = r25801 - r25803; float r25805 = r25802 / r25804; return r25805; } double f_id(double alpha, double beta, double i) { double r25806 = i; double r25807 = alpha; double r25808 = beta; double r25809 = r25807 + r25808; double r25810 = r25809 + r25806; double r25811 = r25806 * r25810; double r25812 = r25808 * r25807; double r25813 = r25812 + r25811; double r25814 = r25811 * r25813; double r25815 = 2; double r25816 = r25815 * r25806; double r25817 = r25809 + r25816; double r25818 = r25817 * r25817; double r25819 = r25814 / r25818; double r25820 = 1.0; double r25821 = r25818 - r25820; double r25822 = r25819 / r25821; return r25822; } double f_of(float alpha, float beta, float i) { float r25823 = alpha; float r25824 = i; float r25825 = r25823 + r25824; float r25826 = beta; float r25827 = r25825 + r25826; float r25828 = r25826 * r25823; float r25829 = fma(r25827, r25824, r25828); float r25830 = 2; float r25831 = fma(r25824, r25830, r25823); float r25832 = r25831 + r25826; float r25833 = 1.0; float r25834 = -r25833; float r25835 = fma(r25832, r25832, r25834); float r25836 = r25829 / r25835; float r25837 = r25827 / r25832; float r25838 = cbrt(r25837); float r25839 = r25838 * r25838; float r25840 = r25839 * r25838; float r25841 = r25824 / r25832; float r25842 = r25840 * r25841; float r25843 = r25836 * r25842; float r25844 = 0.06250026085862662; bool r25845 = r25843 <= r25844; float r25846 = sqrt(r25836); float r25847 = r25846 * r25846; float r25848 = r25837 * r25841; float r25849 = r25847 * r25848; float r25850 = 0; float r25851 = r25845 ? r25849 : r25850; return r25851; } double f_od(double alpha, double beta, double i) { double r25852 = alpha; double r25853 = i; double r25854 = r25852 + r25853; double r25855 = beta; double r25856 = r25854 + r25855; double r25857 = r25855 * r25852; double r25858 = fma(r25856, r25853, r25857); double r25859 = 2; double r25860 = fma(r25853, r25859, r25852); double r25861 = r25860 + r25855; double r25862 = 1.0; double r25863 = -r25862; double r25864 = fma(r25861, r25861, r25863); double r25865 = r25858 / r25864; double r25866 = r25856 / r25861; double r25867 = cbrt(r25866); double r25868 = r25867 * r25867; double r25869 = r25868 * r25867; double r25870 = r25853 / r25861; double r25871 = r25869 * r25870; double r25872 = r25865 * r25871; double r25873 = 0.06250026085862662; bool r25874 = r25872 <= r25873; double r25875 = sqrt(r25865); double r25876 = r25875 * r25875; double r25877 = r25866 * r25870; double r25878 = r25876 * r25877; double r25879 = 0; double r25880 = r25874 ? r25878 : r25879; return r25880; } 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 r25881, r25882, r25883, r25884, r25885, r25886, r25887, r25888, r25889, r25890, r25891, r25892, r25893, r25894, r25895, r25896, r25897; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25881); mpfr_init(r25882); mpfr_init(r25883); mpfr_init(r25884); mpfr_init(r25885); mpfr_init(r25886); mpfr_init(r25887); mpfr_init(r25888); mpfr_init(r25889); mpfr_init_set_str(r25890, "2", 10, MPFR_RNDN); mpfr_init(r25891); mpfr_init(r25892); mpfr_init(r25893); mpfr_init(r25894); mpfr_init_set_str(r25895, "1.0", 10, MPFR_RNDN); mpfr_init(r25896); mpfr_init(r25897); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r25881, i, MPFR_RNDN); mpfr_set_d(r25882, alpha, MPFR_RNDN); mpfr_set_d(r25883, beta, MPFR_RNDN); mpfr_add(r25884, r25882, r25883, MPFR_RNDN); mpfr_add(r25885, r25884, r25881, MPFR_RNDN); mpfr_mul(r25886, r25881, r25885, MPFR_RNDN); mpfr_mul(r25887, r25883, r25882, MPFR_RNDN); mpfr_add(r25888, r25887, r25886, MPFR_RNDN); mpfr_mul(r25889, r25886, r25888, MPFR_RNDN); ; mpfr_mul(r25891, r25890, r25881, MPFR_RNDN); mpfr_add(r25892, r25884, r25891, MPFR_RNDN); mpfr_mul(r25893, r25892, r25892, MPFR_RNDN); mpfr_div(r25894, r25889, r25893, MPFR_RNDN); ; mpfr_sub(r25896, r25893, r25895, MPFR_RNDN); mpfr_div(r25897, r25894, r25896, MPFR_RNDN); return mpfr_get_d(r25897, MPFR_RNDN); } static mpfr_t r25898, r25899, r25900, r25901, r25902, r25903, r25904, r25905, r25906, r25907, r25908, r25909, r25910, r25911, r25912, r25913, r25914, r25915, r25916, r25917, r25918, r25919, r25920, r25921, r25922, r25923, r25924, r25925, r25926; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25898); mpfr_init(r25899); mpfr_init(r25900); mpfr_init(r25901); mpfr_init(r25902); mpfr_init(r25903); mpfr_init(r25904); mpfr_init_set_str(r25905, "2", 10, MPFR_RNDN); mpfr_init(r25906); mpfr_init(r25907); mpfr_init_set_str(r25908, "1.0", 10, MPFR_RNDN); mpfr_init(r25909); mpfr_init(r25910); mpfr_init(r25911); mpfr_init(r25912); mpfr_init(r25913); mpfr_init(r25914); mpfr_init(r25915); mpfr_init(r25916); mpfr_init(r25917); mpfr_init(r25918); mpfr_init_set_str(r25919, "0.06250026085862662", 10, MPFR_RNDN); mpfr_init(r25920); mpfr_init(r25921); mpfr_init(r25922); mpfr_init(r25923); mpfr_init(r25924); mpfr_init_set_str(r25925, "0", 10, MPFR_RNDN); mpfr_init(r25926); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r25898, alpha, MPFR_RNDN); mpfr_set_d(r25899, i, MPFR_RNDN); mpfr_add(r25900, r25898, r25899, MPFR_RNDN); mpfr_set_d(r25901, beta, MPFR_RNDN); mpfr_add(r25902, r25900, r25901, MPFR_RNDN); mpfr_mul(r25903, r25901, r25898, MPFR_RNDN); mpfr_fma(r25904, r25902, r25899, r25903, MPFR_RNDN); ; mpfr_fma(r25906, r25899, r25905, r25898, MPFR_RNDN); mpfr_add(r25907, r25906, r25901, MPFR_RNDN); ; mpfr_neg(r25909, r25908, MPFR_RNDN); mpfr_fma(r25910, r25907, r25907, r25909, MPFR_RNDN); mpfr_div(r25911, r25904, r25910, MPFR_RNDN); mpfr_div(r25912, r25902, r25907, MPFR_RNDN); mpfr_cbrt(r25913, r25912, MPFR_RNDN); mpfr_mul(r25914, r25913, r25913, MPFR_RNDN); mpfr_mul(r25915, r25914, r25913, MPFR_RNDN); mpfr_div(r25916, r25899, r25907, MPFR_RNDN); mpfr_mul(r25917, r25915, r25916, MPFR_RNDN); mpfr_mul(r25918, r25911, r25917, MPFR_RNDN); ; mpfr_set_si(r25920, mpfr_cmp(r25918, r25919) <= 0, MPFR_RNDN); mpfr_sqrt(r25921, r25911, MPFR_RNDN); mpfr_mul(r25922, r25921, r25921, MPFR_RNDN); mpfr_mul(r25923, r25912, r25916, MPFR_RNDN); mpfr_mul(r25924, r25922, r25923, MPFR_RNDN); ; if (mpfr_get_si(r25920, MPFR_RNDN)) { mpfr_set(r25926, r25924, MPFR_RNDN); } else { mpfr_set(r25926, r25925, MPFR_RNDN); }; return mpfr_get_d(r25926, MPFR_RNDN); } static mpfr_t r25927, r25928, r25929, r25930, r25931, r25932, r25933, r25934, r25935, r25936, r25937, r25938, r25939, r25940, r25941, r25942, r25943, r25944, r25945, r25946, r25947, r25948, r25949, r25950, r25951, r25952, r25953, r25954, r25955; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25927); mpfr_init(r25928); mpfr_init(r25929); mpfr_init(r25930); mpfr_init(r25931); mpfr_init(r25932); mpfr_init(r25933); mpfr_init_set_str(r25934, "2", 10, MPFR_RNDN); mpfr_init(r25935); mpfr_init(r25936); mpfr_init_set_str(r25937, "1.0", 10, MPFR_RNDN); mpfr_init(r25938); mpfr_init(r25939); mpfr_init(r25940); mpfr_init(r25941); mpfr_init(r25942); mpfr_init(r25943); mpfr_init(r25944); mpfr_init(r25945); mpfr_init(r25946); mpfr_init(r25947); mpfr_init_set_str(r25948, "0.06250026085862662", 10, MPFR_RNDN); mpfr_init(r25949); mpfr_init(r25950); mpfr_init(r25951); mpfr_init(r25952); mpfr_init(r25953); mpfr_init_set_str(r25954, "0", 10, MPFR_RNDN); mpfr_init(r25955); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r25927, alpha, MPFR_RNDN); mpfr_set_d(r25928, i, MPFR_RNDN); mpfr_add(r25929, r25927, r25928, MPFR_RNDN); mpfr_set_d(r25930, beta, MPFR_RNDN); mpfr_add(r25931, r25929, r25930, MPFR_RNDN); mpfr_mul(r25932, r25930, r25927, MPFR_RNDN); mpfr_fma(r25933, r25931, r25928, r25932, MPFR_RNDN); ; mpfr_fma(r25935, r25928, r25934, r25927, MPFR_RNDN); mpfr_add(r25936, r25935, r25930, MPFR_RNDN); ; mpfr_neg(r25938, r25937, MPFR_RNDN); mpfr_fma(r25939, r25936, r25936, r25938, MPFR_RNDN); mpfr_div(r25940, r25933, r25939, MPFR_RNDN); mpfr_div(r25941, r25931, r25936, MPFR_RNDN); mpfr_cbrt(r25942, r25941, MPFR_RNDN); mpfr_mul(r25943, r25942, r25942, MPFR_RNDN); mpfr_mul(r25944, r25943, r25942, MPFR_RNDN); mpfr_div(r25945, r25928, r25936, MPFR_RNDN); mpfr_mul(r25946, r25944, r25945, MPFR_RNDN); mpfr_mul(r25947, r25940, r25946, MPFR_RNDN); ; mpfr_set_si(r25949, mpfr_cmp(r25947, r25948) <= 0, MPFR_RNDN); mpfr_sqrt(r25950, r25940, MPFR_RNDN); mpfr_mul(r25951, r25950, r25950, MPFR_RNDN); mpfr_mul(r25952, r25941, r25945, MPFR_RNDN); mpfr_mul(r25953, r25951, r25952, MPFR_RNDN); ; if (mpfr_get_si(r25949, MPFR_RNDN)) { mpfr_set(r25955, r25953, MPFR_RNDN); } else { mpfr_set(r25955, r25954, MPFR_RNDN); }; return mpfr_get_d(r25955, MPFR_RNDN); }