#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 r16359 = x; float r16360 = y; float r16361 = 2.0f; float r16362 = z; float r16363 = t; float r16364 = a; float r16365 = r16363 + r16364; float r16366 = sqrt(r16365); float r16367 = r16362 * r16366; float r16368 = r16367 / r16363; float r16369 = b; float r16370 = c; float r16371 = r16369 - r16370; float r16372 = 5.0f; float r16373 = 6.0f; float r16374 = r16372 / r16373; float r16375 = r16364 + r16374; float r16376 = 3.0f; float r16377 = r16363 * r16376; float r16378 = r16361 / r16377; float r16379 = r16375 - r16378; float r16380 = r16371 * r16379; float r16381 = r16368 - r16380; float r16382 = r16361 * r16381; float r16383 = exp(r16382); float r16384 = r16360 * r16383; float r16385 = r16359 + r16384; float r16386 = r16359 / r16385; return r16386; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r16387 = x; double r16388 = y; double r16389 = 2.0; double r16390 = z; double r16391 = t; double r16392 = a; double r16393 = r16391 + r16392; double r16394 = sqrt(r16393); double r16395 = r16390 * r16394; double r16396 = r16395 / r16391; double r16397 = b; double r16398 = c; double r16399 = r16397 - r16398; double r16400 = 5.0; double r16401 = 6.0; double r16402 = r16400 / r16401; double r16403 = r16392 + r16402; double r16404 = 3.0; double r16405 = r16391 * r16404; double r16406 = r16389 / r16405; double r16407 = r16403 - r16406; double r16408 = r16399 * r16407; double r16409 = r16396 - r16408; double r16410 = r16389 * r16409; double r16411 = exp(r16410); double r16412 = r16388 * r16411; double r16413 = r16387 + r16412; double r16414 = r16387 / r16413; return r16414; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r16415 = x; float r16416 = y; float r16417 = 2.0f; float r16418 = z; float r16419 = t; float r16420 = a; float r16421 = r16419 + r16420; float r16422 = sqrt(r16421); float r16423 = r16418 * r16422; float r16424 = r16423 / r16419; float r16425 = b; float r16426 = c; float r16427 = r16425 - r16426; float r16428 = 5.0f; float r16429 = 6.0f; float r16430 = r16428 / r16429; float r16431 = r16420 + r16430; float r16432 = 3.0f; float r16433 = r16419 * r16432; float r16434 = r16417 / r16433; float r16435 = r16431 - r16434; float r16436 = r16427 * r16435; float r16437 = r16424 - r16436; float r16438 = r16417 * r16437; float r16439 = exp(r16438); float r16440 = r16416 * r16439; float r16441 = r16415 + r16440; float r16442 = r16415 / r16441; return r16442; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r16443 = x; double r16444 = y; double r16445 = 2.0; double r16446 = z; double r16447 = t; double r16448 = a; double r16449 = r16447 + r16448; double r16450 = sqrt(r16449); double r16451 = r16446 * r16450; double r16452 = r16451 / r16447; double r16453 = b; double r16454 = c; double r16455 = r16453 - r16454; double r16456 = 5.0; double r16457 = 6.0; double r16458 = r16456 / r16457; double r16459 = r16448 + r16458; double r16460 = 3.0; double r16461 = r16447 * r16460; double r16462 = r16445 / r16461; double r16463 = r16459 - r16462; double r16464 = r16455 * r16463; double r16465 = r16452 - r16464; double r16466 = r16445 * r16465; double r16467 = exp(r16466); double r16468 = r16444 * r16467; double r16469 = r16443 + r16468; double r16470 = r16443 / r16469; return r16470; } 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 r16471, r16472, r16473, r16474, r16475, r16476, r16477, r16478, r16479, r16480, r16481, r16482, r16483, r16484, r16485, r16486, r16487, r16488, r16489, r16490, r16491, r16492, r16493, r16494, r16495, r16496, r16497, r16498; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r16471); mpfr_init(r16472); mpfr_init_set_str(r16473, "2.0", 10, MPFR_RNDN); mpfr_init(r16474); mpfr_init(r16475); mpfr_init(r16476); mpfr_init(r16477); mpfr_init(r16478); mpfr_init(r16479); mpfr_init(r16480); mpfr_init(r16481); mpfr_init(r16482); mpfr_init(r16483); mpfr_init_set_str(r16484, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r16485, "6.0", 10, MPFR_RNDN); mpfr_init(r16486); mpfr_init(r16487); mpfr_init_set_str(r16488, "3.0", 10, MPFR_RNDN); mpfr_init(r16489); mpfr_init(r16490); mpfr_init(r16491); mpfr_init(r16492); mpfr_init(r16493); mpfr_init(r16494); mpfr_init(r16495); mpfr_init(r16496); mpfr_init(r16497); mpfr_init(r16498); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r16471, x, MPFR_RNDN); mpfr_set_d(r16472, y, MPFR_RNDN); ; mpfr_set_d(r16474, z, MPFR_RNDN); mpfr_set_d(r16475, t, MPFR_RNDN); mpfr_set_d(r16476, a, MPFR_RNDN); mpfr_add(r16477, r16475, r16476, MPFR_RNDN); mpfr_sqrt(r16478, r16477, MPFR_RNDN); mpfr_mul(r16479, r16474, r16478, MPFR_RNDN); mpfr_div(r16480, r16479, r16475, MPFR_RNDN); mpfr_set_d(r16481, b, MPFR_RNDN); mpfr_set_d(r16482, c, MPFR_RNDN); mpfr_sub(r16483, r16481, r16482, MPFR_RNDN); ; ; mpfr_div(r16486, r16484, r16485, MPFR_RNDN); mpfr_add(r16487, r16476, r16486, MPFR_RNDN); ; mpfr_mul(r16489, r16475, r16488, MPFR_RNDN); mpfr_div(r16490, r16473, r16489, MPFR_RNDN); mpfr_sub(r16491, r16487, r16490, MPFR_RNDN); mpfr_mul(r16492, r16483, r16491, MPFR_RNDN); mpfr_sub(r16493, r16480, r16492, MPFR_RNDN); mpfr_mul(r16494, r16473, r16493, MPFR_RNDN); mpfr_exp(r16495, r16494, MPFR_RNDN); mpfr_mul(r16496, r16472, r16495, MPFR_RNDN); mpfr_add(r16497, r16471, r16496, MPFR_RNDN); mpfr_div(r16498, r16471, r16497, MPFR_RNDN); return mpfr_get_d(r16498, MPFR_RNDN); } static mpfr_t r16499, r16500, r16501, r16502, r16503, r16504, r16505, r16506, r16507, r16508, r16509, r16510, r16511, r16512, r16513, r16514, r16515, r16516, r16517, r16518, r16519, r16520, r16521, r16522, r16523, r16524, r16525, r16526; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r16499); mpfr_init(r16500); mpfr_init_set_str(r16501, "2.0", 10, MPFR_RNDN); mpfr_init(r16502); mpfr_init(r16503); mpfr_init(r16504); mpfr_init(r16505); mpfr_init(r16506); mpfr_init(r16507); mpfr_init(r16508); mpfr_init(r16509); mpfr_init(r16510); mpfr_init(r16511); mpfr_init_set_str(r16512, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r16513, "6.0", 10, MPFR_RNDN); mpfr_init(r16514); mpfr_init(r16515); mpfr_init_set_str(r16516, "3.0", 10, MPFR_RNDN); mpfr_init(r16517); mpfr_init(r16518); mpfr_init(r16519); mpfr_init(r16520); mpfr_init(r16521); mpfr_init(r16522); mpfr_init(r16523); mpfr_init(r16524); mpfr_init(r16525); mpfr_init(r16526); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r16499, x, MPFR_RNDN); mpfr_set_d(r16500, y, MPFR_RNDN); ; mpfr_set_d(r16502, z, MPFR_RNDN); mpfr_set_d(r16503, t, MPFR_RNDN); mpfr_set_d(r16504, a, MPFR_RNDN); mpfr_add(r16505, r16503, r16504, MPFR_RNDN); mpfr_sqrt(r16506, r16505, MPFR_RNDN); mpfr_mul(r16507, r16502, r16506, MPFR_RNDN); mpfr_div(r16508, r16507, r16503, MPFR_RNDN); mpfr_set_d(r16509, b, MPFR_RNDN); mpfr_set_d(r16510, c, MPFR_RNDN); mpfr_sub(r16511, r16509, r16510, MPFR_RNDN); ; ; mpfr_div(r16514, r16512, r16513, MPFR_RNDN); mpfr_add(r16515, r16504, r16514, MPFR_RNDN); ; mpfr_mul(r16517, r16503, r16516, MPFR_RNDN); mpfr_div(r16518, r16501, r16517, MPFR_RNDN); mpfr_sub(r16519, r16515, r16518, MPFR_RNDN); mpfr_mul(r16520, r16511, r16519, MPFR_RNDN); mpfr_sub(r16521, r16508, r16520, MPFR_RNDN); mpfr_mul(r16522, r16501, r16521, MPFR_RNDN); mpfr_exp(r16523, r16522, MPFR_RNDN); mpfr_mul(r16524, r16500, r16523, MPFR_RNDN); mpfr_add(r16525, r16499, r16524, MPFR_RNDN); mpfr_div(r16526, r16499, r16525, MPFR_RNDN); return mpfr_get_d(r16526, MPFR_RNDN); } static mpfr_t r16527, r16528, r16529, r16530, r16531, r16532, r16533, r16534, r16535, r16536, r16537, r16538, r16539, r16540, r16541, r16542, r16543, r16544, r16545, r16546, r16547, r16548, r16549, r16550, r16551, r16552, r16553, r16554; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r16527); mpfr_init(r16528); mpfr_init_set_str(r16529, "2.0", 10, MPFR_RNDN); mpfr_init(r16530); mpfr_init(r16531); mpfr_init(r16532); mpfr_init(r16533); mpfr_init(r16534); mpfr_init(r16535); mpfr_init(r16536); mpfr_init(r16537); mpfr_init(r16538); mpfr_init(r16539); mpfr_init_set_str(r16540, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r16541, "6.0", 10, MPFR_RNDN); mpfr_init(r16542); mpfr_init(r16543); mpfr_init_set_str(r16544, "3.0", 10, MPFR_RNDN); mpfr_init(r16545); mpfr_init(r16546); mpfr_init(r16547); mpfr_init(r16548); mpfr_init(r16549); mpfr_init(r16550); mpfr_init(r16551); mpfr_init(r16552); mpfr_init(r16553); mpfr_init(r16554); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r16527, x, MPFR_RNDN); mpfr_set_d(r16528, y, MPFR_RNDN); ; mpfr_set_d(r16530, z, MPFR_RNDN); mpfr_set_d(r16531, t, MPFR_RNDN); mpfr_set_d(r16532, a, MPFR_RNDN); mpfr_add(r16533, r16531, r16532, MPFR_RNDN); mpfr_sqrt(r16534, r16533, MPFR_RNDN); mpfr_mul(r16535, r16530, r16534, MPFR_RNDN); mpfr_div(r16536, r16535, r16531, MPFR_RNDN); mpfr_set_d(r16537, b, MPFR_RNDN); mpfr_set_d(r16538, c, MPFR_RNDN); mpfr_sub(r16539, r16537, r16538, MPFR_RNDN); ; ; mpfr_div(r16542, r16540, r16541, MPFR_RNDN); mpfr_add(r16543, r16532, r16542, MPFR_RNDN); ; mpfr_mul(r16545, r16531, r16544, MPFR_RNDN); mpfr_div(r16546, r16529, r16545, MPFR_RNDN); mpfr_sub(r16547, r16543, r16546, MPFR_RNDN); mpfr_mul(r16548, r16539, r16547, MPFR_RNDN); mpfr_sub(r16549, r16536, r16548, MPFR_RNDN); mpfr_mul(r16550, r16529, r16549, MPFR_RNDN); mpfr_exp(r16551, r16550, MPFR_RNDN); mpfr_mul(r16552, r16528, r16551, MPFR_RNDN); mpfr_add(r16553, r16527, r16552, MPFR_RNDN); mpfr_div(r16554, r16527, r16553, MPFR_RNDN); return mpfr_get_d(r16554, MPFR_RNDN); }