#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 r25464 = x; float r25465 = y; float r25466 = 2.0; float r25467 = z; float r25468 = t; float r25469 = a; float r25470 = r25468 + r25469; float r25471 = sqrt(r25470); float r25472 = r25467 * r25471; float r25473 = r25472 / r25468; float r25474 = b; float r25475 = c; float r25476 = r25474 - r25475; float r25477 = 5.0; float r25478 = 6.0; float r25479 = r25477 / r25478; float r25480 = r25469 + r25479; float r25481 = 3.0; float r25482 = r25468 * r25481; float r25483 = r25466 / r25482; float r25484 = r25480 - r25483; float r25485 = r25476 * r25484; float r25486 = r25473 - r25485; float r25487 = r25466 * r25486; float r25488 = exp(r25487); float r25489 = r25465 * r25488; float r25490 = r25464 + r25489; float r25491 = r25464 / r25490; return r25491; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r25492 = x; double r25493 = y; double r25494 = 2.0; double r25495 = z; double r25496 = t; double r25497 = a; double r25498 = r25496 + r25497; double r25499 = sqrt(r25498); double r25500 = r25495 * r25499; double r25501 = r25500 / r25496; double r25502 = b; double r25503 = c; double r25504 = r25502 - r25503; double r25505 = 5.0; double r25506 = 6.0; double r25507 = r25505 / r25506; double r25508 = r25497 + r25507; double r25509 = 3.0; double r25510 = r25496 * r25509; double r25511 = r25494 / r25510; double r25512 = r25508 - r25511; double r25513 = r25504 * r25512; double r25514 = r25501 - r25513; double r25515 = r25494 * r25514; double r25516 = exp(r25515); double r25517 = r25493 * r25516; double r25518 = r25492 + r25517; double r25519 = r25492 / r25518; return r25519; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r25520 = x; float r25521 = y; float r25522 = 2.0; float r25523 = z; float r25524 = t; float r25525 = a; float r25526 = r25524 + r25525; float r25527 = sqrt(r25526); float r25528 = r25524 / r25527; float r25529 = r25523 / r25528; float r25530 = b; float r25531 = c; float r25532 = r25530 - r25531; float r25533 = 5.0; float r25534 = 6.0; float r25535 = r25533 / r25534; float r25536 = r25525 + r25535; float r25537 = 3.0; float r25538 = r25524 * r25537; float r25539 = r25522 / r25538; float r25540 = r25536 - r25539; float r25541 = r25532 * r25540; float r25542 = r25529 - r25541; float r25543 = r25522 * r25542; float r25544 = exp(r25543); float r25545 = r25521 * r25544; float r25546 = r25520 + r25545; float r25547 = r25520 / r25546; return r25547; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r25548 = x; double r25549 = y; double r25550 = 2.0; double r25551 = z; double r25552 = t; double r25553 = a; double r25554 = r25552 + r25553; double r25555 = sqrt(r25554); double r25556 = r25552 / r25555; double r25557 = r25551 / r25556; double r25558 = b; double r25559 = c; double r25560 = r25558 - r25559; double r25561 = 5.0; double r25562 = 6.0; double r25563 = r25561 / r25562; double r25564 = r25553 + r25563; double r25565 = 3.0; double r25566 = r25552 * r25565; double r25567 = r25550 / r25566; double r25568 = r25564 - r25567; double r25569 = r25560 * r25568; double r25570 = r25557 - r25569; double r25571 = r25550 * r25570; double r25572 = exp(r25571); double r25573 = r25549 * r25572; double r25574 = r25548 + r25573; double r25575 = r25548 / r25574; return r25575; } 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 r25576, r25577, r25578, r25579, r25580, r25581, r25582, r25583, r25584, r25585, r25586, r25587, r25588, r25589, r25590, r25591, r25592, r25593, r25594, r25595, r25596, r25597, r25598, r25599, r25600, r25601, r25602, r25603; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25576); mpfr_init(r25577); mpfr_init_set_str(r25578, "2.0", 10, MPFR_RNDN); mpfr_init(r25579); mpfr_init(r25580); mpfr_init(r25581); mpfr_init(r25582); mpfr_init(r25583); mpfr_init(r25584); mpfr_init(r25585); mpfr_init(r25586); mpfr_init(r25587); mpfr_init(r25588); mpfr_init_set_str(r25589, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25590, "6.0", 10, MPFR_RNDN); mpfr_init(r25591); mpfr_init(r25592); mpfr_init_set_str(r25593, "3.0", 10, MPFR_RNDN); mpfr_init(r25594); mpfr_init(r25595); mpfr_init(r25596); mpfr_init(r25597); mpfr_init(r25598); mpfr_init(r25599); mpfr_init(r25600); mpfr_init(r25601); mpfr_init(r25602); mpfr_init(r25603); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25576, x, MPFR_RNDN); mpfr_set_d(r25577, y, MPFR_RNDN); ; mpfr_set_d(r25579, z, MPFR_RNDN); mpfr_set_d(r25580, t, MPFR_RNDN); mpfr_set_d(r25581, a, MPFR_RNDN); mpfr_add(r25582, r25580, r25581, MPFR_RNDN); mpfr_sqrt(r25583, r25582, MPFR_RNDN); mpfr_mul(r25584, r25579, r25583, MPFR_RNDN); mpfr_div(r25585, r25584, r25580, MPFR_RNDN); mpfr_set_d(r25586, b, MPFR_RNDN); mpfr_set_d(r25587, c, MPFR_RNDN); mpfr_sub(r25588, r25586, r25587, MPFR_RNDN); ; ; mpfr_div(r25591, r25589, r25590, MPFR_RNDN); mpfr_add(r25592, r25581, r25591, MPFR_RNDN); ; mpfr_mul(r25594, r25580, r25593, MPFR_RNDN); mpfr_div(r25595, r25578, r25594, MPFR_RNDN); mpfr_sub(r25596, r25592, r25595, MPFR_RNDN); mpfr_mul(r25597, r25588, r25596, MPFR_RNDN); mpfr_sub(r25598, r25585, r25597, MPFR_RNDN); mpfr_mul(r25599, r25578, r25598, MPFR_RNDN); mpfr_exp(r25600, r25599, MPFR_RNDN); mpfr_mul(r25601, r25577, r25600, MPFR_RNDN); mpfr_add(r25602, r25576, r25601, MPFR_RNDN); mpfr_div(r25603, r25576, r25602, MPFR_RNDN); return mpfr_get_d(r25603, MPFR_RNDN); } static mpfr_t r25604, r25605, r25606, r25607, r25608, r25609, r25610, r25611, r25612, r25613, r25614, r25615, r25616, r25617, r25618, r25619, r25620, r25621, r25622, r25623, r25624, r25625, r25626, r25627, r25628, r25629, r25630, r25631; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25604); mpfr_init(r25605); mpfr_init_set_str(r25606, "2.0", 10, MPFR_RNDN); mpfr_init(r25607); mpfr_init(r25608); mpfr_init(r25609); mpfr_init(r25610); mpfr_init(r25611); mpfr_init(r25612); mpfr_init(r25613); mpfr_init(r25614); mpfr_init(r25615); mpfr_init(r25616); mpfr_init_set_str(r25617, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25618, "6.0", 10, MPFR_RNDN); mpfr_init(r25619); mpfr_init(r25620); mpfr_init_set_str(r25621, "3.0", 10, MPFR_RNDN); mpfr_init(r25622); mpfr_init(r25623); mpfr_init(r25624); mpfr_init(r25625); mpfr_init(r25626); mpfr_init(r25627); mpfr_init(r25628); mpfr_init(r25629); mpfr_init(r25630); mpfr_init(r25631); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25604, x, MPFR_RNDN); mpfr_set_d(r25605, y, MPFR_RNDN); ; mpfr_set_d(r25607, z, MPFR_RNDN); mpfr_set_d(r25608, t, MPFR_RNDN); mpfr_set_d(r25609, a, MPFR_RNDN); mpfr_add(r25610, r25608, r25609, MPFR_RNDN); mpfr_sqrt(r25611, r25610, MPFR_RNDN); mpfr_div(r25612, r25608, r25611, MPFR_RNDN); mpfr_div(r25613, r25607, r25612, MPFR_RNDN); mpfr_set_d(r25614, b, MPFR_RNDN); mpfr_set_d(r25615, c, MPFR_RNDN); mpfr_sub(r25616, r25614, r25615, MPFR_RNDN); ; ; mpfr_div(r25619, r25617, r25618, MPFR_RNDN); mpfr_add(r25620, r25609, r25619, MPFR_RNDN); ; mpfr_mul(r25622, r25608, r25621, MPFR_RNDN); mpfr_div(r25623, r25606, r25622, MPFR_RNDN); mpfr_sub(r25624, r25620, r25623, MPFR_RNDN); mpfr_mul(r25625, r25616, r25624, MPFR_RNDN); mpfr_sub(r25626, r25613, r25625, MPFR_RNDN); mpfr_mul(r25627, r25606, r25626, MPFR_RNDN); mpfr_exp(r25628, r25627, MPFR_RNDN); mpfr_mul(r25629, r25605, r25628, MPFR_RNDN); mpfr_add(r25630, r25604, r25629, MPFR_RNDN); mpfr_div(r25631, r25604, r25630, MPFR_RNDN); return mpfr_get_d(r25631, MPFR_RNDN); } static mpfr_t r25632, r25633, r25634, r25635, r25636, r25637, r25638, r25639, r25640, r25641, r25642, r25643, r25644, r25645, r25646, r25647, r25648, r25649, r25650, r25651, r25652, r25653, r25654, r25655, r25656, r25657, r25658, r25659; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25632); mpfr_init(r25633); mpfr_init_set_str(r25634, "2.0", 10, MPFR_RNDN); mpfr_init(r25635); mpfr_init(r25636); mpfr_init(r25637); mpfr_init(r25638); mpfr_init(r25639); mpfr_init(r25640); mpfr_init(r25641); mpfr_init(r25642); mpfr_init(r25643); mpfr_init(r25644); mpfr_init_set_str(r25645, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25646, "6.0", 10, MPFR_RNDN); mpfr_init(r25647); mpfr_init(r25648); mpfr_init_set_str(r25649, "3.0", 10, MPFR_RNDN); mpfr_init(r25650); mpfr_init(r25651); mpfr_init(r25652); mpfr_init(r25653); mpfr_init(r25654); mpfr_init(r25655); mpfr_init(r25656); mpfr_init(r25657); mpfr_init(r25658); mpfr_init(r25659); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25632, x, MPFR_RNDN); mpfr_set_d(r25633, y, MPFR_RNDN); ; mpfr_set_d(r25635, z, MPFR_RNDN); mpfr_set_d(r25636, t, MPFR_RNDN); mpfr_set_d(r25637, a, MPFR_RNDN); mpfr_add(r25638, r25636, r25637, MPFR_RNDN); mpfr_sqrt(r25639, r25638, MPFR_RNDN); mpfr_div(r25640, r25636, r25639, MPFR_RNDN); mpfr_div(r25641, r25635, r25640, MPFR_RNDN); mpfr_set_d(r25642, b, MPFR_RNDN); mpfr_set_d(r25643, c, MPFR_RNDN); mpfr_sub(r25644, r25642, r25643, MPFR_RNDN); ; ; mpfr_div(r25647, r25645, r25646, MPFR_RNDN); mpfr_add(r25648, r25637, r25647, MPFR_RNDN); ; mpfr_mul(r25650, r25636, r25649, MPFR_RNDN); mpfr_div(r25651, r25634, r25650, MPFR_RNDN); mpfr_sub(r25652, r25648, r25651, MPFR_RNDN); mpfr_mul(r25653, r25644, r25652, MPFR_RNDN); mpfr_sub(r25654, r25641, r25653, MPFR_RNDN); mpfr_mul(r25655, r25634, r25654, MPFR_RNDN); mpfr_exp(r25656, r25655, MPFR_RNDN); mpfr_mul(r25657, r25633, r25656, MPFR_RNDN); mpfr_add(r25658, r25632, r25657, MPFR_RNDN); mpfr_div(r25659, r25632, r25658, MPFR_RNDN); return mpfr_get_d(r25659, MPFR_RNDN); }