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