\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)double f(double t, double l, double Om, double Omc) {
double r66813 = 1.0;
double r66814 = Om;
double r66815 = Omc;
double r66816 = r66814 / r66815;
double r66817 = 2.0;
double r66818 = pow(r66816, r66817);
double r66819 = r66813 - r66818;
double r66820 = t;
double r66821 = l;
double r66822 = r66820 / r66821;
double r66823 = pow(r66822, r66817);
double r66824 = r66817 * r66823;
double r66825 = r66813 + r66824;
double r66826 = r66819 / r66825;
double r66827 = sqrt(r66826);
double r66828 = asin(r66827);
return r66828;
}
double f(double t, double l, double Om, double Omc) {
double r66829 = 1.0;
double r66830 = Om;
double r66831 = Omc;
double r66832 = r66830 / r66831;
double r66833 = 2.0;
double r66834 = pow(r66832, r66833);
double r66835 = r66829 - r66834;
double r66836 = t;
double r66837 = l;
double r66838 = r66836 / r66837;
double r66839 = pow(r66838, r66833);
double r66840 = r66833 * r66839;
double r66841 = r66829 + r66840;
double r66842 = r66835 / r66841;
double r66843 = sqrt(r66842);
double r66844 = asin(r66843);
return r66844;
}



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus Omc
Results
Initial program 10.2
Final simplification10.2
herbie shell --seed 2020047
(FPCore (t l Om Omc)
:name "Toniolo and Linder, Equation (2)"
:precision binary64
(asin (sqrt (/ (- 1 (pow (/ Om Omc) 2)) (+ 1 (* 2 (pow (/ t l) 2)))))))