#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 r25024 = x; float r25025 = y; float r25026 = r25024 + r25025; float r25027 = log(r25026); float r25028 = z; float r25029 = log(r25028); float r25030 = r25027 + r25029; float r25031 = t; float r25032 = r25030 - r25031; float r25033 = a; float r25034 = 0.5; float r25035 = r25033 - r25034; float r25036 = log(r25031); float r25037 = r25035 * r25036; float r25038 = r25032 + r25037; return r25038; } double f_id(double x, double y, double z, double t, double a) { double r25039 = x; double r25040 = y; double r25041 = r25039 + r25040; double r25042 = log(r25041); double r25043 = z; double r25044 = log(r25043); double r25045 = r25042 + r25044; double r25046 = t; double r25047 = r25045 - r25046; double r25048 = a; double r25049 = 0.5; double r25050 = r25048 - r25049; double r25051 = log(r25046); double r25052 = r25050 * r25051; double r25053 = r25047 + r25052; return r25053; } double f_of(float x, float y, float z, float t, float a) { float r25054 = x; float r25055 = y; float r25056 = r25054 + r25055; float r25057 = log(r25056); float r25058 = z; float r25059 = log(r25058); float r25060 = r25057 + r25059; float r25061 = t; float r25062 = r25060 - r25061; float r25063 = a; float r25064 = 0.5; float r25065 = r25063 - r25064; float r25066 = log(r25061); float r25067 = r25065 * r25066; float r25068 = r25062 + r25067; return r25068; } double f_od(double x, double y, double z, double t, double a) { double r25069 = x; double r25070 = y; double r25071 = r25069 + r25070; double r25072 = log(r25071); double r25073 = z; double r25074 = log(r25073); double r25075 = r25072 + r25074; double r25076 = t; double r25077 = r25075 - r25076; double r25078 = a; double r25079 = 0.5; double r25080 = r25078 - r25079; double r25081 = log(r25076); double r25082 = r25080 * r25081; double r25083 = r25077 + r25082; return r25083; } 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 r25084, r25085, r25086, r25087, r25088, r25089, r25090, r25091, r25092, r25093, r25094, r25095, r25096, r25097, r25098; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25084); mpfr_init(r25085); mpfr_init(r25086); mpfr_init(r25087); mpfr_init(r25088); mpfr_init(r25089); mpfr_init(r25090); mpfr_init(r25091); mpfr_init(r25092); mpfr_init(r25093); mpfr_init_set_str(r25094, "0.5", 10, MPFR_RNDN); mpfr_init(r25095); mpfr_init(r25096); mpfr_init(r25097); mpfr_init(r25098); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r25084, x, MPFR_RNDN); mpfr_set_d(r25085, y, MPFR_RNDN); mpfr_add(r25086, r25084, r25085, MPFR_RNDN); mpfr_log(r25087, r25086, MPFR_RNDN); mpfr_set_d(r25088, z, MPFR_RNDN); mpfr_log(r25089, r25088, MPFR_RNDN); mpfr_add(r25090, r25087, r25089, MPFR_RNDN); mpfr_set_d(r25091, t, MPFR_RNDN); mpfr_sub(r25092, r25090, r25091, MPFR_RNDN); mpfr_set_d(r25093, a, MPFR_RNDN); ; mpfr_sub(r25095, r25093, r25094, MPFR_RNDN); mpfr_log(r25096, r25091, MPFR_RNDN); mpfr_mul(r25097, r25095, r25096, MPFR_RNDN); mpfr_add(r25098, r25092, r25097, MPFR_RNDN); return mpfr_get_d(r25098, MPFR_RNDN); } static mpfr_t r25099, r25100, r25101, r25102, r25103, r25104, r25105, r25106, r25107, r25108, r25109, r25110, r25111, r25112, r25113; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25099); mpfr_init(r25100); mpfr_init(r25101); mpfr_init(r25102); mpfr_init(r25103); mpfr_init(r25104); mpfr_init(r25105); mpfr_init(r25106); mpfr_init(r25107); mpfr_init(r25108); mpfr_init_set_str(r25109, "0.5", 10, MPFR_RNDN); mpfr_init(r25110); mpfr_init(r25111); mpfr_init(r25112); mpfr_init(r25113); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r25099, x, MPFR_RNDN); mpfr_set_d(r25100, y, MPFR_RNDN); mpfr_add(r25101, r25099, r25100, MPFR_RNDN); mpfr_log(r25102, r25101, MPFR_RNDN); mpfr_set_d(r25103, z, MPFR_RNDN); mpfr_log(r25104, r25103, MPFR_RNDN); mpfr_add(r25105, r25102, r25104, MPFR_RNDN); mpfr_set_d(r25106, t, MPFR_RNDN); mpfr_sub(r25107, r25105, r25106, MPFR_RNDN); mpfr_set_d(r25108, a, MPFR_RNDN); ; mpfr_sub(r25110, r25108, r25109, MPFR_RNDN); mpfr_log(r25111, r25106, MPFR_RNDN); mpfr_mul(r25112, r25110, r25111, MPFR_RNDN); mpfr_add(r25113, r25107, r25112, MPFR_RNDN); return mpfr_get_d(r25113, MPFR_RNDN); } static mpfr_t r25114, r25115, r25116, r25117, r25118, r25119, r25120, r25121, r25122, r25123, r25124, r25125, r25126, r25127, r25128; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25114); mpfr_init(r25115); mpfr_init(r25116); mpfr_init(r25117); mpfr_init(r25118); mpfr_init(r25119); mpfr_init(r25120); mpfr_init(r25121); mpfr_init(r25122); mpfr_init(r25123); mpfr_init_set_str(r25124, "0.5", 10, MPFR_RNDN); mpfr_init(r25125); mpfr_init(r25126); mpfr_init(r25127); mpfr_init(r25128); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r25114, x, MPFR_RNDN); mpfr_set_d(r25115, y, MPFR_RNDN); mpfr_add(r25116, r25114, r25115, MPFR_RNDN); mpfr_log(r25117, r25116, MPFR_RNDN); mpfr_set_d(r25118, z, MPFR_RNDN); mpfr_log(r25119, r25118, MPFR_RNDN); mpfr_add(r25120, r25117, r25119, MPFR_RNDN); mpfr_set_d(r25121, t, MPFR_RNDN); mpfr_sub(r25122, r25120, r25121, MPFR_RNDN); mpfr_set_d(r25123, a, MPFR_RNDN); ; mpfr_sub(r25125, r25123, r25124, MPFR_RNDN); mpfr_log(r25126, r25121, MPFR_RNDN); mpfr_mul(r25127, r25125, r25126, MPFR_RNDN); mpfr_add(r25128, r25122, r25127, MPFR_RNDN); return mpfr_get_d(r25128, MPFR_RNDN); }