#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r20570 = 1; float r20571 = Om; float r20572 = Omc; float r20573 = r20571 / r20572; float r20574 = 2; float r20575 = pow(r20573, r20574); float r20576 = r20570 - r20575; float r20577 = t; float r20578 = l; float r20579 = r20577 / r20578; float r20580 = pow(r20579, r20574); float r20581 = r20574 * r20580; float r20582 = r20570 + r20581; float r20583 = r20576 / r20582; float r20584 = sqrt(r20583); float r20585 = asin(r20584); return r20585; } double f_id(double t, double l, double Om, double Omc) { double r20586 = 1; double r20587 = Om; double r20588 = Omc; double r20589 = r20587 / r20588; double r20590 = 2; double r20591 = pow(r20589, r20590); double r20592 = r20586 - r20591; double r20593 = t; double r20594 = l; double r20595 = r20593 / r20594; double r20596 = pow(r20595, r20590); double r20597 = r20590 * r20596; double r20598 = r20586 + r20597; double r20599 = r20592 / r20598; double r20600 = sqrt(r20599); double r20601 = asin(r20600); return r20601; } double f_of(float t, float l, float Om, float Omc) { float r20602 = 1; float r20603 = Om; float r20604 = Omc; float r20605 = r20603 / r20604; float r20606 = 2; float r20607 = pow(r20605, r20606); float r20608 = r20602 - r20607; float r20609 = t; float r20610 = l; float r20611 = r20609 / r20610; float r20612 = pow(r20611, r20606); float r20613 = r20606 * r20612; float r20614 = r20602 + r20613; float r20615 = r20608 / r20614; float r20616 = sqrt(r20615); float r20617 = asin(r20616); return r20617; } double f_od(double t, double l, double Om, double Omc) { double r20618 = 1; double r20619 = Om; double r20620 = Omc; double r20621 = r20619 / r20620; double r20622 = 2; double r20623 = pow(r20621, r20622); double r20624 = r20618 - r20623; double r20625 = t; double r20626 = l; double r20627 = r20625 / r20626; double r20628 = pow(r20627, r20622); double r20629 = r20622 * r20628; double r20630 = r20618 + r20629; double r20631 = r20624 / r20630; double r20632 = sqrt(r20631); double r20633 = asin(r20632); return r20633; } 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 r20634, r20635, r20636, r20637, r20638, r20639, r20640, r20641, r20642, r20643, r20644, r20645, r20646, r20647, r20648, r20649; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20634, "1", 10, MPFR_RNDN); mpfr_init(r20635); mpfr_init(r20636); mpfr_init(r20637); mpfr_init_set_str(r20638, "2", 10, MPFR_RNDN); mpfr_init(r20639); mpfr_init(r20640); mpfr_init(r20641); mpfr_init(r20642); mpfr_init(r20643); mpfr_init(r20644); mpfr_init(r20645); mpfr_init(r20646); mpfr_init(r20647); mpfr_init(r20648); mpfr_init(r20649); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20635, Om, MPFR_RNDN); mpfr_set_d(r20636, Omc, MPFR_RNDN); mpfr_div(r20637, r20635, r20636, MPFR_RNDN); ; mpfr_pow(r20639, r20637, r20638, MPFR_RNDN); mpfr_sub(r20640, r20634, r20639, MPFR_RNDN); mpfr_set_d(r20641, t, MPFR_RNDN); mpfr_set_d(r20642, l, MPFR_RNDN); mpfr_div(r20643, r20641, r20642, MPFR_RNDN); mpfr_pow(r20644, r20643, r20638, MPFR_RNDN); mpfr_mul(r20645, r20638, r20644, MPFR_RNDN); mpfr_add(r20646, r20634, r20645, MPFR_RNDN); mpfr_div(r20647, r20640, r20646, MPFR_RNDN); mpfr_sqrt(r20648, r20647, MPFR_RNDN); mpfr_asin(r20649, r20648, MPFR_RNDN); return mpfr_get_d(r20649, MPFR_RNDN); } static mpfr_t r20650, r20651, r20652, r20653, r20654, r20655, r20656, r20657, r20658, r20659, r20660, r20661, r20662, r20663, r20664, r20665; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20650, "1", 10, MPFR_RNDN); mpfr_init(r20651); mpfr_init(r20652); mpfr_init(r20653); mpfr_init_set_str(r20654, "2", 10, MPFR_RNDN); mpfr_init(r20655); mpfr_init(r20656); mpfr_init(r20657); mpfr_init(r20658); mpfr_init(r20659); mpfr_init(r20660); mpfr_init(r20661); mpfr_init(r20662); mpfr_init(r20663); mpfr_init(r20664); mpfr_init(r20665); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20651, Om, MPFR_RNDN); mpfr_set_d(r20652, Omc, MPFR_RNDN); mpfr_div(r20653, r20651, r20652, MPFR_RNDN); ; mpfr_pow(r20655, r20653, r20654, MPFR_RNDN); mpfr_sub(r20656, r20650, r20655, MPFR_RNDN); mpfr_set_d(r20657, t, MPFR_RNDN); mpfr_set_d(r20658, l, MPFR_RNDN); mpfr_div(r20659, r20657, r20658, MPFR_RNDN); mpfr_pow(r20660, r20659, r20654, MPFR_RNDN); mpfr_mul(r20661, r20654, r20660, MPFR_RNDN); mpfr_add(r20662, r20650, r20661, MPFR_RNDN); mpfr_div(r20663, r20656, r20662, MPFR_RNDN); mpfr_sqrt(r20664, r20663, MPFR_RNDN); mpfr_asin(r20665, r20664, MPFR_RNDN); return mpfr_get_d(r20665, MPFR_RNDN); } static mpfr_t r20666, r20667, r20668, r20669, r20670, r20671, r20672, r20673, r20674, r20675, r20676, r20677, r20678, r20679, r20680, r20681; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20666, "1", 10, MPFR_RNDN); mpfr_init(r20667); mpfr_init(r20668); mpfr_init(r20669); mpfr_init_set_str(r20670, "2", 10, MPFR_RNDN); mpfr_init(r20671); mpfr_init(r20672); mpfr_init(r20673); mpfr_init(r20674); mpfr_init(r20675); mpfr_init(r20676); mpfr_init(r20677); mpfr_init(r20678); mpfr_init(r20679); mpfr_init(r20680); mpfr_init(r20681); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20667, Om, MPFR_RNDN); mpfr_set_d(r20668, Omc, MPFR_RNDN); mpfr_div(r20669, r20667, r20668, MPFR_RNDN); ; mpfr_pow(r20671, r20669, r20670, MPFR_RNDN); mpfr_sub(r20672, r20666, r20671, MPFR_RNDN); mpfr_set_d(r20673, t, MPFR_RNDN); mpfr_set_d(r20674, l, MPFR_RNDN); mpfr_div(r20675, r20673, r20674, MPFR_RNDN); mpfr_pow(r20676, r20675, r20670, MPFR_RNDN); mpfr_mul(r20677, r20670, r20676, MPFR_RNDN); mpfr_add(r20678, r20666, r20677, MPFR_RNDN); mpfr_div(r20679, r20672, r20678, MPFR_RNDN); mpfr_sqrt(r20680, r20679, MPFR_RNDN); mpfr_asin(r20681, r20680, MPFR_RNDN); return mpfr_get_d(r20681, MPFR_RNDN); }