#include #include #include #include #include char *name = "Octave 3.8, jcobi/4"; double f_if(float alpha, float beta, float i) { float r22834 = i; float r22835 = alpha; float r22836 = beta; float r22837 = r22835 + r22836; float r22838 = r22837 + r22834; float r22839 = r22834 * r22838; float r22840 = r22836 * r22835; float r22841 = r22840 + r22839; float r22842 = r22839 * r22841; float r22843 = 2; float r22844 = r22843 * r22834; float r22845 = r22837 + r22844; float r22846 = r22845 * r22845; float r22847 = r22842 / r22846; float r22848 = 1.0; float r22849 = r22846 - r22848; float r22850 = r22847 / r22849; return r22850; } double f_id(double alpha, double beta, double i) { double r22851 = i; double r22852 = alpha; double r22853 = beta; double r22854 = r22852 + r22853; double r22855 = r22854 + r22851; double r22856 = r22851 * r22855; double r22857 = r22853 * r22852; double r22858 = r22857 + r22856; double r22859 = r22856 * r22858; double r22860 = 2; double r22861 = r22860 * r22851; double r22862 = r22854 + r22861; double r22863 = r22862 * r22862; double r22864 = r22859 / r22863; double r22865 = 1.0; double r22866 = r22863 - r22865; double r22867 = r22864 / r22866; return r22867; } double f_of(float alpha, float beta, float i) { float r22868 = beta; float r22869 = 2.368086067055097e+152; bool r22870 = r22868 <= r22869; float r22871 = i; float r22872 = alpha; float r22873 = r22872 + r22868; float r22874 = 2; float r22875 = r22874 * r22871; float r22876 = r22873 + r22875; float r22877 = r22876 * r22876; float r22878 = r22868 * r22872; float r22879 = r22873 + r22871; float r22880 = r22871 * r22879; float r22881 = r22878 + r22880; float r22882 = r22877 / r22881; float r22883 = sqrt(r22882); float r22884 = r22871 / r22883; float r22885 = 1.0; float r22886 = r22877 - r22885; float r22887 = sqrt(r22886); float r22888 = r22884 / r22887; float r22889 = r22879 / r22883; float r22890 = r22889 / r22887; float r22891 = r22888 * r22890; float r22892 = 0; float r22893 = r22870 ? r22891 : r22892; return r22893; } double f_od(double alpha, double beta, double i) { double r22894 = beta; double r22895 = 2.368086067055097e+152; bool r22896 = r22894 <= r22895; double r22897 = i; double r22898 = alpha; double r22899 = r22898 + r22894; double r22900 = 2; double r22901 = r22900 * r22897; double r22902 = r22899 + r22901; double r22903 = r22902 * r22902; double r22904 = r22894 * r22898; double r22905 = r22899 + r22897; double r22906 = r22897 * r22905; double r22907 = r22904 + r22906; double r22908 = r22903 / r22907; double r22909 = sqrt(r22908); double r22910 = r22897 / r22909; double r22911 = 1.0; double r22912 = r22903 - r22911; double r22913 = sqrt(r22912); double r22914 = r22910 / r22913; double r22915 = r22905 / r22909; double r22916 = r22915 / r22913; double r22917 = r22914 * r22916; double r22918 = 0; double r22919 = r22896 ? r22917 : r22918; return r22919; } 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 r22920, r22921, r22922, r22923, r22924, r22925, r22926, r22927, r22928, r22929, r22930, r22931, r22932, r22933, r22934, r22935, r22936; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22920); mpfr_init(r22921); mpfr_init(r22922); mpfr_init(r22923); mpfr_init(r22924); mpfr_init(r22925); mpfr_init(r22926); mpfr_init(r22927); mpfr_init(r22928); mpfr_init_set_str(r22929, "2", 10, MPFR_RNDN); mpfr_init(r22930); mpfr_init(r22931); mpfr_init(r22932); mpfr_init(r22933); mpfr_init_set_str(r22934, "1.0", 10, MPFR_RNDN); mpfr_init(r22935); mpfr_init(r22936); } double f_im(double alpha, double beta, double i) { mpfr_set_d(r22920, i, MPFR_RNDN); mpfr_set_d(r22921, alpha, MPFR_RNDN); mpfr_set_d(r22922, beta, MPFR_RNDN); mpfr_add(r22923, r22921, r22922, MPFR_RNDN); mpfr_add(r22924, r22923, r22920, MPFR_RNDN); mpfr_mul(r22925, r22920, r22924, MPFR_RNDN); mpfr_mul(r22926, r22922, r22921, MPFR_RNDN); mpfr_add(r22927, r22926, r22925, MPFR_RNDN); mpfr_mul(r22928, r22925, r22927, MPFR_RNDN); ; mpfr_mul(r22930, r22929, r22920, MPFR_RNDN); mpfr_add(r22931, r22923, r22930, MPFR_RNDN); mpfr_mul(r22932, r22931, r22931, MPFR_RNDN); mpfr_div(r22933, r22928, r22932, MPFR_RNDN); ; mpfr_sub(r22935, r22932, r22934, MPFR_RNDN); mpfr_div(r22936, r22933, r22935, MPFR_RNDN); return mpfr_get_d(r22936, MPFR_RNDN); } static mpfr_t r22937, r22938, r22939, r22940, r22941, r22942, r22943, r22944, r22945, r22946, r22947, r22948, r22949, r22950, r22951, r22952, r22953, r22954, r22955, r22956, r22957, r22958, r22959, r22960, r22961, r22962; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22937); mpfr_init_set_str(r22938, "2.368086067055097e+152", 10, MPFR_RNDN); mpfr_init(r22939); mpfr_init(r22940); mpfr_init(r22941); mpfr_init(r22942); mpfr_init_set_str(r22943, "2", 10, MPFR_RNDN); mpfr_init(r22944); mpfr_init(r22945); mpfr_init(r22946); mpfr_init(r22947); mpfr_init(r22948); mpfr_init(r22949); mpfr_init(r22950); mpfr_init(r22951); mpfr_init(r22952); mpfr_init(r22953); mpfr_init_set_str(r22954, "1.0", 10, MPFR_RNDN); mpfr_init(r22955); mpfr_init(r22956); mpfr_init(r22957); mpfr_init(r22958); mpfr_init(r22959); mpfr_init(r22960); mpfr_init_set_str(r22961, "0", 10, MPFR_RNDN); mpfr_init(r22962); } double f_fm(double alpha, double beta, double i) { mpfr_set_d(r22937, beta, MPFR_RNDN); ; mpfr_set_si(r22939, mpfr_cmp(r22937, r22938) <= 0, MPFR_RNDN); mpfr_set_d(r22940, i, MPFR_RNDN); mpfr_set_d(r22941, alpha, MPFR_RNDN); mpfr_add(r22942, r22941, r22937, MPFR_RNDN); ; mpfr_mul(r22944, r22943, r22940, MPFR_RNDN); mpfr_add(r22945, r22942, r22944, MPFR_RNDN); mpfr_mul(r22946, r22945, r22945, MPFR_RNDN); mpfr_mul(r22947, r22937, r22941, MPFR_RNDN); mpfr_add(r22948, r22942, r22940, MPFR_RNDN); mpfr_mul(r22949, r22940, r22948, MPFR_RNDN); mpfr_add(r22950, r22947, r22949, MPFR_RNDN); mpfr_div(r22951, r22946, r22950, MPFR_RNDN); mpfr_sqrt(r22952, r22951, MPFR_RNDN); mpfr_div(r22953, r22940, r22952, MPFR_RNDN); ; mpfr_sub(r22955, r22946, r22954, MPFR_RNDN); mpfr_sqrt(r22956, r22955, MPFR_RNDN); mpfr_div(r22957, r22953, r22956, MPFR_RNDN); mpfr_div(r22958, r22948, r22952, MPFR_RNDN); mpfr_div(r22959, r22958, r22956, MPFR_RNDN); mpfr_mul(r22960, r22957, r22959, MPFR_RNDN); ; if (mpfr_get_si(r22939, MPFR_RNDN)) { mpfr_set(r22962, r22960, MPFR_RNDN); } else { mpfr_set(r22962, r22961, MPFR_RNDN); }; return mpfr_get_d(r22962, MPFR_RNDN); } static mpfr_t r22963, r22964, r22965, r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975, r22976, r22977, r22978, r22979, r22980, r22981, r22982, r22983, r22984, r22985, r22986, r22987, r22988; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22963); mpfr_init_set_str(r22964, "2.368086067055097e+152", 10, MPFR_RNDN); mpfr_init(r22965); mpfr_init(r22966); mpfr_init(r22967); mpfr_init(r22968); mpfr_init_set_str(r22969, "2", 10, MPFR_RNDN); mpfr_init(r22970); mpfr_init(r22971); mpfr_init(r22972); mpfr_init(r22973); mpfr_init(r22974); mpfr_init(r22975); mpfr_init(r22976); mpfr_init(r22977); mpfr_init(r22978); mpfr_init(r22979); mpfr_init_set_str(r22980, "1.0", 10, MPFR_RNDN); mpfr_init(r22981); mpfr_init(r22982); mpfr_init(r22983); mpfr_init(r22984); mpfr_init(r22985); mpfr_init(r22986); mpfr_init_set_str(r22987, "0", 10, MPFR_RNDN); mpfr_init(r22988); } double f_dm(double alpha, double beta, double i) { mpfr_set_d(r22963, beta, MPFR_RNDN); ; mpfr_set_si(r22965, mpfr_cmp(r22963, r22964) <= 0, MPFR_RNDN); mpfr_set_d(r22966, i, MPFR_RNDN); mpfr_set_d(r22967, alpha, MPFR_RNDN); mpfr_add(r22968, r22967, r22963, MPFR_RNDN); ; mpfr_mul(r22970, r22969, r22966, MPFR_RNDN); mpfr_add(r22971, r22968, r22970, MPFR_RNDN); mpfr_mul(r22972, r22971, r22971, MPFR_RNDN); mpfr_mul(r22973, r22963, r22967, MPFR_RNDN); mpfr_add(r22974, r22968, r22966, MPFR_RNDN); mpfr_mul(r22975, r22966, r22974, MPFR_RNDN); mpfr_add(r22976, r22973, r22975, MPFR_RNDN); mpfr_div(r22977, r22972, r22976, MPFR_RNDN); mpfr_sqrt(r22978, r22977, MPFR_RNDN); mpfr_div(r22979, r22966, r22978, MPFR_RNDN); ; mpfr_sub(r22981, r22972, r22980, MPFR_RNDN); mpfr_sqrt(r22982, r22981, MPFR_RNDN); mpfr_div(r22983, r22979, r22982, MPFR_RNDN); mpfr_div(r22984, r22974, r22978, MPFR_RNDN); mpfr_div(r22985, r22984, r22982, MPFR_RNDN); mpfr_mul(r22986, r22983, r22985, MPFR_RNDN); ; if (mpfr_get_si(r22965, MPFR_RNDN)) { mpfr_set(r22988, r22986, MPFR_RNDN); } else { mpfr_set(r22988, r22987, MPFR_RNDN); }; return mpfr_get_d(r22988, MPFR_RNDN); }