#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 r19357 = x; float r19358 = y; float r19359 = r19357 + r19358; float r19360 = log(r19359); float r19361 = z; float r19362 = log(r19361); float r19363 = r19360 + r19362; float r19364 = t; float r19365 = r19363 - r19364; float r19366 = a; float r19367 = 0.5f; float r19368 = r19366 - r19367; float r19369 = log(r19364); float r19370 = r19368 * r19369; float r19371 = r19365 + r19370; return r19371; } double f_id(double x, double y, double z, double t, double a) { double r19372 = x; double r19373 = y; double r19374 = r19372 + r19373; double r19375 = log(r19374); double r19376 = z; double r19377 = log(r19376); double r19378 = r19375 + r19377; double r19379 = t; double r19380 = r19378 - r19379; double r19381 = a; double r19382 = 0.5; double r19383 = r19381 - r19382; double r19384 = log(r19379); double r19385 = r19383 * r19384; double r19386 = r19380 + r19385; return r19386; } double f_of(float x, float y, float z, float t, float a) { float r19387 = y; float r19388 = x; float r19389 = r19387 + r19388; float r19390 = log(r19389); float r19391 = r19390 * (r19390 * r19390); float r19392 = z; float r19393 = log(r19392); float r19394 = r19393 * (r19393 * r19393); float r19395 = r19391 + r19394; float r19396 = r19393 - r19390; float r19397 = r19393 * r19396; float r19398 = r19390 * r19390; float r19399 = r19397 + r19398; float r19400 = r19395 / r19399; float r19401 = t; float r19402 = r19400 - r19401; float r19403 = a; float r19404 = 0.5f; float r19405 = r19403 - r19404; float r19406 = log(r19401); float r19407 = r19405 * r19406; float r19408 = r19402 + r19407; return r19408; } double f_od(double x, double y, double z, double t, double a) { double r19409 = y; double r19410 = x; double r19411 = r19409 + r19410; double r19412 = log(r19411); double r19413 = r19412 * (r19412 * r19412); double r19414 = z; double r19415 = log(r19414); double r19416 = r19415 * (r19415 * r19415); double r19417 = r19413 + r19416; double r19418 = r19415 - r19412; double r19419 = r19415 * r19418; double r19420 = r19412 * r19412; double r19421 = r19419 + r19420; double r19422 = r19417 / r19421; double r19423 = t; double r19424 = r19422 - r19423; double r19425 = a; double r19426 = 0.5; double r19427 = r19425 - r19426; double r19428 = log(r19423); double r19429 = r19427 * r19428; double r19430 = r19424 + r19429; return r19430; } 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 r19431, r19432, r19433, r19434, r19435, r19436, r19437, r19438, r19439, r19440, r19441, r19442, r19443, r19444, r19445; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r19431); mpfr_init(r19432); mpfr_init(r19433); mpfr_init(r19434); mpfr_init(r19435); mpfr_init(r19436); mpfr_init(r19437); mpfr_init(r19438); mpfr_init(r19439); mpfr_init(r19440); mpfr_init_set_str(r19441, "0.5", 10, MPFR_RNDN); mpfr_init(r19442); mpfr_init(r19443); mpfr_init(r19444); mpfr_init(r19445); } double f_im(double x, double y, double z, double t, double a) { mpfr_set_d(r19431, x, MPFR_RNDN); mpfr_set_d(r19432, y, MPFR_RNDN); mpfr_add(r19433, r19431, r19432, MPFR_RNDN); mpfr_log(r19434, r19433, MPFR_RNDN); mpfr_set_d(r19435, z, MPFR_RNDN); mpfr_log(r19436, r19435, MPFR_RNDN); mpfr_add(r19437, r19434, r19436, MPFR_RNDN); mpfr_set_d(r19438, t, MPFR_RNDN); mpfr_sub(r19439, r19437, r19438, MPFR_RNDN); mpfr_set_d(r19440, a, MPFR_RNDN); ; mpfr_sub(r19442, r19440, r19441, MPFR_RNDN); mpfr_log(r19443, r19438, MPFR_RNDN); mpfr_mul(r19444, r19442, r19443, MPFR_RNDN); mpfr_add(r19445, r19439, r19444, MPFR_RNDN); return mpfr_get_d(r19445, MPFR_RNDN); } static mpfr_t r19446, r19447, r19448, r19449, r19450, r19451, r19452, r19453, r19454, r19455, r19456, r19457, r19458, r19459, r19460, r19461, r19462, r19463, r19464, r19465, r19466, r19467; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r19446); mpfr_init(r19447); mpfr_init(r19448); mpfr_init(r19449); mpfr_init(r19450); mpfr_init(r19451); mpfr_init(r19452); mpfr_init(r19453); mpfr_init(r19454); mpfr_init(r19455); mpfr_init(r19456); mpfr_init(r19457); mpfr_init(r19458); mpfr_init(r19459); mpfr_init(r19460); mpfr_init(r19461); mpfr_init(r19462); mpfr_init_set_str(r19463, "0.5", 10, MPFR_RNDN); mpfr_init(r19464); mpfr_init(r19465); mpfr_init(r19466); mpfr_init(r19467); } double f_fm(double x, double y, double z, double t, double a) { mpfr_set_d(r19446, y, MPFR_RNDN); mpfr_set_d(r19447, x, MPFR_RNDN); mpfr_add(r19448, r19446, r19447, MPFR_RNDN); mpfr_log(r19449, r19448, MPFR_RNDN); mpfr_mul(r19450, r19449, r19449, MPFR_RNDN); mpfr_mul(r19450, r19450, r19449, MPFR_RNDN); mpfr_set_d(r19451, z, MPFR_RNDN); mpfr_log(r19452, r19451, MPFR_RNDN); mpfr_mul(r19453, r19452, r19452, MPFR_RNDN); mpfr_mul(r19453, r19453, r19452, MPFR_RNDN); mpfr_add(r19454, r19450, r19453, MPFR_RNDN); mpfr_sub(r19455, r19452, r19449, MPFR_RNDN); mpfr_mul(r19456, r19452, r19455, MPFR_RNDN); mpfr_sqr(r19457, r19449, MPFR_RNDN); mpfr_add(r19458, r19456, r19457, MPFR_RNDN); mpfr_div(r19459, r19454, r19458, MPFR_RNDN); mpfr_set_d(r19460, t, MPFR_RNDN); mpfr_sub(r19461, r19459, r19460, MPFR_RNDN); mpfr_set_d(r19462, a, MPFR_RNDN); ; mpfr_sub(r19464, r19462, r19463, MPFR_RNDN); mpfr_log(r19465, r19460, MPFR_RNDN); mpfr_mul(r19466, r19464, r19465, MPFR_RNDN); mpfr_add(r19467, r19461, r19466, MPFR_RNDN); return mpfr_get_d(r19467, MPFR_RNDN); } static mpfr_t r19468, r19469, r19470, r19471, r19472, r19473, r19474, r19475, r19476, r19477, r19478, r19479, r19480, r19481, r19482, r19483, r19484, r19485, r19486, r19487, r19488, r19489; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r19468); mpfr_init(r19469); mpfr_init(r19470); mpfr_init(r19471); mpfr_init(r19472); mpfr_init(r19473); mpfr_init(r19474); mpfr_init(r19475); mpfr_init(r19476); mpfr_init(r19477); mpfr_init(r19478); mpfr_init(r19479); mpfr_init(r19480); mpfr_init(r19481); mpfr_init(r19482); mpfr_init(r19483); mpfr_init(r19484); mpfr_init_set_str(r19485, "0.5", 10, MPFR_RNDN); mpfr_init(r19486); mpfr_init(r19487); mpfr_init(r19488); mpfr_init(r19489); } double f_dm(double x, double y, double z, double t, double a) { mpfr_set_d(r19468, y, MPFR_RNDN); mpfr_set_d(r19469, x, MPFR_RNDN); mpfr_add(r19470, r19468, r19469, MPFR_RNDN); mpfr_log(r19471, r19470, MPFR_RNDN); mpfr_mul(r19472, r19471, r19471, MPFR_RNDN); mpfr_mul(r19472, r19472, r19471, MPFR_RNDN); mpfr_set_d(r19473, z, MPFR_RNDN); mpfr_log(r19474, r19473, MPFR_RNDN); mpfr_mul(r19475, r19474, r19474, MPFR_RNDN); mpfr_mul(r19475, r19475, r19474, MPFR_RNDN); mpfr_add(r19476, r19472, r19475, MPFR_RNDN); mpfr_sub(r19477, r19474, r19471, MPFR_RNDN); mpfr_mul(r19478, r19474, r19477, MPFR_RNDN); mpfr_sqr(r19479, r19471, MPFR_RNDN); mpfr_add(r19480, r19478, r19479, MPFR_RNDN); mpfr_div(r19481, r19476, r19480, MPFR_RNDN); mpfr_set_d(r19482, t, MPFR_RNDN); mpfr_sub(r19483, r19481, r19482, MPFR_RNDN); mpfr_set_d(r19484, a, MPFR_RNDN); ; mpfr_sub(r19486, r19484, r19485, MPFR_RNDN); mpfr_log(r19487, r19482, MPFR_RNDN); mpfr_mul(r19488, r19486, r19487, MPFR_RNDN); mpfr_add(r19489, r19483, r19488, MPFR_RNDN); return mpfr_get_d(r19489, MPFR_RNDN); }