#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r24793 = 1; float r24794 = Om; float r24795 = Omc; float r24796 = r24794 / r24795; float r24797 = 2; float r24798 = pow(r24796, r24797); float r24799 = r24793 - r24798; float r24800 = t; float r24801 = l; float r24802 = r24800 / r24801; float r24803 = pow(r24802, r24797); float r24804 = r24797 * r24803; float r24805 = r24793 + r24804; float r24806 = r24799 / r24805; float r24807 = sqrt(r24806); float r24808 = asin(r24807); return r24808; } double f_id(double t, double l, double Om, double Omc) { double r24809 = 1; double r24810 = Om; double r24811 = Omc; double r24812 = r24810 / r24811; double r24813 = 2; double r24814 = pow(r24812, r24813); double r24815 = r24809 - r24814; double r24816 = t; double r24817 = l; double r24818 = r24816 / r24817; double r24819 = pow(r24818, r24813); double r24820 = r24813 * r24819; double r24821 = r24809 + r24820; double r24822 = r24815 / r24821; double r24823 = sqrt(r24822); double r24824 = asin(r24823); return r24824; } double f_of(float t, float l, float Om, float Omc) { float r24825 = 1; float r24826 = Om; float r24827 = Omc; float r24828 = r24826 / r24827; float r24829 = 2; float r24830 = pow(r24828, r24829); float r24831 = r24825 - r24830; float r24832 = t; float r24833 = l; float r24834 = r24832 / r24833; float r24835 = pow(r24834, r24829); float r24836 = r24829 * r24835; float r24837 = r24825 + r24836; float r24838 = r24831 / r24837; float r24839 = sqrt(r24838); float r24840 = asin(r24839); return r24840; } double f_od(double t, double l, double Om, double Omc) { double r24841 = 1; double r24842 = Om; double r24843 = Omc; double r24844 = r24842 / r24843; double r24845 = 2; double r24846 = pow(r24844, r24845); double r24847 = r24841 - r24846; double r24848 = t; double r24849 = l; double r24850 = r24848 / r24849; double r24851 = pow(r24850, r24845); double r24852 = r24845 * r24851; double r24853 = r24841 + r24852; double r24854 = r24847 / r24853; double r24855 = sqrt(r24854); double r24856 = asin(r24855); return r24856; } 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 r24857, r24858, r24859, r24860, r24861, r24862, r24863, r24864, r24865, r24866, r24867, r24868, r24869, r24870, r24871, r24872; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24857, "1", 10, MPFR_RNDN); mpfr_init(r24858); mpfr_init(r24859); mpfr_init(r24860); mpfr_init_set_str(r24861, "2", 10, MPFR_RNDN); mpfr_init(r24862); mpfr_init(r24863); mpfr_init(r24864); mpfr_init(r24865); mpfr_init(r24866); mpfr_init(r24867); mpfr_init(r24868); mpfr_init(r24869); mpfr_init(r24870); mpfr_init(r24871); mpfr_init(r24872); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24858, Om, MPFR_RNDN); mpfr_set_d(r24859, Omc, MPFR_RNDN); mpfr_div(r24860, r24858, r24859, MPFR_RNDN); ; mpfr_pow(r24862, r24860, r24861, MPFR_RNDN); mpfr_sub(r24863, r24857, r24862, MPFR_RNDN); mpfr_set_d(r24864, t, MPFR_RNDN); mpfr_set_d(r24865, l, MPFR_RNDN); mpfr_div(r24866, r24864, r24865, MPFR_RNDN); mpfr_pow(r24867, r24866, r24861, MPFR_RNDN); mpfr_mul(r24868, r24861, r24867, MPFR_RNDN); mpfr_add(r24869, r24857, r24868, MPFR_RNDN); mpfr_div(r24870, r24863, r24869, MPFR_RNDN); mpfr_sqrt(r24871, r24870, MPFR_RNDN); mpfr_asin(r24872, r24871, MPFR_RNDN); return mpfr_get_d(r24872, MPFR_RNDN); } static mpfr_t r24873, r24874, r24875, r24876, r24877, r24878, r24879, r24880, r24881, r24882, r24883, r24884, r24885, r24886, r24887, r24888; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24873, "1", 10, MPFR_RNDN); mpfr_init(r24874); mpfr_init(r24875); mpfr_init(r24876); mpfr_init_set_str(r24877, "2", 10, MPFR_RNDN); mpfr_init(r24878); mpfr_init(r24879); mpfr_init(r24880); mpfr_init(r24881); mpfr_init(r24882); mpfr_init(r24883); mpfr_init(r24884); mpfr_init(r24885); mpfr_init(r24886); mpfr_init(r24887); mpfr_init(r24888); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24874, Om, MPFR_RNDN); mpfr_set_d(r24875, Omc, MPFR_RNDN); mpfr_div(r24876, r24874, r24875, MPFR_RNDN); ; mpfr_pow(r24878, r24876, r24877, MPFR_RNDN); mpfr_sub(r24879, r24873, r24878, MPFR_RNDN); mpfr_set_d(r24880, t, MPFR_RNDN); mpfr_set_d(r24881, l, MPFR_RNDN); mpfr_div(r24882, r24880, r24881, MPFR_RNDN); mpfr_pow(r24883, r24882, r24877, MPFR_RNDN); mpfr_mul(r24884, r24877, r24883, MPFR_RNDN); mpfr_add(r24885, r24873, r24884, MPFR_RNDN); mpfr_div(r24886, r24879, r24885, MPFR_RNDN); mpfr_sqrt(r24887, r24886, MPFR_RNDN); mpfr_asin(r24888, r24887, MPFR_RNDN); return mpfr_get_d(r24888, MPFR_RNDN); } static mpfr_t r24889, r24890, r24891, r24892, r24893, r24894, r24895, r24896, r24897, r24898, r24899, r24900, r24901, r24902, r24903, r24904; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24889, "1", 10, MPFR_RNDN); mpfr_init(r24890); mpfr_init(r24891); mpfr_init(r24892); mpfr_init_set_str(r24893, "2", 10, MPFR_RNDN); mpfr_init(r24894); mpfr_init(r24895); mpfr_init(r24896); mpfr_init(r24897); mpfr_init(r24898); mpfr_init(r24899); mpfr_init(r24900); mpfr_init(r24901); mpfr_init(r24902); mpfr_init(r24903); mpfr_init(r24904); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24890, Om, MPFR_RNDN); mpfr_set_d(r24891, Omc, MPFR_RNDN); mpfr_div(r24892, r24890, r24891, MPFR_RNDN); ; mpfr_pow(r24894, r24892, r24893, MPFR_RNDN); mpfr_sub(r24895, r24889, r24894, MPFR_RNDN); mpfr_set_d(r24896, t, MPFR_RNDN); mpfr_set_d(r24897, l, MPFR_RNDN); mpfr_div(r24898, r24896, r24897, MPFR_RNDN); mpfr_pow(r24899, r24898, r24893, MPFR_RNDN); mpfr_mul(r24900, r24893, r24899, MPFR_RNDN); mpfr_add(r24901, r24889, r24900, MPFR_RNDN); mpfr_div(r24902, r24895, r24901, MPFR_RNDN); mpfr_sqrt(r24903, r24902, MPFR_RNDN); mpfr_asin(r24904, r24903, MPFR_RNDN); return mpfr_get_d(r24904, MPFR_RNDN); }