#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 r22067 = x; float r22068 = y; float r22069 = r22067 + r22068; float r22070 = log(r22069); float r22071 = z; float r22072 = log(r22071); float r22073 = r22070 + r22072; float r22074 = t; float r22075 = r22073 - r22074; float r22076 = a; float r22077 = 0.5; float r22078 = r22076 - r22077; float r22079 = log(r22074); float r22080 = r22078 * r22079; float r22081 = r22075 + r22080; return r22081; } double f_id(double x, double y, double z, double t, double a) { double r22082 = x; double r22083 = y; double r22084 = r22082 + r22083; double r22085 = log(r22084); double r22086 = z; double r22087 = log(r22086); double r22088 = r22085 + r22087; double r22089 = t; double r22090 = r22088 - r22089; double r22091 = a; double r22092 = 0.5; double r22093 = r22091 - r22092; double r22094 = log(r22089); double r22095 = r22093 * r22094; double r22096 = r22090 + r22095; return r22096; } double f_of(float x, float y, float z, float t, float a) { float r22097 = a; float r22098 = 0.5; float r22099 = r22097 - r22098; float r22100 = t; float r22101 = log(r22100); float r22102 = y; float r22103 = x; float r22104 = r22102 + r22103; float r22105 = log(r22104); float r22106 = cbrt(r22105); float r22107 = r22106 * r22106; float r22108 = z; float r22109 = log(r22108); float r22110 = r22109 - r22100; float r22111 = fma(r22107, r22106, r22110); float r22112 = fma(r22099, r22101, r22111); return r22112; } double f_od(double x, double y, double z, double t, double a) { double r22113 = a; double r22114 = 0.5; double r22115 = r22113 - r22114; double r22116 = t; double r22117 = log(r22116); double r22118 = y; double r22119 = x; double r22120 = r22118 + r22119; double r22121 = log(r22120); double r22122 = cbrt(r22121); double r22123 = r22122 * r22122; double r22124 = z; double r22125 = log(r22124); double r22126 = r22125 - r22116; double r22127 = fma(r22123, r22122, r22126); double r22128 = fma(r22115, r22117, r22127); return r22128; } 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 r22129, r22130, r22131, r22132, r22133, r22134, r22135, r22136, r22137, r22138, r22139, r22140, r22141, r22142, r22143; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r22129); mpfr_init(r22130); mpfr_init(r22131); mpfr_init(r22132); mpfr_init(r22133); mpfr_init(r22134); mpfr_init(r22135); mpfr_init(r22136); mpfr_init(r22137); mpfr_init(r22138); mpfr_init_set_str(r22139, "0.5", 10, MPFR_RNDN); mpfr_init(r22140); mpfr_init(r22141); mpfr_init(r22142); mpfr_init(r22143); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r22129, x, MPFR_RNDN); mpfr_set_d(r22130, y, MPFR_RNDN); mpfr_add(r22131, r22129, r22130, MPFR_RNDN); mpfr_log(r22132, r22131, MPFR_RNDN); mpfr_set_d(r22133, z, MPFR_RNDN); mpfr_log(r22134, r22133, MPFR_RNDN); mpfr_add(r22135, r22132, r22134, MPFR_RNDN); mpfr_set_d(r22136, t, MPFR_RNDN); mpfr_sub(r22137, r22135, r22136, MPFR_RNDN); mpfr_set_d(r22138, a, MPFR_RNDN); ; mpfr_sub(r22140, r22138, r22139, MPFR_RNDN); mpfr_log(r22141, r22136, MPFR_RNDN); mpfr_mul(r22142, r22140, r22141, MPFR_RNDN); mpfr_add(r22143, r22137, r22142, MPFR_RNDN); return mpfr_get_d(r22143, MPFR_RNDN); } static mpfr_t r22144, r22145, r22146, r22147, r22148, r22149, r22150, r22151, r22152, r22153, r22154, r22155, r22156, r22157, r22158, r22159; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r22144); mpfr_init_set_str(r22145, "0.5", 10, MPFR_RNDN); mpfr_init(r22146); mpfr_init(r22147); mpfr_init(r22148); mpfr_init(r22149); mpfr_init(r22150); mpfr_init(r22151); mpfr_init(r22152); mpfr_init(r22153); mpfr_init(r22154); mpfr_init(r22155); mpfr_init(r22156); mpfr_init(r22157); mpfr_init(r22158); mpfr_init(r22159); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r22144, a, MPFR_RNDN); ; mpfr_sub(r22146, r22144, r22145, MPFR_RNDN); mpfr_set_d(r22147, t, MPFR_RNDN); mpfr_log(r22148, r22147, MPFR_RNDN); mpfr_set_d(r22149, y, MPFR_RNDN); mpfr_set_d(r22150, x, MPFR_RNDN); mpfr_add(r22151, r22149, r22150, MPFR_RNDN); mpfr_log(r22152, r22151, MPFR_RNDN); mpfr_cbrt(r22153, r22152, MPFR_RNDN); mpfr_mul(r22154, r22153, r22153, MPFR_RNDN); mpfr_set_d(r22155, z, MPFR_RNDN); mpfr_log(r22156, r22155, MPFR_RNDN); mpfr_sub(r22157, r22156, r22147, MPFR_RNDN); mpfr_fma(r22158, r22154, r22153, r22157, MPFR_RNDN); mpfr_fma(r22159, r22146, r22148, r22158, MPFR_RNDN); return mpfr_get_d(r22159, MPFR_RNDN); } static mpfr_t r22160, r22161, r22162, r22163, r22164, r22165, r22166, r22167, r22168, r22169, r22170, r22171, r22172, r22173, r22174, r22175; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r22160); mpfr_init_set_str(r22161, "0.5", 10, MPFR_RNDN); mpfr_init(r22162); mpfr_init(r22163); mpfr_init(r22164); mpfr_init(r22165); mpfr_init(r22166); mpfr_init(r22167); mpfr_init(r22168); mpfr_init(r22169); mpfr_init(r22170); mpfr_init(r22171); mpfr_init(r22172); mpfr_init(r22173); mpfr_init(r22174); mpfr_init(r22175); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r22160, a, MPFR_RNDN); ; mpfr_sub(r22162, r22160, r22161, MPFR_RNDN); mpfr_set_d(r22163, t, MPFR_RNDN); mpfr_log(r22164, r22163, MPFR_RNDN); mpfr_set_d(r22165, y, MPFR_RNDN); mpfr_set_d(r22166, x, MPFR_RNDN); mpfr_add(r22167, r22165, r22166, MPFR_RNDN); mpfr_log(r22168, r22167, MPFR_RNDN); mpfr_cbrt(r22169, r22168, MPFR_RNDN); mpfr_mul(r22170, r22169, r22169, MPFR_RNDN); mpfr_set_d(r22171, z, MPFR_RNDN); mpfr_log(r22172, r22171, MPFR_RNDN); mpfr_sub(r22173, r22172, r22163, MPFR_RNDN); mpfr_fma(r22174, r22170, r22169, r22173, MPFR_RNDN); mpfr_fma(r22175, r22162, r22164, r22174, MPFR_RNDN); return mpfr_get_d(r22175, MPFR_RNDN); }