#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 r15957 = x; float r15958 = y; float r15959 = r15957 + r15958; float r15960 = log(r15959); float r15961 = z; float r15962 = log(r15961); float r15963 = r15960 + r15962; float r15964 = t; float r15965 = r15963 - r15964; float r15966 = a; float r15967 = 0.5f; float r15968 = r15966 - r15967; float r15969 = log(r15964); float r15970 = r15968 * r15969; float r15971 = r15965 + r15970; return r15971; } double f_id(double x, double y, double z, double t, double a) { double r15972 = x; double r15973 = y; double r15974 = r15972 + r15973; double r15975 = log(r15974); double r15976 = z; double r15977 = log(r15976); double r15978 = r15975 + r15977; double r15979 = t; double r15980 = r15978 - r15979; double r15981 = a; double r15982 = 0.5; double r15983 = r15981 - r15982; double r15984 = log(r15979); double r15985 = r15983 * r15984; double r15986 = r15980 + r15985; return r15986; } double f_of(float x, float y, float z, float t, float a) { float r15987 = a; float r15988 = 0.5f; float r15989 = r15987 - r15988; float r15990 = t; float r15991 = log(r15990); float r15992 = x; float r15993 = y; float r15994 = r15992 + r15993; float r15995 = log(r15994); float r15996 = fma(r15989, r15991, r15995); float r15997 = z; float r15998 = log(r15997); float r15999 = r15998 - r15990; float r16000 = r15996 + r15999; return r16000; } double f_od(double x, double y, double z, double t, double a) { double r16001 = a; double r16002 = 0.5; double r16003 = r16001 - r16002; double r16004 = t; double r16005 = log(r16004); double r16006 = x; double r16007 = y; double r16008 = r16006 + r16007; double r16009 = log(r16008); double r16010 = fma(r16003, r16005, r16009); double r16011 = z; double r16012 = log(r16011); double r16013 = r16012 - r16004; double r16014 = r16010 + r16013; return r16014; } 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 r16015, r16016, r16017, r16018, r16019, r16020, r16021, r16022, r16023, r16024, r16025, r16026, r16027, r16028, r16029; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16015); mpfr_init(r16016); mpfr_init(r16017); mpfr_init(r16018); mpfr_init(r16019); mpfr_init(r16020); mpfr_init(r16021); mpfr_init(r16022); mpfr_init(r16023); mpfr_init(r16024); mpfr_init_set_str(r16025, "0.5", 10, MPFR_RNDN); mpfr_init(r16026); mpfr_init(r16027); mpfr_init(r16028); mpfr_init(r16029); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r16015, x, MPFR_RNDN); mpfr_set_d(r16016, y, MPFR_RNDN); mpfr_add(r16017, r16015, r16016, MPFR_RNDN); mpfr_log(r16018, r16017, MPFR_RNDN); mpfr_set_d(r16019, z, MPFR_RNDN); mpfr_log(r16020, r16019, MPFR_RNDN); mpfr_add(r16021, r16018, r16020, MPFR_RNDN); mpfr_set_d(r16022, t, MPFR_RNDN); mpfr_sub(r16023, r16021, r16022, MPFR_RNDN); mpfr_set_d(r16024, a, MPFR_RNDN); ; mpfr_sub(r16026, r16024, r16025, MPFR_RNDN); mpfr_log(r16027, r16022, MPFR_RNDN); mpfr_mul(r16028, r16026, r16027, MPFR_RNDN); mpfr_add(r16029, r16023, r16028, MPFR_RNDN); return mpfr_get_d(r16029, MPFR_RNDN); } static mpfr_t r16030, r16031, r16032, r16033, r16034, r16035, r16036, r16037, r16038, r16039, r16040, r16041, r16042, r16043; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16030); mpfr_init_set_str(r16031, "0.5", 10, MPFR_RNDN); mpfr_init(r16032); mpfr_init(r16033); mpfr_init(r16034); mpfr_init(r16035); mpfr_init(r16036); mpfr_init(r16037); mpfr_init(r16038); mpfr_init(r16039); mpfr_init(r16040); mpfr_init(r16041); mpfr_init(r16042); mpfr_init(r16043); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r16030, a, MPFR_RNDN); ; mpfr_sub(r16032, r16030, r16031, MPFR_RNDN); mpfr_set_d(r16033, t, MPFR_RNDN); mpfr_log(r16034, r16033, MPFR_RNDN); mpfr_set_d(r16035, x, MPFR_RNDN); mpfr_set_d(r16036, y, MPFR_RNDN); mpfr_add(r16037, r16035, r16036, MPFR_RNDN); mpfr_log(r16038, r16037, MPFR_RNDN); mpfr_fma(r16039, r16032, r16034, r16038, MPFR_RNDN); mpfr_set_d(r16040, z, MPFR_RNDN); mpfr_log(r16041, r16040, MPFR_RNDN); mpfr_sub(r16042, r16041, r16033, MPFR_RNDN); mpfr_add(r16043, r16039, r16042, MPFR_RNDN); return mpfr_get_d(r16043, MPFR_RNDN); } static mpfr_t r16044, r16045, r16046, r16047, r16048, r16049, r16050, r16051, r16052, r16053, r16054, r16055, r16056, r16057; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16044); mpfr_init_set_str(r16045, "0.5", 10, MPFR_RNDN); mpfr_init(r16046); mpfr_init(r16047); mpfr_init(r16048); mpfr_init(r16049); mpfr_init(r16050); mpfr_init(r16051); mpfr_init(r16052); mpfr_init(r16053); mpfr_init(r16054); mpfr_init(r16055); mpfr_init(r16056); mpfr_init(r16057); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r16044, a, MPFR_RNDN); ; mpfr_sub(r16046, r16044, r16045, MPFR_RNDN); mpfr_set_d(r16047, t, MPFR_RNDN); mpfr_log(r16048, r16047, MPFR_RNDN); mpfr_set_d(r16049, x, MPFR_RNDN); mpfr_set_d(r16050, y, MPFR_RNDN); mpfr_add(r16051, r16049, r16050, MPFR_RNDN); mpfr_log(r16052, r16051, MPFR_RNDN); mpfr_fma(r16053, r16046, r16048, r16052, MPFR_RNDN); mpfr_set_d(r16054, z, MPFR_RNDN); mpfr_log(r16055, r16054, MPFR_RNDN); mpfr_sub(r16056, r16055, r16047, MPFR_RNDN); mpfr_add(r16057, r16053, r16056, MPFR_RNDN); return mpfr_get_d(r16057, MPFR_RNDN); }