#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r20493 = 1; float r20494 = Om; float r20495 = Omc; float r20496 = r20494 / r20495; float r20497 = 2; float r20498 = pow(r20496, r20497); float r20499 = r20493 - r20498; float r20500 = t; float r20501 = l; float r20502 = r20500 / r20501; float r20503 = pow(r20502, r20497); float r20504 = r20497 * r20503; float r20505 = r20493 + r20504; float r20506 = r20499 / r20505; float r20507 = sqrt(r20506); float r20508 = asin(r20507); return r20508; } double f_id(double t, double l, double Om, double Omc) { double r20509 = 1; double r20510 = Om; double r20511 = Omc; double r20512 = r20510 / r20511; double r20513 = 2; double r20514 = pow(r20512, r20513); double r20515 = r20509 - r20514; double r20516 = t; double r20517 = l; double r20518 = r20516 / r20517; double r20519 = pow(r20518, r20513); double r20520 = r20513 * r20519; double r20521 = r20509 + r20520; double r20522 = r20515 / r20521; double r20523 = sqrt(r20522); double r20524 = asin(r20523); return r20524; } double f_of(float t, float l, float Om, float Omc) { float r20525 = 1; float r20526 = Om; float r20527 = Omc; float r20528 = r20526 / r20527; float r20529 = 2; float r20530 = pow(r20528, r20529); float r20531 = r20525 - r20530; float r20532 = t; float r20533 = l; float r20534 = r20532 / r20533; float r20535 = pow(r20534, r20529); float r20536 = r20529 * r20535; float r20537 = r20525 + r20536; float r20538 = r20531 / r20537; float r20539 = sqrt(r20538); float r20540 = asin(r20539); return r20540; } double f_od(double t, double l, double Om, double Omc) { double r20541 = 1; double r20542 = Om; double r20543 = Omc; double r20544 = r20542 / r20543; double r20545 = 2; double r20546 = pow(r20544, r20545); double r20547 = r20541 - r20546; double r20548 = t; double r20549 = l; double r20550 = r20548 / r20549; double r20551 = pow(r20550, r20545); double r20552 = r20545 * r20551; double r20553 = r20541 + r20552; double r20554 = r20547 / r20553; double r20555 = sqrt(r20554); double r20556 = asin(r20555); return r20556; } 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 r20557, r20558, r20559, r20560, r20561, r20562, r20563, r20564, r20565, r20566, r20567, r20568, r20569, r20570, r20571, r20572; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20557, "1", 10, MPFR_RNDN); mpfr_init(r20558); mpfr_init(r20559); mpfr_init(r20560); mpfr_init_set_str(r20561, "2", 10, MPFR_RNDN); mpfr_init(r20562); mpfr_init(r20563); mpfr_init(r20564); mpfr_init(r20565); mpfr_init(r20566); mpfr_init(r20567); mpfr_init(r20568); mpfr_init(r20569); mpfr_init(r20570); mpfr_init(r20571); mpfr_init(r20572); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20558, Om, MPFR_RNDN); mpfr_set_d(r20559, Omc, MPFR_RNDN); mpfr_div(r20560, r20558, r20559, MPFR_RNDN); ; mpfr_pow(r20562, r20560, r20561, MPFR_RNDN); mpfr_sub(r20563, r20557, r20562, MPFR_RNDN); mpfr_set_d(r20564, t, MPFR_RNDN); mpfr_set_d(r20565, l, MPFR_RNDN); mpfr_div(r20566, r20564, r20565, MPFR_RNDN); mpfr_pow(r20567, r20566, r20561, MPFR_RNDN); mpfr_mul(r20568, r20561, r20567, MPFR_RNDN); mpfr_add(r20569, r20557, r20568, MPFR_RNDN); mpfr_div(r20570, r20563, r20569, MPFR_RNDN); mpfr_sqrt(r20571, r20570, MPFR_RNDN); mpfr_asin(r20572, r20571, MPFR_RNDN); return mpfr_get_d(r20572, MPFR_RNDN); } static mpfr_t r20573, r20574, r20575, r20576, r20577, r20578, r20579, r20580, r20581, r20582, r20583, r20584, r20585, r20586, r20587, r20588; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20573, "1", 10, MPFR_RNDN); mpfr_init(r20574); mpfr_init(r20575); mpfr_init(r20576); mpfr_init_set_str(r20577, "2", 10, MPFR_RNDN); mpfr_init(r20578); mpfr_init(r20579); mpfr_init(r20580); mpfr_init(r20581); mpfr_init(r20582); mpfr_init(r20583); mpfr_init(r20584); mpfr_init(r20585); mpfr_init(r20586); mpfr_init(r20587); mpfr_init(r20588); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20574, Om, MPFR_RNDN); mpfr_set_d(r20575, Omc, MPFR_RNDN); mpfr_div(r20576, r20574, r20575, MPFR_RNDN); ; mpfr_pow(r20578, r20576, r20577, MPFR_RNDN); mpfr_sub(r20579, r20573, r20578, MPFR_RNDN); mpfr_set_d(r20580, t, MPFR_RNDN); mpfr_set_d(r20581, l, MPFR_RNDN); mpfr_div(r20582, r20580, r20581, MPFR_RNDN); mpfr_pow(r20583, r20582, r20577, MPFR_RNDN); mpfr_mul(r20584, r20577, r20583, MPFR_RNDN); mpfr_add(r20585, r20573, r20584, MPFR_RNDN); mpfr_div(r20586, r20579, r20585, MPFR_RNDN); mpfr_sqrt(r20587, r20586, MPFR_RNDN); mpfr_asin(r20588, r20587, MPFR_RNDN); return mpfr_get_d(r20588, MPFR_RNDN); } static mpfr_t r20589, r20590, r20591, r20592, r20593, r20594, r20595, r20596, r20597, r20598, r20599, r20600, r20601, r20602, r20603, r20604; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20589, "1", 10, MPFR_RNDN); mpfr_init(r20590); mpfr_init(r20591); mpfr_init(r20592); mpfr_init_set_str(r20593, "2", 10, MPFR_RNDN); mpfr_init(r20594); mpfr_init(r20595); mpfr_init(r20596); mpfr_init(r20597); mpfr_init(r20598); mpfr_init(r20599); mpfr_init(r20600); mpfr_init(r20601); mpfr_init(r20602); mpfr_init(r20603); mpfr_init(r20604); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20590, Om, MPFR_RNDN); mpfr_set_d(r20591, Omc, MPFR_RNDN); mpfr_div(r20592, r20590, r20591, MPFR_RNDN); ; mpfr_pow(r20594, r20592, r20593, MPFR_RNDN); mpfr_sub(r20595, r20589, r20594, MPFR_RNDN); mpfr_set_d(r20596, t, MPFR_RNDN); mpfr_set_d(r20597, l, MPFR_RNDN); mpfr_div(r20598, r20596, r20597, MPFR_RNDN); mpfr_pow(r20599, r20598, r20593, MPFR_RNDN); mpfr_mul(r20600, r20593, r20599, MPFR_RNDN); mpfr_add(r20601, r20589, r20600, MPFR_RNDN); mpfr_div(r20602, r20595, r20601, MPFR_RNDN); mpfr_sqrt(r20603, r20602, MPFR_RNDN); mpfr_asin(r20604, r20603, MPFR_RNDN); return mpfr_get_d(r20604, MPFR_RNDN); }