#include #include #include #include #include char *name = "Numeric.SpecFunctions:incompleteBetaWorker from math-functions-0.1.5.2"; double f_if(float x, float y, float z, float t, float a, float b) { float r25406 = x; float r25407 = y; float r25408 = z; float r25409 = log(r25408); float r25410 = r25407 * r25409; float r25411 = t; float r25412 = 1.0; float r25413 = r25411 - r25412; float r25414 = a; float r25415 = log(r25414); float r25416 = r25413 * r25415; float r25417 = r25410 + r25416; float r25418 = b; float r25419 = r25417 - r25418; float r25420 = exp(r25419); float r25421 = r25406 * r25420; float r25422 = r25421 / r25407; return r25422; } double f_id(double x, double y, double z, double t, double a, double b) { double r25423 = x; double r25424 = y; double r25425 = z; double r25426 = log(r25425); double r25427 = r25424 * r25426; double r25428 = t; double r25429 = 1.0; double r25430 = r25428 - r25429; double r25431 = a; double r25432 = log(r25431); double r25433 = r25430 * r25432; double r25434 = r25427 + r25433; double r25435 = b; double r25436 = r25434 - r25435; double r25437 = exp(r25436); double r25438 = r25423 * r25437; double r25439 = r25438 / r25424; return r25439; } double f_of(float x, float y, float z, float t, float a, float b) { float r25440 = z; float r25441 = y; float r25442 = pow(r25440, r25441); float r25443 = x; float r25444 = r25441 / r25443; float r25445 = r25442 / r25444; float r25446 = a; float r25447 = t; float r25448 = 1.0; float r25449 = r25447 - r25448; float r25450 = pow(r25446, r25449); float r25451 = b; float r25452 = exp(r25451); float r25453 = r25450 / r25452; float r25454 = r25445 * r25453; float r25455 = 2.079865453438447e+280; bool r25456 = r25454 <= r25455; float r25457 = log(r25440); float r25458 = r25441 * r25457; float r25459 = log(r25446); float r25460 = r25449 * r25459; float r25461 = r25458 + r25460; float r25462 = r25461 - r25451; float r25463 = exp(r25462); float r25464 = r25443 * r25463; float r25465 = r25464 / r25441; float r25466 = r25456 ? r25454 : r25465; return r25466; } double f_od(double x, double y, double z, double t, double a, double b) { double r25467 = z; double r25468 = y; double r25469 = pow(r25467, r25468); double r25470 = x; double r25471 = r25468 / r25470; double r25472 = r25469 / r25471; double r25473 = a; double r25474 = t; double r25475 = 1.0; double r25476 = r25474 - r25475; double r25477 = pow(r25473, r25476); double r25478 = b; double r25479 = exp(r25478); double r25480 = r25477 / r25479; double r25481 = r25472 * r25480; double r25482 = 2.079865453438447e+280; bool r25483 = r25481 <= r25482; double r25484 = log(r25467); double r25485 = r25468 * r25484; double r25486 = log(r25473); double r25487 = r25476 * r25486; double r25488 = r25485 + r25487; double r25489 = r25488 - r25478; double r25490 = exp(r25489); double r25491 = r25470 * r25490; double r25492 = r25491 / r25468; double r25493 = r25483 ? r25481 : r25492; return r25493; } 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 r25494, r25495, r25496, r25497, r25498, r25499, r25500, r25501, r25502, r25503, r25504, r25505, r25506, r25507, r25508, r25509, r25510; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init(r25494); mpfr_init(r25495); mpfr_init(r25496); mpfr_init(r25497); mpfr_init(r25498); mpfr_init(r25499); mpfr_init_set_str(r25500, "1.0", 10, MPFR_RNDN); mpfr_init(r25501); mpfr_init(r25502); mpfr_init(r25503); mpfr_init(r25504); mpfr_init(r25505); mpfr_init(r25506); mpfr_init(r25507); mpfr_init(r25508); mpfr_init(r25509); mpfr_init(r25510); } double f_im(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r25494, x, MPFR_RNDN); mpfr_set_d(r25495, y, MPFR_RNDN); mpfr_set_d(r25496, z, MPFR_RNDN); mpfr_log(r25497, r25496, MPFR_RNDN); mpfr_mul(r25498, r25495, r25497, MPFR_RNDN); mpfr_set_d(r25499, t, MPFR_RNDN); ; mpfr_sub(r25501, r25499, r25500, MPFR_RNDN); mpfr_set_d(r25502, a, MPFR_RNDN); mpfr_log(r25503, r25502, MPFR_RNDN); mpfr_mul(r25504, r25501, r25503, MPFR_RNDN); mpfr_add(r25505, r25498, r25504, MPFR_RNDN); mpfr_set_d(r25506, b, MPFR_RNDN); mpfr_sub(r25507, r25505, r25506, MPFR_RNDN); mpfr_exp(r25508, r25507, MPFR_RNDN); mpfr_mul(r25509, r25494, r25508, MPFR_RNDN); mpfr_div(r25510, r25509, r25495, MPFR_RNDN); return mpfr_get_d(r25510, MPFR_RNDN); } static mpfr_t r25511, r25512, r25513, r25514, r25515, r25516, r25517, r25518, r25519, r25520, r25521, r25522, r25523, r25524, r25525, r25526, r25527, r25528, r25529, r25530, r25531, r25532, r25533, r25534, r25535, r25536, r25537; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init(r25514); mpfr_init(r25515); mpfr_init(r25516); mpfr_init(r25517); mpfr_init(r25518); mpfr_init_set_str(r25519, "1.0", 10, MPFR_RNDN); mpfr_init(r25520); mpfr_init(r25521); mpfr_init(r25522); mpfr_init(r25523); mpfr_init(r25524); mpfr_init(r25525); mpfr_init_set_str(r25526, "2.079865453438447e+280", 10, MPFR_RNDN); mpfr_init(r25527); mpfr_init(r25528); mpfr_init(r25529); mpfr_init(r25530); mpfr_init(r25531); mpfr_init(r25532); mpfr_init(r25533); mpfr_init(r25534); mpfr_init(r25535); mpfr_init(r25536); mpfr_init(r25537); } double f_fm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r25511, z, MPFR_RNDN); mpfr_set_d(r25512, y, MPFR_RNDN); mpfr_pow(r25513, r25511, r25512, MPFR_RNDN); mpfr_set_d(r25514, x, MPFR_RNDN); mpfr_div(r25515, r25512, r25514, MPFR_RNDN); mpfr_div(r25516, r25513, r25515, MPFR_RNDN); mpfr_set_d(r25517, a, MPFR_RNDN); mpfr_set_d(r25518, t, MPFR_RNDN); ; mpfr_sub(r25520, r25518, r25519, MPFR_RNDN); mpfr_pow(r25521, r25517, r25520, MPFR_RNDN); mpfr_set_d(r25522, b, MPFR_RNDN); mpfr_exp(r25523, r25522, MPFR_RNDN); mpfr_div(r25524, r25521, r25523, MPFR_RNDN); mpfr_mul(r25525, r25516, r25524, MPFR_RNDN); ; mpfr_set_si(r25527, mpfr_cmp(r25525, r25526) <= 0, MPFR_RNDN); mpfr_log(r25528, r25511, MPFR_RNDN); mpfr_mul(r25529, r25512, r25528, MPFR_RNDN); mpfr_log(r25530, r25517, MPFR_RNDN); mpfr_mul(r25531, r25520, r25530, MPFR_RNDN); mpfr_add(r25532, r25529, r25531, MPFR_RNDN); mpfr_sub(r25533, r25532, r25522, MPFR_RNDN); mpfr_exp(r25534, r25533, MPFR_RNDN); mpfr_mul(r25535, r25514, r25534, MPFR_RNDN); mpfr_div(r25536, r25535, r25512, MPFR_RNDN); if (mpfr_get_si(r25527, MPFR_RNDN)) { mpfr_set(r25537, r25525, MPFR_RNDN); } else { mpfr_set(r25537, r25536, MPFR_RNDN); }; return mpfr_get_d(r25537, MPFR_RNDN); } static mpfr_t r25538, r25539, r25540, r25541, r25542, r25543, r25544, r25545, r25546, r25547, r25548, r25549, r25550, r25551, r25552, r25553, r25554, r25555, r25556, r25557, r25558, r25559, r25560, r25561, r25562, r25563, r25564; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init(r25538); mpfr_init(r25539); mpfr_init(r25540); mpfr_init(r25541); mpfr_init(r25542); mpfr_init(r25543); mpfr_init(r25544); mpfr_init(r25545); mpfr_init_set_str(r25546, "1.0", 10, MPFR_RNDN); mpfr_init(r25547); mpfr_init(r25548); mpfr_init(r25549); mpfr_init(r25550); mpfr_init(r25551); mpfr_init(r25552); mpfr_init_set_str(r25553, "2.079865453438447e+280", 10, MPFR_RNDN); mpfr_init(r25554); mpfr_init(r25555); mpfr_init(r25556); mpfr_init(r25557); mpfr_init(r25558); mpfr_init(r25559); mpfr_init(r25560); mpfr_init(r25561); mpfr_init(r25562); mpfr_init(r25563); mpfr_init(r25564); } double f_dm(double x, double y, double z, double t, double a, double b) { mpfr_set_d(r25538, z, MPFR_RNDN); mpfr_set_d(r25539, y, MPFR_RNDN); mpfr_pow(r25540, r25538, r25539, MPFR_RNDN); mpfr_set_d(r25541, x, MPFR_RNDN); mpfr_div(r25542, r25539, r25541, MPFR_RNDN); mpfr_div(r25543, r25540, r25542, MPFR_RNDN); mpfr_set_d(r25544, a, MPFR_RNDN); mpfr_set_d(r25545, t, MPFR_RNDN); ; mpfr_sub(r25547, r25545, r25546, MPFR_RNDN); mpfr_pow(r25548, r25544, r25547, MPFR_RNDN); mpfr_set_d(r25549, b, MPFR_RNDN); mpfr_exp(r25550, r25549, MPFR_RNDN); mpfr_div(r25551, r25548, r25550, MPFR_RNDN); mpfr_mul(r25552, r25543, r25551, MPFR_RNDN); ; mpfr_set_si(r25554, mpfr_cmp(r25552, r25553) <= 0, MPFR_RNDN); mpfr_log(r25555, r25538, MPFR_RNDN); mpfr_mul(r25556, r25539, r25555, MPFR_RNDN); mpfr_log(r25557, r25544, MPFR_RNDN); mpfr_mul(r25558, r25547, r25557, MPFR_RNDN); mpfr_add(r25559, r25556, r25558, MPFR_RNDN); mpfr_sub(r25560, r25559, r25549, MPFR_RNDN); mpfr_exp(r25561, r25560, MPFR_RNDN); mpfr_mul(r25562, r25541, r25561, MPFR_RNDN); mpfr_div(r25563, r25562, r25539, MPFR_RNDN); if (mpfr_get_si(r25554, MPFR_RNDN)) { mpfr_set(r25564, r25552, MPFR_RNDN); } else { mpfr_set(r25564, r25563, MPFR_RNDN); }; return mpfr_get_d(r25564, MPFR_RNDN); }