Average Error: 10.4 → 10.4
Time: 10.7s
Precision: 64
\[\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)\]
\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 r79847 = 1.0;
        double r79848 = Om;
        double r79849 = Omc;
        double r79850 = r79848 / r79849;
        double r79851 = 2.0;
        double r79852 = pow(r79850, r79851);
        double r79853 = r79847 - r79852;
        double r79854 = t;
        double r79855 = l;
        double r79856 = r79854 / r79855;
        double r79857 = pow(r79856, r79851);
        double r79858 = r79851 * r79857;
        double r79859 = r79847 + r79858;
        double r79860 = r79853 / r79859;
        double r79861 = sqrt(r79860);
        double r79862 = asin(r79861);
        return r79862;
}
↓
double f(double t, double l, double Om, double Omc) {
        double r79863 = 1.0;
        double r79864 = Om;
        double r79865 = Omc;
        double r79866 = r79864 / r79865;
        double r79867 = 2.0;
        double r79868 = pow(r79866, r79867);
        double r79869 = r79863 - r79868;
        double r79870 = t;
        double r79871 = l;
        double r79872 = r79870 / r79871;
        double r79873 = pow(r79872, r79867);
        double r79874 = r79867 * r79873;
        double r79875 = r79863 + r79874;
        double r79876 = r79869 / r79875;
        double r79877 = sqrt(r79876);
        double r79878 = asin(r79877);
        return r79878;
}

Error

Bits error versus t

Bits error versus l

Bits error versus Om

Bits error versus Omc

Try it out

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation

  1. Initial program 10.4

    \[\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\]

  2. Final simplification10.4

    \[\leadsto \sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\]

Reproduce

herbie shell --seed 2019346 
(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)))))))