#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r8615 = 1; float r8616 = Om; float r8617 = Omc; float r8618 = r8616 / r8617; float r8619 = 2; float r8620 = pow(r8618, r8619); float r8621 = r8615 - r8620; float r8622 = t; float r8623 = l; float r8624 = r8622 / r8623; float r8625 = pow(r8624, r8619); float r8626 = r8619 * r8625; float r8627 = r8615 + r8626; float r8628 = r8621 / r8627; float r8629 = sqrt(r8628); float r8630 = asin(r8629); return r8630; } double f_id(double t, double l, double Om, double Omc) { double r8631 = 1; double r8632 = Om; double r8633 = Omc; double r8634 = r8632 / r8633; double r8635 = 2; double r8636 = pow(r8634, r8635); double r8637 = r8631 - r8636; double r8638 = t; double r8639 = l; double r8640 = r8638 / r8639; double r8641 = pow(r8640, r8635); double r8642 = r8635 * r8641; double r8643 = r8631 + r8642; double r8644 = r8637 / r8643; double r8645 = sqrt(r8644); double r8646 = asin(r8645); return r8646; } double f_of(float t, float l, float Om, float Omc) { float r8647 = 1; float r8648 = Om; float r8649 = Omc; float r8650 = r8648 / r8649; float r8651 = 2; float r8652 = pow(r8650, r8651); float r8653 = r8647 - r8652; float r8654 = t; float r8655 = l; float r8656 = r8654 / r8655; float r8657 = pow(r8656, r8651); float r8658 = r8651 * r8657; float r8659 = r8647 + r8658; float r8660 = r8653 / r8659; float r8661 = sqrt(r8660); float r8662 = asin(r8661); return r8662; } double f_od(double t, double l, double Om, double Omc) { double r8663 = 1; double r8664 = Om; double r8665 = Omc; double r8666 = r8664 / r8665; double r8667 = 2; double r8668 = pow(r8666, r8667); double r8669 = r8663 - r8668; double r8670 = t; double r8671 = l; double r8672 = r8670 / r8671; double r8673 = pow(r8672, r8667); double r8674 = r8667 * r8673; double r8675 = r8663 + r8674; double r8676 = r8669 / r8675; double r8677 = sqrt(r8676); double r8678 = asin(r8677); return r8678; } 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 r8679, r8680, r8681, r8682, r8683, r8684, r8685, r8686, r8687, r8688, r8689, r8690, r8691, r8692, r8693, r8694; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r8679, "1", 10, MPFR_RNDN); mpfr_init(r8680); mpfr_init(r8681); mpfr_init(r8682); mpfr_init_set_str(r8683, "2", 10, MPFR_RNDN); mpfr_init(r8684); mpfr_init(r8685); mpfr_init(r8686); mpfr_init(r8687); mpfr_init(r8688); mpfr_init(r8689); mpfr_init(r8690); mpfr_init(r8691); mpfr_init(r8692); mpfr_init(r8693); mpfr_init(r8694); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8680, Om, MPFR_RNDN); mpfr_set_d(r8681, Omc, MPFR_RNDN); mpfr_div(r8682, r8680, r8681, MPFR_RNDN); ; mpfr_pow(r8684, r8682, r8683, MPFR_RNDN); mpfr_sub(r8685, r8679, r8684, MPFR_RNDN); mpfr_set_d(r8686, t, MPFR_RNDN); mpfr_set_d(r8687, l, MPFR_RNDN); mpfr_div(r8688, r8686, r8687, MPFR_RNDN); mpfr_pow(r8689, r8688, r8683, MPFR_RNDN); mpfr_mul(r8690, r8683, r8689, MPFR_RNDN); mpfr_add(r8691, r8679, r8690, MPFR_RNDN); mpfr_div(r8692, r8685, r8691, MPFR_RNDN); mpfr_sqrt(r8693, r8692, MPFR_RNDN); mpfr_asin(r8694, r8693, MPFR_RNDN); return mpfr_get_d(r8694, MPFR_RNDN); } static mpfr_t r8695, r8696, r8697, r8698, r8699, r8700, r8701, r8702, r8703, r8704, r8705, r8706, r8707, r8708, r8709, r8710; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8695, "1", 10, MPFR_RNDN); mpfr_init(r8696); mpfr_init(r8697); mpfr_init(r8698); mpfr_init_set_str(r8699, "2", 10, MPFR_RNDN); mpfr_init(r8700); mpfr_init(r8701); mpfr_init(r8702); mpfr_init(r8703); mpfr_init(r8704); mpfr_init(r8705); mpfr_init(r8706); mpfr_init(r8707); mpfr_init(r8708); mpfr_init(r8709); mpfr_init(r8710); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8696, Om, MPFR_RNDN); mpfr_set_d(r8697, Omc, MPFR_RNDN); mpfr_div(r8698, r8696, r8697, MPFR_RNDN); ; mpfr_pow(r8700, r8698, r8699, MPFR_RNDN); mpfr_sub(r8701, r8695, r8700, MPFR_RNDN); mpfr_set_d(r8702, t, MPFR_RNDN); mpfr_set_d(r8703, l, MPFR_RNDN); mpfr_div(r8704, r8702, r8703, MPFR_RNDN); mpfr_pow(r8705, r8704, r8699, MPFR_RNDN); mpfr_mul(r8706, r8699, r8705, MPFR_RNDN); mpfr_add(r8707, r8695, r8706, MPFR_RNDN); mpfr_div(r8708, r8701, r8707, MPFR_RNDN); mpfr_sqrt(r8709, r8708, MPFR_RNDN); mpfr_asin(r8710, r8709, MPFR_RNDN); return mpfr_get_d(r8710, MPFR_RNDN); } static mpfr_t r8711, r8712, r8713, r8714, r8715, r8716, r8717, r8718, r8719, r8720, r8721, r8722, r8723, r8724, r8725, r8726; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8711, "1", 10, MPFR_RNDN); mpfr_init(r8712); mpfr_init(r8713); mpfr_init(r8714); mpfr_init_set_str(r8715, "2", 10, MPFR_RNDN); mpfr_init(r8716); mpfr_init(r8717); mpfr_init(r8718); mpfr_init(r8719); mpfr_init(r8720); mpfr_init(r8721); mpfr_init(r8722); mpfr_init(r8723); mpfr_init(r8724); mpfr_init(r8725); mpfr_init(r8726); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8712, Om, MPFR_RNDN); mpfr_set_d(r8713, Omc, MPFR_RNDN); mpfr_div(r8714, r8712, r8713, MPFR_RNDN); ; mpfr_pow(r8716, r8714, r8715, MPFR_RNDN); mpfr_sub(r8717, r8711, r8716, MPFR_RNDN); mpfr_set_d(r8718, t, MPFR_RNDN); mpfr_set_d(r8719, l, MPFR_RNDN); mpfr_div(r8720, r8718, r8719, MPFR_RNDN); mpfr_pow(r8721, r8720, r8715, MPFR_RNDN); mpfr_mul(r8722, r8715, r8721, MPFR_RNDN); mpfr_add(r8723, r8711, r8722, MPFR_RNDN); mpfr_div(r8724, r8717, r8723, MPFR_RNDN); mpfr_sqrt(r8725, r8724, MPFR_RNDN); mpfr_asin(r8726, r8725, MPFR_RNDN); return mpfr_get_d(r8726, MPFR_RNDN); }