#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 r25322 = x; float r25323 = y; float r25324 = 2.0; float r25325 = z; float r25326 = t; float r25327 = a; float r25328 = r25326 + r25327; float r25329 = sqrt(r25328); float r25330 = r25325 * r25329; float r25331 = r25330 / r25326; float r25332 = b; float r25333 = c; float r25334 = r25332 - r25333; float r25335 = 5.0; float r25336 = 6.0; float r25337 = r25335 / r25336; float r25338 = r25327 + r25337; float r25339 = 3.0; float r25340 = r25326 * r25339; float r25341 = r25324 / r25340; float r25342 = r25338 - r25341; float r25343 = r25334 * r25342; float r25344 = r25331 - r25343; float r25345 = r25324 * r25344; float r25346 = exp(r25345); float r25347 = r25323 * r25346; float r25348 = r25322 + r25347; float r25349 = r25322 / r25348; return r25349; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r25350 = x; double r25351 = y; double r25352 = 2.0; double r25353 = z; double r25354 = t; double r25355 = a; double r25356 = r25354 + r25355; double r25357 = sqrt(r25356); double r25358 = r25353 * r25357; double r25359 = r25358 / r25354; double r25360 = b; double r25361 = c; double r25362 = r25360 - r25361; double r25363 = 5.0; double r25364 = 6.0; double r25365 = r25363 / r25364; double r25366 = r25355 + r25365; double r25367 = 3.0; double r25368 = r25354 * r25367; double r25369 = r25352 / r25368; double r25370 = r25366 - r25369; double r25371 = r25362 * r25370; double r25372 = r25359 - r25371; double r25373 = r25352 * r25372; double r25374 = exp(r25373); double r25375 = r25351 * r25374; double r25376 = r25350 + r25375; double r25377 = r25350 / r25376; return r25377; } double f_of(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 = a; float r25382 = 5.0; float r25383 = 6.0; float r25384 = r25382 / r25383; float r25385 = r25381 + r25384; float r25386 = 3.0; float r25387 = t; float r25388 = r25386 * r25387; float r25389 = r25380 / r25388; float r25390 = r25385 - r25389; float r25391 = c; float r25392 = b; float r25393 = r25391 - r25392; float r25394 = r25387 + r25381; float r25395 = sqrt(r25394); float r25396 = z; float r25397 = r25396 / r25387; float r25398 = r25395 * r25397; float r25399 = fma(r25390, r25393, r25398); float r25400 = 0; float r25401 = r25399 + r25400; float r25402 = r25380 * r25401; float r25403 = exp(r25402); float r25404 = r25379 * r25403; float r25405 = r25378 + r25404; float r25406 = r25378 / r25405; return r25406; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r25407 = x; double r25408 = y; double r25409 = 2.0; double r25410 = a; double r25411 = 5.0; double r25412 = 6.0; double r25413 = r25411 / r25412; double r25414 = r25410 + r25413; double r25415 = 3.0; double r25416 = t; double r25417 = r25415 * r25416; double r25418 = r25409 / r25417; double r25419 = r25414 - r25418; double r25420 = c; double r25421 = b; double r25422 = r25420 - r25421; double r25423 = r25416 + r25410; double r25424 = sqrt(r25423); double r25425 = z; double r25426 = r25425 / r25416; double r25427 = r25424 * r25426; double r25428 = fma(r25419, r25422, r25427); double r25429 = 0; double r25430 = r25428 + r25429; double r25431 = r25409 * r25430; double r25432 = exp(r25431); double r25433 = r25408 * r25432; double r25434 = r25407 + r25433; double r25435 = r25407 / r25434; return r25435; } 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 r25436, r25437, r25438, r25439, r25440, r25441, r25442, r25443, r25444, r25445, r25446, r25447, r25448, r25449, r25450, r25451, r25452, r25453, r25454, r25455, r25456, r25457, r25458, r25459, r25460, r25461, r25462, r25463; void setup_mpfr_f_im() { mpfr_set_default_prec(336); mpfr_init(r25436); mpfr_init(r25437); mpfr_init_set_str(r25438, "2.0", 10, MPFR_RNDN); mpfr_init(r25439); mpfr_init(r25440); mpfr_init(r25441); mpfr_init(r25442); mpfr_init(r25443); mpfr_init(r25444); mpfr_init(r25445); mpfr_init(r25446); mpfr_init(r25447); mpfr_init(r25448); mpfr_init_set_str(r25449, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25450, "6.0", 10, MPFR_RNDN); mpfr_init(r25451); mpfr_init(r25452); mpfr_init_set_str(r25453, "3.0", 10, MPFR_RNDN); mpfr_init(r25454); mpfr_init(r25455); mpfr_init(r25456); mpfr_init(r25457); mpfr_init(r25458); mpfr_init(r25459); mpfr_init(r25460); mpfr_init(r25461); mpfr_init(r25462); mpfr_init(r25463); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25436, x, MPFR_RNDN); mpfr_set_d(r25437, y, MPFR_RNDN); ; mpfr_set_d(r25439, z, MPFR_RNDN); mpfr_set_d(r25440, t, MPFR_RNDN); mpfr_set_d(r25441, a, MPFR_RNDN); mpfr_add(r25442, r25440, r25441, MPFR_RNDN); mpfr_sqrt(r25443, r25442, MPFR_RNDN); mpfr_mul(r25444, r25439, r25443, MPFR_RNDN); mpfr_div(r25445, r25444, r25440, MPFR_RNDN); mpfr_set_d(r25446, b, MPFR_RNDN); mpfr_set_d(r25447, c, MPFR_RNDN); mpfr_sub(r25448, r25446, r25447, MPFR_RNDN); ; ; mpfr_div(r25451, r25449, r25450, MPFR_RNDN); mpfr_add(r25452, r25441, r25451, MPFR_RNDN); ; mpfr_mul(r25454, r25440, r25453, MPFR_RNDN); mpfr_div(r25455, r25438, r25454, MPFR_RNDN); mpfr_sub(r25456, r25452, r25455, MPFR_RNDN); mpfr_mul(r25457, r25448, r25456, MPFR_RNDN); mpfr_sub(r25458, r25445, r25457, MPFR_RNDN); mpfr_mul(r25459, r25438, r25458, MPFR_RNDN); mpfr_exp(r25460, r25459, MPFR_RNDN); mpfr_mul(r25461, r25437, r25460, MPFR_RNDN); mpfr_add(r25462, r25436, r25461, MPFR_RNDN); mpfr_div(r25463, r25436, r25462, MPFR_RNDN); return mpfr_get_d(r25463, MPFR_RNDN); } static mpfr_t r25464, r25465, r25466, r25467, r25468, r25469, r25470, r25471, r25472, r25473, r25474, r25475, r25476, r25477, r25478, r25479, r25480, r25481, r25482, r25483, r25484, r25485, r25486, r25487, r25488, r25489, r25490, r25491, r25492; void setup_mpfr_f_fm() { mpfr_set_default_prec(336); mpfr_init(r25464); mpfr_init(r25465); mpfr_init_set_str(r25466, "2.0", 10, MPFR_RNDN); mpfr_init(r25467); mpfr_init_set_str(r25468, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25469, "6.0", 10, MPFR_RNDN); mpfr_init(r25470); mpfr_init(r25471); mpfr_init_set_str(r25472, "3.0", 10, MPFR_RNDN); mpfr_init(r25473); mpfr_init(r25474); mpfr_init(r25475); 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, "0", 10, MPFR_RNDN); mpfr_init(r25487); mpfr_init(r25488); mpfr_init(r25489); mpfr_init(r25490); mpfr_init(r25491); mpfr_init(r25492); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25464, x, MPFR_RNDN); mpfr_set_d(r25465, y, MPFR_RNDN); ; mpfr_set_d(r25467, a, MPFR_RNDN); ; ; mpfr_div(r25470, r25468, r25469, MPFR_RNDN); mpfr_add(r25471, r25467, r25470, MPFR_RNDN); ; mpfr_set_d(r25473, t, MPFR_RNDN); mpfr_mul(r25474, r25472, r25473, MPFR_RNDN); mpfr_div(r25475, r25466, r25474, MPFR_RNDN); mpfr_sub(r25476, r25471, r25475, MPFR_RNDN); mpfr_set_d(r25477, c, MPFR_RNDN); mpfr_set_d(r25478, b, MPFR_RNDN); mpfr_sub(r25479, r25477, r25478, MPFR_RNDN); mpfr_add(r25480, r25473, r25467, MPFR_RNDN); mpfr_sqrt(r25481, r25480, MPFR_RNDN); mpfr_set_d(r25482, z, MPFR_RNDN); mpfr_div(r25483, r25482, r25473, MPFR_RNDN); mpfr_mul(r25484, r25481, r25483, MPFR_RNDN); mpfr_fma(r25485, r25476, r25479, r25484, MPFR_RNDN); ; mpfr_add(r25487, r25485, r25486, MPFR_RNDN); mpfr_mul(r25488, r25466, r25487, MPFR_RNDN); mpfr_exp(r25489, r25488, MPFR_RNDN); mpfr_mul(r25490, r25465, r25489, MPFR_RNDN); mpfr_add(r25491, r25464, r25490, MPFR_RNDN); mpfr_div(r25492, r25464, r25491, MPFR_RNDN); return mpfr_get_d(r25492, MPFR_RNDN); } static mpfr_t 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, r25518, r25519, r25520, r25521; void setup_mpfr_f_dm() { mpfr_set_default_prec(336); mpfr_init(r25493); mpfr_init(r25494); mpfr_init_set_str(r25495, "2.0", 10, MPFR_RNDN); mpfr_init(r25496); mpfr_init_set_str(r25497, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25498, "6.0", 10, MPFR_RNDN); mpfr_init(r25499); mpfr_init(r25500); mpfr_init_set_str(r25501, "3.0", 10, MPFR_RNDN); 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); mpfr_init(r25511); mpfr_init(r25512); mpfr_init(r25513); mpfr_init(r25514); mpfr_init_set_str(r25515, "0", 10, MPFR_RNDN); mpfr_init(r25516); mpfr_init(r25517); mpfr_init(r25518); mpfr_init(r25519); mpfr_init(r25520); mpfr_init(r25521); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25493, x, MPFR_RNDN); mpfr_set_d(r25494, y, MPFR_RNDN); ; mpfr_set_d(r25496, a, MPFR_RNDN); ; ; mpfr_div(r25499, r25497, r25498, MPFR_RNDN); mpfr_add(r25500, r25496, r25499, MPFR_RNDN); ; mpfr_set_d(r25502, t, MPFR_RNDN); mpfr_mul(r25503, r25501, r25502, MPFR_RNDN); mpfr_div(r25504, r25495, r25503, MPFR_RNDN); mpfr_sub(r25505, r25500, r25504, MPFR_RNDN); mpfr_set_d(r25506, c, MPFR_RNDN); mpfr_set_d(r25507, b, MPFR_RNDN); mpfr_sub(r25508, r25506, r25507, MPFR_RNDN); mpfr_add(r25509, r25502, r25496, MPFR_RNDN); mpfr_sqrt(r25510, r25509, MPFR_RNDN); mpfr_set_d(r25511, z, MPFR_RNDN); mpfr_div(r25512, r25511, r25502, MPFR_RNDN); mpfr_mul(r25513, r25510, r25512, MPFR_RNDN); mpfr_fma(r25514, r25505, r25508, r25513, MPFR_RNDN); ; mpfr_add(r25516, r25514, r25515, MPFR_RNDN); mpfr_mul(r25517, r25495, r25516, MPFR_RNDN); mpfr_exp(r25518, r25517, MPFR_RNDN); mpfr_mul(r25519, r25494, r25518, MPFR_RNDN); mpfr_add(r25520, r25493, r25519, MPFR_RNDN); mpfr_div(r25521, r25493, r25520, MPFR_RNDN); return mpfr_get_d(r25521, MPFR_RNDN); }