#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r8809 = 1; float r8810 = Om; float r8811 = Omc; float r8812 = r8810 / r8811; float r8813 = 2; float r8814 = pow(r8812, r8813); float r8815 = r8809 - r8814; float r8816 = t; float r8817 = l; float r8818 = r8816 / r8817; float r8819 = pow(r8818, r8813); float r8820 = r8813 * r8819; float r8821 = r8809 + r8820; float r8822 = r8815 / r8821; float r8823 = sqrt(r8822); float r8824 = asin(r8823); return r8824; } double f_id(double t, double l, double Om, double Omc) { double r8825 = 1; double r8826 = Om; double r8827 = Omc; double r8828 = r8826 / r8827; double r8829 = 2; double r8830 = pow(r8828, r8829); double r8831 = r8825 - r8830; double r8832 = t; double r8833 = l; double r8834 = r8832 / r8833; double r8835 = pow(r8834, r8829); double r8836 = r8829 * r8835; double r8837 = r8825 + r8836; double r8838 = r8831 / r8837; double r8839 = sqrt(r8838); double r8840 = asin(r8839); return r8840; } double f_of(float t, float l, float Om, float Omc) { float r8841 = 1; float r8842 = Om; float r8843 = Omc; float r8844 = r8842 / r8843; float r8845 = 2; float r8846 = pow(r8844, r8845); float r8847 = r8841 - r8846; float r8848 = t; float r8849 = l; float r8850 = r8848 / r8849; float r8851 = pow(r8850, r8845); float r8852 = r8845 * r8851; float r8853 = r8841 + r8852; float r8854 = r8847 / r8853; float r8855 = sqrt(r8854); float r8856 = asin(r8855); float r8857 = log1p(r8856); float r8858 = expm1(r8857); return r8858; } double f_od(double t, double l, double Om, double Omc) { double r8859 = 1; double r8860 = Om; double r8861 = Omc; double r8862 = r8860 / r8861; double r8863 = 2; double r8864 = pow(r8862, r8863); double r8865 = r8859 - r8864; double r8866 = t; double r8867 = l; double r8868 = r8866 / r8867; double r8869 = pow(r8868, r8863); double r8870 = r8863 * r8869; double r8871 = r8859 + r8870; double r8872 = r8865 / r8871; double r8873 = sqrt(r8872); double r8874 = asin(r8873); double r8875 = log1p(r8874); double r8876 = expm1(r8875); return r8876; } 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 r8877, r8878, r8879, r8880, r8881, r8882, r8883, r8884, r8885, r8886, r8887, r8888, r8889, r8890, r8891, r8892; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r8877, "1", 10, MPFR_RNDN); mpfr_init(r8878); mpfr_init(r8879); mpfr_init(r8880); mpfr_init_set_str(r8881, "2", 10, MPFR_RNDN); mpfr_init(r8882); mpfr_init(r8883); mpfr_init(r8884); mpfr_init(r8885); mpfr_init(r8886); mpfr_init(r8887); mpfr_init(r8888); mpfr_init(r8889); mpfr_init(r8890); mpfr_init(r8891); mpfr_init(r8892); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8878, Om, MPFR_RNDN); mpfr_set_d(r8879, Omc, MPFR_RNDN); mpfr_div(r8880, r8878, r8879, MPFR_RNDN); ; mpfr_pow(r8882, r8880, r8881, MPFR_RNDN); mpfr_sub(r8883, r8877, r8882, MPFR_RNDN); mpfr_set_d(r8884, t, MPFR_RNDN); mpfr_set_d(r8885, l, MPFR_RNDN); mpfr_div(r8886, r8884, r8885, MPFR_RNDN); mpfr_pow(r8887, r8886, r8881, MPFR_RNDN); mpfr_mul(r8888, r8881, r8887, MPFR_RNDN); mpfr_add(r8889, r8877, r8888, MPFR_RNDN); mpfr_div(r8890, r8883, r8889, MPFR_RNDN); mpfr_sqrt(r8891, r8890, MPFR_RNDN); mpfr_asin(r8892, r8891, MPFR_RNDN); return mpfr_get_d(r8892, MPFR_RNDN); } static mpfr_t r8893, r8894, r8895, r8896, r8897, r8898, r8899, r8900, r8901, r8902, r8903, r8904, r8905, r8906, r8907, r8908, r8909, r8910; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8893, "1", 10, MPFR_RNDN); mpfr_init(r8894); mpfr_init(r8895); mpfr_init(r8896); mpfr_init_set_str(r8897, "2", 10, MPFR_RNDN); mpfr_init(r8898); mpfr_init(r8899); mpfr_init(r8900); mpfr_init(r8901); mpfr_init(r8902); mpfr_init(r8903); mpfr_init(r8904); mpfr_init(r8905); mpfr_init(r8906); mpfr_init(r8907); mpfr_init(r8908); mpfr_init(r8909); mpfr_init(r8910); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8894, Om, MPFR_RNDN); mpfr_set_d(r8895, Omc, MPFR_RNDN); mpfr_div(r8896, r8894, r8895, MPFR_RNDN); ; mpfr_pow(r8898, r8896, r8897, MPFR_RNDN); mpfr_sub(r8899, r8893, r8898, MPFR_RNDN); mpfr_set_d(r8900, t, MPFR_RNDN); mpfr_set_d(r8901, l, MPFR_RNDN); mpfr_div(r8902, r8900, r8901, MPFR_RNDN); mpfr_pow(r8903, r8902, r8897, MPFR_RNDN); mpfr_mul(r8904, r8897, r8903, MPFR_RNDN); mpfr_add(r8905, r8893, r8904, MPFR_RNDN); mpfr_div(r8906, r8899, r8905, MPFR_RNDN); mpfr_sqrt(r8907, r8906, MPFR_RNDN); mpfr_asin(r8908, r8907, MPFR_RNDN); mpfr_log1p(r8909, r8908, MPFR_RNDN); mpfr_expm1(r8910, r8909, MPFR_RNDN); return mpfr_get_d(r8910, MPFR_RNDN); } static mpfr_t r8911, r8912, r8913, r8914, r8915, r8916, r8917, r8918, r8919, r8920, r8921, r8922, r8923, r8924, r8925, r8926, r8927, r8928; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r8911, "1", 10, MPFR_RNDN); mpfr_init(r8912); mpfr_init(r8913); mpfr_init(r8914); mpfr_init_set_str(r8915, "2", 10, MPFR_RNDN); mpfr_init(r8916); mpfr_init(r8917); mpfr_init(r8918); mpfr_init(r8919); mpfr_init(r8920); mpfr_init(r8921); mpfr_init(r8922); mpfr_init(r8923); mpfr_init(r8924); mpfr_init(r8925); mpfr_init(r8926); mpfr_init(r8927); mpfr_init(r8928); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r8912, Om, MPFR_RNDN); mpfr_set_d(r8913, Omc, MPFR_RNDN); mpfr_div(r8914, r8912, r8913, MPFR_RNDN); ; mpfr_pow(r8916, r8914, r8915, MPFR_RNDN); mpfr_sub(r8917, r8911, r8916, MPFR_RNDN); mpfr_set_d(r8918, t, MPFR_RNDN); mpfr_set_d(r8919, l, MPFR_RNDN); mpfr_div(r8920, r8918, r8919, MPFR_RNDN); mpfr_pow(r8921, r8920, r8915, MPFR_RNDN); mpfr_mul(r8922, r8915, r8921, MPFR_RNDN); mpfr_add(r8923, r8911, r8922, MPFR_RNDN); mpfr_div(r8924, r8917, r8923, MPFR_RNDN); mpfr_sqrt(r8925, r8924, MPFR_RNDN); mpfr_asin(r8926, r8925, MPFR_RNDN); mpfr_log1p(r8927, r8926, MPFR_RNDN); mpfr_expm1(r8928, r8927, MPFR_RNDN); return mpfr_get_d(r8928, MPFR_RNDN); }