#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 r21907 = x; float r21908 = y; float r21909 = r21907 + r21908; float r21910 = log(r21909); float r21911 = z; float r21912 = log(r21911); float r21913 = r21910 + r21912; float r21914 = t; float r21915 = r21913 - r21914; float r21916 = a; float r21917 = 0.5; float r21918 = r21916 - r21917; float r21919 = log(r21914); float r21920 = r21918 * r21919; float r21921 = r21915 + r21920; return r21921; } double f_id(double x, double y, double z, double t, double a) { double r21922 = x; double r21923 = y; double r21924 = r21922 + r21923; double r21925 = log(r21924); double r21926 = z; double r21927 = log(r21926); double r21928 = r21925 + r21927; double r21929 = t; double r21930 = r21928 - r21929; double r21931 = a; double r21932 = 0.5; double r21933 = r21931 - r21932; double r21934 = log(r21929); double r21935 = r21933 * r21934; double r21936 = r21930 + r21935; return r21936; } double f_of(float x, float y, float z, float t, float a) { float r21937 = x; float r21938 = y; float r21939 = r21937 + r21938; float r21940 = log(r21939); float r21941 = 3; float r21942 = pow(r21940, r21941); float r21943 = z; float r21944 = log(r21943); float r21945 = pow(r21944, r21941); float r21946 = r21942 + r21945; float r21947 = r21938 + r21937; float r21948 = log(r21947); float r21949 = r21944 - r21948; float r21950 = r21949 * r21944; float r21951 = r21948 * r21948; float r21952 = r21950 + r21951; float r21953 = r21946 / r21952; float r21954 = t; float r21955 = r21953 - r21954; float r21956 = a; float r21957 = 0.5; float r21958 = r21956 - r21957; float r21959 = log(r21954); float r21960 = r21958 * r21959; float r21961 = r21955 + r21960; return r21961; } double f_od(double x, double y, double z, double t, double a) { double r21962 = x; double r21963 = y; double r21964 = r21962 + r21963; double r21965 = log(r21964); double r21966 = 3; double r21967 = pow(r21965, r21966); double r21968 = z; double r21969 = log(r21968); double r21970 = pow(r21969, r21966); double r21971 = r21967 + r21970; double r21972 = r21963 + r21962; double r21973 = log(r21972); double r21974 = r21969 - r21973; double r21975 = r21974 * r21969; double r21976 = r21973 * r21973; double r21977 = r21975 + r21976; double r21978 = r21971 / r21977; double r21979 = t; double r21980 = r21978 - r21979; double r21981 = a; double r21982 = 0.5; double r21983 = r21981 - r21982; double r21984 = log(r21979); double r21985 = r21983 * r21984; double r21986 = r21980 + r21985; return r21986; } 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 r21987, r21988, r21989, r21990, r21991, r21992, r21993, r21994, r21995, r21996, r21997, r21998, r21999, r22000, r22001; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r21987); mpfr_init(r21988); mpfr_init(r21989); mpfr_init(r21990); mpfr_init(r21991); mpfr_init(r21992); mpfr_init(r21993); mpfr_init(r21994); mpfr_init(r21995); mpfr_init(r21996); mpfr_init_set_str(r21997, "0.5", 10, MPFR_RNDN); mpfr_init(r21998); mpfr_init(r21999); mpfr_init(r22000); mpfr_init(r22001); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r21987, x, MPFR_RNDN); mpfr_set_d(r21988, y, MPFR_RNDN); mpfr_add(r21989, r21987, r21988, MPFR_RNDN); mpfr_log(r21990, r21989, MPFR_RNDN); mpfr_set_d(r21991, z, MPFR_RNDN); mpfr_log(r21992, r21991, MPFR_RNDN); mpfr_add(r21993, r21990, r21992, MPFR_RNDN); mpfr_set_d(r21994, t, MPFR_RNDN); mpfr_sub(r21995, r21993, r21994, MPFR_RNDN); mpfr_set_d(r21996, a, MPFR_RNDN); ; mpfr_sub(r21998, r21996, r21997, MPFR_RNDN); mpfr_log(r21999, r21994, MPFR_RNDN); mpfr_mul(r22000, r21998, r21999, MPFR_RNDN); mpfr_add(r22001, r21995, r22000, MPFR_RNDN); return mpfr_get_d(r22001, MPFR_RNDN); } static mpfr_t r22002, r22003, r22004, r22005, r22006, r22007, r22008, r22009, r22010, r22011, r22012, r22013, r22014, r22015, r22016, r22017, r22018, r22019, r22020, r22021, r22022, r22023, r22024, r22025, r22026; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22002); mpfr_init(r22003); mpfr_init(r22004); mpfr_init(r22005); mpfr_init_set_str(r22006, "3", 10, MPFR_RNDN); mpfr_init(r22007); mpfr_init(r22008); mpfr_init(r22009); mpfr_init(r22010); mpfr_init(r22011); mpfr_init(r22012); mpfr_init(r22013); mpfr_init(r22014); mpfr_init(r22015); mpfr_init(r22016); mpfr_init(r22017); mpfr_init(r22018); mpfr_init(r22019); mpfr_init(r22020); mpfr_init(r22021); mpfr_init_set_str(r22022, "0.5", 10, MPFR_RNDN); mpfr_init(r22023); mpfr_init(r22024); mpfr_init(r22025); mpfr_init(r22026); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r22002, x, MPFR_RNDN); mpfr_set_d(r22003, y, MPFR_RNDN); mpfr_add(r22004, r22002, r22003, MPFR_RNDN); mpfr_log(r22005, r22004, MPFR_RNDN); ; mpfr_pow(r22007, r22005, r22006, MPFR_RNDN); mpfr_set_d(r22008, z, MPFR_RNDN); mpfr_log(r22009, r22008, MPFR_RNDN); mpfr_pow(r22010, r22009, r22006, MPFR_RNDN); mpfr_add(r22011, r22007, r22010, MPFR_RNDN); mpfr_add(r22012, r22003, r22002, MPFR_RNDN); mpfr_log(r22013, r22012, MPFR_RNDN); mpfr_sub(r22014, r22009, r22013, MPFR_RNDN); mpfr_mul(r22015, r22014, r22009, MPFR_RNDN); mpfr_mul(r22016, r22013, r22013, MPFR_RNDN); mpfr_add(r22017, r22015, r22016, MPFR_RNDN); mpfr_div(r22018, r22011, r22017, MPFR_RNDN); mpfr_set_d(r22019, t, MPFR_RNDN); mpfr_sub(r22020, r22018, r22019, MPFR_RNDN); mpfr_set_d(r22021, a, MPFR_RNDN); ; mpfr_sub(r22023, r22021, r22022, MPFR_RNDN); mpfr_log(r22024, r22019, MPFR_RNDN); mpfr_mul(r22025, r22023, r22024, MPFR_RNDN); mpfr_add(r22026, r22020, r22025, MPFR_RNDN); return mpfr_get_d(r22026, MPFR_RNDN); } static mpfr_t r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036, r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22027); mpfr_init(r22028); mpfr_init(r22029); mpfr_init(r22030); mpfr_init_set_str(r22031, "3", 10, MPFR_RNDN); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); mpfr_init(r22036); mpfr_init(r22037); mpfr_init(r22038); mpfr_init(r22039); mpfr_init(r22040); mpfr_init(r22041); mpfr_init(r22042); mpfr_init(r22043); mpfr_init(r22044); mpfr_init(r22045); mpfr_init(r22046); mpfr_init_set_str(r22047, "0.5", 10, MPFR_RNDN); mpfr_init(r22048); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r22027, x, MPFR_RNDN); mpfr_set_d(r22028, y, MPFR_RNDN); mpfr_add(r22029, r22027, r22028, MPFR_RNDN); mpfr_log(r22030, r22029, MPFR_RNDN); ; mpfr_pow(r22032, r22030, r22031, MPFR_RNDN); mpfr_set_d(r22033, z, MPFR_RNDN); mpfr_log(r22034, r22033, MPFR_RNDN); mpfr_pow(r22035, r22034, r22031, MPFR_RNDN); mpfr_add(r22036, r22032, r22035, MPFR_RNDN); mpfr_add(r22037, r22028, r22027, MPFR_RNDN); mpfr_log(r22038, r22037, MPFR_RNDN); mpfr_sub(r22039, r22034, r22038, MPFR_RNDN); mpfr_mul(r22040, r22039, r22034, MPFR_RNDN); mpfr_mul(r22041, r22038, r22038, MPFR_RNDN); mpfr_add(r22042, r22040, r22041, MPFR_RNDN); mpfr_div(r22043, r22036, r22042, MPFR_RNDN); mpfr_set_d(r22044, t, MPFR_RNDN); mpfr_sub(r22045, r22043, r22044, MPFR_RNDN); mpfr_set_d(r22046, a, MPFR_RNDN); ; mpfr_sub(r22048, r22046, r22047, MPFR_RNDN); mpfr_log(r22049, r22044, MPFR_RNDN); mpfr_mul(r22050, r22048, r22049, MPFR_RNDN); mpfr_add(r22051, r22045, r22050, MPFR_RNDN); return mpfr_get_d(r22051, MPFR_RNDN); }