#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 r25378 = x; float r25379 = y; float r25380 = 2.0; float r25381 = z; float r25382 = t; float r25383 = a; float r25384 = r25382 + r25383; float r25385 = sqrt(r25384); float r25386 = r25381 * r25385; float r25387 = r25386 / r25382; float r25388 = b; float r25389 = c; float r25390 = r25388 - r25389; float r25391 = 5.0; float r25392 = 6.0; float r25393 = r25391 / r25392; float r25394 = r25383 + r25393; float r25395 = 3.0; float r25396 = r25382 * r25395; float r25397 = r25380 / r25396; float r25398 = r25394 - r25397; float r25399 = r25390 * r25398; float r25400 = r25387 - r25399; float r25401 = r25380 * r25400; float r25402 = exp(r25401); float r25403 = r25379 * r25402; float r25404 = r25378 + r25403; float r25405 = r25378 / r25404; return r25405; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r25406 = x; double r25407 = y; double r25408 = 2.0; double r25409 = z; double r25410 = t; double r25411 = a; double r25412 = r25410 + r25411; double r25413 = sqrt(r25412); double r25414 = r25409 * r25413; double r25415 = r25414 / r25410; double r25416 = b; double r25417 = c; double r25418 = r25416 - r25417; double r25419 = 5.0; double r25420 = 6.0; double r25421 = r25419 / r25420; double r25422 = r25411 + r25421; double r25423 = 3.0; double r25424 = r25410 * r25423; double r25425 = r25408 / r25424; double r25426 = r25422 - r25425; double r25427 = r25418 * r25426; double r25428 = r25415 - r25427; double r25429 = r25408 * r25428; double r25430 = exp(r25429); double r25431 = r25407 * r25430; double r25432 = r25406 + r25431; double r25433 = r25406 / r25432; return r25433; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r25434 = x; float r25435 = y; float r25436 = 2.0; float r25437 = z; float r25438 = t; float r25439 = a; float r25440 = r25438 + r25439; float r25441 = sqrt(r25440); float r25442 = r25438 / r25441; float r25443 = r25437 / r25442; float r25444 = b; float r25445 = c; float r25446 = r25444 - r25445; float r25447 = 5.0; float r25448 = 6.0; float r25449 = r25447 / r25448; float r25450 = r25439 + r25449; float r25451 = 3.0; float r25452 = r25438 * r25451; float r25453 = r25436 / r25452; float r25454 = r25450 - r25453; float r25455 = r25446 * r25454; float r25456 = r25443 - r25455; float r25457 = r25436 * r25456; float r25458 = exp(r25457); float r25459 = r25435 * r25458; float r25460 = r25434 + r25459; float r25461 = r25434 / r25460; return r25461; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r25462 = x; double r25463 = y; double r25464 = 2.0; double r25465 = z; double r25466 = t; double r25467 = a; double r25468 = r25466 + r25467; double r25469 = sqrt(r25468); double r25470 = r25466 / r25469; double r25471 = r25465 / r25470; double r25472 = b; double r25473 = c; double r25474 = r25472 - r25473; double r25475 = 5.0; double r25476 = 6.0; double r25477 = r25475 / r25476; double r25478 = r25467 + r25477; double r25479 = 3.0; double r25480 = r25466 * r25479; double r25481 = r25464 / r25480; double r25482 = r25478 - r25481; double r25483 = r25474 * r25482; double r25484 = r25471 - r25483; double r25485 = r25464 * r25484; double r25486 = exp(r25485); double r25487 = r25463 * r25486; double r25488 = r25462 + r25487; double r25489 = r25462 / r25488; return r25489; } 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 r25490, r25491, r25492, r25493, r25494, r25495, r25496, r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504, r25505, r25506, r25507, r25508, r25509, r25510, r25511, r25512, r25513, r25514, r25515, r25516, r25517; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25490); mpfr_init(r25491); mpfr_init_set_str(r25492, "2.0", 10, MPFR_RNDN); mpfr_init(r25493); mpfr_init(r25494); mpfr_init(r25495); mpfr_init(r25496); mpfr_init(r25497); mpfr_init(r25498); mpfr_init(r25499); mpfr_init(r25500); mpfr_init(r25501); mpfr_init(r25502); mpfr_init_set_str(r25503, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25504, "6.0", 10, MPFR_RNDN); mpfr_init(r25505); mpfr_init(r25506); mpfr_init_set_str(r25507, "3.0", 10, MPFR_RNDN); mpfr_init(r25508); mpfr_init(r25509); mpfr_init(r25510); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init(r25514); mpfr_init(r25515); mpfr_init(r25516); mpfr_init(r25517); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25490, x, MPFR_RNDN); mpfr_set_d(r25491, y, MPFR_RNDN); ; mpfr_set_d(r25493, z, MPFR_RNDN); mpfr_set_d(r25494, t, MPFR_RNDN); mpfr_set_d(r25495, a, MPFR_RNDN); mpfr_add(r25496, r25494, r25495, MPFR_RNDN); mpfr_sqrt(r25497, r25496, MPFR_RNDN); mpfr_mul(r25498, r25493, r25497, MPFR_RNDN); mpfr_div(r25499, r25498, r25494, MPFR_RNDN); mpfr_set_d(r25500, b, MPFR_RNDN); mpfr_set_d(r25501, c, MPFR_RNDN); mpfr_sub(r25502, r25500, r25501, MPFR_RNDN); ; ; mpfr_div(r25505, r25503, r25504, MPFR_RNDN); mpfr_add(r25506, r25495, r25505, MPFR_RNDN); ; mpfr_mul(r25508, r25494, r25507, MPFR_RNDN); mpfr_div(r25509, r25492, r25508, MPFR_RNDN); mpfr_sub(r25510, r25506, r25509, MPFR_RNDN); mpfr_mul(r25511, r25502, r25510, MPFR_RNDN); mpfr_sub(r25512, r25499, r25511, MPFR_RNDN); mpfr_mul(r25513, r25492, r25512, MPFR_RNDN); mpfr_exp(r25514, r25513, MPFR_RNDN); mpfr_mul(r25515, r25491, r25514, MPFR_RNDN); mpfr_add(r25516, r25490, r25515, MPFR_RNDN); mpfr_div(r25517, r25490, r25516, MPFR_RNDN); return mpfr_get_d(r25517, MPFR_RNDN); } static mpfr_t r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537, r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25518); mpfr_init(r25519); mpfr_init_set_str(r25520, "2.0", 10, MPFR_RNDN); mpfr_init(r25521); mpfr_init(r25522); mpfr_init(r25523); mpfr_init(r25524); mpfr_init(r25525); mpfr_init(r25526); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); mpfr_init(r25530); mpfr_init_set_str(r25531, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25532, "6.0", 10, MPFR_RNDN); mpfr_init(r25533); mpfr_init(r25534); mpfr_init_set_str(r25535, "3.0", 10, MPFR_RNDN); mpfr_init(r25536); mpfr_init(r25537); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25518, x, MPFR_RNDN); mpfr_set_d(r25519, y, MPFR_RNDN); ; mpfr_set_d(r25521, z, MPFR_RNDN); mpfr_set_d(r25522, t, MPFR_RNDN); mpfr_set_d(r25523, a, MPFR_RNDN); mpfr_add(r25524, r25522, r25523, MPFR_RNDN); mpfr_sqrt(r25525, r25524, MPFR_RNDN); mpfr_div(r25526, r25522, r25525, MPFR_RNDN); mpfr_div(r25527, r25521, r25526, MPFR_RNDN); mpfr_set_d(r25528, b, MPFR_RNDN); mpfr_set_d(r25529, c, MPFR_RNDN); mpfr_sub(r25530, r25528, r25529, MPFR_RNDN); ; ; mpfr_div(r25533, r25531, r25532, MPFR_RNDN); mpfr_add(r25534, r25523, r25533, MPFR_RNDN); ; mpfr_mul(r25536, r25522, r25535, MPFR_RNDN); mpfr_div(r25537, r25520, r25536, MPFR_RNDN); mpfr_sub(r25538, r25534, r25537, MPFR_RNDN); mpfr_mul(r25539, r25530, r25538, MPFR_RNDN); mpfr_sub(r25540, r25527, r25539, MPFR_RNDN); mpfr_mul(r25541, r25520, r25540, MPFR_RNDN); mpfr_exp(r25542, r25541, MPFR_RNDN); mpfr_mul(r25543, r25519, r25542, MPFR_RNDN); mpfr_add(r25544, r25518, r25543, MPFR_RNDN); mpfr_div(r25545, r25518, r25544, MPFR_RNDN); return mpfr_get_d(r25545, MPFR_RNDN); } static mpfr_t 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, r25572, r25573; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25546); mpfr_init(r25547); mpfr_init_set_str(r25548, "2.0", 10, MPFR_RNDN); 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(r25557); mpfr_init(r25558); mpfr_init_set_str(r25559, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25560, "6.0", 10, MPFR_RNDN); mpfr_init(r25561); mpfr_init(r25562); mpfr_init_set_str(r25563, "3.0", 10, MPFR_RNDN); mpfr_init(r25564); mpfr_init(r25565); mpfr_init(r25566); mpfr_init(r25567); mpfr_init(r25568); mpfr_init(r25569); mpfr_init(r25570); mpfr_init(r25571); mpfr_init(r25572); mpfr_init(r25573); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25546, x, MPFR_RNDN); mpfr_set_d(r25547, y, MPFR_RNDN); ; mpfr_set_d(r25549, z, MPFR_RNDN); mpfr_set_d(r25550, t, MPFR_RNDN); mpfr_set_d(r25551, a, MPFR_RNDN); mpfr_add(r25552, r25550, r25551, MPFR_RNDN); mpfr_sqrt(r25553, r25552, MPFR_RNDN); mpfr_div(r25554, r25550, r25553, MPFR_RNDN); mpfr_div(r25555, r25549, r25554, MPFR_RNDN); mpfr_set_d(r25556, b, MPFR_RNDN); mpfr_set_d(r25557, c, MPFR_RNDN); mpfr_sub(r25558, r25556, r25557, MPFR_RNDN); ; ; mpfr_div(r25561, r25559, r25560, MPFR_RNDN); mpfr_add(r25562, r25551, r25561, MPFR_RNDN); ; mpfr_mul(r25564, r25550, r25563, MPFR_RNDN); mpfr_div(r25565, r25548, r25564, MPFR_RNDN); mpfr_sub(r25566, r25562, r25565, MPFR_RNDN); mpfr_mul(r25567, r25558, r25566, MPFR_RNDN); mpfr_sub(r25568, r25555, r25567, MPFR_RNDN); mpfr_mul(r25569, r25548, r25568, MPFR_RNDN); mpfr_exp(r25570, r25569, MPFR_RNDN); mpfr_mul(r25571, r25547, r25570, MPFR_RNDN); mpfr_add(r25572, r25546, r25571, MPFR_RNDN); mpfr_div(r25573, r25546, r25572, MPFR_RNDN); return mpfr_get_d(r25573, MPFR_RNDN); }