#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r22498 = 1; float r22499 = Om; float r22500 = Omc; float r22501 = r22499 / r22500; float r22502 = 2; float r22503 = pow(r22501, r22502); float r22504 = r22498 - r22503; float r22505 = t; float r22506 = l; float r22507 = r22505 / r22506; float r22508 = pow(r22507, r22502); float r22509 = r22502 * r22508; float r22510 = r22498 + r22509; float r22511 = r22504 / r22510; float r22512 = sqrt(r22511); float r22513 = asin(r22512); return r22513; } double f_id(double t, double l, double Om, double Omc) { double r22514 = 1; double r22515 = Om; double r22516 = Omc; double r22517 = r22515 / r22516; double r22518 = 2; double r22519 = pow(r22517, r22518); double r22520 = r22514 - r22519; double r22521 = t; double r22522 = l; double r22523 = r22521 / r22522; double r22524 = pow(r22523, r22518); double r22525 = r22518 * r22524; double r22526 = r22514 + r22525; double r22527 = r22520 / r22526; double r22528 = sqrt(r22527); double r22529 = asin(r22528); return r22529; } double f_of(float t, float l, float Om, float Omc) { float r22530 = 1; float r22531 = Om; float r22532 = Omc; float r22533 = r22531 / r22532; float r22534 = 2; float r22535 = pow(r22533, r22534); float r22536 = r22530 - r22535; float r22537 = t; float r22538 = l; float r22539 = r22537 / r22538; float r22540 = pow(r22539, r22534); float r22541 = r22534 * r22540; float r22542 = r22530 + r22541; float r22543 = r22536 / r22542; float r22544 = sqrt(r22543); float r22545 = asin(r22544); return r22545; } double f_od(double t, double l, double Om, double Omc) { double r22546 = 1; double r22547 = Om; double r22548 = Omc; double r22549 = r22547 / r22548; double r22550 = 2; double r22551 = pow(r22549, r22550); double r22552 = r22546 - r22551; double r22553 = t; double r22554 = l; double r22555 = r22553 / r22554; double r22556 = pow(r22555, r22550); double r22557 = r22550 * r22556; double r22558 = r22546 + r22557; double r22559 = r22552 / r22558; double r22560 = sqrt(r22559); double r22561 = asin(r22560); return r22561; } 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 r22562, r22563, r22564, r22565, r22566, r22567, r22568, r22569, r22570, r22571, r22572, r22573, r22574, r22575, r22576, r22577; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r22562, "1", 10, MPFR_RNDN); mpfr_init(r22563); mpfr_init(r22564); mpfr_init(r22565); mpfr_init_set_str(r22566, "2", 10, MPFR_RNDN); mpfr_init(r22567); mpfr_init(r22568); mpfr_init(r22569); mpfr_init(r22570); mpfr_init(r22571); mpfr_init(r22572); mpfr_init(r22573); mpfr_init(r22574); mpfr_init(r22575); mpfr_init(r22576); mpfr_init(r22577); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r22563, Om, MPFR_RNDN); mpfr_set_d(r22564, Omc, MPFR_RNDN); mpfr_div(r22565, r22563, r22564, MPFR_RNDN); ; mpfr_pow(r22567, r22565, r22566, MPFR_RNDN); mpfr_sub(r22568, r22562, r22567, MPFR_RNDN); mpfr_set_d(r22569, t, MPFR_RNDN); mpfr_set_d(r22570, l, MPFR_RNDN); mpfr_div(r22571, r22569, r22570, MPFR_RNDN); mpfr_pow(r22572, r22571, r22566, MPFR_RNDN); mpfr_mul(r22573, r22566, r22572, MPFR_RNDN); mpfr_add(r22574, r22562, r22573, MPFR_RNDN); mpfr_div(r22575, r22568, r22574, MPFR_RNDN); mpfr_sqrt(r22576, r22575, MPFR_RNDN); mpfr_asin(r22577, r22576, MPFR_RNDN); return mpfr_get_d(r22577, MPFR_RNDN); } static mpfr_t r22578, r22579, r22580, r22581, r22582, r22583, r22584, r22585, r22586, r22587, r22588, r22589, r22590, r22591, r22592, r22593; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22578, "1", 10, MPFR_RNDN); mpfr_init(r22579); mpfr_init(r22580); mpfr_init(r22581); mpfr_init_set_str(r22582, "2", 10, MPFR_RNDN); mpfr_init(r22583); mpfr_init(r22584); mpfr_init(r22585); mpfr_init(r22586); mpfr_init(r22587); mpfr_init(r22588); mpfr_init(r22589); mpfr_init(r22590); mpfr_init(r22591); mpfr_init(r22592); mpfr_init(r22593); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r22579, Om, MPFR_RNDN); mpfr_set_d(r22580, Omc, MPFR_RNDN); mpfr_div(r22581, r22579, r22580, MPFR_RNDN); ; mpfr_pow(r22583, r22581, r22582, MPFR_RNDN); mpfr_sub(r22584, r22578, r22583, MPFR_RNDN); mpfr_set_d(r22585, t, MPFR_RNDN); mpfr_set_d(r22586, l, MPFR_RNDN); mpfr_div(r22587, r22585, r22586, MPFR_RNDN); mpfr_pow(r22588, r22587, r22582, MPFR_RNDN); mpfr_mul(r22589, r22582, r22588, MPFR_RNDN); mpfr_add(r22590, r22578, r22589, MPFR_RNDN); mpfr_div(r22591, r22584, r22590, MPFR_RNDN); mpfr_sqrt(r22592, r22591, MPFR_RNDN); mpfr_asin(r22593, r22592, MPFR_RNDN); return mpfr_get_d(r22593, MPFR_RNDN); } static mpfr_t r22594, r22595, r22596, r22597, r22598, r22599, r22600, r22601, r22602, r22603, r22604, r22605, r22606, r22607, r22608, r22609; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r22594, "1", 10, MPFR_RNDN); mpfr_init(r22595); mpfr_init(r22596); mpfr_init(r22597); mpfr_init_set_str(r22598, "2", 10, MPFR_RNDN); mpfr_init(r22599); mpfr_init(r22600); mpfr_init(r22601); mpfr_init(r22602); mpfr_init(r22603); mpfr_init(r22604); mpfr_init(r22605); mpfr_init(r22606); mpfr_init(r22607); mpfr_init(r22608); mpfr_init(r22609); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r22595, Om, MPFR_RNDN); mpfr_set_d(r22596, Omc, MPFR_RNDN); mpfr_div(r22597, r22595, r22596, MPFR_RNDN); ; mpfr_pow(r22599, r22597, r22598, MPFR_RNDN); mpfr_sub(r22600, r22594, r22599, MPFR_RNDN); mpfr_set_d(r22601, t, MPFR_RNDN); mpfr_set_d(r22602, l, MPFR_RNDN); mpfr_div(r22603, r22601, r22602, MPFR_RNDN); mpfr_pow(r22604, r22603, r22598, MPFR_RNDN); mpfr_mul(r22605, r22598, r22604, MPFR_RNDN); mpfr_add(r22606, r22594, r22605, MPFR_RNDN); mpfr_div(r22607, r22600, r22606, MPFR_RNDN); mpfr_sqrt(r22608, r22607, MPFR_RNDN); mpfr_asin(r22609, r22608, MPFR_RNDN); return mpfr_get_d(r22609, MPFR_RNDN); }