#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r21957 = 1; float r21958 = Om; float r21959 = Omc; float r21960 = r21958 / r21959; float r21961 = 2; float r21962 = pow(r21960, r21961); float r21963 = r21957 - r21962; float r21964 = t; float r21965 = l; float r21966 = r21964 / r21965; float r21967 = pow(r21966, r21961); float r21968 = r21961 * r21967; float r21969 = r21957 + r21968; float r21970 = r21963 / r21969; float r21971 = sqrt(r21970); float r21972 = asin(r21971); return r21972; } double f_id(double t, double l, double Om, double Omc) { double r21973 = 1; double r21974 = Om; double r21975 = Omc; double r21976 = r21974 / r21975; double r21977 = 2; double r21978 = pow(r21976, r21977); double r21979 = r21973 - r21978; double r21980 = t; double r21981 = l; double r21982 = r21980 / r21981; double r21983 = pow(r21982, r21977); double r21984 = r21977 * r21983; double r21985 = r21973 + r21984; double r21986 = r21979 / r21985; double r21987 = sqrt(r21986); double r21988 = asin(r21987); return r21988; } double f_of(float t, float l, float Om, float Omc) { float r21989 = 1; float r21990 = Om; float r21991 = Omc; float r21992 = r21990 / r21991; float r21993 = 2; float r21994 = pow(r21992, r21993); float r21995 = r21989 - r21994; float r21996 = t; float r21997 = l; float r21998 = r21996 / r21997; float r21999 = pow(r21998, r21993); float r22000 = r21993 * r21999; float r22001 = r21989 + r22000; float r22002 = r21995 / r22001; float r22003 = sqrt(r22002); float r22004 = asin(r22003); return r22004; } double f_od(double t, double l, double Om, double Omc) { double r22005 = 1; double r22006 = Om; double r22007 = Omc; double r22008 = r22006 / r22007; double r22009 = 2; double r22010 = pow(r22008, r22009); double r22011 = r22005 - r22010; double r22012 = t; double r22013 = l; double r22014 = r22012 / r22013; double r22015 = pow(r22014, r22009); double r22016 = r22009 * r22015; double r22017 = r22005 + r22016; double r22018 = r22011 / r22017; double r22019 = sqrt(r22018); double r22020 = asin(r22019); return r22020; } 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 r22021, r22022, r22023, r22024, r22025, r22026, r22027, r22028, r22029, r22030, r22031, r22032, r22033, r22034, r22035, r22036; void setup_mpfr_f_im() { mpfr_set_default_prec(400); mpfr_init_set_str(r22021, "1", 10, MPFR_RNDN); mpfr_init(r22022); mpfr_init(r22023); mpfr_init(r22024); mpfr_init_set_str(r22025, "2", 10, MPFR_RNDN); mpfr_init(r22026); mpfr_init(r22027); mpfr_init(r22028); mpfr_init(r22029); mpfr_init(r22030); mpfr_init(r22031); mpfr_init(r22032); mpfr_init(r22033); mpfr_init(r22034); mpfr_init(r22035); mpfr_init(r22036); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r22022, Om, MPFR_RNDN); mpfr_set_d(r22023, Omc, MPFR_RNDN); mpfr_div(r22024, r22022, r22023, MPFR_RNDN); ; mpfr_pow(r22026, r22024, r22025, MPFR_RNDN); mpfr_sub(r22027, r22021, r22026, MPFR_RNDN); mpfr_set_d(r22028, t, MPFR_RNDN); mpfr_set_d(r22029, l, MPFR_RNDN); mpfr_div(r22030, r22028, r22029, MPFR_RNDN); mpfr_pow(r22031, r22030, r22025, MPFR_RNDN); mpfr_mul(r22032, r22025, r22031, MPFR_RNDN); mpfr_add(r22033, r22021, r22032, MPFR_RNDN); mpfr_div(r22034, r22027, r22033, MPFR_RNDN); mpfr_sqrt(r22035, r22034, MPFR_RNDN); mpfr_asin(r22036, r22035, MPFR_RNDN); return mpfr_get_d(r22036, MPFR_RNDN); } static mpfr_t r22037, r22038, r22039, r22040, r22041, r22042, r22043, r22044, r22045, r22046, r22047, r22048, r22049, r22050, r22051, r22052; void setup_mpfr_f_fm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22037, "1", 10, MPFR_RNDN); mpfr_init(r22038); mpfr_init(r22039); mpfr_init(r22040); mpfr_init_set_str(r22041, "2", 10, MPFR_RNDN); mpfr_init(r22042); mpfr_init(r22043); mpfr_init(r22044); mpfr_init(r22045); mpfr_init(r22046); mpfr_init(r22047); mpfr_init(r22048); mpfr_init(r22049); mpfr_init(r22050); mpfr_init(r22051); mpfr_init(r22052); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r22038, Om, MPFR_RNDN); mpfr_set_d(r22039, Omc, MPFR_RNDN); mpfr_div(r22040, r22038, r22039, MPFR_RNDN); ; mpfr_pow(r22042, r22040, r22041, MPFR_RNDN); mpfr_sub(r22043, r22037, r22042, MPFR_RNDN); mpfr_set_d(r22044, t, MPFR_RNDN); mpfr_set_d(r22045, l, MPFR_RNDN); mpfr_div(r22046, r22044, r22045, MPFR_RNDN); mpfr_pow(r22047, r22046, r22041, MPFR_RNDN); mpfr_mul(r22048, r22041, r22047, MPFR_RNDN); mpfr_add(r22049, r22037, r22048, MPFR_RNDN); mpfr_div(r22050, r22043, r22049, MPFR_RNDN); mpfr_sqrt(r22051, r22050, MPFR_RNDN); mpfr_asin(r22052, r22051, MPFR_RNDN); return mpfr_get_d(r22052, MPFR_RNDN); } static mpfr_t r22053, r22054, r22055, r22056, r22057, r22058, r22059, r22060, r22061, r22062, r22063, r22064, r22065, r22066, r22067, r22068; void setup_mpfr_f_dm() { mpfr_set_default_prec(400); mpfr_init_set_str(r22053, "1", 10, MPFR_RNDN); mpfr_init(r22054); mpfr_init(r22055); mpfr_init(r22056); mpfr_init_set_str(r22057, "2", 10, MPFR_RNDN); mpfr_init(r22058); mpfr_init(r22059); mpfr_init(r22060); mpfr_init(r22061); mpfr_init(r22062); mpfr_init(r22063); mpfr_init(r22064); mpfr_init(r22065); mpfr_init(r22066); mpfr_init(r22067); mpfr_init(r22068); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r22054, Om, MPFR_RNDN); mpfr_set_d(r22055, Omc, MPFR_RNDN); mpfr_div(r22056, r22054, r22055, MPFR_RNDN); ; mpfr_pow(r22058, r22056, r22057, MPFR_RNDN); mpfr_sub(r22059, r22053, r22058, MPFR_RNDN); mpfr_set_d(r22060, t, MPFR_RNDN); mpfr_set_d(r22061, l, MPFR_RNDN); mpfr_div(r22062, r22060, r22061, MPFR_RNDN); mpfr_pow(r22063, r22062, r22057, MPFR_RNDN); mpfr_mul(r22064, r22057, r22063, MPFR_RNDN); mpfr_add(r22065, r22053, r22064, MPFR_RNDN); mpfr_div(r22066, r22059, r22065, MPFR_RNDN); mpfr_sqrt(r22067, r22066, MPFR_RNDN); mpfr_asin(r22068, r22067, MPFR_RNDN); return mpfr_get_d(r22068, MPFR_RNDN); }