#include #include #include #include #include char *name = "Octave 3.8, jcobi/3"; double f_if(float alpha, float beta) { float r23884 = alpha; float r23885 = beta; float r23886 = r23884 + r23885; float r23887 = r23885 * r23884; float r23888 = r23886 + r23887; float r23889 = 1.0; float r23890 = r23888 + r23889; float r23891 = 2; float r23892 = 1; float r23893 = r23891 * r23892; float r23894 = r23886 + r23893; float r23895 = r23890 / r23894; float r23896 = r23895 / r23894; float r23897 = r23894 + r23889; float r23898 = r23896 / r23897; return r23898; } double f_id(double alpha, double beta) { double r23899 = alpha; double r23900 = beta; double r23901 = r23899 + r23900; double r23902 = r23900 * r23899; double r23903 = r23901 + r23902; double r23904 = 1.0; double r23905 = r23903 + r23904; double r23906 = 2; double r23907 = 1; double r23908 = r23906 * r23907; double r23909 = r23901 + r23908; double r23910 = r23905 / r23909; double r23911 = r23910 / r23909; double r23912 = r23909 + r23904; double r23913 = r23911 / r23912; return r23913; } double f_of(float alpha, float beta) { float r23914 = beta; float r23915 = 1.0; float r23916 = r23914 + r23915; float r23917 = alpha; float r23918 = r23914 * r23917; float r23919 = r23916 + r23918; float r23920 = r23919 + r23917; float r23921 = 2; float r23922 = r23914 + r23921; float r23923 = r23917 + r23922; float r23924 = r23920 / r23923; float r23925 = sqrt(r23924); float r23926 = sqrt(r23923); float r23927 = r23925 / r23926; float r23928 = r23917 + r23919; float r23929 = r23914 + r23917; float r23930 = r23921 + r23929; float r23931 = r23928 / r23930; float r23932 = sqrt(r23931); float r23933 = r23917 + r23921; float r23934 = r23916 + r23933; float r23935 = sqrt(r23930); float r23936 = r23934 * r23935; float r23937 = r23932 / r23936; float r23938 = r23927 * r23937; return r23938; } double f_od(double alpha, double beta) { double r23939 = beta; double r23940 = 1.0; double r23941 = r23939 + r23940; double r23942 = alpha; double r23943 = r23939 * r23942; double r23944 = r23941 + r23943; double r23945 = r23944 + r23942; double r23946 = 2; double r23947 = r23939 + r23946; double r23948 = r23942 + r23947; double r23949 = r23945 / r23948; double r23950 = sqrt(r23949); double r23951 = sqrt(r23948); double r23952 = r23950 / r23951; double r23953 = r23942 + r23944; double r23954 = r23939 + r23942; double r23955 = r23946 + r23954; double r23956 = r23953 / r23955; double r23957 = sqrt(r23956); double r23958 = r23942 + r23946; double r23959 = r23941 + r23958; double r23960 = sqrt(r23955); double r23961 = r23959 * r23960; double r23962 = r23957 / r23961; double r23963 = r23952 * r23962; return r23963; } 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 r23964, r23965, r23966, r23967, r23968, r23969, r23970, r23971, r23972, r23973, r23974, r23975, r23976, r23977, r23978; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r23964); mpfr_init(r23965); mpfr_init(r23966); mpfr_init(r23967); mpfr_init(r23968); mpfr_init_set_str(r23969, "1.0", 10, MPFR_RNDN); mpfr_init(r23970); mpfr_init_set_str(r23971, "2", 10, MPFR_RNDN); mpfr_init_set_str(r23972, "1", 10, MPFR_RNDN); mpfr_init(r23973); mpfr_init(r23974); mpfr_init(r23975); mpfr_init(r23976); mpfr_init(r23977); mpfr_init(r23978); } double f_im(double alpha, double beta) { mpfr_set_d(r23964, alpha, MPFR_RNDN); mpfr_set_d(r23965, beta, MPFR_RNDN); mpfr_add(r23966, r23964, r23965, MPFR_RNDN); mpfr_mul(r23967, r23965, r23964, MPFR_RNDN); mpfr_add(r23968, r23966, r23967, MPFR_RNDN); ; mpfr_add(r23970, r23968, r23969, MPFR_RNDN); ; ; mpfr_mul(r23973, r23971, r23972, MPFR_RNDN); mpfr_add(r23974, r23966, r23973, MPFR_RNDN); mpfr_div(r23975, r23970, r23974, MPFR_RNDN); mpfr_div(r23976, r23975, r23974, MPFR_RNDN); mpfr_add(r23977, r23974, r23969, MPFR_RNDN); mpfr_div(r23978, r23976, r23977, MPFR_RNDN); return mpfr_get_d(r23978, MPFR_RNDN); } static mpfr_t r23979, r23980, r23981, r23982, r23983, r23984, r23985, r23986, r23987, r23988, r23989, r23990, r23991, r23992, r23993, r23994, r23995, r23996, r23997, r23998, r23999, r24000, r24001, r24002, r24003; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r23979); mpfr_init_set_str(r23980, "1.0", 10, MPFR_RNDN); mpfr_init(r23981); mpfr_init(r23982); mpfr_init(r23983); mpfr_init(r23984); mpfr_init(r23985); mpfr_init_set_str(r23986, "2", 10, MPFR_RNDN); mpfr_init(r23987); mpfr_init(r23988); mpfr_init(r23989); mpfr_init(r23990); mpfr_init(r23991); mpfr_init(r23992); mpfr_init(r23993); mpfr_init(r23994); mpfr_init(r23995); mpfr_init(r23996); mpfr_init(r23997); mpfr_init(r23998); mpfr_init(r23999); mpfr_init(r24000); mpfr_init(r24001); mpfr_init(r24002); mpfr_init(r24003); } double f_fm(double alpha, double beta) { mpfr_set_d(r23979, beta, MPFR_RNDN); ; mpfr_add(r23981, r23979, r23980, MPFR_RNDN); mpfr_set_d(r23982, alpha, MPFR_RNDN); mpfr_mul(r23983, r23979, r23982, MPFR_RNDN); mpfr_add(r23984, r23981, r23983, MPFR_RNDN); mpfr_add(r23985, r23984, r23982, MPFR_RNDN); ; mpfr_add(r23987, r23979, r23986, MPFR_RNDN); mpfr_add(r23988, r23982, r23987, MPFR_RNDN); mpfr_div(r23989, r23985, r23988, MPFR_RNDN); mpfr_sqrt(r23990, r23989, MPFR_RNDN); mpfr_sqrt(r23991, r23988, MPFR_RNDN); mpfr_div(r23992, r23990, r23991, MPFR_RNDN); mpfr_add(r23993, r23982, r23984, MPFR_RNDN); mpfr_add(r23994, r23979, r23982, MPFR_RNDN); mpfr_add(r23995, r23986, r23994, MPFR_RNDN); mpfr_div(r23996, r23993, r23995, MPFR_RNDN); mpfr_sqrt(r23997, r23996, MPFR_RNDN); mpfr_add(r23998, r23982, r23986, MPFR_RNDN); mpfr_add(r23999, r23981, r23998, MPFR_RNDN); mpfr_sqrt(r24000, r23995, MPFR_RNDN); mpfr_mul(r24001, r23999, r24000, MPFR_RNDN); mpfr_div(r24002, r23997, r24001, MPFR_RNDN); mpfr_mul(r24003, r23992, r24002, MPFR_RNDN); return mpfr_get_d(r24003, MPFR_RNDN); } static mpfr_t r24004, r24005, r24006, r24007, r24008, r24009, r24010, r24011, r24012, r24013, r24014, r24015, r24016, r24017, r24018, r24019, r24020, r24021, r24022, r24023, r24024, r24025, r24026, r24027, r24028; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r24004); mpfr_init_set_str(r24005, "1.0", 10, MPFR_RNDN); mpfr_init(r24006); mpfr_init(r24007); mpfr_init(r24008); mpfr_init(r24009); mpfr_init(r24010); mpfr_init_set_str(r24011, "2", 10, MPFR_RNDN); mpfr_init(r24012); mpfr_init(r24013); mpfr_init(r24014); mpfr_init(r24015); mpfr_init(r24016); mpfr_init(r24017); mpfr_init(r24018); mpfr_init(r24019); mpfr_init(r24020); mpfr_init(r24021); mpfr_init(r24022); mpfr_init(r24023); mpfr_init(r24024); mpfr_init(r24025); mpfr_init(r24026); mpfr_init(r24027); mpfr_init(r24028); } double f_dm(double alpha, double beta) { mpfr_set_d(r24004, beta, MPFR_RNDN); ; mpfr_add(r24006, r24004, r24005, MPFR_RNDN); mpfr_set_d(r24007, alpha, MPFR_RNDN); mpfr_mul(r24008, r24004, r24007, MPFR_RNDN); mpfr_add(r24009, r24006, r24008, MPFR_RNDN); mpfr_add(r24010, r24009, r24007, MPFR_RNDN); ; mpfr_add(r24012, r24004, r24011, MPFR_RNDN); mpfr_add(r24013, r24007, r24012, MPFR_RNDN); mpfr_div(r24014, r24010, r24013, MPFR_RNDN); mpfr_sqrt(r24015, r24014, MPFR_RNDN); mpfr_sqrt(r24016, r24013, MPFR_RNDN); mpfr_div(r24017, r24015, r24016, MPFR_RNDN); mpfr_add(r24018, r24007, r24009, MPFR_RNDN); mpfr_add(r24019, r24004, r24007, MPFR_RNDN); mpfr_add(r24020, r24011, r24019, MPFR_RNDN); mpfr_div(r24021, r24018, r24020, MPFR_RNDN); mpfr_sqrt(r24022, r24021, MPFR_RNDN); mpfr_add(r24023, r24007, r24011, MPFR_RNDN); mpfr_add(r24024, r24006, r24023, MPFR_RNDN); mpfr_sqrt(r24025, r24020, MPFR_RNDN); mpfr_mul(r24026, r24024, r24025, MPFR_RNDN); mpfr_div(r24027, r24022, r24026, MPFR_RNDN); mpfr_mul(r24028, r24017, r24027, MPFR_RNDN); return mpfr_get_d(r24028, MPFR_RNDN); }