#include #include #include #include #include char *name = "Numeric.SpecFunctions:logGammaL from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a) { float r22898 = x; float r22899 = y; float r22900 = r22898 + r22899; float r22901 = log(r22900); float r22902 = z; float r22903 = log(r22902); float r22904 = r22901 + r22903; float r22905 = t; float r22906 = r22904 - r22905; float r22907 = a; float r22908 = 0.5; float r22909 = r22907 - r22908; float r22910 = log(r22905); float r22911 = r22909 * r22910; float r22912 = r22906 + r22911; return r22912; } double f_id(double x, double y, double z, double t, double a) { double r22913 = x; double r22914 = y; double r22915 = r22913 + r22914; double r22916 = log(r22915); double r22917 = z; double r22918 = log(r22917); double r22919 = r22916 + r22918; double r22920 = t; double r22921 = r22919 - r22920; double r22922 = a; double r22923 = 0.5; double r22924 = r22922 - r22923; double r22925 = log(r22920); double r22926 = r22924 * r22925; double r22927 = r22921 + r22926; return r22927; } double f_of(float x, float y, float z, float t, float a) { float r22928 = x; float r22929 = y; float r22930 = r22928 + r22929; float r22931 = log(r22930); float r22932 = 3; float r22933 = pow(r22931, r22932); float r22934 = z; float r22935 = log(r22934); float r22936 = pow(r22935, r22932); float r22937 = r22933 + r22936; float r22938 = r22929 + r22928; float r22939 = log(r22938); float r22940 = r22939 * r22939; float r22941 = r22939 - r22935; float r22942 = r22935 * r22941; float r22943 = r22940 - r22942; float r22944 = r22937 / r22943; float r22945 = t; float r22946 = r22944 - r22945; float r22947 = a; float r22948 = 0.5; float r22949 = r22947 - r22948; float r22950 = log(r22945); float r22951 = r22949 * r22950; float r22952 = r22946 + r22951; return r22952; } double f_od(double x, double y, double z, double t, double a) { double r22953 = x; double r22954 = y; double r22955 = r22953 + r22954; double r22956 = log(r22955); double r22957 = 3; double r22958 = pow(r22956, r22957); double r22959 = z; double r22960 = log(r22959); double r22961 = pow(r22960, r22957); double r22962 = r22958 + r22961; double r22963 = r22954 + r22953; double r22964 = log(r22963); double r22965 = r22964 * r22964; double r22966 = r22964 - r22960; double r22967 = r22960 * r22966; double r22968 = r22965 - r22967; double r22969 = r22962 / r22968; double r22970 = t; double r22971 = r22969 - r22970; double r22972 = a; double r22973 = 0.5; double r22974 = r22972 - r22973; double r22975 = log(r22970); double r22976 = r22974 * r22975; double r22977 = r22971 + r22976; return r22977; } 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 r22978, r22979, r22980, r22981, r22982, r22983, r22984, r22985, r22986, r22987, r22988, r22989, r22990, r22991, r22992; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22978); mpfr_init(r22979); mpfr_init(r22980); mpfr_init(r22981); mpfr_init(r22982); mpfr_init(r22983); mpfr_init(r22984); mpfr_init(r22985); mpfr_init(r22986); mpfr_init(r22987); mpfr_init_set_str(r22988, "0.5", 10, MPFR_RNDN); mpfr_init(r22989); mpfr_init(r22990); mpfr_init(r22991); mpfr_init(r22992); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r22978, x, MPFR_RNDN); mpfr_set_d(r22979, y, MPFR_RNDN); mpfr_add(r22980, r22978, r22979, MPFR_RNDN); mpfr_log(r22981, r22980, MPFR_RNDN); mpfr_set_d(r22982, z, MPFR_RNDN); mpfr_log(r22983, r22982, MPFR_RNDN); mpfr_add(r22984, r22981, r22983, MPFR_RNDN); mpfr_set_d(r22985, t, MPFR_RNDN); mpfr_sub(r22986, r22984, r22985, MPFR_RNDN); mpfr_set_d(r22987, a, MPFR_RNDN); ; mpfr_sub(r22989, r22987, r22988, MPFR_RNDN); mpfr_log(r22990, r22985, MPFR_RNDN); mpfr_mul(r22991, r22989, r22990, MPFR_RNDN); mpfr_add(r22992, r22986, r22991, MPFR_RNDN); return mpfr_get_d(r22992, MPFR_RNDN); } static mpfr_t r22993, r22994, r22995, r22996, r22997, r22998, r22999, r23000, r23001, r23002, r23003, r23004, r23005, r23006, r23007, r23008, r23009, r23010, r23011, r23012, r23013, r23014, r23015, r23016, r23017; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22993); mpfr_init(r22994); mpfr_init(r22995); mpfr_init(r22996); mpfr_init_set_str(r22997, "3", 10, MPFR_RNDN); mpfr_init(r22998); mpfr_init(r22999); mpfr_init(r23000); mpfr_init(r23001); mpfr_init(r23002); mpfr_init(r23003); mpfr_init(r23004); mpfr_init(r23005); mpfr_init(r23006); mpfr_init(r23007); mpfr_init(r23008); mpfr_init(r23009); mpfr_init(r23010); mpfr_init(r23011); mpfr_init(r23012); mpfr_init_set_str(r23013, "0.5", 10, MPFR_RNDN); mpfr_init(r23014); mpfr_init(r23015); mpfr_init(r23016); mpfr_init(r23017); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r22993, x, MPFR_RNDN); mpfr_set_d(r22994, y, MPFR_RNDN); mpfr_add(r22995, r22993, r22994, MPFR_RNDN); mpfr_log(r22996, r22995, MPFR_RNDN); ; mpfr_pow(r22998, r22996, r22997, MPFR_RNDN); mpfr_set_d(r22999, z, MPFR_RNDN); mpfr_log(r23000, r22999, MPFR_RNDN); mpfr_pow(r23001, r23000, r22997, MPFR_RNDN); mpfr_add(r23002, r22998, r23001, MPFR_RNDN); mpfr_add(r23003, r22994, r22993, MPFR_RNDN); mpfr_log(r23004, r23003, MPFR_RNDN); mpfr_mul(r23005, r23004, r23004, MPFR_RNDN); mpfr_sub(r23006, r23004, r23000, MPFR_RNDN); mpfr_mul(r23007, r23000, r23006, MPFR_RNDN); mpfr_sub(r23008, r23005, r23007, MPFR_RNDN); mpfr_div(r23009, r23002, r23008, MPFR_RNDN); mpfr_set_d(r23010, t, MPFR_RNDN); mpfr_sub(r23011, r23009, r23010, MPFR_RNDN); mpfr_set_d(r23012, a, MPFR_RNDN); ; mpfr_sub(r23014, r23012, r23013, MPFR_RNDN); mpfr_log(r23015, r23010, MPFR_RNDN); mpfr_mul(r23016, r23014, r23015, MPFR_RNDN); mpfr_add(r23017, r23011, r23016, MPFR_RNDN); return mpfr_get_d(r23017, MPFR_RNDN); } static mpfr_t r23018, r23019, r23020, r23021, r23022, r23023, r23024, r23025, r23026, r23027, r23028, r23029, r23030, r23031, r23032, r23033, r23034, r23035, r23036, r23037, r23038, r23039, r23040, r23041, r23042; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r23018); mpfr_init(r23019); mpfr_init(r23020); mpfr_init(r23021); mpfr_init_set_str(r23022, "3", 10, MPFR_RNDN); mpfr_init(r23023); mpfr_init(r23024); mpfr_init(r23025); mpfr_init(r23026); mpfr_init(r23027); mpfr_init(r23028); mpfr_init(r23029); mpfr_init(r23030); mpfr_init(r23031); mpfr_init(r23032); mpfr_init(r23033); mpfr_init(r23034); mpfr_init(r23035); mpfr_init(r23036); mpfr_init(r23037); mpfr_init_set_str(r23038, "0.5", 10, MPFR_RNDN); mpfr_init(r23039); mpfr_init(r23040); mpfr_init(r23041); mpfr_init(r23042); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r23018, x, MPFR_RNDN); mpfr_set_d(r23019, y, MPFR_RNDN); mpfr_add(r23020, r23018, r23019, MPFR_RNDN); mpfr_log(r23021, r23020, MPFR_RNDN); ; mpfr_pow(r23023, r23021, r23022, MPFR_RNDN); mpfr_set_d(r23024, z, MPFR_RNDN); mpfr_log(r23025, r23024, MPFR_RNDN); mpfr_pow(r23026, r23025, r23022, MPFR_RNDN); mpfr_add(r23027, r23023, r23026, MPFR_RNDN); mpfr_add(r23028, r23019, r23018, MPFR_RNDN); mpfr_log(r23029, r23028, MPFR_RNDN); mpfr_mul(r23030, r23029, r23029, MPFR_RNDN); mpfr_sub(r23031, r23029, r23025, MPFR_RNDN); mpfr_mul(r23032, r23025, r23031, MPFR_RNDN); mpfr_sub(r23033, r23030, r23032, MPFR_RNDN); mpfr_div(r23034, r23027, r23033, MPFR_RNDN); mpfr_set_d(r23035, t, MPFR_RNDN); mpfr_sub(r23036, r23034, r23035, MPFR_RNDN); mpfr_set_d(r23037, a, MPFR_RNDN); ; mpfr_sub(r23039, r23037, r23038, MPFR_RNDN); mpfr_log(r23040, r23035, MPFR_RNDN); mpfr_mul(r23041, r23039, r23040, MPFR_RNDN); mpfr_add(r23042, r23036, r23041, MPFR_RNDN); return mpfr_get_d(r23042, MPFR_RNDN); }