#include #include #include #include #include char *name = "Toniolo and Linder, Equation (2)"; double f_if(float t, float l, float Om, float Omc) { float r24902 = 1; float r24903 = Om; float r24904 = Omc; float r24905 = r24903 / r24904; float r24906 = 2; float r24907 = pow(r24905, r24906); float r24908 = r24902 - r24907; float r24909 = t; float r24910 = l; float r24911 = r24909 / r24910; float r24912 = pow(r24911, r24906); float r24913 = r24906 * r24912; float r24914 = r24902 + r24913; float r24915 = r24908 / r24914; float r24916 = sqrt(r24915); float r24917 = asin(r24916); return r24917; } double f_id(double t, double l, double Om, double Omc) { double r24918 = 1; double r24919 = Om; double r24920 = Omc; double r24921 = r24919 / r24920; double r24922 = 2; double r24923 = pow(r24921, r24922); double r24924 = r24918 - r24923; double r24925 = t; double r24926 = l; double r24927 = r24925 / r24926; double r24928 = pow(r24927, r24922); double r24929 = r24922 * r24928; double r24930 = r24918 + r24929; double r24931 = r24924 / r24930; double r24932 = sqrt(r24931); double r24933 = asin(r24932); return r24933; } double f_of(float t, float l, float Om, float Omc) { float r24934 = 1; float r24935 = Om; float r24936 = Omc; float r24937 = r24935 / r24936; float r24938 = r24937 * r24937; float r24939 = r24934 - r24938; float r24940 = 2; float r24941 = l; float r24942 = t; float r24943 = r24941 / r24942; float r24944 = r24940 / r24943; float r24945 = r24944 / r24943; float r24946 = r24945 + r24934; float r24947 = r24939 / r24946; float r24948 = sqrt(r24947); float r24949 = asin(r24948); return r24949; } double f_od(double t, double l, double Om, double Omc) { double r24950 = 1; double r24951 = Om; double r24952 = Omc; double r24953 = r24951 / r24952; double r24954 = r24953 * r24953; double r24955 = r24950 - r24954; double r24956 = 2; double r24957 = l; double r24958 = t; double r24959 = r24957 / r24958; double r24960 = r24956 / r24959; double r24961 = r24960 / r24959; double r24962 = r24961 + r24950; double r24963 = r24955 / r24962; double r24964 = sqrt(r24963); double r24965 = asin(r24964); return r24965; } 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 r24966, r24967, r24968, r24969, r24970, r24971, r24972, r24973, r24974, r24975, r24976, r24977, r24978, r24979, r24980, r24981; void setup_mpfr_f_im() { mpfr_set_default_prec(592); mpfr_init_set_str(r24966, "1", 10, MPFR_RNDN); mpfr_init(r24967); mpfr_init(r24968); mpfr_init(r24969); mpfr_init_set_str(r24970, "2", 10, MPFR_RNDN); mpfr_init(r24971); mpfr_init(r24972); mpfr_init(r24973); mpfr_init(r24974); mpfr_init(r24975); mpfr_init(r24976); mpfr_init(r24977); mpfr_init(r24978); mpfr_init(r24979); mpfr_init(r24980); mpfr_init(r24981); } double f_im(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24967, Om, MPFR_RNDN); mpfr_set_d(r24968, Omc, MPFR_RNDN); mpfr_div(r24969, r24967, r24968, MPFR_RNDN); ; mpfr_pow(r24971, r24969, r24970, MPFR_RNDN); mpfr_sub(r24972, r24966, r24971, MPFR_RNDN); mpfr_set_d(r24973, t, MPFR_RNDN); mpfr_set_d(r24974, l, MPFR_RNDN); mpfr_div(r24975, r24973, r24974, MPFR_RNDN); mpfr_pow(r24976, r24975, r24970, MPFR_RNDN); mpfr_mul(r24977, r24970, r24976, MPFR_RNDN); mpfr_add(r24978, r24966, r24977, MPFR_RNDN); mpfr_div(r24979, r24972, r24978, MPFR_RNDN); mpfr_sqrt(r24980, r24979, MPFR_RNDN); mpfr_asin(r24981, r24980, MPFR_RNDN); return mpfr_get_d(r24981, MPFR_RNDN); } static mpfr_t r24982, r24983, r24984, r24985, r24986, r24987, r24988, r24989, r24990, r24991, r24992, r24993, r24994, r24995, r24996, r24997; void setup_mpfr_f_fm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24982, "1", 10, MPFR_RNDN); mpfr_init(r24983); mpfr_init(r24984); mpfr_init(r24985); mpfr_init(r24986); mpfr_init(r24987); mpfr_init_set_str(r24988, "2", 10, MPFR_RNDN); mpfr_init(r24989); mpfr_init(r24990); mpfr_init(r24991); mpfr_init(r24992); mpfr_init(r24993); mpfr_init(r24994); mpfr_init(r24995); mpfr_init(r24996); mpfr_init(r24997); } double f_fm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24983, Om, MPFR_RNDN); mpfr_set_d(r24984, Omc, MPFR_RNDN); mpfr_div(r24985, r24983, r24984, MPFR_RNDN); mpfr_mul(r24986, r24985, r24985, MPFR_RNDN); mpfr_sub(r24987, r24982, r24986, MPFR_RNDN); ; mpfr_set_d(r24989, l, MPFR_RNDN); mpfr_set_d(r24990, t, MPFR_RNDN); mpfr_div(r24991, r24989, r24990, MPFR_RNDN); mpfr_div(r24992, r24988, r24991, MPFR_RNDN); mpfr_div(r24993, r24992, r24991, MPFR_RNDN); mpfr_add(r24994, r24993, r24982, MPFR_RNDN); mpfr_div(r24995, r24987, r24994, MPFR_RNDN); mpfr_sqrt(r24996, r24995, MPFR_RNDN); mpfr_asin(r24997, r24996, MPFR_RNDN); return mpfr_get_d(r24997, MPFR_RNDN); } static mpfr_t r24998, r24999, r25000, r25001, r25002, r25003, r25004, r25005, r25006, r25007, r25008, r25009, r25010, r25011, r25012, r25013; void setup_mpfr_f_dm() { mpfr_set_default_prec(592); mpfr_init_set_str(r24998, "1", 10, MPFR_RNDN); mpfr_init(r24999); mpfr_init(r25000); mpfr_init(r25001); mpfr_init(r25002); mpfr_init(r25003); mpfr_init_set_str(r25004, "2", 10, MPFR_RNDN); mpfr_init(r25005); mpfr_init(r25006); mpfr_init(r25007); mpfr_init(r25008); mpfr_init(r25009); mpfr_init(r25010); mpfr_init(r25011); mpfr_init(r25012); mpfr_init(r25013); } double f_dm(double t, double l, double Om, double Omc) { ; mpfr_set_d(r24999, Om, MPFR_RNDN); mpfr_set_d(r25000, Omc, MPFR_RNDN); mpfr_div(r25001, r24999, r25000, MPFR_RNDN); mpfr_mul(r25002, r25001, r25001, MPFR_RNDN); mpfr_sub(r25003, r24998, r25002, MPFR_RNDN); ; mpfr_set_d(r25005, l, MPFR_RNDN); mpfr_set_d(r25006, t, MPFR_RNDN); mpfr_div(r25007, r25005, r25006, MPFR_RNDN); mpfr_div(r25008, r25004, r25007, MPFR_RNDN); mpfr_div(r25009, r25008, r25007, MPFR_RNDN); mpfr_add(r25010, r25009, r24998, MPFR_RNDN); mpfr_div(r25011, r25003, r25010, MPFR_RNDN); mpfr_sqrt(r25012, r25011, MPFR_RNDN); mpfr_asin(r25013, r25012, MPFR_RNDN); return mpfr_get_d(r25013, MPFR_RNDN); }