#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 r25051 = x; float r25052 = y; float r25053 = r25051 + r25052; float r25054 = log(r25053); float r25055 = z; float r25056 = log(r25055); float r25057 = r25054 + r25056; float r25058 = t; float r25059 = r25057 - r25058; float r25060 = a; float r25061 = 0.5; float r25062 = r25060 - r25061; float r25063 = log(r25058); float r25064 = r25062 * r25063; float r25065 = r25059 + r25064; return r25065; } double f_id(double x, double y, double z, double t, double a) { double r25066 = x; double r25067 = y; double r25068 = r25066 + r25067; double r25069 = log(r25068); double r25070 = z; double r25071 = log(r25070); double r25072 = r25069 + r25071; double r25073 = t; double r25074 = r25072 - r25073; double r25075 = a; double r25076 = 0.5; double r25077 = r25075 - r25076; double r25078 = log(r25073); double r25079 = r25077 * r25078; double r25080 = r25074 + r25079; return r25080; } double f_of(float x, float y, float z, float t, float a) { float r25081 = y; float r25082 = x; float r25083 = r25081 + r25082; float r25084 = log(r25083); float r25085 = t; float r25086 = a; float r25087 = 0.5; float r25088 = r25086 - r25087; float r25089 = log(r25085); float r25090 = z; float r25091 = log(r25090); float r25092 = fma(r25088, r25089, r25091); float r25093 = r25085 - r25092; float r25094 = r25084 - r25093; return r25094; } double f_od(double x, double y, double z, double t, double a) { double r25095 = y; double r25096 = x; double r25097 = r25095 + r25096; double r25098 = log(r25097); double r25099 = t; double r25100 = a; double r25101 = 0.5; double r25102 = r25100 - r25101; double r25103 = log(r25099); double r25104 = z; double r25105 = log(r25104); double r25106 = fma(r25102, r25103, r25105); double r25107 = r25099 - r25106; double r25108 = r25098 - r25107; return r25108; } 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 r25109, r25110, r25111, r25112, r25113, r25114, r25115, r25116, r25117, r25118, r25119, r25120, r25121, r25122, r25123; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25109); mpfr_init(r25110); mpfr_init(r25111); mpfr_init(r25112); mpfr_init(r25113); mpfr_init(r25114); mpfr_init(r25115); mpfr_init(r25116); mpfr_init(r25117); mpfr_init(r25118); mpfr_init_set_str(r25119, "0.5", 10, MPFR_RNDN); mpfr_init(r25120); mpfr_init(r25121); mpfr_init(r25122); mpfr_init(r25123); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r25109, x, MPFR_RNDN); mpfr_set_d(r25110, y, MPFR_RNDN); mpfr_add(r25111, r25109, r25110, MPFR_RNDN); mpfr_log(r25112, r25111, MPFR_RNDN); mpfr_set_d(r25113, z, MPFR_RNDN); mpfr_log(r25114, r25113, MPFR_RNDN); mpfr_add(r25115, r25112, r25114, MPFR_RNDN); mpfr_set_d(r25116, t, MPFR_RNDN); mpfr_sub(r25117, r25115, r25116, MPFR_RNDN); mpfr_set_d(r25118, a, MPFR_RNDN); ; mpfr_sub(r25120, r25118, r25119, MPFR_RNDN); mpfr_log(r25121, r25116, MPFR_RNDN); mpfr_mul(r25122, r25120, r25121, MPFR_RNDN); mpfr_add(r25123, r25117, r25122, MPFR_RNDN); return mpfr_get_d(r25123, MPFR_RNDN); } static mpfr_t r25124, r25125, r25126, r25127, r25128, r25129, r25130, r25131, r25132, r25133, r25134, r25135, r25136, r25137; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25124); mpfr_init(r25125); mpfr_init(r25126); mpfr_init(r25127); mpfr_init(r25128); mpfr_init(r25129); mpfr_init_set_str(r25130, "0.5", 10, MPFR_RNDN); mpfr_init(r25131); mpfr_init(r25132); mpfr_init(r25133); mpfr_init(r25134); mpfr_init(r25135); mpfr_init(r25136); mpfr_init(r25137); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r25124, y, MPFR_RNDN); mpfr_set_d(r25125, x, MPFR_RNDN); mpfr_add(r25126, r25124, r25125, MPFR_RNDN); mpfr_log(r25127, r25126, MPFR_RNDN); mpfr_set_d(r25128, t, MPFR_RNDN); mpfr_set_d(r25129, a, MPFR_RNDN); ; mpfr_sub(r25131, r25129, r25130, MPFR_RNDN); mpfr_log(r25132, r25128, MPFR_RNDN); mpfr_set_d(r25133, z, MPFR_RNDN); mpfr_log(r25134, r25133, MPFR_RNDN); mpfr_fma(r25135, r25131, r25132, r25134, MPFR_RNDN); mpfr_sub(r25136, r25128, r25135, MPFR_RNDN); mpfr_sub(r25137, r25127, r25136, MPFR_RNDN); return mpfr_get_d(r25137, MPFR_RNDN); } static mpfr_t r25138, r25139, r25140, r25141, r25142, r25143, r25144, r25145, r25146, r25147, r25148, r25149, r25150, r25151; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25138); mpfr_init(r25139); mpfr_init(r25140); mpfr_init(r25141); mpfr_init(r25142); mpfr_init(r25143); mpfr_init_set_str(r25144, "0.5", 10, MPFR_RNDN); mpfr_init(r25145); mpfr_init(r25146); mpfr_init(r25147); mpfr_init(r25148); mpfr_init(r25149); mpfr_init(r25150); mpfr_init(r25151); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r25138, y, MPFR_RNDN); mpfr_set_d(r25139, x, MPFR_RNDN); mpfr_add(r25140, r25138, r25139, MPFR_RNDN); mpfr_log(r25141, r25140, MPFR_RNDN); mpfr_set_d(r25142, t, MPFR_RNDN); mpfr_set_d(r25143, a, MPFR_RNDN); ; mpfr_sub(r25145, r25143, r25144, MPFR_RNDN); mpfr_log(r25146, r25142, MPFR_RNDN); mpfr_set_d(r25147, z, MPFR_RNDN); mpfr_log(r25148, r25147, MPFR_RNDN); mpfr_fma(r25149, r25145, r25146, r25148, MPFR_RNDN); mpfr_sub(r25150, r25142, r25149, MPFR_RNDN); mpfr_sub(r25151, r25141, r25150, MPFR_RNDN); return mpfr_get_d(r25151, MPFR_RNDN); }