#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 r25430 = x; float r25431 = y; float r25432 = 2.0; float r25433 = z; float r25434 = t; float r25435 = a; float r25436 = r25434 + r25435; float r25437 = sqrt(r25436); float r25438 = r25433 * r25437; float r25439 = r25438 / r25434; float r25440 = b; float r25441 = c; float r25442 = r25440 - r25441; float r25443 = 5.0; float r25444 = 6.0; float r25445 = r25443 / r25444; float r25446 = r25435 + r25445; float r25447 = 3.0; float r25448 = r25434 * r25447; float r25449 = r25432 / r25448; float r25450 = r25446 - r25449; float r25451 = r25442 * r25450; float r25452 = r25439 - r25451; float r25453 = r25432 * r25452; float r25454 = exp(r25453); float r25455 = r25431 * r25454; float r25456 = r25430 + r25455; float r25457 = r25430 / r25456; return r25457; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r25458 = x; double r25459 = y; double r25460 = 2.0; double r25461 = z; double r25462 = t; double r25463 = a; double r25464 = r25462 + r25463; double r25465 = sqrt(r25464); double r25466 = r25461 * r25465; double r25467 = r25466 / r25462; double r25468 = b; double r25469 = c; double r25470 = r25468 - r25469; double r25471 = 5.0; double r25472 = 6.0; double r25473 = r25471 / r25472; double r25474 = r25463 + r25473; double r25475 = 3.0; double r25476 = r25462 * r25475; double r25477 = r25460 / r25476; double r25478 = r25474 - r25477; double r25479 = r25470 * r25478; double r25480 = r25467 - r25479; double r25481 = r25460 * r25480; double r25482 = exp(r25481); double r25483 = r25459 * r25482; double r25484 = r25458 + r25483; double r25485 = r25458 / r25484; return r25485; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r25486 = x; float r25487 = y; float r25488 = 2.0; float r25489 = a; float r25490 = 5.0; float r25491 = 6.0; float r25492 = r25490 / r25491; float r25493 = r25489 + r25492; float r25494 = 3.0; float r25495 = t; float r25496 = r25494 * r25495; float r25497 = r25488 / r25496; float r25498 = r25493 - r25497; float r25499 = c; float r25500 = b; float r25501 = r25499 - r25500; float r25502 = r25495 + r25489; float r25503 = sqrt(r25502); float r25504 = z; float r25505 = r25504 / r25495; float r25506 = r25503 * r25505; float r25507 = fma(r25498, r25501, r25506); float r25508 = 0; float r25509 = r25507 + r25508; float r25510 = r25488 * r25509; float r25511 = exp(r25510); float r25512 = r25487 * r25511; float r25513 = r25486 + r25512; float r25514 = r25486 / r25513; return r25514; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r25515 = x; double r25516 = y; double r25517 = 2.0; double r25518 = a; double r25519 = 5.0; double r25520 = 6.0; double r25521 = r25519 / r25520; double r25522 = r25518 + r25521; double r25523 = 3.0; double r25524 = t; double r25525 = r25523 * r25524; double r25526 = r25517 / r25525; double r25527 = r25522 - r25526; double r25528 = c; double r25529 = b; double r25530 = r25528 - r25529; double r25531 = r25524 + r25518; double r25532 = sqrt(r25531); double r25533 = z; double r25534 = r25533 / r25524; double r25535 = r25532 * r25534; double r25536 = fma(r25527, r25530, r25535); double r25537 = 0; double r25538 = r25536 + r25537; double r25539 = r25517 * r25538; double r25540 = exp(r25539); double r25541 = r25516 * r25540; double r25542 = r25515 + r25541; double r25543 = r25515 / r25542; return r25543; } 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 r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564, r25565, r25566, r25567, r25568, r25569, r25570, r25571; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25544); mpfr_init(r25545); mpfr_init_set_str(r25546, "2.0", 10, MPFR_RNDN); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); mpfr_init(r25552); mpfr_init(r25553); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init_set_str(r25557, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25558, "6.0", 10, MPFR_RNDN); mpfr_init(r25559); mpfr_init(r25560); mpfr_init_set_str(r25561, "3.0", 10, MPFR_RNDN); mpfr_init(r25562); mpfr_init(r25563); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); mpfr_init(r25567); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); mpfr_init(r25571); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25544, x, MPFR_RNDN); mpfr_set_d(r25545, y, MPFR_RNDN); ; mpfr_set_d(r25547, z, MPFR_RNDN); mpfr_set_d(r25548, t, MPFR_RNDN); mpfr_set_d(r25549, a, MPFR_RNDN); mpfr_add(r25550, r25548, r25549, MPFR_RNDN); mpfr_sqrt(r25551, r25550, MPFR_RNDN); mpfr_mul(r25552, r25547, r25551, MPFR_RNDN); mpfr_div(r25553, r25552, r25548, MPFR_RNDN); mpfr_set_d(r25554, b, MPFR_RNDN); mpfr_set_d(r25555, c, MPFR_RNDN); mpfr_sub(r25556, r25554, r25555, MPFR_RNDN); ; ; mpfr_div(r25559, r25557, r25558, MPFR_RNDN); mpfr_add(r25560, r25549, r25559, MPFR_RNDN); ; mpfr_mul(r25562, r25548, r25561, MPFR_RNDN); mpfr_div(r25563, r25546, r25562, MPFR_RNDN); mpfr_sub(r25564, r25560, r25563, MPFR_RNDN); mpfr_mul(r25565, r25556, r25564, MPFR_RNDN); mpfr_sub(r25566, r25553, r25565, MPFR_RNDN); mpfr_mul(r25567, r25546, r25566, MPFR_RNDN); mpfr_exp(r25568, r25567, MPFR_RNDN); mpfr_mul(r25569, r25545, r25568, MPFR_RNDN); mpfr_add(r25570, r25544, r25569, MPFR_RNDN); mpfr_div(r25571, r25544, r25570, MPFR_RNDN); return mpfr_get_d(r25571, MPFR_RNDN); } static mpfr_t r25572, r25573, r25574, r25575, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25572); mpfr_init(r25573); mpfr_init_set_str(r25574, "2.0", 10, MPFR_RNDN); mpfr_init(r25575); mpfr_init_set_str(r25576, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25577, "6.0", 10, MPFR_RNDN); mpfr_init(r25578); mpfr_init(r25579); mpfr_init_set_str(r25580, "3.0", 10, MPFR_RNDN); 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(r25589); mpfr_init(r25590); mpfr_init(r25591); mpfr_init(r25592); mpfr_init(r25593); mpfr_init_set_str(r25594, "0", 10, MPFR_RNDN); mpfr_init(r25595); mpfr_init(r25596); mpfr_init(r25597); mpfr_init(r25598); mpfr_init(r25599); mpfr_init(r25600); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25572, x, MPFR_RNDN); mpfr_set_d(r25573, y, MPFR_RNDN); ; mpfr_set_d(r25575, a, MPFR_RNDN); ; ; mpfr_div(r25578, r25576, r25577, MPFR_RNDN); mpfr_add(r25579, r25575, r25578, MPFR_RNDN); ; mpfr_set_d(r25581, t, MPFR_RNDN); mpfr_mul(r25582, r25580, r25581, MPFR_RNDN); mpfr_div(r25583, r25574, r25582, MPFR_RNDN); mpfr_sub(r25584, r25579, r25583, MPFR_RNDN); mpfr_set_d(r25585, c, MPFR_RNDN); mpfr_set_d(r25586, b, MPFR_RNDN); mpfr_sub(r25587, r25585, r25586, MPFR_RNDN); mpfr_add(r25588, r25581, r25575, MPFR_RNDN); mpfr_sqrt(r25589, r25588, MPFR_RNDN); mpfr_set_d(r25590, z, MPFR_RNDN); mpfr_div(r25591, r25590, r25581, MPFR_RNDN); mpfr_mul(r25592, r25589, r25591, MPFR_RNDN); mpfr_fma(r25593, r25584, r25587, r25592, MPFR_RNDN); ; mpfr_add(r25595, r25593, r25594, MPFR_RNDN); mpfr_mul(r25596, r25574, r25595, MPFR_RNDN); mpfr_exp(r25597, r25596, MPFR_RNDN); mpfr_mul(r25598, r25573, r25597, MPFR_RNDN); mpfr_add(r25599, r25572, r25598, MPFR_RNDN); mpfr_div(r25600, r25572, r25599, MPFR_RNDN); return mpfr_get_d(r25600, MPFR_RNDN); } static mpfr_t r25601, r25602, r25603, 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; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25601); mpfr_init(r25602); mpfr_init_set_str(r25603, "2.0", 10, MPFR_RNDN); mpfr_init(r25604); mpfr_init_set_str(r25605, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25606, "6.0", 10, MPFR_RNDN); mpfr_init(r25607); mpfr_init(r25608); mpfr_init_set_str(r25609, "3.0", 10, MPFR_RNDN); mpfr_init(r25610); mpfr_init(r25611); mpfr_init(r25612); mpfr_init(r25613); mpfr_init(r25614); mpfr_init(r25615); mpfr_init(r25616); mpfr_init(r25617); mpfr_init(r25618); mpfr_init(r25619); mpfr_init(r25620); mpfr_init(r25621); mpfr_init(r25622); mpfr_init_set_str(r25623, "0", 10, MPFR_RNDN); mpfr_init(r25624); mpfr_init(r25625); mpfr_init(r25626); mpfr_init(r25627); mpfr_init(r25628); mpfr_init(r25629); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25601, x, MPFR_RNDN); mpfr_set_d(r25602, y, MPFR_RNDN); ; mpfr_set_d(r25604, a, MPFR_RNDN); ; ; mpfr_div(r25607, r25605, r25606, MPFR_RNDN); mpfr_add(r25608, r25604, r25607, MPFR_RNDN); ; mpfr_set_d(r25610, t, MPFR_RNDN); mpfr_mul(r25611, r25609, r25610, MPFR_RNDN); mpfr_div(r25612, r25603, r25611, MPFR_RNDN); mpfr_sub(r25613, r25608, r25612, MPFR_RNDN); mpfr_set_d(r25614, c, MPFR_RNDN); mpfr_set_d(r25615, b, MPFR_RNDN); mpfr_sub(r25616, r25614, r25615, MPFR_RNDN); mpfr_add(r25617, r25610, r25604, MPFR_RNDN); mpfr_sqrt(r25618, r25617, MPFR_RNDN); mpfr_set_d(r25619, z, MPFR_RNDN); mpfr_div(r25620, r25619, r25610, MPFR_RNDN); mpfr_mul(r25621, r25618, r25620, MPFR_RNDN); mpfr_fma(r25622, r25613, r25616, r25621, MPFR_RNDN); ; mpfr_add(r25624, r25622, r25623, MPFR_RNDN); mpfr_mul(r25625, r25603, r25624, MPFR_RNDN); mpfr_exp(r25626, r25625, MPFR_RNDN); mpfr_mul(r25627, r25602, r25626, MPFR_RNDN); mpfr_add(r25628, r25601, r25627, MPFR_RNDN); mpfr_div(r25629, r25601, r25628, MPFR_RNDN); return mpfr_get_d(r25629, MPFR_RNDN); }