#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 r22873 = x; float r22874 = y; float r22875 = r22873 + r22874; float r22876 = log(r22875); float r22877 = z; float r22878 = log(r22877); float r22879 = r22876 + r22878; float r22880 = t; float r22881 = r22879 - r22880; float r22882 = a; float r22883 = 0.5; float r22884 = r22882 - r22883; float r22885 = log(r22880); float r22886 = r22884 * r22885; float r22887 = r22881 + r22886; return r22887; } double f_id(double x, double y, double z, double t, double a) { double r22888 = x; double r22889 = y; double r22890 = r22888 + r22889; double r22891 = log(r22890); double r22892 = z; double r22893 = log(r22892); double r22894 = r22891 + r22893; double r22895 = t; double r22896 = r22894 - r22895; double r22897 = a; double r22898 = 0.5; double r22899 = r22897 - r22898; double r22900 = log(r22895); double r22901 = r22899 * r22900; double r22902 = r22896 + r22901; return r22902; } double f_of(float x, float y, float z, float t, float a) { float r22903 = x; float r22904 = y; float r22905 = r22903 + r22904; float r22906 = log(r22905); float r22907 = 3; float r22908 = pow(r22906, r22907); float r22909 = z; float r22910 = log(r22909); float r22911 = pow(r22910, r22907); float r22912 = r22908 + r22911; float r22913 = r22904 + r22903; float r22914 = log(r22913); float r22915 = r22914 - r22910; float r22916 = r22910 * r22910; float r22917 = fma(r22915, r22914, r22916); float r22918 = r22912 / r22917; float r22919 = t; float r22920 = r22918 - r22919; float r22921 = a; float r22922 = 0.5; float r22923 = r22921 - r22922; float r22924 = log(r22919); float r22925 = r22923 * r22924; float r22926 = r22920 + r22925; return r22926; } double f_od(double x, double y, double z, double t, double a) { double r22927 = x; double r22928 = y; double r22929 = r22927 + r22928; double r22930 = log(r22929); double r22931 = 3; double r22932 = pow(r22930, r22931); double r22933 = z; double r22934 = log(r22933); double r22935 = pow(r22934, r22931); double r22936 = r22932 + r22935; double r22937 = r22928 + r22927; double r22938 = log(r22937); double r22939 = r22938 - r22934; double r22940 = r22934 * r22934; double r22941 = fma(r22939, r22938, r22940); double r22942 = r22936 / r22941; double r22943 = t; double r22944 = r22942 - r22943; double r22945 = a; double r22946 = 0.5; double r22947 = r22945 - r22946; double r22948 = log(r22943); double r22949 = r22947 * r22948; double r22950 = r22944 + r22949; return r22950; } 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 r22951, r22952, r22953, r22954, r22955, r22956, r22957, r22958, r22959, r22960, r22961, r22962, r22963, r22964, r22965; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r22951); mpfr_init(r22952); mpfr_init(r22953); mpfr_init(r22954); mpfr_init(r22955); mpfr_init(r22956); mpfr_init(r22957); mpfr_init(r22958); mpfr_init(r22959); mpfr_init(r22960); mpfr_init_set_str(r22961, "0.5", 10, MPFR_RNDN); mpfr_init(r22962); mpfr_init(r22963); mpfr_init(r22964); mpfr_init(r22965); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r22951, x, MPFR_RNDN); mpfr_set_d(r22952, y, MPFR_RNDN); mpfr_add(r22953, r22951, r22952, MPFR_RNDN); mpfr_log(r22954, r22953, MPFR_RNDN); mpfr_set_d(r22955, z, MPFR_RNDN); mpfr_log(r22956, r22955, MPFR_RNDN); mpfr_add(r22957, r22954, r22956, MPFR_RNDN); mpfr_set_d(r22958, t, MPFR_RNDN); mpfr_sub(r22959, r22957, r22958, MPFR_RNDN); mpfr_set_d(r22960, a, MPFR_RNDN); ; mpfr_sub(r22962, r22960, r22961, MPFR_RNDN); mpfr_log(r22963, r22958, MPFR_RNDN); mpfr_mul(r22964, r22962, r22963, MPFR_RNDN); mpfr_add(r22965, r22959, r22964, MPFR_RNDN); return mpfr_get_d(r22965, MPFR_RNDN); } static mpfr_t r22966, r22967, r22968, r22969, r22970, r22971, r22972, r22973, r22974, r22975, r22976, r22977, r22978, r22979, r22980, r22981, r22982, r22983, r22984, r22985, r22986, r22987, r22988, r22989; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r22966); mpfr_init(r22967); mpfr_init(r22968); mpfr_init(r22969); mpfr_init_set_str(r22970, "3", 10, MPFR_RNDN); 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(r22980); mpfr_init(r22981); mpfr_init(r22982); mpfr_init(r22983); mpfr_init(r22984); mpfr_init_set_str(r22985, "0.5", 10, MPFR_RNDN); mpfr_init(r22986); mpfr_init(r22987); mpfr_init(r22988); mpfr_init(r22989); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r22966, x, MPFR_RNDN); mpfr_set_d(r22967, y, MPFR_RNDN); mpfr_add(r22968, r22966, r22967, MPFR_RNDN); mpfr_log(r22969, r22968, MPFR_RNDN); ; mpfr_pow(r22971, r22969, r22970, MPFR_RNDN); mpfr_set_d(r22972, z, MPFR_RNDN); mpfr_log(r22973, r22972, MPFR_RNDN); mpfr_pow(r22974, r22973, r22970, MPFR_RNDN); mpfr_add(r22975, r22971, r22974, MPFR_RNDN); mpfr_add(r22976, r22967, r22966, MPFR_RNDN); mpfr_log(r22977, r22976, MPFR_RNDN); mpfr_sub(r22978, r22977, r22973, MPFR_RNDN); mpfr_mul(r22979, r22973, r22973, MPFR_RNDN); mpfr_fma(r22980, r22978, r22977, r22979, MPFR_RNDN); mpfr_div(r22981, r22975, r22980, MPFR_RNDN); mpfr_set_d(r22982, t, MPFR_RNDN); mpfr_sub(r22983, r22981, r22982, MPFR_RNDN); mpfr_set_d(r22984, a, MPFR_RNDN); ; mpfr_sub(r22986, r22984, r22985, MPFR_RNDN); mpfr_log(r22987, r22982, MPFR_RNDN); mpfr_mul(r22988, r22986, r22987, MPFR_RNDN); mpfr_add(r22989, r22983, r22988, MPFR_RNDN); return mpfr_get_d(r22989, MPFR_RNDN); } static mpfr_t r22990, r22991, r22992, r22993, r22994, r22995, r22996, r22997, r22998, r22999, r23000, r23001, r23002, r23003, r23004, r23005, r23006, r23007, r23008, r23009, r23010, r23011, r23012, r23013; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r22990); mpfr_init(r22991); mpfr_init(r22992); mpfr_init(r22993); mpfr_init_set_str(r22994, "3", 10, MPFR_RNDN); mpfr_init(r22995); mpfr_init(r22996); mpfr_init(r22997); 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_set_str(r23009, "0.5", 10, MPFR_RNDN); mpfr_init(r23010); mpfr_init(r23011); mpfr_init(r23012); mpfr_init(r23013); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r22990, x, MPFR_RNDN); mpfr_set_d(r22991, y, MPFR_RNDN); mpfr_add(r22992, r22990, r22991, MPFR_RNDN); mpfr_log(r22993, r22992, MPFR_RNDN); ; mpfr_pow(r22995, r22993, r22994, MPFR_RNDN); mpfr_set_d(r22996, z, MPFR_RNDN); mpfr_log(r22997, r22996, MPFR_RNDN); mpfr_pow(r22998, r22997, r22994, MPFR_RNDN); mpfr_add(r22999, r22995, r22998, MPFR_RNDN); mpfr_add(r23000, r22991, r22990, MPFR_RNDN); mpfr_log(r23001, r23000, MPFR_RNDN); mpfr_sub(r23002, r23001, r22997, MPFR_RNDN); mpfr_mul(r23003, r22997, r22997, MPFR_RNDN); mpfr_fma(r23004, r23002, r23001, r23003, MPFR_RNDN); mpfr_div(r23005, r22999, r23004, MPFR_RNDN); mpfr_set_d(r23006, t, MPFR_RNDN); mpfr_sub(r23007, r23005, r23006, MPFR_RNDN); mpfr_set_d(r23008, a, MPFR_RNDN); ; mpfr_sub(r23010, r23008, r23009, MPFR_RNDN); mpfr_log(r23011, r23006, MPFR_RNDN); mpfr_mul(r23012, r23010, r23011, MPFR_RNDN); mpfr_add(r23013, r23007, r23012, MPFR_RNDN); return mpfr_get_d(r23013, MPFR_RNDN); }