#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r24713 = 1; float r24714 = Om; float r24715 = Omc; float r24716 = r24714 / r24715; float r24717 = 2; float r24718 = pow(r24716, r24717); float r24719 = r24713 - r24718; float r24720 = t; float r24721 = l; float r24722 = r24720 / r24721; float r24723 = pow(r24722, r24717); float r24724 = r24717 * r24723; float r24725 = r24713 + r24724; float r24726 = r24719 / r24725; float r24727 = sqrt(r24726); float r24728 = asin(r24727); return r24728; } double f_id(double t, double l, double Om, double Omc) { double r24729 = 1; double r24730 = Om; double r24731 = Omc; double r24732 = r24730 / r24731; double r24733 = 2; double r24734 = pow(r24732, r24733); double r24735 = r24729 - r24734; double r24736 = t; double r24737 = l; double r24738 = r24736 / r24737; double r24739 = pow(r24738, r24733); double r24740 = r24733 * r24739; double r24741 = r24729 + r24740; double r24742 = r24735 / r24741; double r24743 = sqrt(r24742); double r24744 = asin(r24743); return r24744; } double f_of(float t, float l, float Om, float Omc) { float r24745 = 1; float r24746 = Om; float r24747 = Omc; float r24748 = r24746 / r24747; float r24749 = 2; float r24750 = pow(r24748, r24749); float r24751 = r24745 - r24750; float r24752 = t; float r24753 = l; float r24754 = r24752 / r24753; float r24755 = pow(r24754, r24749); float r24756 = r24749 * r24755; float r24757 = r24745 + r24756; float r24758 = r24751 / r24757; float r24759 = sqrt(r24758); float r24760 = asin(r24759); return r24760; } double f_od(double t, double l, double Om, double Omc) { double r24761 = 1; double r24762 = Om; double r24763 = Omc; double r24764 = r24762 / r24763; double r24765 = 2; double r24766 = pow(r24764, r24765); double r24767 = r24761 - r24766; double r24768 = t; double r24769 = l; double r24770 = r24768 / r24769; double r24771 = pow(r24770, r24765); double r24772 = r24765 * r24771; double r24773 = r24761 + r24772; double r24774 = r24767 / r24773; double r24775 = sqrt(r24774); double r24776 = asin(r24775); return r24776; } 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 r24777, r24778, r24779, r24780, r24781, r24782, r24783, r24784, r24785, r24786, r24787, r24788, r24789, r24790, r24791, r24792; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24777, "1", 10, MPFR_RNDN); mpfr_init(r24778); mpfr_init(r24779); mpfr_init(r24780); mpfr_init_set_str(r24781, "2", 10, MPFR_RNDN); mpfr_init(r24782); mpfr_init(r24783); mpfr_init(r24784); mpfr_init(r24785); mpfr_init(r24786); mpfr_init(r24787); mpfr_init(r24788); mpfr_init(r24789); mpfr_init(r24790); mpfr_init(r24791); mpfr_init(r24792); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24778, Om, MPFR_RNDN); mpfr_set_d(r24779, Omc, MPFR_RNDN); mpfr_div(r24780, r24778, r24779, MPFR_RNDN); ; mpfr_pow(r24782, r24780, r24781, MPFR_RNDN); mpfr_sub(r24783, r24777, r24782, MPFR_RNDN); mpfr_set_d(r24784, t, MPFR_RNDN); mpfr_set_d(r24785, l, MPFR_RNDN); mpfr_div(r24786, r24784, r24785, MPFR_RNDN); mpfr_pow(r24787, r24786, r24781, MPFR_RNDN); mpfr_mul(r24788, r24781, r24787, MPFR_RNDN); mpfr_add(r24789, r24777, r24788, MPFR_RNDN); mpfr_div(r24790, r24783, r24789, MPFR_RNDN); mpfr_sqrt(r24791, r24790, MPFR_RNDN); mpfr_asin(r24792, r24791, MPFR_RNDN); return mpfr_get_d(r24792, MPFR_RNDN); } static mpfr_t r24793, r24794, r24795, r24796, r24797, r24798, r24799, r24800, r24801, r24802, r24803, r24804, r24805, r24806, r24807, r24808; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24793, "1", 10, MPFR_RNDN); mpfr_init(r24794); mpfr_init(r24795); mpfr_init(r24796); mpfr_init_set_str(r24797, "2", 10, MPFR_RNDN); mpfr_init(r24798); mpfr_init(r24799); mpfr_init(r24800); mpfr_init(r24801); mpfr_init(r24802); mpfr_init(r24803); mpfr_init(r24804); mpfr_init(r24805); mpfr_init(r24806); mpfr_init(r24807); mpfr_init(r24808); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24794, Om, MPFR_RNDN); mpfr_set_d(r24795, Omc, MPFR_RNDN); mpfr_div(r24796, r24794, r24795, MPFR_RNDN); ; mpfr_pow(r24798, r24796, r24797, MPFR_RNDN); mpfr_sub(r24799, r24793, r24798, MPFR_RNDN); mpfr_set_d(r24800, t, MPFR_RNDN); mpfr_set_d(r24801, l, MPFR_RNDN); mpfr_div(r24802, r24800, r24801, MPFR_RNDN); mpfr_pow(r24803, r24802, r24797, MPFR_RNDN); mpfr_mul(r24804, r24797, r24803, MPFR_RNDN); mpfr_add(r24805, r24793, r24804, MPFR_RNDN); mpfr_div(r24806, r24799, r24805, MPFR_RNDN); mpfr_sqrt(r24807, r24806, MPFR_RNDN); mpfr_asin(r24808, r24807, MPFR_RNDN); return mpfr_get_d(r24808, MPFR_RNDN); } static mpfr_t r24809, r24810, r24811, r24812, r24813, r24814, r24815, r24816, r24817, r24818, r24819, r24820, r24821, r24822, r24823, r24824; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24809, "1", 10, MPFR_RNDN); mpfr_init(r24810); mpfr_init(r24811); mpfr_init(r24812); mpfr_init_set_str(r24813, "2", 10, MPFR_RNDN); mpfr_init(r24814); mpfr_init(r24815); mpfr_init(r24816); mpfr_init(r24817); mpfr_init(r24818); mpfr_init(r24819); mpfr_init(r24820); mpfr_init(r24821); mpfr_init(r24822); mpfr_init(r24823); mpfr_init(r24824); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24810, Om, MPFR_RNDN); mpfr_set_d(r24811, Omc, MPFR_RNDN); mpfr_div(r24812, r24810, r24811, MPFR_RNDN); ; mpfr_pow(r24814, r24812, r24813, MPFR_RNDN); mpfr_sub(r24815, r24809, r24814, MPFR_RNDN); mpfr_set_d(r24816, t, MPFR_RNDN); mpfr_set_d(r24817, l, MPFR_RNDN); mpfr_div(r24818, r24816, r24817, MPFR_RNDN); mpfr_pow(r24819, r24818, r24813, MPFR_RNDN); mpfr_mul(r24820, r24813, r24819, MPFR_RNDN); mpfr_add(r24821, r24809, r24820, MPFR_RNDN); mpfr_div(r24822, r24815, r24821, MPFR_RNDN); mpfr_sqrt(r24823, r24822, MPFR_RNDN); mpfr_asin(r24824, r24823, MPFR_RNDN); return mpfr_get_d(r24824, MPFR_RNDN); }