#include #include #include #include #include char *name = "Numeric.SpecFunctions:logGamma from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b, float c, float i) { float r19490 = x; float r19491 = y; float r19492 = r19490 * r19491; float r19493 = z; float r19494 = r19492 + r19493; float r19495 = r19494 * r19491; float r19496 = 27464.7644705f; float r19497 = r19495 + r19496; float r19498 = r19497 * r19491; float r19499 = 230661.510616f; float r19500 = r19498 + r19499; float r19501 = r19500 * r19491; float r19502 = t; float r19503 = r19501 + r19502; float r19504 = a; float r19505 = r19491 + r19504; float r19506 = r19505 * r19491; float r19507 = b; float r19508 = r19506 + r19507; float r19509 = r19508 * r19491; float r19510 = c; float r19511 = r19509 + r19510; float r19512 = r19511 * r19491; float r19513 = i; float r19514 = r19512 + r19513; float r19515 = r19503 / r19514; return r19515; } double f_id(double x, double y, double z, double t, double a, double b, double c, double i) { double r19516 = x; double r19517 = y; double r19518 = r19516 * r19517; double r19519 = z; double r19520 = r19518 + r19519; double r19521 = r19520 * r19517; double r19522 = 27464.7644705; double r19523 = r19521 + r19522; double r19524 = r19523 * r19517; double r19525 = 230661.510616; double r19526 = r19524 + r19525; double r19527 = r19526 * r19517; double r19528 = t; double r19529 = r19527 + r19528; double r19530 = a; double r19531 = r19517 + r19530; double r19532 = r19531 * r19517; double r19533 = b; double r19534 = r19532 + r19533; double r19535 = r19534 * r19517; double r19536 = c; double r19537 = r19535 + r19536; double r19538 = r19537 * r19517; double r19539 = i; double r19540 = r19538 + r19539; double r19541 = r19529 / r19540; return r19541; } double f_of(float x, float y, float z, float t, float a, float b, float c, float i) { float r19542 = 1.0f; float r19543 = y; float r19544 = a; float r19545 = r19543 + r19544; float r19546 = r19545 * r19543; float r19547 = b; float r19548 = r19546 + r19547; float r19549 = r19548 * r19543; float r19550 = c; float r19551 = r19549 + r19550; float r19552 = r19551 * r19543; float r19553 = i; float r19554 = r19552 + r19553; float r19555 = x; float r19556 = r19555 * r19543; float r19557 = z; float r19558 = r19556 + r19557; float r19559 = r19558 * r19543; float r19560 = 27464.7644705f; float r19561 = r19559 + r19560; float r19562 = r19561 * r19543; float r19563 = 230661.510616f; float r19564 = r19562 + r19563; float r19565 = r19564 * r19543; float r19566 = t; float r19567 = r19565 + r19566; float r19568 = r19554 / r19567; float r19569 = r19542 / r19568; return r19569; } double f_od(double x, double y, double z, double t, double a, double b, double c, double i) { double r19570 = 1.0; double r19571 = y; double r19572 = a; double r19573 = r19571 + r19572; double r19574 = r19573 * r19571; double r19575 = b; double r19576 = r19574 + r19575; double r19577 = r19576 * r19571; double r19578 = c; double r19579 = r19577 + r19578; double r19580 = r19579 * r19571; double r19581 = i; double r19582 = r19580 + r19581; double r19583 = x; double r19584 = r19583 * r19571; double r19585 = z; double r19586 = r19584 + r19585; double r19587 = r19586 * r19571; double r19588 = 27464.7644705; double r19589 = r19587 + r19588; double r19590 = r19589 * r19571; double r19591 = 230661.510616; double r19592 = r19590 + r19591; double r19593 = r19592 * r19571; double r19594 = t; double r19595 = r19593 + r19594; double r19596 = r19582 / r19595; double r19597 = r19570 / r19596; return r19597; } 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 r19598, r19599, r19600, r19601, r19602, r19603, r19604, r19605, r19606, r19607, r19608, r19609, r19610, r19611, r19612, r19613, r19614, r19615, r19616, r19617, r19618, r19619, r19620, r19621, r19622, r19623; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init(r19598); mpfr_init(r19599); mpfr_init(r19600); mpfr_init(r19601); mpfr_init(r19602); mpfr_init(r19603); mpfr_init_set_str(r19604, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r19605); mpfr_init(r19606); mpfr_init_set_str(r19607, "230661.510616", 10, MPFR_RNDN); mpfr_init(r19608); mpfr_init(r19609); mpfr_init(r19610); mpfr_init(r19611); mpfr_init(r19612); mpfr_init(r19613); mpfr_init(r19614); mpfr_init(r19615); mpfr_init(r19616); mpfr_init(r19617); mpfr_init(r19618); mpfr_init(r19619); mpfr_init(r19620); mpfr_init(r19621); mpfr_init(r19622); mpfr_init(r19623); } double f_im(double x, double y, double z, double t, double a, double b, double c, double i) { mpfr_set_d(r19598, x, MPFR_RNDN); mpfr_set_d(r19599, y, MPFR_RNDN); mpfr_mul(r19600, r19598, r19599, MPFR_RNDN); mpfr_set_d(r19601, z, MPFR_RNDN); mpfr_add(r19602, r19600, r19601, MPFR_RNDN); mpfr_mul(r19603, r19602, r19599, MPFR_RNDN); ; mpfr_add(r19605, r19603, r19604, MPFR_RNDN); mpfr_mul(r19606, r19605, r19599, MPFR_RNDN); ; mpfr_add(r19608, r19606, r19607, MPFR_RNDN); mpfr_mul(r19609, r19608, r19599, MPFR_RNDN); mpfr_set_d(r19610, t, MPFR_RNDN); mpfr_add(r19611, r19609, r19610, MPFR_RNDN); mpfr_set_d(r19612, a, MPFR_RNDN); mpfr_add(r19613, r19599, r19612, MPFR_RNDN); mpfr_mul(r19614, r19613, r19599, MPFR_RNDN); mpfr_set_d(r19615, b, MPFR_RNDN); mpfr_add(r19616, r19614, r19615, MPFR_RNDN); mpfr_mul(r19617, r19616, r19599, MPFR_RNDN); mpfr_set_d(r19618, c, MPFR_RNDN); mpfr_add(r19619, r19617, r19618, MPFR_RNDN); mpfr_mul(r19620, r19619, r19599, MPFR_RNDN); mpfr_set_d(r19621, i, MPFR_RNDN); mpfr_add(r19622, r19620, r19621, MPFR_RNDN); mpfr_div(r19623, r19611, r19622, MPFR_RNDN); return mpfr_get_d(r19623, MPFR_RNDN); } static mpfr_t r19624, r19625, r19626, r19627, r19628, r19629, r19630, r19631, r19632, r19633, r19634, r19635, r19636, r19637, r19638, r19639, r19640, r19641, r19642, r19643, r19644, r19645, r19646, r19647, r19648, r19649, r19650, r19651; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r19624, "1", 10, MPFR_RNDN); mpfr_init(r19625); mpfr_init(r19626); mpfr_init(r19627); mpfr_init(r19628); mpfr_init(r19629); mpfr_init(r19630); mpfr_init(r19631); mpfr_init(r19632); mpfr_init(r19633); mpfr_init(r19634); mpfr_init(r19635); mpfr_init(r19636); mpfr_init(r19637); mpfr_init(r19638); mpfr_init(r19639); mpfr_init(r19640); mpfr_init(r19641); mpfr_init_set_str(r19642, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r19643); mpfr_init(r19644); mpfr_init_set_str(r19645, "230661.510616", 10, MPFR_RNDN); mpfr_init(r19646); mpfr_init(r19647); mpfr_init(r19648); mpfr_init(r19649); mpfr_init(r19650); mpfr_init(r19651); } double f_fm(double x, double y, double z, double t, double a, double b, double c, double i) { ; mpfr_set_d(r19625, y, MPFR_RNDN); mpfr_set_d(r19626, a, MPFR_RNDN); mpfr_add(r19627, r19625, r19626, MPFR_RNDN); mpfr_mul(r19628, r19627, r19625, MPFR_RNDN); mpfr_set_d(r19629, b, MPFR_RNDN); mpfr_add(r19630, r19628, r19629, MPFR_RNDN); mpfr_mul(r19631, r19630, r19625, MPFR_RNDN); mpfr_set_d(r19632, c, MPFR_RNDN); mpfr_add(r19633, r19631, r19632, MPFR_RNDN); mpfr_mul(r19634, r19633, r19625, MPFR_RNDN); mpfr_set_d(r19635, i, MPFR_RNDN); mpfr_add(r19636, r19634, r19635, MPFR_RNDN); mpfr_set_d(r19637, x, MPFR_RNDN); mpfr_mul(r19638, r19637, r19625, MPFR_RNDN); mpfr_set_d(r19639, z, MPFR_RNDN); mpfr_add(r19640, r19638, r19639, MPFR_RNDN); mpfr_mul(r19641, r19640, r19625, MPFR_RNDN); ; mpfr_add(r19643, r19641, r19642, MPFR_RNDN); mpfr_mul(r19644, r19643, r19625, MPFR_RNDN); ; mpfr_add(r19646, r19644, r19645, MPFR_RNDN); mpfr_mul(r19647, r19646, r19625, MPFR_RNDN); mpfr_set_d(r19648, t, MPFR_RNDN); mpfr_add(r19649, r19647, r19648, MPFR_RNDN); mpfr_div(r19650, r19636, r19649, MPFR_RNDN); mpfr_div(r19651, r19624, r19650, MPFR_RNDN); return mpfr_get_d(r19651, MPFR_RNDN); } static mpfr_t r19652, r19653, r19654, r19655, r19656, r19657, r19658, r19659, r19660, r19661, r19662, r19663, r19664, r19665, r19666, r19667, r19668, r19669, r19670, r19671, r19672, r19673, r19674, r19675, r19676, r19677, r19678, r19679; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r19652, "1", 10, MPFR_RNDN); mpfr_init(r19653); mpfr_init(r19654); mpfr_init(r19655); mpfr_init(r19656); mpfr_init(r19657); mpfr_init(r19658); mpfr_init(r19659); mpfr_init(r19660); mpfr_init(r19661); mpfr_init(r19662); mpfr_init(r19663); mpfr_init(r19664); mpfr_init(r19665); mpfr_init(r19666); mpfr_init(r19667); mpfr_init(r19668); mpfr_init(r19669); mpfr_init_set_str(r19670, "27464.7644705", 10, MPFR_RNDN); mpfr_init(r19671); mpfr_init(r19672); mpfr_init_set_str(r19673, "230661.510616", 10, MPFR_RNDN); mpfr_init(r19674); mpfr_init(r19675); mpfr_init(r19676); mpfr_init(r19677); mpfr_init(r19678); mpfr_init(r19679); } double f_dm(double x, double y, double z, double t, double a, double b, double c, double i) { ; mpfr_set_d(r19653, y, MPFR_RNDN); mpfr_set_d(r19654, a, MPFR_RNDN); mpfr_add(r19655, r19653, r19654, MPFR_RNDN); mpfr_mul(r19656, r19655, r19653, MPFR_RNDN); mpfr_set_d(r19657, b, MPFR_RNDN); mpfr_add(r19658, r19656, r19657, MPFR_RNDN); mpfr_mul(r19659, r19658, r19653, MPFR_RNDN); mpfr_set_d(r19660, c, MPFR_RNDN); mpfr_add(r19661, r19659, r19660, MPFR_RNDN); mpfr_mul(r19662, r19661, r19653, MPFR_RNDN); mpfr_set_d(r19663, i, MPFR_RNDN); mpfr_add(r19664, r19662, r19663, MPFR_RNDN); mpfr_set_d(r19665, x, MPFR_RNDN); mpfr_mul(r19666, r19665, r19653, MPFR_RNDN); mpfr_set_d(r19667, z, MPFR_RNDN); mpfr_add(r19668, r19666, r19667, MPFR_RNDN); mpfr_mul(r19669, r19668, r19653, MPFR_RNDN); ; mpfr_add(r19671, r19669, r19670, MPFR_RNDN); mpfr_mul(r19672, r19671, r19653, MPFR_RNDN); ; mpfr_add(r19674, r19672, r19673, MPFR_RNDN); mpfr_mul(r19675, r19674, r19653, MPFR_RNDN); mpfr_set_d(r19676, t, MPFR_RNDN); mpfr_add(r19677, r19675, r19676, MPFR_RNDN); mpfr_div(r19678, r19664, r19677, MPFR_RNDN); mpfr_div(r19679, r19652, r19678, MPFR_RNDN); return mpfr_get_d(r19679, MPFR_RNDN); }