#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 r25413 = x; float r25414 = y; float r25415 = 2.0f; float r25416 = z; float r25417 = t; float r25418 = a; float r25419 = r25417 + r25418; float r25420 = sqrt(r25419); float r25421 = r25416 * r25420; float r25422 = r25421 / r25417; float r25423 = b; float r25424 = c; float r25425 = r25423 - r25424; float r25426 = 5.0f; float r25427 = 6.0f; float r25428 = r25426 / r25427; float r25429 = r25418 + r25428; float r25430 = 3.0f; float r25431 = r25417 * r25430; float r25432 = r25415 / r25431; float r25433 = r25429 - r25432; float r25434 = r25425 * r25433; float r25435 = r25422 - r25434; float r25436 = r25415 * r25435; float r25437 = exp(r25436); float r25438 = r25414 * r25437; float r25439 = r25413 + r25438; float r25440 = r25413 / r25439; return r25440; } double f_id(double x, double y, double z, double t, double a, double b, double c) { double r25441 = x; double r25442 = y; double r25443 = 2.0; double r25444 = z; double r25445 = t; double r25446 = a; double r25447 = r25445 + r25446; double r25448 = sqrt(r25447); double r25449 = r25444 * r25448; double r25450 = r25449 / r25445; double r25451 = b; double r25452 = c; double r25453 = r25451 - r25452; double r25454 = 5.0; double r25455 = 6.0; double r25456 = r25454 / r25455; double r25457 = r25446 + r25456; double r25458 = 3.0; double r25459 = r25445 * r25458; double r25460 = r25443 / r25459; double r25461 = r25457 - r25460; double r25462 = r25453 * r25461; double r25463 = r25450 - r25462; double r25464 = r25443 * r25463; double r25465 = exp(r25464); double r25466 = r25442 * r25465; double r25467 = r25441 + r25466; double r25468 = r25441 / r25467; return r25468; } double f_of(float x, float y, float z, float t, float a, float b, float c) { float r25469 = t; float r25470 = -2.4368284028755833e+34f; bool r25471 = r25469 <= r25470; float r25472 = x; float r25473 = y; float r25474 = 2.0f; float r25475 = a; float r25476 = c; float r25477 = r25475 * r25476; float r25478 = 0.8333333333333334f; float r25479 = r25478 * r25476; float r25480 = r25477 + r25479; float r25481 = b; float r25482 = r25481 * r25475; float r25483 = r25480 - r25482; float r25484 = r25474 * r25483; float r25485 = exp(r25484); float r25486 = r25473 * r25485; float r25487 = r25472 + r25486; float r25488 = r25472 / r25487; float r25489 = 2.3351005256458113e-142f; bool r25490 = r25469 <= r25489; float r25491 = r25474 / r25469; float r25492 = 3.0f; float r25493 = r25491 / r25492; float r25494 = 5.0f; float r25495 = 6.0f; float r25496 = r25494 / r25495; float r25497 = r25496 + r25475; float r25498 = r25493 + r25497; float r25499 = z; float r25500 = r25475 + r25469; float r25501 = sqrt(r25500); float r25502 = r25499 * r25501; float r25503 = r25481 - r25476; float r25504 = r25503 * r25469; float r25505 = r25497 - r25493; float r25506 = r25504 * r25505; float r25507 = r25502 - r25506; float r25508 = r25498 * r25507; float r25509 = r25475 + r25496; float r25510 = r25469 * r25492; float r25511 = r25474 / r25510; float r25512 = r25509 + r25511; float r25513 = r25469 * r25512; float r25514 = r25508 / r25513; float r25515 = r25474 * r25514; float r25516 = exp(r25515); float r25517 = r25473 * r25516; float r25518 = r25472 + r25517; float r25519 = r25472 / r25518; float r25520 = r25469 + r25475; float r25521 = sqrt(r25520); float r25522 = r25499 * r25521; float r25523 = r25522 / r25469; float r25524 = r25509 - r25511; float r25525 = r25503 * r25524; float r25526 = r25523 - r25525; float r25527 = r25474 * r25526; float r25528 = exp(r25527); float r25529 = r25473 * r25528; float r25530 = r25472 + r25529; float r25531 = r25472 / r25530; float r25532 = r25490 ? r25519 : r25531; float r25533 = r25471 ? r25488 : r25532; return r25533; } double f_od(double x, double y, double z, double t, double a, double b, double c) { double r25534 = t; double r25535 = -2.4368284028755833e+34; bool r25536 = r25534 <= r25535; double r25537 = x; double r25538 = y; double r25539 = 2.0; double r25540 = a; double r25541 = c; double r25542 = r25540 * r25541; double r25543 = 0.8333333333333334; double r25544 = r25543 * r25541; double r25545 = r25542 + r25544; double r25546 = b; double r25547 = r25546 * r25540; double r25548 = r25545 - r25547; double r25549 = r25539 * r25548; double r25550 = exp(r25549); double r25551 = r25538 * r25550; double r25552 = r25537 + r25551; double r25553 = r25537 / r25552; double r25554 = 2.3351005256458113e-142; bool r25555 = r25534 <= r25554; double r25556 = r25539 / r25534; double r25557 = 3.0; double r25558 = r25556 / r25557; double r25559 = 5.0; double r25560 = 6.0; double r25561 = r25559 / r25560; double r25562 = r25561 + r25540; double r25563 = r25558 + r25562; double r25564 = z; double r25565 = r25540 + r25534; double r25566 = sqrt(r25565); double r25567 = r25564 * r25566; double r25568 = r25546 - r25541; double r25569 = r25568 * r25534; double r25570 = r25562 - r25558; double r25571 = r25569 * r25570; double r25572 = r25567 - r25571; double r25573 = r25563 * r25572; double r25574 = r25540 + r25561; double r25575 = r25534 * r25557; double r25576 = r25539 / r25575; double r25577 = r25574 + r25576; double r25578 = r25534 * r25577; double r25579 = r25573 / r25578; double r25580 = r25539 * r25579; double r25581 = exp(r25580); double r25582 = r25538 * r25581; double r25583 = r25537 + r25582; double r25584 = r25537 / r25583; double r25585 = r25534 + r25540; double r25586 = sqrt(r25585); double r25587 = r25564 * r25586; double r25588 = r25587 / r25534; double r25589 = r25574 - r25576; double r25590 = r25568 * r25589; double r25591 = r25588 - r25590; double r25592 = r25539 * r25591; double r25593 = exp(r25592); double r25594 = r25538 * r25593; double r25595 = r25537 + r25594; double r25596 = r25537 / r25595; double r25597 = r25555 ? r25584 : r25596; double r25598 = r25536 ? r25553 : r25597; return r25598; } 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 r25599, r25600, r25601, r25602, r25603, r25604, r25605, r25606, r25607, r25608, r25609, r25610, r25611, r25612, r25613, r25614, r25615, r25616, r25617, r25618, r25619, r25620, r25621, r25622, r25623, r25624, r25625, r25626; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init(r25599); mpfr_init(r25600); mpfr_init_set_str(r25601, "2.0", 10, MPFR_RNDN); mpfr_init(r25602); mpfr_init(r25603); mpfr_init(r25604); mpfr_init(r25605); mpfr_init(r25606); mpfr_init(r25607); mpfr_init(r25608); mpfr_init(r25609); mpfr_init(r25610); mpfr_init(r25611); mpfr_init_set_str(r25612, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25613, "6.0", 10, MPFR_RNDN); mpfr_init(r25614); mpfr_init(r25615); mpfr_init_set_str(r25616, "3.0", 10, MPFR_RNDN); mpfr_init(r25617); mpfr_init(r25618); mpfr_init(r25619); mpfr_init(r25620); mpfr_init(r25621); mpfr_init(r25622); mpfr_init(r25623); mpfr_init(r25624); mpfr_init(r25625); mpfr_init(r25626); } double f_im(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25599, x, MPFR_RNDN); mpfr_set_d(r25600, y, MPFR_RNDN); ; mpfr_set_d(r25602, z, MPFR_RNDN); mpfr_set_d(r25603, t, MPFR_RNDN); mpfr_set_d(r25604, a, MPFR_RNDN); mpfr_add(r25605, r25603, r25604, MPFR_RNDN); mpfr_sqrt(r25606, r25605, MPFR_RNDN); mpfr_mul(r25607, r25602, r25606, MPFR_RNDN); mpfr_div(r25608, r25607, r25603, MPFR_RNDN); mpfr_set_d(r25609, b, MPFR_RNDN); mpfr_set_d(r25610, c, MPFR_RNDN); mpfr_sub(r25611, r25609, r25610, MPFR_RNDN); ; ; mpfr_div(r25614, r25612, r25613, MPFR_RNDN); mpfr_add(r25615, r25604, r25614, MPFR_RNDN); ; mpfr_mul(r25617, r25603, r25616, MPFR_RNDN); mpfr_div(r25618, r25601, r25617, MPFR_RNDN); mpfr_sub(r25619, r25615, r25618, MPFR_RNDN); mpfr_mul(r25620, r25611, r25619, MPFR_RNDN); mpfr_sub(r25621, r25608, r25620, MPFR_RNDN); mpfr_mul(r25622, r25601, r25621, MPFR_RNDN); mpfr_exp(r25623, r25622, MPFR_RNDN); mpfr_mul(r25624, r25600, r25623, MPFR_RNDN); mpfr_add(r25625, r25599, r25624, MPFR_RNDN); mpfr_div(r25626, r25599, r25625, MPFR_RNDN); return mpfr_get_d(r25626, MPFR_RNDN); } static mpfr_t r25627, r25628, r25629, r25630, r25631, r25632, r25633, r25634, r25635, r25636, r25637, r25638, r25639, r25640, r25641, r25642, r25643, r25644, r25645, r25646, r25647, r25648, r25649, r25650, r25651, r25652, r25653, r25654, r25655, r25656, r25657, r25658, r25659, r25660, r25661, r25662, r25663, r25664, r25665, r25666, r25667, r25668, r25669, r25670, r25671, r25672, r25673, r25674, r25675, r25676, r25677, r25678, r25679, r25680, r25681, r25682, r25683, r25684, r25685, r25686, r25687, r25688, r25689, r25690, r25691; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init(r25627); mpfr_init_set_str(r25628, "-2.4368284028755833e+34", 10, MPFR_RNDN); mpfr_init(r25629); mpfr_init(r25630); mpfr_init(r25631); mpfr_init_set_str(r25632, "2.0", 10, MPFR_RNDN); mpfr_init(r25633); mpfr_init(r25634); mpfr_init(r25635); mpfr_init_set_str(r25636, "0.8333333333333334", 10, MPFR_RNDN); mpfr_init(r25637); mpfr_init(r25638); mpfr_init(r25639); mpfr_init(r25640); mpfr_init(r25641); mpfr_init(r25642); mpfr_init(r25643); mpfr_init(r25644); mpfr_init(r25645); mpfr_init(r25646); mpfr_init_set_str(r25647, "2.3351005256458113e-142", 10, MPFR_RNDN); mpfr_init(r25648); mpfr_init(r25649); mpfr_init_set_str(r25650, "3.0", 10, MPFR_RNDN); mpfr_init(r25651); mpfr_init_set_str(r25652, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25653, "6.0", 10, MPFR_RNDN); mpfr_init(r25654); mpfr_init(r25655); mpfr_init(r25656); mpfr_init(r25657); mpfr_init(r25658); mpfr_init(r25659); mpfr_init(r25660); mpfr_init(r25661); mpfr_init(r25662); mpfr_init(r25663); mpfr_init(r25664); mpfr_init(r25665); mpfr_init(r25666); mpfr_init(r25667); mpfr_init(r25668); mpfr_init(r25669); mpfr_init(r25670); mpfr_init(r25671); mpfr_init(r25672); mpfr_init(r25673); mpfr_init(r25674); mpfr_init(r25675); mpfr_init(r25676); mpfr_init(r25677); mpfr_init(r25678); mpfr_init(r25679); mpfr_init(r25680); mpfr_init(r25681); mpfr_init(r25682); mpfr_init(r25683); mpfr_init(r25684); mpfr_init(r25685); mpfr_init(r25686); mpfr_init(r25687); mpfr_init(r25688); mpfr_init(r25689); mpfr_init(r25690); mpfr_init(r25691); } double f_fm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25627, t, MPFR_RNDN); ; mpfr_set_si(r25629, mpfr_cmp(r25627, r25628) <= 0, MPFR_RNDN); mpfr_set_d(r25630, x, MPFR_RNDN); mpfr_set_d(r25631, y, MPFR_RNDN); ; mpfr_set_d(r25633, a, MPFR_RNDN); mpfr_set_d(r25634, c, MPFR_RNDN); mpfr_mul(r25635, r25633, r25634, MPFR_RNDN); ; mpfr_mul(r25637, r25636, r25634, MPFR_RNDN); mpfr_add(r25638, r25635, r25637, MPFR_RNDN); mpfr_set_d(r25639, b, MPFR_RNDN); mpfr_mul(r25640, r25639, r25633, MPFR_RNDN); mpfr_sub(r25641, r25638, r25640, MPFR_RNDN); mpfr_mul(r25642, r25632, r25641, MPFR_RNDN); mpfr_exp(r25643, r25642, MPFR_RNDN); mpfr_mul(r25644, r25631, r25643, MPFR_RNDN); mpfr_add(r25645, r25630, r25644, MPFR_RNDN); mpfr_div(r25646, r25630, r25645, MPFR_RNDN); ; mpfr_set_si(r25648, mpfr_cmp(r25627, r25647) <= 0, MPFR_RNDN); mpfr_div(r25649, r25632, r25627, MPFR_RNDN); ; mpfr_div(r25651, r25649, r25650, MPFR_RNDN); ; ; mpfr_div(r25654, r25652, r25653, MPFR_RNDN); mpfr_add(r25655, r25654, r25633, MPFR_RNDN); mpfr_add(r25656, r25651, r25655, MPFR_RNDN); mpfr_set_d(r25657, z, MPFR_RNDN); mpfr_add(r25658, r25633, r25627, MPFR_RNDN); mpfr_sqrt(r25659, r25658, MPFR_RNDN); mpfr_mul(r25660, r25657, r25659, MPFR_RNDN); mpfr_sub(r25661, r25639, r25634, MPFR_RNDN); mpfr_mul(r25662, r25661, r25627, MPFR_RNDN); mpfr_sub(r25663, r25655, r25651, MPFR_RNDN); mpfr_mul(r25664, r25662, r25663, MPFR_RNDN); mpfr_sub(r25665, r25660, r25664, MPFR_RNDN); mpfr_mul(r25666, r25656, r25665, MPFR_RNDN); mpfr_add(r25667, r25633, r25654, MPFR_RNDN); mpfr_mul(r25668, r25627, r25650, MPFR_RNDN); mpfr_div(r25669, r25632, r25668, MPFR_RNDN); mpfr_add(r25670, r25667, r25669, MPFR_RNDN); mpfr_mul(r25671, r25627, r25670, MPFR_RNDN); mpfr_div(r25672, r25666, r25671, MPFR_RNDN); mpfr_mul(r25673, r25632, r25672, MPFR_RNDN); mpfr_exp(r25674, r25673, MPFR_RNDN); mpfr_mul(r25675, r25631, r25674, MPFR_RNDN); mpfr_add(r25676, r25630, r25675, MPFR_RNDN); mpfr_div(r25677, r25630, r25676, MPFR_RNDN); mpfr_add(r25678, r25627, r25633, MPFR_RNDN); mpfr_sqrt(r25679, r25678, MPFR_RNDN); mpfr_mul(r25680, r25657, r25679, MPFR_RNDN); mpfr_div(r25681, r25680, r25627, MPFR_RNDN); mpfr_sub(r25682, r25667, r25669, MPFR_RNDN); mpfr_mul(r25683, r25661, r25682, MPFR_RNDN); mpfr_sub(r25684, r25681, r25683, MPFR_RNDN); mpfr_mul(r25685, r25632, r25684, MPFR_RNDN); mpfr_exp(r25686, r25685, MPFR_RNDN); mpfr_mul(r25687, r25631, r25686, MPFR_RNDN); mpfr_add(r25688, r25630, r25687, MPFR_RNDN); mpfr_div(r25689, r25630, r25688, MPFR_RNDN); if (mpfr_get_si(r25648, MPFR_RNDN)) { mpfr_set(r25690, r25677, MPFR_RNDN); } else { mpfr_set(r25690, r25689, MPFR_RNDN); }; if (mpfr_get_si(r25629, MPFR_RNDN)) { mpfr_set(r25691, r25646, MPFR_RNDN); } else { mpfr_set(r25691, r25690, MPFR_RNDN); }; return mpfr_get_d(r25691, MPFR_RNDN); } static mpfr_t r25692, r25693, r25694, r25695, r25696, r25697, r25698, r25699, r25700, r25701, r25702, r25703, r25704, r25705, r25706, r25707, r25708, r25709, r25710, r25711, r25712, r25713, r25714, r25715, r25716, r25717, r25718, r25719, r25720, r25721, r25722, r25723, r25724, r25725, r25726, r25727, r25728, r25729, r25730, r25731, r25732, r25733, r25734, r25735, r25736, r25737, r25738, r25739, r25740, r25741, r25742, r25743, r25744, r25745, r25746, r25747, r25748, r25749, r25750, r25751, r25752, r25753, r25754, r25755, r25756; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init(r25692); mpfr_init_set_str(r25693, "-2.4368284028755833e+34", 10, MPFR_RNDN); mpfr_init(r25694); mpfr_init(r25695); mpfr_init(r25696); mpfr_init_set_str(r25697, "2.0", 10, MPFR_RNDN); mpfr_init(r25698); mpfr_init(r25699); mpfr_init(r25700); mpfr_init_set_str(r25701, "0.8333333333333334", 10, MPFR_RNDN); mpfr_init(r25702); mpfr_init(r25703); mpfr_init(r25704); mpfr_init(r25705); mpfr_init(r25706); mpfr_init(r25707); mpfr_init(r25708); mpfr_init(r25709); mpfr_init(r25710); mpfr_init(r25711); mpfr_init_set_str(r25712, "2.3351005256458113e-142", 10, MPFR_RNDN); mpfr_init(r25713); mpfr_init(r25714); mpfr_init_set_str(r25715, "3.0", 10, MPFR_RNDN); mpfr_init(r25716); mpfr_init_set_str(r25717, "5.0", 10, MPFR_RNDN); mpfr_init_set_str(r25718, "6.0", 10, MPFR_RNDN); mpfr_init(r25719); mpfr_init(r25720); mpfr_init(r25721); mpfr_init(r25722); mpfr_init(r25723); mpfr_init(r25724); mpfr_init(r25725); mpfr_init(r25726); mpfr_init(r25727); mpfr_init(r25728); mpfr_init(r25729); mpfr_init(r25730); mpfr_init(r25731); mpfr_init(r25732); mpfr_init(r25733); mpfr_init(r25734); mpfr_init(r25735); mpfr_init(r25736); mpfr_init(r25737); mpfr_init(r25738); mpfr_init(r25739); mpfr_init(r25740); mpfr_init(r25741); mpfr_init(r25742); mpfr_init(r25743); mpfr_init(r25744); mpfr_init(r25745); mpfr_init(r25746); mpfr_init(r25747); mpfr_init(r25748); mpfr_init(r25749); mpfr_init(r25750); mpfr_init(r25751); mpfr_init(r25752); mpfr_init(r25753); mpfr_init(r25754); mpfr_init(r25755); mpfr_init(r25756); } double f_dm(double x, double y, double z, double t, double a, double b, double c) { mpfr_set_d(r25692, t, MPFR_RNDN); ; mpfr_set_si(r25694, mpfr_cmp(r25692, r25693) <= 0, MPFR_RNDN); mpfr_set_d(r25695, x, MPFR_RNDN); mpfr_set_d(r25696, y, MPFR_RNDN); ; mpfr_set_d(r25698, a, MPFR_RNDN); mpfr_set_d(r25699, c, MPFR_RNDN); mpfr_mul(r25700, r25698, r25699, MPFR_RNDN); ; mpfr_mul(r25702, r25701, r25699, MPFR_RNDN); mpfr_add(r25703, r25700, r25702, MPFR_RNDN); mpfr_set_d(r25704, b, MPFR_RNDN); mpfr_mul(r25705, r25704, r25698, MPFR_RNDN); mpfr_sub(r25706, r25703, r25705, MPFR_RNDN); mpfr_mul(r25707, r25697, r25706, MPFR_RNDN); mpfr_exp(r25708, r25707, MPFR_RNDN); mpfr_mul(r25709, r25696, r25708, MPFR_RNDN); mpfr_add(r25710, r25695, r25709, MPFR_RNDN); mpfr_div(r25711, r25695, r25710, MPFR_RNDN); ; mpfr_set_si(r25713, mpfr_cmp(r25692, r25712) <= 0, MPFR_RNDN); mpfr_div(r25714, r25697, r25692, MPFR_RNDN); ; mpfr_div(r25716, r25714, r25715, MPFR_RNDN); ; ; mpfr_div(r25719, r25717, r25718, MPFR_RNDN); mpfr_add(r25720, r25719, r25698, MPFR_RNDN); mpfr_add(r25721, r25716, r25720, MPFR_RNDN); mpfr_set_d(r25722, z, MPFR_RNDN); mpfr_add(r25723, r25698, r25692, MPFR_RNDN); mpfr_sqrt(r25724, r25723, MPFR_RNDN); mpfr_mul(r25725, r25722, r25724, MPFR_RNDN); mpfr_sub(r25726, r25704, r25699, MPFR_RNDN); mpfr_mul(r25727, r25726, r25692, MPFR_RNDN); mpfr_sub(r25728, r25720, r25716, MPFR_RNDN); mpfr_mul(r25729, r25727, r25728, MPFR_RNDN); mpfr_sub(r25730, r25725, r25729, MPFR_RNDN); mpfr_mul(r25731, r25721, r25730, MPFR_RNDN); mpfr_add(r25732, r25698, r25719, MPFR_RNDN); mpfr_mul(r25733, r25692, r25715, MPFR_RNDN); mpfr_div(r25734, r25697, r25733, MPFR_RNDN); mpfr_add(r25735, r25732, r25734, MPFR_RNDN); mpfr_mul(r25736, r25692, r25735, MPFR_RNDN); mpfr_div(r25737, r25731, r25736, MPFR_RNDN); mpfr_mul(r25738, r25697, r25737, MPFR_RNDN); mpfr_exp(r25739, r25738, MPFR_RNDN); mpfr_mul(r25740, r25696, r25739, MPFR_RNDN); mpfr_add(r25741, r25695, r25740, MPFR_RNDN); mpfr_div(r25742, r25695, r25741, MPFR_RNDN); mpfr_add(r25743, r25692, r25698, MPFR_RNDN); mpfr_sqrt(r25744, r25743, MPFR_RNDN); mpfr_mul(r25745, r25722, r25744, MPFR_RNDN); mpfr_div(r25746, r25745, r25692, MPFR_RNDN); mpfr_sub(r25747, r25732, r25734, MPFR_RNDN); mpfr_mul(r25748, r25726, r25747, MPFR_RNDN); mpfr_sub(r25749, r25746, r25748, MPFR_RNDN); mpfr_mul(r25750, r25697, r25749, MPFR_RNDN); mpfr_exp(r25751, r25750, MPFR_RNDN); mpfr_mul(r25752, r25696, r25751, MPFR_RNDN); mpfr_add(r25753, r25695, r25752, MPFR_RNDN); mpfr_div(r25754, r25695, r25753, MPFR_RNDN); if (mpfr_get_si(r25713, MPFR_RNDN)) { mpfr_set(r25755, r25742, MPFR_RNDN); } else { mpfr_set(r25755, r25754, MPFR_RNDN); }; if (mpfr_get_si(r25694, MPFR_RNDN)) { mpfr_set(r25756, r25711, MPFR_RNDN); } else { mpfr_set(r25756, r25755, MPFR_RNDN); }; return mpfr_get_d(r25756, MPFR_RNDN); }