#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r20466 = 1; float r20467 = Om; float r20468 = Omc; float r20469 = r20467 / r20468; float r20470 = 2; float r20471 = pow(r20469, r20470); float r20472 = r20466 - r20471; float r20473 = t; float r20474 = l; float r20475 = r20473 / r20474; float r20476 = pow(r20475, r20470); float r20477 = r20470 * r20476; float r20478 = r20466 + r20477; float r20479 = r20472 / r20478; float r20480 = sqrt(r20479); float r20481 = asin(r20480); return r20481; } double f_id(double t, double l, double Om, double Omc) { double r20482 = 1; double r20483 = Om; double r20484 = Omc; double r20485 = r20483 / r20484; double r20486 = 2; double r20487 = pow(r20485, r20486); double r20488 = r20482 - r20487; double r20489 = t; double r20490 = l; double r20491 = r20489 / r20490; double r20492 = pow(r20491, r20486); double r20493 = r20486 * r20492; double r20494 = r20482 + r20493; double r20495 = r20488 / r20494; double r20496 = sqrt(r20495); double r20497 = asin(r20496); return r20497; } double f_of(float t, float l, float Om, float Omc) { float r20498 = 1; float r20499 = Om; float r20500 = Omc; float r20501 = r20499 / r20500; float r20502 = 2; float r20503 = pow(r20501, r20502); float r20504 = r20498 - r20503; float r20505 = t; float r20506 = l; float r20507 = r20505 / r20506; float r20508 = pow(r20507, r20502); float r20509 = r20502 * r20508; float r20510 = r20498 + r20509; float r20511 = r20504 / r20510; float r20512 = sqrt(r20511); float r20513 = asin(r20512); return r20513; } double f_od(double t, double l, double Om, double Omc) { double r20514 = 1; double r20515 = Om; double r20516 = Omc; double r20517 = r20515 / r20516; double r20518 = 2; double r20519 = pow(r20517, r20518); double r20520 = r20514 - r20519; double r20521 = t; double r20522 = l; double r20523 = r20521 / r20522; double r20524 = pow(r20523, r20518); double r20525 = r20518 * r20524; double r20526 = r20514 + r20525; double r20527 = r20520 / r20526; double r20528 = sqrt(r20527); double r20529 = asin(r20528); return r20529; } 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 r20530, r20531, r20532, r20533, r20534, r20535, r20536, r20537, r20538, r20539, r20540, r20541, r20542, r20543, r20544, r20545; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r20530, "1", 10, MPFR_RNDN); mpfr_init(r20531); mpfr_init(r20532); mpfr_init(r20533); mpfr_init_set_str(r20534, "2", 10, MPFR_RNDN); mpfr_init(r20535); mpfr_init(r20536); mpfr_init(r20537); mpfr_init(r20538); mpfr_init(r20539); mpfr_init(r20540); mpfr_init(r20541); mpfr_init(r20542); mpfr_init(r20543); mpfr_init(r20544); mpfr_init(r20545); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20531, Om, MPFR_RNDN); mpfr_set_d(r20532, Omc, MPFR_RNDN); mpfr_div(r20533, r20531, r20532, MPFR_RNDN); ; mpfr_pow(r20535, r20533, r20534, MPFR_RNDN); mpfr_sub(r20536, r20530, r20535, MPFR_RNDN); mpfr_set_d(r20537, t, MPFR_RNDN); mpfr_set_d(r20538, l, MPFR_RNDN); mpfr_div(r20539, r20537, r20538, MPFR_RNDN); mpfr_pow(r20540, r20539, r20534, MPFR_RNDN); mpfr_mul(r20541, r20534, r20540, MPFR_RNDN); mpfr_add(r20542, r20530, r20541, MPFR_RNDN); mpfr_div(r20543, r20536, r20542, MPFR_RNDN); mpfr_sqrt(r20544, r20543, MPFR_RNDN); mpfr_asin(r20545, r20544, MPFR_RNDN); return mpfr_get_d(r20545, MPFR_RNDN); } static mpfr_t r20546, r20547, r20548, r20549, r20550, r20551, r20552, r20553, r20554, r20555, r20556, r20557, r20558, r20559, r20560, r20561; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20546, "1", 10, MPFR_RNDN); mpfr_init(r20547); mpfr_init(r20548); mpfr_init(r20549); mpfr_init_set_str(r20550, "2", 10, MPFR_RNDN); mpfr_init(r20551); mpfr_init(r20552); mpfr_init(r20553); mpfr_init(r20554); mpfr_init(r20555); mpfr_init(r20556); mpfr_init(r20557); mpfr_init(r20558); mpfr_init(r20559); mpfr_init(r20560); mpfr_init(r20561); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20547, Om, MPFR_RNDN); mpfr_set_d(r20548, Omc, MPFR_RNDN); mpfr_div(r20549, r20547, r20548, MPFR_RNDN); ; mpfr_pow(r20551, r20549, r20550, MPFR_RNDN); mpfr_sub(r20552, r20546, r20551, MPFR_RNDN); mpfr_set_d(r20553, t, MPFR_RNDN); mpfr_set_d(r20554, l, MPFR_RNDN); mpfr_div(r20555, r20553, r20554, MPFR_RNDN); mpfr_pow(r20556, r20555, r20550, MPFR_RNDN); mpfr_mul(r20557, r20550, r20556, MPFR_RNDN); mpfr_add(r20558, r20546, r20557, MPFR_RNDN); mpfr_div(r20559, r20552, r20558, MPFR_RNDN); mpfr_sqrt(r20560, r20559, MPFR_RNDN); mpfr_asin(r20561, r20560, MPFR_RNDN); return mpfr_get_d(r20561, MPFR_RNDN); } static mpfr_t r20562, r20563, r20564, r20565, r20566, r20567, r20568, r20569, r20570, r20571, r20572, r20573, r20574, r20575, r20576, r20577; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r20562, "1", 10, MPFR_RNDN); mpfr_init(r20563); mpfr_init(r20564); mpfr_init(r20565); mpfr_init_set_str(r20566, "2", 10, MPFR_RNDN); mpfr_init(r20567); mpfr_init(r20568); mpfr_init(r20569); mpfr_init(r20570); mpfr_init(r20571); mpfr_init(r20572); mpfr_init(r20573); mpfr_init(r20574); mpfr_init(r20575); mpfr_init(r20576); mpfr_init(r20577); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r20563, Om, MPFR_RNDN); mpfr_set_d(r20564, Omc, MPFR_RNDN); mpfr_div(r20565, r20563, r20564, MPFR_RNDN); ; mpfr_pow(r20567, r20565, r20566, MPFR_RNDN); mpfr_sub(r20568, r20562, r20567, MPFR_RNDN); mpfr_set_d(r20569, t, MPFR_RNDN); mpfr_set_d(r20570, l, MPFR_RNDN); mpfr_div(r20571, r20569, r20570, MPFR_RNDN); mpfr_pow(r20572, r20571, r20566, MPFR_RNDN); mpfr_mul(r20573, r20566, r20572, MPFR_RNDN); mpfr_add(r20574, r20562, r20573, MPFR_RNDN); mpfr_div(r20575, r20568, r20574, MPFR_RNDN); mpfr_sqrt(r20576, r20575, MPFR_RNDN); mpfr_asin(r20577, r20576, MPFR_RNDN); return mpfr_get_d(r20577, MPFR_RNDN); }