#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r21459 = 1; float r21460 = Om; float r21461 = Omc; float r21462 = r21460 / r21461; float r21463 = 2; float r21464 = pow(r21462, r21463); float r21465 = r21459 - r21464; float r21466 = t; float r21467 = l; float r21468 = r21466 / r21467; float r21469 = pow(r21468, r21463); float r21470 = r21463 * r21469; float r21471 = r21459 + r21470; float r21472 = r21465 / r21471; float r21473 = sqrt(r21472); float r21474 = asin(r21473); return r21474; } double f_id(double t, double l, double Om, double Omc) { double r21475 = 1; double r21476 = Om; double r21477 = Omc; double r21478 = r21476 / r21477; double r21479 = 2; double r21480 = pow(r21478, r21479); double r21481 = r21475 - r21480; double r21482 = t; double r21483 = l; double r21484 = r21482 / r21483; double r21485 = pow(r21484, r21479); double r21486 = r21479 * r21485; double r21487 = r21475 + r21486; double r21488 = r21481 / r21487; double r21489 = sqrt(r21488); double r21490 = asin(r21489); return r21490; } double f_of(float t, float l, float Om, float Omc) { float r21491 = 1; float r21492 = Om; float r21493 = Omc; float r21494 = r21492 / r21493; float r21495 = 2; float r21496 = pow(r21494, r21495); float r21497 = r21491 - r21496; float r21498 = t; float r21499 = l; float r21500 = r21498 / r21499; float r21501 = pow(r21500, r21495); float r21502 = r21495 * r21501; float r21503 = r21491 + r21502; float r21504 = r21497 / r21503; float r21505 = sqrt(r21504); float r21506 = asin(r21505); return r21506; } double f_od(double t, double l, double Om, double Omc) { double r21507 = 1; double r21508 = Om; double r21509 = Omc; double r21510 = r21508 / r21509; double r21511 = 2; double r21512 = pow(r21510, r21511); double r21513 = r21507 - r21512; double r21514 = t; double r21515 = l; double r21516 = r21514 / r21515; double r21517 = pow(r21516, r21511); double r21518 = r21511 * r21517; double r21519 = r21507 + r21518; double r21520 = r21513 / r21519; double r21521 = sqrt(r21520); double r21522 = asin(r21521); return r21522; } 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 r21523, r21524, r21525, r21526, r21527, r21528, r21529, r21530, r21531, r21532, r21533, r21534, r21535, r21536, r21537, r21538; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r21523, "1", 10, MPFR_RNDN); mpfr_init(r21524); mpfr_init(r21525); mpfr_init(r21526); mpfr_init_set_str(r21527, "2", 10, MPFR_RNDN); mpfr_init(r21528); mpfr_init(r21529); mpfr_init(r21530); mpfr_init(r21531); mpfr_init(r21532); mpfr_init(r21533); mpfr_init(r21534); mpfr_init(r21535); mpfr_init(r21536); mpfr_init(r21537); mpfr_init(r21538); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r21524, Om, MPFR_RNDN); mpfr_set_d(r21525, Omc, MPFR_RNDN); mpfr_div(r21526, r21524, r21525, MPFR_RNDN); ; mpfr_pow(r21528, r21526, r21527, MPFR_RNDN); mpfr_sub(r21529, r21523, r21528, MPFR_RNDN); mpfr_set_d(r21530, t, MPFR_RNDN); mpfr_set_d(r21531, l, MPFR_RNDN); mpfr_div(r21532, r21530, r21531, MPFR_RNDN); mpfr_pow(r21533, r21532, r21527, MPFR_RNDN); mpfr_mul(r21534, r21527, r21533, MPFR_RNDN); mpfr_add(r21535, r21523, r21534, MPFR_RNDN); mpfr_div(r21536, r21529, r21535, MPFR_RNDN); mpfr_sqrt(r21537, r21536, MPFR_RNDN); mpfr_asin(r21538, r21537, MPFR_RNDN); return mpfr_get_d(r21538, MPFR_RNDN); } static mpfr_t r21539, r21540, r21541, r21542, r21543, r21544, r21545, r21546, r21547, r21548, r21549, r21550, r21551, r21552, r21553, r21554; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21539, "1", 10, MPFR_RNDN); mpfr_init(r21540); mpfr_init(r21541); mpfr_init(r21542); mpfr_init_set_str(r21543, "2", 10, MPFR_RNDN); mpfr_init(r21544); mpfr_init(r21545); mpfr_init(r21546); mpfr_init(r21547); mpfr_init(r21548); mpfr_init(r21549); mpfr_init(r21550); mpfr_init(r21551); mpfr_init(r21552); mpfr_init(r21553); mpfr_init(r21554); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r21540, Om, MPFR_RNDN); mpfr_set_d(r21541, Omc, MPFR_RNDN); mpfr_div(r21542, r21540, r21541, MPFR_RNDN); ; mpfr_pow(r21544, r21542, r21543, MPFR_RNDN); mpfr_sub(r21545, r21539, r21544, MPFR_RNDN); mpfr_set_d(r21546, t, MPFR_RNDN); mpfr_set_d(r21547, l, MPFR_RNDN); mpfr_div(r21548, r21546, r21547, MPFR_RNDN); mpfr_pow(r21549, r21548, r21543, MPFR_RNDN); mpfr_mul(r21550, r21543, r21549, MPFR_RNDN); mpfr_add(r21551, r21539, r21550, MPFR_RNDN); mpfr_div(r21552, r21545, r21551, MPFR_RNDN); mpfr_sqrt(r21553, r21552, MPFR_RNDN); mpfr_asin(r21554, r21553, MPFR_RNDN); return mpfr_get_d(r21554, MPFR_RNDN); } static mpfr_t r21555, r21556, r21557, r21558, r21559, r21560, r21561, r21562, r21563, r21564, r21565, r21566, r21567, r21568, r21569, r21570; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r21555, "1", 10, MPFR_RNDN); mpfr_init(r21556); mpfr_init(r21557); mpfr_init(r21558); mpfr_init_set_str(r21559, "2", 10, MPFR_RNDN); mpfr_init(r21560); mpfr_init(r21561); mpfr_init(r21562); mpfr_init(r21563); mpfr_init(r21564); mpfr_init(r21565); mpfr_init(r21566); mpfr_init(r21567); mpfr_init(r21568); mpfr_init(r21569); mpfr_init(r21570); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r21556, Om, MPFR_RNDN); mpfr_set_d(r21557, Omc, MPFR_RNDN); mpfr_div(r21558, r21556, r21557, MPFR_RNDN); ; mpfr_pow(r21560, r21558, r21559, MPFR_RNDN); mpfr_sub(r21561, r21555, r21560, MPFR_RNDN); mpfr_set_d(r21562, t, MPFR_RNDN); mpfr_set_d(r21563, l, MPFR_RNDN); mpfr_div(r21564, r21562, r21563, MPFR_RNDN); mpfr_pow(r21565, r21564, r21559, MPFR_RNDN); mpfr_mul(r21566, r21559, r21565, MPFR_RNDN); mpfr_add(r21567, r21555, r21566, MPFR_RNDN); mpfr_div(r21568, r21561, r21567, MPFR_RNDN); mpfr_sqrt(r21569, r21568, MPFR_RNDN); mpfr_asin(r21570, r21569, MPFR_RNDN); return mpfr_get_d(r21570, MPFR_RNDN); }