#include #include #include #include #include char *name = "Toniolo and Linder, Equation (13)"; double f_if(float n, float U, float t, float l, float Om, float U_) { float r30462522 = 2; float r30462523 = n; float r30462524 = r30462522 * r30462523; float r30462525 = U; float r30462526 = r30462524 * r30462525; float r30462527 = t; float r30462528 = l; float r30462529 = r30462528 * r30462528; float r30462530 = Om; float r30462531 = r30462529 / r30462530; float r30462532 = r30462522 * r30462531; float r30462533 = r30462527 - r30462532; float r30462534 = r30462528 / r30462530; float r30462535 = pow(r30462534, r30462522); float r30462536 = r30462523 * r30462535; float r30462537 = U_; float r30462538 = r30462525 - r30462537; float r30462539 = r30462536 * r30462538; float r30462540 = r30462533 - r30462539; float r30462541 = r30462526 * r30462540; float r30462542 = sqrt(r30462541); return r30462542; } double f_id(double n, double U, double t, double l, double Om, double U_) { double r30462543 = 2; double r30462544 = n; double r30462545 = r30462543 * r30462544; double r30462546 = U; double r30462547 = r30462545 * r30462546; double r30462548 = t; double r30462549 = l; double r30462550 = r30462549 * r30462549; double r30462551 = Om; double r30462552 = r30462550 / r30462551; double r30462553 = r30462543 * r30462552; double r30462554 = r30462548 - r30462553; double r30462555 = r30462549 / r30462551; double r30462556 = pow(r30462555, r30462543); double r30462557 = r30462544 * r30462556; double r30462558 = U_; double r30462559 = r30462546 - r30462558; double r30462560 = r30462557 * r30462559; double r30462561 = r30462554 - r30462560; double r30462562 = r30462547 * r30462561; double r30462563 = sqrt(r30462562); return r30462563; } double f_of(float n, float U, float t, float l, float Om, float U_) { float r30462564 = t; float r30462565 = l; float r30462566 = Om; float r30462567 = r30462565 / r30462566; float r30462568 = r30462565 + r30462565; float r30462569 = r30462567 * r30462568; float r30462570 = r30462564 - r30462569; float r30462571 = -inf.0; bool r30462572 = r30462570 <= r30462571; float r30462573 = n; float r30462574 = r30462564 * r30462573; float r30462575 = U; float r30462576 = r30462575 + r30462575; float r30462577 = r30462574 * r30462576; float r30462578 = r30462573 / r30462566; float r30462579 = r30462575 * r30462565; float r30462580 = r30462578 * r30462579; float r30462581 = 4; float r30462582 = r30462581 * r30462565; float r30462583 = r30462580 * r30462582; float r30462584 = r30462577 - r30462583; float r30462585 = sqrt(r30462584); float r30462586 = 3.437860867674477e-207; bool r30462587 = r30462570 <= r30462586; float r30462588 = 2; float r30462589 = r30462588 * r30462573; float r30462590 = r30462566 / r30462565; float r30462591 = r30462565 / r30462590; float r30462592 = r30462588 * r30462591; float r30462593 = r30462564 - r30462592; float r30462594 = pow(r30462567, r30462588); float r30462595 = r30462573 * r30462594; float r30462596 = U_; float r30462597 = r30462575 - r30462596; float r30462598 = r30462595 * r30462597; float r30462599 = r30462593 - r30462598; float r30462600 = r30462575 * r30462599; float r30462601 = r30462589 * r30462600; float r30462602 = sqrt(r30462601); float r30462603 = 9.025178386511885e+265; bool r30462604 = r30462570 <= r30462603; float r30462605 = r30462589 * r30462575; float r30462606 = sqrt(r30462605); float r30462607 = r30462565 * r30462565; float r30462608 = r30462607 / r30462566; float r30462609 = r30462588 * r30462608; float r30462610 = r30462564 - r30462609; float r30462611 = r30462610 - r30462598; float r30462612 = sqrt(r30462611); float r30462613 = r30462606 * r30462612; float r30462614 = r30462604 ? r30462613 : r30462585; float r30462615 = r30462587 ? r30462602 : r30462614; float r30462616 = r30462572 ? r30462585 : r30462615; return r30462616; } double f_od(double n, double U, double t, double l, double Om, double U_) { double r30462617 = t; double r30462618 = l; double r30462619 = Om; double r30462620 = r30462618 / r30462619; double r30462621 = r30462618 + r30462618; double r30462622 = r30462620 * r30462621; double r30462623 = r30462617 - r30462622; double r30462624 = -inf.0; bool r30462625 = r30462623 <= r30462624; double r30462626 = n; double r30462627 = r30462617 * r30462626; double r30462628 = U; double r30462629 = r30462628 + r30462628; double r30462630 = r30462627 * r30462629; double r30462631 = r30462626 / r30462619; double r30462632 = r30462628 * r30462618; double r30462633 = r30462631 * r30462632; double r30462634 = 4; double r30462635 = r30462634 * r30462618; double r30462636 = r30462633 * r30462635; double r30462637 = r30462630 - r30462636; double r30462638 = sqrt(r30462637); double r30462639 = 3.437860867674477e-207; bool r30462640 = r30462623 <= r30462639; double r30462641 = 2; double r30462642 = r30462641 * r30462626; double r30462643 = r30462619 / r30462618; double r30462644 = r30462618 / r30462643; double r30462645 = r30462641 * r30462644; double r30462646 = r30462617 - r30462645; double r30462647 = pow(r30462620, r30462641); double r30462648 = r30462626 * r30462647; double r30462649 = U_; double r30462650 = r30462628 - r30462649; double r30462651 = r30462648 * r30462650; double r30462652 = r30462646 - r30462651; double r30462653 = r30462628 * r30462652; double r30462654 = r30462642 * r30462653; double r30462655 = sqrt(r30462654); double r30462656 = 9.025178386511885e+265; bool r30462657 = r30462623 <= r30462656; double r30462658 = r30462642 * r30462628; double r30462659 = sqrt(r30462658); double r30462660 = r30462618 * r30462618; double r30462661 = r30462660 / r30462619; double r30462662 = r30462641 * r30462661; double r30462663 = r30462617 - r30462662; double r30462664 = r30462663 - r30462651; double r30462665 = sqrt(r30462664); double r30462666 = r30462659 * r30462665; double r30462667 = r30462657 ? r30462666 : r30462638; double r30462668 = r30462640 ? r30462655 : r30462667; double r30462669 = r30462625 ? r30462638 : r30462668; return r30462669; } 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 r30462670, r30462671, r30462672, r30462673, r30462674, r30462675, r30462676, r30462677, r30462678, r30462679, r30462680, r30462681, r30462682, r30462683, r30462684, r30462685, r30462686, r30462687, r30462688, r30462689, r30462690; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r30462670, "2", 10, MPFR_RNDN); mpfr_init(r30462671); mpfr_init(r30462672); mpfr_init(r30462673); mpfr_init(r30462674); mpfr_init(r30462675); mpfr_init(r30462676); mpfr_init(r30462677); mpfr_init(r30462678); mpfr_init(r30462679); mpfr_init(r30462680); mpfr_init(r30462681); mpfr_init(r30462682); mpfr_init(r30462683); mpfr_init(r30462684); mpfr_init(r30462685); mpfr_init(r30462686); mpfr_init(r30462687); mpfr_init(r30462688); mpfr_init(r30462689); mpfr_init(r30462690); } double f_im(double n, double U, double t, double l, double Om, double U_) { ; mpfr_set_d(r30462671, n, MPFR_RNDN); mpfr_mul(r30462672, r30462670, r30462671, MPFR_RNDN); mpfr_set_d(r30462673, U, MPFR_RNDN); mpfr_mul(r30462674, r30462672, r30462673, MPFR_RNDN); mpfr_set_d(r30462675, t, MPFR_RNDN); mpfr_set_d(r30462676, l, MPFR_RNDN); mpfr_mul(r30462677, r30462676, r30462676, MPFR_RNDN); mpfr_set_d(r30462678, Om, MPFR_RNDN); mpfr_div(r30462679, r30462677, r30462678, MPFR_RNDN); mpfr_mul(r30462680, r30462670, r30462679, MPFR_RNDN); mpfr_sub(r30462681, r30462675, r30462680, MPFR_RNDN); mpfr_div(r30462682, r30462676, r30462678, MPFR_RNDN); mpfr_pow(r30462683, r30462682, r30462670, MPFR_RNDN); mpfr_mul(r30462684, r30462671, r30462683, MPFR_RNDN); mpfr_set_d(r30462685, U_, MPFR_RNDN); mpfr_sub(r30462686, r30462673, r30462685, MPFR_RNDN); mpfr_mul(r30462687, r30462684, r30462686, MPFR_RNDN); mpfr_sub(r30462688, r30462681, r30462687, MPFR_RNDN); mpfr_mul(r30462689, r30462674, r30462688, MPFR_RNDN); mpfr_sqrt(r30462690, r30462689, MPFR_RNDN); return mpfr_get_d(r30462690, MPFR_RNDN); } static mpfr_t r30462691, r30462692, r30462693, r30462694, r30462695, r30462696, r30462697, r30462698, r30462699, r30462700, r30462701, r30462702, r30462703, r30462704, r30462705, r30462706, r30462707, r30462708, r30462709, r30462710, r30462711, r30462712, r30462713, r30462714, r30462715, r30462716, r30462717, r30462718, r30462719, r30462720, r30462721, r30462722, r30462723, r30462724, r30462725, r30462726, r30462727, r30462728, r30462729, r30462730, r30462731, r30462732, r30462733, r30462734, r30462735, r30462736, r30462737, r30462738, r30462739, r30462740, r30462741, r30462742, r30462743; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init(r30462691); mpfr_init(r30462692); mpfr_init(r30462693); mpfr_init(r30462694); mpfr_init(r30462695); mpfr_init(r30462696); mpfr_init(r30462697); mpfr_init_set_str(r30462698, "-inf.0", 10, MPFR_RNDN); mpfr_init(r30462699); mpfr_init(r30462700); mpfr_init(r30462701); mpfr_init(r30462702); mpfr_init(r30462703); mpfr_init(r30462704); mpfr_init(r30462705); mpfr_init(r30462706); mpfr_init(r30462707); mpfr_init_set_str(r30462708, "4", 10, MPFR_RNDN); mpfr_init(r30462709); mpfr_init(r30462710); mpfr_init(r30462711); mpfr_init(r30462712); mpfr_init_set_str(r30462713, "3.437860867674477e-207", 10, MPFR_RNDN); mpfr_init(r30462714); mpfr_init_set_str(r30462715, "2", 10, MPFR_RNDN); mpfr_init(r30462716); mpfr_init(r30462717); mpfr_init(r30462718); mpfr_init(r30462719); mpfr_init(r30462720); mpfr_init(r30462721); mpfr_init(r30462722); mpfr_init(r30462723); mpfr_init(r30462724); mpfr_init(r30462725); mpfr_init(r30462726); mpfr_init(r30462727); mpfr_init(r30462728); mpfr_init(r30462729); mpfr_init_set_str(r30462730, "9.025178386511885e+265", 10, MPFR_RNDN); mpfr_init(r30462731); mpfr_init(r30462732); mpfr_init(r30462733); mpfr_init(r30462734); mpfr_init(r30462735); mpfr_init(r30462736); mpfr_init(r30462737); mpfr_init(r30462738); mpfr_init(r30462739); mpfr_init(r30462740); mpfr_init(r30462741); mpfr_init(r30462742); mpfr_init(r30462743); } double f_fm(double n, double U, double t, double l, double Om, double U_) { mpfr_set_d(r30462691, t, MPFR_RNDN); mpfr_set_d(r30462692, l, MPFR_RNDN); mpfr_set_d(r30462693, Om, MPFR_RNDN); mpfr_div(r30462694, r30462692, r30462693, MPFR_RNDN); mpfr_add(r30462695, r30462692, r30462692, MPFR_RNDN); mpfr_mul(r30462696, r30462694, r30462695, MPFR_RNDN); mpfr_sub(r30462697, r30462691, r30462696, MPFR_RNDN); ; mpfr_set_si(r30462699, mpfr_cmp(r30462697, r30462698) <= 0, MPFR_RNDN); mpfr_set_d(r30462700, n, MPFR_RNDN); mpfr_mul(r30462701, r30462691, r30462700, MPFR_RNDN); mpfr_set_d(r30462702, U, MPFR_RNDN); mpfr_add(r30462703, r30462702, r30462702, MPFR_RNDN); mpfr_mul(r30462704, r30462701, r30462703, MPFR_RNDN); mpfr_div(r30462705, r30462700, r30462693, MPFR_RNDN); mpfr_mul(r30462706, r30462702, r30462692, MPFR_RNDN); mpfr_mul(r30462707, r30462705, r30462706, MPFR_RNDN); ; mpfr_mul(r30462709, r30462708, r30462692, MPFR_RNDN); mpfr_mul(r30462710, r30462707, r30462709, MPFR_RNDN); mpfr_sub(r30462711, r30462704, r30462710, MPFR_RNDN); mpfr_sqrt(r30462712, r30462711, MPFR_RNDN); ; mpfr_set_si(r30462714, mpfr_cmp(r30462697, r30462713) <= 0, MPFR_RNDN); ; mpfr_mul(r30462716, r30462715, r30462700, MPFR_RNDN); mpfr_div(r30462717, r30462693, r30462692, MPFR_RNDN); mpfr_div(r30462718, r30462692, r30462717, MPFR_RNDN); mpfr_mul(r30462719, r30462715, r30462718, MPFR_RNDN); mpfr_sub(r30462720, r30462691, r30462719, MPFR_RNDN); mpfr_pow(r30462721, r30462694, r30462715, MPFR_RNDN); mpfr_mul(r30462722, r30462700, r30462721, MPFR_RNDN); mpfr_set_d(r30462723, U_, MPFR_RNDN); mpfr_sub(r30462724, r30462702, r30462723, MPFR_RNDN); mpfr_mul(r30462725, r30462722, r30462724, MPFR_RNDN); mpfr_sub(r30462726, r30462720, r30462725, MPFR_RNDN); mpfr_mul(r30462727, r30462702, r30462726, MPFR_RNDN); mpfr_mul(r30462728, r30462716, r30462727, MPFR_RNDN); mpfr_sqrt(r30462729, r30462728, MPFR_RNDN); ; mpfr_set_si(r30462731, mpfr_cmp(r30462697, r30462730) <= 0, MPFR_RNDN); mpfr_mul(r30462732, r30462716, r30462702, MPFR_RNDN); mpfr_sqrt(r30462733, r30462732, MPFR_RNDN); mpfr_mul(r30462734, r30462692, r30462692, MPFR_RNDN); mpfr_div(r30462735, r30462734, r30462693, MPFR_RNDN); mpfr_mul(r30462736, r30462715, r30462735, MPFR_RNDN); mpfr_sub(r30462737, r30462691, r30462736, MPFR_RNDN); mpfr_sub(r30462738, r30462737, r30462725, MPFR_RNDN); mpfr_sqrt(r30462739, r30462738, MPFR_RNDN); mpfr_mul(r30462740, r30462733, r30462739, MPFR_RNDN); if (mpfr_get_si(r30462731, MPFR_RNDN)) { mpfr_set(r30462741, r30462740, MPFR_RNDN); } else { mpfr_set(r30462741, r30462712, MPFR_RNDN); }; if (mpfr_get_si(r30462714, MPFR_RNDN)) { mpfr_set(r30462742, r30462729, MPFR_RNDN); } else { mpfr_set(r30462742, r30462741, MPFR_RNDN); }; if (mpfr_get_si(r30462699, MPFR_RNDN)) { mpfr_set(r30462743, r30462712, MPFR_RNDN); } else { mpfr_set(r30462743, r30462742, MPFR_RNDN); }; return mpfr_get_d(r30462743, MPFR_RNDN); } static mpfr_t r30462744, r30462745, r30462746, r30462747, r30462748, r30462749, r30462750, r30462751, r30462752, r30462753, r30462754, r30462755, r30462756, r30462757, r30462758, r30462759, r30462760, r30462761, r30462762, r30462763, r30462764, r30462765, r30462766, r30462767, r30462768, r30462769, r30462770, r30462771, r30462772, r30462773, r30462774, r30462775, r30462776, r30462777, r30462778, r30462779, r30462780, r30462781, r30462782, r30462783, r30462784, r30462785, r30462786, r30462787, r30462788, r30462789, r30462790, r30462791, r30462792, r30462793, r30462794, r30462795, r30462796; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init(r30462744); mpfr_init(r30462745); mpfr_init(r30462746); mpfr_init(r30462747); mpfr_init(r30462748); mpfr_init(r30462749); mpfr_init(r30462750); mpfr_init_set_str(r30462751, "-inf.0", 10, MPFR_RNDN); mpfr_init(r30462752); mpfr_init(r30462753); mpfr_init(r30462754); mpfr_init(r30462755); mpfr_init(r30462756); mpfr_init(r30462757); mpfr_init(r30462758); mpfr_init(r30462759); mpfr_init(r30462760); mpfr_init_set_str(r30462761, "4", 10, MPFR_RNDN); mpfr_init(r30462762); mpfr_init(r30462763); mpfr_init(r30462764); mpfr_init(r30462765); mpfr_init_set_str(r30462766, "3.437860867674477e-207", 10, MPFR_RNDN); mpfr_init(r30462767); mpfr_init_set_str(r30462768, "2", 10, MPFR_RNDN); mpfr_init(r30462769); mpfr_init(r30462770); mpfr_init(r30462771); mpfr_init(r30462772); mpfr_init(r30462773); mpfr_init(r30462774); mpfr_init(r30462775); mpfr_init(r30462776); mpfr_init(r30462777); mpfr_init(r30462778); mpfr_init(r30462779); mpfr_init(r30462780); mpfr_init(r30462781); mpfr_init(r30462782); mpfr_init_set_str(r30462783, "9.025178386511885e+265", 10, MPFR_RNDN); mpfr_init(r30462784); mpfr_init(r30462785); mpfr_init(r30462786); mpfr_init(r30462787); mpfr_init(r30462788); mpfr_init(r30462789); mpfr_init(r30462790); mpfr_init(r30462791); mpfr_init(r30462792); mpfr_init(r30462793); mpfr_init(r30462794); mpfr_init(r30462795); mpfr_init(r30462796); } double f_dm(double n, double U, double t, double l, double Om, double U_) { mpfr_set_d(r30462744, t, MPFR_RNDN); mpfr_set_d(r30462745, l, MPFR_RNDN); mpfr_set_d(r30462746, Om, MPFR_RNDN); mpfr_div(r30462747, r30462745, r30462746, MPFR_RNDN); mpfr_add(r30462748, r30462745, r30462745, MPFR_RNDN); mpfr_mul(r30462749, r30462747, r30462748, MPFR_RNDN); mpfr_sub(r30462750, r30462744, r30462749, MPFR_RNDN); ; mpfr_set_si(r30462752, mpfr_cmp(r30462750, r30462751) <= 0, MPFR_RNDN); mpfr_set_d(r30462753, n, MPFR_RNDN); mpfr_mul(r30462754, r30462744, r30462753, MPFR_RNDN); mpfr_set_d(r30462755, U, MPFR_RNDN); mpfr_add(r30462756, r30462755, r30462755, MPFR_RNDN); mpfr_mul(r30462757, r30462754, r30462756, MPFR_RNDN); mpfr_div(r30462758, r30462753, r30462746, MPFR_RNDN); mpfr_mul(r30462759, r30462755, r30462745, MPFR_RNDN); mpfr_mul(r30462760, r30462758, r30462759, MPFR_RNDN); ; mpfr_mul(r30462762, r30462761, r30462745, MPFR_RNDN); mpfr_mul(r30462763, r30462760, r30462762, MPFR_RNDN); mpfr_sub(r30462764, r30462757, r30462763, MPFR_RNDN); mpfr_sqrt(r30462765, r30462764, MPFR_RNDN); ; mpfr_set_si(r30462767, mpfr_cmp(r30462750, r30462766) <= 0, MPFR_RNDN); ; mpfr_mul(r30462769, r30462768, r30462753, MPFR_RNDN); mpfr_div(r30462770, r30462746, r30462745, MPFR_RNDN); mpfr_div(r30462771, r30462745, r30462770, MPFR_RNDN); mpfr_mul(r30462772, r30462768, r30462771, MPFR_RNDN); mpfr_sub(r30462773, r30462744, r30462772, MPFR_RNDN); mpfr_pow(r30462774, r30462747, r30462768, MPFR_RNDN); mpfr_mul(r30462775, r30462753, r30462774, MPFR_RNDN); mpfr_set_d(r30462776, U_, MPFR_RNDN); mpfr_sub(r30462777, r30462755, r30462776, MPFR_RNDN); mpfr_mul(r30462778, r30462775, r30462777, MPFR_RNDN); mpfr_sub(r30462779, r30462773, r30462778, MPFR_RNDN); mpfr_mul(r30462780, r30462755, r30462779, MPFR_RNDN); mpfr_mul(r30462781, r30462769, r30462780, MPFR_RNDN); mpfr_sqrt(r30462782, r30462781, MPFR_RNDN); ; mpfr_set_si(r30462784, mpfr_cmp(r30462750, r30462783) <= 0, MPFR_RNDN); mpfr_mul(r30462785, r30462769, r30462755, MPFR_RNDN); mpfr_sqrt(r30462786, r30462785, MPFR_RNDN); mpfr_mul(r30462787, r30462745, r30462745, MPFR_RNDN); mpfr_div(r30462788, r30462787, r30462746, MPFR_RNDN); mpfr_mul(r30462789, r30462768, r30462788, MPFR_RNDN); mpfr_sub(r30462790, r30462744, r30462789, MPFR_RNDN); mpfr_sub(r30462791, r30462790, r30462778, MPFR_RNDN); mpfr_sqrt(r30462792, r30462791, MPFR_RNDN); mpfr_mul(r30462793, r30462786, r30462792, MPFR_RNDN); if (mpfr_get_si(r30462784, MPFR_RNDN)) { mpfr_set(r30462794, r30462793, MPFR_RNDN); } else { mpfr_set(r30462794, r30462765, MPFR_RNDN); }; if (mpfr_get_si(r30462767, MPFR_RNDN)) { mpfr_set(r30462795, r30462782, MPFR_RNDN); } else { mpfr_set(r30462795, r30462794, MPFR_RNDN); }; if (mpfr_get_si(r30462752, MPFR_RNDN)) { mpfr_set(r30462796, r30462765, MPFR_RNDN); } else { mpfr_set(r30462796, r30462795, MPFR_RNDN); }; return mpfr_get_d(r30462796, MPFR_RNDN); }