#include #include #include #include #include char *name = "Numeric.SpecFunctions:invIncompleteBetaWorker 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 r23432 = x; float r23433 = y; float r23434 = 2.0; float r23435 = z; float r23436 = t; float r23437 = a; float r23438 = r23436 + r23437; float r23439 = sqrt(r23438); float r23440 = r23435 * r23439; float r23441 = r23440 / r23436; float r23442 = b; float r23443 = c; float r23444 = r23442 - r23443; float r23445 = 5.0; float r23446 = 6.0; float r23447 = r23445 / r23446; float r23448 = r23437 + r23447; float r23449 = 3.0; float r23450 = r23436 * r23449; float r23451 = r23434 / r23450; float r23452 = r23448 - r23451; float r23453 = r23444 * r23452; float r23454 = r23441 - r23453; float r23455 = r23434 * r23454; float r23456 = exp(r23455); float r23457 = r23433 * r23456; float r23458 = r23432 + r23457; float r23459 = r23432 / r23458; return r23459; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r23460 = x; double r23461 = y; double r23462 = 2.0; double r23463 = z; double r23464 = t; double r23465 = a; double r23466 = r23464 + r23465; double r23467 = sqrt(r23466); double r23468 = r23463 * r23467; double r23469 = r23468 / r23464; double r23470 = b; double r23471 = c; double r23472 = r23470 - r23471; double r23473 = 5.0; double r23474 = 6.0; double r23475 = r23473 / r23474; double r23476 = r23465 + r23475; double r23477 = 3.0; double r23478 = r23464 * r23477; double r23479 = r23462 / r23478; double r23480 = r23476 - r23479; double r23481 = r23472 * r23480; double r23482 = r23469 - r23481; double r23483 = r23462 * r23482; double r23484 = exp(r23483); double r23485 = r23461 * r23484; double r23486 = r23460 + r23485; double r23487 = r23460 / r23486; return r23487; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r23488 = x; float r23489 = y; float r23490 = 2.0; float r23491 = z; float r23492 = t; float r23493 = a; float r23494 = r23492 + r23493; float r23495 = sqrt(r23494); float r23496 = r23492 / r23495; float r23497 = r23491 / r23496; float r23498 = b; float r23499 = c; float r23500 = r23498 - r23499; float r23501 = 5.0; float r23502 = 6.0; float r23503 = r23501 / r23502; float r23504 = r23493 + r23503; float r23505 = 3.0; float r23506 = r23492 * r23505; float r23507 = r23490 / r23506; float r23508 = r23504 - r23507; float r23509 = r23500 * r23508; float r23510 = r23497 - r23509; float r23511 = r23490 * r23510; float r23512 = exp(r23511); float r23513 = r23489 * r23512; float r23514 = r23488 + r23513; float r23515 = r23488 / r23514; return r23515; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r23516 = x; double r23517 = y; double r23518 = 2.0; double r23519 = z; double r23520 = t; double r23521 = a; double r23522 = r23520 + r23521; double r23523 = sqrt(r23522); double r23524 = r23520 / r23523; double r23525 = r23519 / r23524; double r23526 = b; double r23527 = c; double r23528 = r23526 - r23527; double r23529 = 5.0; double r23530 = 6.0; double r23531 = r23529 / r23530; double r23532 = r23521 + r23531; double r23533 = 3.0; double r23534 = r23520 * r23533; double r23535 = r23518 / r23534; double r23536 = r23532 - r23535; double r23537 = r23528 * r23536; double r23538 = r23525 - r23537; double r23539 = r23518 * r23538; double r23540 = exp(r23539); double r23541 = r23517 * r23540; double r23542 = r23516 + r23541; double r23543 = r23516 / r23542; return r23543; } 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 r23544, r23545, r23546, r23547, r23548, r23549, r23550, r23551, r23552, r23553, r23554, r23555, r23556, r23557, r23558, r23559, r23560, r23561, r23562, r23563, r23564, r23565, r23566, r23567, r23568, r23569, r23570, r23571; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r23544); mpfr_init(r23545); mpfr_init_set_str(r23546, "2.0", 10, MPFR_RNDN); mpfr_init(r23547); mpfr_init(r23548); mpfr_init(r23549); mpfr_init(r23550); mpfr_init(r23551); mpfr_init(r23552); mpfr_init(r23553); mpfr_init(r23554); mpfr_init(r23555); mpfr_init(r23556); mpfr_init_set_str(r23557, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23558, "6.0", 10, MPFR_RNDN); mpfr_init(r23559); mpfr_init(r23560); mpfr_init_set_str(r23561, "3.0", 10, MPFR_RNDN); mpfr_init(r23562); mpfr_init(r23563); mpfr_init(r23564); mpfr_init(r23565); mpfr_init(r23566); mpfr_init(r23567); mpfr_init(r23568); mpfr_init(r23569); mpfr_init(r23570); mpfr_init(r23571); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23544, x, MPFR_RNDN); mpfr_set_d(r23545, y, MPFR_RNDN); ; mpfr_set_d(r23547, z, MPFR_RNDN); mpfr_set_d(r23548, t, MPFR_RNDN); mpfr_set_d(r23549, a, MPFR_RNDN); mpfr_add(r23550, r23548, r23549, MPFR_RNDN); mpfr_sqrt(r23551, r23550, MPFR_RNDN); mpfr_mul(r23552, r23547, r23551, MPFR_RNDN); mpfr_div(r23553, r23552, r23548, MPFR_RNDN); mpfr_set_d(r23554, b, MPFR_RNDN); mpfr_set_d(r23555, c, MPFR_RNDN); mpfr_sub(r23556, r23554, r23555, MPFR_RNDN); ; ; mpfr_div(r23559, r23557, r23558, MPFR_RNDN); mpfr_add(r23560, r23549, r23559, MPFR_RNDN); ; mpfr_mul(r23562, r23548, r23561, MPFR_RNDN); mpfr_div(r23563, r23546, r23562, MPFR_RNDN); mpfr_sub(r23564, r23560, r23563, MPFR_RNDN); mpfr_mul(r23565, r23556, r23564, MPFR_RNDN); mpfr_sub(r23566, r23553, r23565, MPFR_RNDN); mpfr_mul(r23567, r23546, r23566, MPFR_RNDN); mpfr_exp(r23568, r23567, MPFR_RNDN); mpfr_mul(r23569, r23545, r23568, MPFR_RNDN); mpfr_add(r23570, r23544, r23569, MPFR_RNDN); mpfr_div(r23571, r23544, r23570, MPFR_RNDN); return mpfr_get_d(r23571, MPFR_RNDN); } static mpfr_t r23572, r23573, r23574, r23575, r23576, r23577, r23578, r23579, r23580, r23581, r23582, r23583, r23584, r23585, r23586, r23587, r23588, r23589, r23590, r23591, r23592, r23593, r23594, r23595, r23596, r23597, r23598, r23599; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r23572); mpfr_init(r23573); mpfr_init_set_str(r23574, "2.0", 10, MPFR_RNDN); mpfr_init(r23575); mpfr_init(r23576); mpfr_init(r23577); mpfr_init(r23578); mpfr_init(r23579); mpfr_init(r23580); mpfr_init(r23581); mpfr_init(r23582); mpfr_init(r23583); mpfr_init(r23584); mpfr_init_set_str(r23585, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23586, "6.0", 10, MPFR_RNDN); mpfr_init(r23587); mpfr_init(r23588); mpfr_init_set_str(r23589, "3.0", 10, MPFR_RNDN); mpfr_init(r23590); mpfr_init(r23591); mpfr_init(r23592); mpfr_init(r23593); mpfr_init(r23594); mpfr_init(r23595); mpfr_init(r23596); mpfr_init(r23597); mpfr_init(r23598); mpfr_init(r23599); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23572, x, MPFR_RNDN); mpfr_set_d(r23573, y, MPFR_RNDN); ; mpfr_set_d(r23575, z, MPFR_RNDN); mpfr_set_d(r23576, t, MPFR_RNDN); mpfr_set_d(r23577, a, MPFR_RNDN); mpfr_add(r23578, r23576, r23577, MPFR_RNDN); mpfr_sqrt(r23579, r23578, MPFR_RNDN); mpfr_div(r23580, r23576, r23579, MPFR_RNDN); mpfr_div(r23581, r23575, r23580, MPFR_RNDN); mpfr_set_d(r23582, b, MPFR_RNDN); mpfr_set_d(r23583, c, MPFR_RNDN); mpfr_sub(r23584, r23582, r23583, MPFR_RNDN); ; ; mpfr_div(r23587, r23585, r23586, MPFR_RNDN); mpfr_add(r23588, r23577, r23587, MPFR_RNDN); ; mpfr_mul(r23590, r23576, r23589, MPFR_RNDN); mpfr_div(r23591, r23574, r23590, MPFR_RNDN); mpfr_sub(r23592, r23588, r23591, MPFR_RNDN); mpfr_mul(r23593, r23584, r23592, MPFR_RNDN); mpfr_sub(r23594, r23581, r23593, MPFR_RNDN); mpfr_mul(r23595, r23574, r23594, MPFR_RNDN); mpfr_exp(r23596, r23595, MPFR_RNDN); mpfr_mul(r23597, r23573, r23596, MPFR_RNDN); mpfr_add(r23598, r23572, r23597, MPFR_RNDN); mpfr_div(r23599, r23572, r23598, MPFR_RNDN); return mpfr_get_d(r23599, MPFR_RNDN); } static mpfr_t r23600, r23601, r23602, r23603, r23604, r23605, r23606, r23607, r23608, r23609, r23610, r23611, r23612, r23613, r23614, r23615, r23616, r23617, r23618, r23619, r23620, r23621, r23622, r23623, r23624, r23625, r23626, r23627; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r23600); mpfr_init(r23601); mpfr_init_set_str(r23602, "2.0", 10, MPFR_RNDN); mpfr_init(r23603); mpfr_init(r23604); mpfr_init(r23605); mpfr_init(r23606); mpfr_init(r23607); mpfr_init(r23608); mpfr_init(r23609); mpfr_init(r23610); mpfr_init(r23611); mpfr_init(r23612); mpfr_init_set_str(r23613, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r23614, "6.0", 10, MPFR_RNDN); mpfr_init(r23615); mpfr_init(r23616); mpfr_init_set_str(r23617, "3.0", 10, MPFR_RNDN); mpfr_init(r23618); mpfr_init(r23619); mpfr_init(r23620); mpfr_init(r23621); mpfr_init(r23622); mpfr_init(r23623); mpfr_init(r23624); mpfr_init(r23625); mpfr_init(r23626); mpfr_init(r23627); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r23600, x, MPFR_RNDN); mpfr_set_d(r23601, y, MPFR_RNDN); ; mpfr_set_d(r23603, z, MPFR_RNDN); mpfr_set_d(r23604, t, MPFR_RNDN); mpfr_set_d(r23605, a, MPFR_RNDN); mpfr_add(r23606, r23604, r23605, MPFR_RNDN); mpfr_sqrt(r23607, r23606, MPFR_RNDN); mpfr_div(r23608, r23604, r23607, MPFR_RNDN); mpfr_div(r23609, r23603, r23608, MPFR_RNDN); mpfr_set_d(r23610, b, MPFR_RNDN); mpfr_set_d(r23611, c, MPFR_RNDN); mpfr_sub(r23612, r23610, r23611, MPFR_RNDN); ; ; mpfr_div(r23615, r23613, r23614, MPFR_RNDN); mpfr_add(r23616, r23605, r23615, MPFR_RNDN); ; mpfr_mul(r23618, r23604, r23617, MPFR_RNDN); mpfr_div(r23619, r23602, r23618, MPFR_RNDN); mpfr_sub(r23620, r23616, r23619, MPFR_RNDN); mpfr_mul(r23621, r23612, r23620, MPFR_RNDN); mpfr_sub(r23622, r23609, r23621, MPFR_RNDN); mpfr_mul(r23623, r23602, r23622, MPFR_RNDN); mpfr_exp(r23624, r23623, MPFR_RNDN); mpfr_mul(r23625, r23601, r23624, MPFR_RNDN); mpfr_add(r23626, r23600, r23625, MPFR_RNDN); mpfr_div(r23627, r23600, r23626, MPFR_RNDN); return mpfr_get_d(r23627, MPFR_RNDN); }