#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 r16463 = x; float r16464 = y; float r16465 = 2.0f; float r16466 = z; float r16467 = t; float r16468 = a; float r16469 = r16467 + r16468; float r16470 = sqrt(r16469); float r16471 = r16466 * r16470; float r16472 = r16471 / r16467; float r16473 = b; float r16474 = c; float r16475 = r16473 - r16474; float r16476 = 5.0f; float r16477 = 6.0f; float r16478 = r16476 / r16477; float r16479 = r16468 + r16478; float r16480 = 3.0f; float r16481 = r16467 * r16480; float r16482 = r16465 / r16481; float r16483 = r16479 - r16482; float r16484 = r16475 * r16483; float r16485 = r16472 - r16484; float r16486 = r16465 * r16485; float r16487 = exp(r16486); float r16488 = r16464 * r16487; float r16489 = r16463 + r16488; float r16490 = r16463 / r16489; return r16490; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r16491 = x; double r16492 = y; double r16493 = 2.0; double r16494 = z; double r16495 = t; double r16496 = a; double r16497 = r16495 + r16496; double r16498 = sqrt(r16497); double r16499 = r16494 * r16498; double r16500 = r16499 / r16495; double r16501 = b; double r16502 = c; double r16503 = r16501 - r16502; double r16504 = 5.0; double r16505 = 6.0; double r16506 = r16504 / r16505; double r16507 = r16496 + r16506; double r16508 = 3.0; double r16509 = r16495 * r16508; double r16510 = r16493 / r16509; double r16511 = r16507 - r16510; double r16512 = r16503 * r16511; double r16513 = r16500 - r16512; double r16514 = r16493 * r16513; double r16515 = exp(r16514); double r16516 = r16492 * r16515; double r16517 = r16491 + r16516; double r16518 = r16491 / r16517; return r16518; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r16519 = x; float r16520 = y; float r16521 = 2.0f; float r16522 = z; float r16523 = t; float r16524 = a; float r16525 = r16523 + r16524; float r16526 = sqrt(r16525); float r16527 = r16522 * r16526; float r16528 = r16527 / r16523; float r16529 = b; float r16530 = c; float r16531 = r16529 - r16530; float r16532 = 5.0f; float r16533 = 6.0f; float r16534 = r16532 / r16533; float r16535 = r16524 + r16534; float r16536 = 3.0f; float r16537 = r16523 * r16536; float r16538 = r16521 / r16537; float r16539 = r16535 - r16538; float r16540 = r16531 * r16539; float r16541 = r16528 - r16540; float r16542 = r16521 * r16541; float r16543 = exp(r16542); float r16544 = r16520 * r16543; float r16545 = r16519 + r16544; float r16546 = r16519 / r16545; return r16546; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r16547 = x; double r16548 = y; double r16549 = 2.0; double r16550 = z; double r16551 = t; double r16552 = a; double r16553 = r16551 + r16552; double r16554 = sqrt(r16553); double r16555 = r16550 * r16554; double r16556 = r16555 / r16551; double r16557 = b; double r16558 = c; double r16559 = r16557 - r16558; double r16560 = 5.0; double r16561 = 6.0; double r16562 = r16560 / r16561; double r16563 = r16552 + r16562; double r16564 = 3.0; double r16565 = r16551 * r16564; double r16566 = r16549 / r16565; double r16567 = r16563 - r16566; double r16568 = r16559 * r16567; double r16569 = r16556 - r16568; double r16570 = r16549 * r16569; double r16571 = exp(r16570); double r16572 = r16548 * r16571; double r16573 = r16547 + r16572; double r16574 = r16547 / r16573; return r16574; } 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 r16575, r16576, r16577, r16578, r16579, r16580, r16581, r16582, r16583, r16584, r16585, r16586, r16587, r16588, r16589, r16590, r16591, r16592, r16593, r16594, r16595, r16596, r16597, r16598, r16599, r16600, r16601, r16602; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16575); mpfr_init(r16576); mpfr_init_set_str(r16577, "2.0", 10, MPFR_RNDN); mpfr_init(r16578); mpfr_init(r16579); mpfr_init(r16580); mpfr_init(r16581); mpfr_init(r16582); mpfr_init(r16583); mpfr_init(r16584); mpfr_init(r16585); mpfr_init(r16586); mpfr_init(r16587); mpfr_init_set_str(r16588, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r16589, "6.0", 10, MPFR_RNDN); mpfr_init(r16590); mpfr_init(r16591); mpfr_init_set_str(r16592, "3.0", 10, MPFR_RNDN); mpfr_init(r16593); mpfr_init(r16594); mpfr_init(r16595); mpfr_init(r16596); mpfr_init(r16597); mpfr_init(r16598); mpfr_init(r16599); mpfr_init(r16600); mpfr_init(r16601); mpfr_init(r16602); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r16575, x, MPFR_RNDN); mpfr_set_d(r16576, y, MPFR_RNDN); ; mpfr_set_d(r16578, z, MPFR_RNDN); mpfr_set_d(r16579, t, MPFR_RNDN); mpfr_set_d(r16580, a, MPFR_RNDN); mpfr_add(r16581, r16579, r16580, MPFR_RNDN); mpfr_sqrt(r16582, r16581, MPFR_RNDN); mpfr_mul(r16583, r16578, r16582, MPFR_RNDN); mpfr_div(r16584, r16583, r16579, MPFR_RNDN); mpfr_set_d(r16585, b, MPFR_RNDN); mpfr_set_d(r16586, c, MPFR_RNDN); mpfr_sub(r16587, r16585, r16586, MPFR_RNDN); ; ; mpfr_div(r16590, r16588, r16589, MPFR_RNDN); mpfr_add(r16591, r16580, r16590, MPFR_RNDN); ; mpfr_mul(r16593, r16579, r16592, MPFR_RNDN); mpfr_div(r16594, r16577, r16593, MPFR_RNDN); mpfr_sub(r16595, r16591, r16594, MPFR_RNDN); mpfr_mul(r16596, r16587, r16595, MPFR_RNDN); mpfr_sub(r16597, r16584, r16596, MPFR_RNDN); mpfr_mul(r16598, r16577, r16597, MPFR_RNDN); mpfr_exp(r16599, r16598, MPFR_RNDN); mpfr_mul(r16600, r16576, r16599, MPFR_RNDN); mpfr_add(r16601, r16575, r16600, MPFR_RNDN); mpfr_div(r16602, r16575, r16601, MPFR_RNDN); return mpfr_get_d(r16602, MPFR_RNDN); } static mpfr_t r16603, r16604, r16605, r16606, r16607, r16608, r16609, r16610, r16611, r16612, r16613, r16614, r16615, r16616, r16617, r16618, r16619, r16620, r16621, r16622, r16623, r16624, r16625, r16626, r16627, r16628, r16629, r16630; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16603); mpfr_init(r16604); mpfr_init_set_str(r16605, "2.0", 10, MPFR_RNDN); mpfr_init(r16606); mpfr_init(r16607); mpfr_init(r16608); mpfr_init(r16609); mpfr_init(r16610); mpfr_init(r16611); mpfr_init(r16612); mpfr_init(r16613); mpfr_init(r16614); mpfr_init(r16615); mpfr_init_set_str(r16616, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r16617, "6.0", 10, MPFR_RNDN); mpfr_init(r16618); mpfr_init(r16619); mpfr_init_set_str(r16620, "3.0", 10, MPFR_RNDN); mpfr_init(r16621); mpfr_init(r16622); mpfr_init(r16623); mpfr_init(r16624); mpfr_init(r16625); mpfr_init(r16626); mpfr_init(r16627); mpfr_init(r16628); mpfr_init(r16629); mpfr_init(r16630); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r16603, x, MPFR_RNDN); mpfr_set_d(r16604, y, MPFR_RNDN); ; mpfr_set_d(r16606, z, MPFR_RNDN); mpfr_set_d(r16607, t, MPFR_RNDN); mpfr_set_d(r16608, a, MPFR_RNDN); mpfr_add(r16609, r16607, r16608, MPFR_RNDN); mpfr_sqrt(r16610, r16609, MPFR_RNDN); mpfr_mul(r16611, r16606, r16610, MPFR_RNDN); mpfr_div(r16612, r16611, r16607, MPFR_RNDN); mpfr_set_d(r16613, b, MPFR_RNDN); mpfr_set_d(r16614, c, MPFR_RNDN); mpfr_sub(r16615, r16613, r16614, MPFR_RNDN); ; ; mpfr_div(r16618, r16616, r16617, MPFR_RNDN); mpfr_add(r16619, r16608, r16618, MPFR_RNDN); ; mpfr_mul(r16621, r16607, r16620, MPFR_RNDN); mpfr_div(r16622, r16605, r16621, MPFR_RNDN); mpfr_sub(r16623, r16619, r16622, MPFR_RNDN); mpfr_mul(r16624, r16615, r16623, MPFR_RNDN); mpfr_sub(r16625, r16612, r16624, MPFR_RNDN); mpfr_mul(r16626, r16605, r16625, MPFR_RNDN); mpfr_exp(r16627, r16626, MPFR_RNDN); mpfr_mul(r16628, r16604, r16627, MPFR_RNDN); mpfr_add(r16629, r16603, r16628, MPFR_RNDN); mpfr_div(r16630, r16603, r16629, MPFR_RNDN); return mpfr_get_d(r16630, MPFR_RNDN); } static mpfr_t r16631, r16632, r16633, r16634, r16635, r16636, r16637, r16638, r16639, r16640, r16641, r16642, r16643, r16644, r16645, r16646, r16647, r16648, r16649, r16650, r16651, r16652, r16653, r16654, r16655, r16656, r16657, r16658; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16631); mpfr_init(r16632); mpfr_init_set_str(r16633, "2.0", 10, MPFR_RNDN); mpfr_init(r16634); mpfr_init(r16635); mpfr_init(r16636); mpfr_init(r16637); mpfr_init(r16638); mpfr_init(r16639); mpfr_init(r16640); mpfr_init(r16641); mpfr_init(r16642); mpfr_init(r16643); mpfr_init_set_str(r16644, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r16645, "6.0", 10, MPFR_RNDN); mpfr_init(r16646); mpfr_init(r16647); mpfr_init_set_str(r16648, "3.0", 10, MPFR_RNDN); mpfr_init(r16649); mpfr_init(r16650); mpfr_init(r16651); mpfr_init(r16652); mpfr_init(r16653); mpfr_init(r16654); mpfr_init(r16655); mpfr_init(r16656); mpfr_init(r16657); mpfr_init(r16658); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r16631, x, MPFR_RNDN); mpfr_set_d(r16632, y, MPFR_RNDN); ; mpfr_set_d(r16634, z, MPFR_RNDN); mpfr_set_d(r16635, t, MPFR_RNDN); mpfr_set_d(r16636, a, MPFR_RNDN); mpfr_add(r16637, r16635, r16636, MPFR_RNDN); mpfr_sqrt(r16638, r16637, MPFR_RNDN); mpfr_mul(r16639, r16634, r16638, MPFR_RNDN); mpfr_div(r16640, r16639, r16635, MPFR_RNDN); mpfr_set_d(r16641, b, MPFR_RNDN); mpfr_set_d(r16642, c, MPFR_RNDN); mpfr_sub(r16643, r16641, r16642, MPFR_RNDN); ; ; mpfr_div(r16646, r16644, r16645, MPFR_RNDN); mpfr_add(r16647, r16636, r16646, MPFR_RNDN); ; mpfr_mul(r16649, r16635, r16648, MPFR_RNDN); mpfr_div(r16650, r16633, r16649, MPFR_RNDN); mpfr_sub(r16651, r16647, r16650, MPFR_RNDN); mpfr_mul(r16652, r16643, r16651, MPFR_RNDN); mpfr_sub(r16653, r16640, r16652, MPFR_RNDN); mpfr_mul(r16654, r16633, r16653, MPFR_RNDN); mpfr_exp(r16655, r16654, MPFR_RNDN); mpfr_mul(r16656, r16632, r16655, MPFR_RNDN); mpfr_add(r16657, r16631, r16656, MPFR_RNDN); mpfr_div(r16658, r16631, r16657, MPFR_RNDN); return mpfr_get_d(r16658, MPFR_RNDN); }