Average Error: 10.3 → 10.3
Time: 46.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)\]
↓
\[\mathsf{expm1}\left(\mathsf{log1p}\left(\sin^{-1} \left(\sqrt{\frac{{1}^{3} - {\left({\left(\frac{Om}{Omc}\right)}^{2}\right)}^{3}}{\left({\left(\frac{Om}{Omc}\right)}^{2} \cdot \left({\left(\frac{Om}{Omc}\right)}^{2} + 1\right) + 1 \cdot 1\right) \cdot \left(1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}\right)}}\right)\right)\right)\]
\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)
↓
\mathsf{expm1}\left(\mathsf{log1p}\left(\sin^{-1} \left(\sqrt{\frac{{1}^{3} - {\left({\left(\frac{Om}{Omc}\right)}^{2}\right)}^{3}}{\left({\left(\frac{Om}{Omc}\right)}^{2} \cdot \left({\left(\frac{Om}{Omc}\right)}^{2} + 1\right) + 1 \cdot 1\right) \cdot \left(1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}\right)}}\right)\right)\right)
double code(double t, double l, double Om, double Omc) {
	return ((double) asin(((double) sqrt(((double) (((double) (1.0 - ((double) pow(((double) (Om / Omc)), 2.0)))) / ((double) (1.0 + ((double) (2.0 * ((double) pow(((double) (t / l)), 2.0))))))))))));
}
↓
double code(double t, double l, double Om, double Omc) {
	return ((double) expm1(((double) log1p(((double) asin(((double) sqrt(((double) (((double) (((double) pow(1.0, 3.0)) - ((double) pow(((double) pow(((double) (Om / Omc)), 2.0)), 3.0)))) / ((double) (((double) (((double) (((double) pow(((double) (Om / Omc)), 2.0)) * ((double) (((double) pow(((double) (Om / Omc)), 2.0)) + 1.0)))) + ((double) (1.0 * 1.0)))) * ((double) (1.0 + ((double) (2.0 * ((double) pow(((double) (t / l)), 2.0))))))))))))))))));
}

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.3

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

  3. Applied flip3--10.3

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

  4. Applied associate-/l/10.3

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

  5. Simplified10.3

    \[\leadsto \sin^{-1} \left(\sqrt{\frac{{1}^{3} - {\left({\left(\frac{Om}{Omc}\right)}^{2}\right)}^{3}}{\color{blue}{\left({\left(\frac{Om}{Omc}\right)}^{2} \cdot \left({\left(\frac{Om}{Omc}\right)}^{2} + 1\right) + 1 \cdot 1\right) \cdot \left(1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}\right)}}}\right)\]
  6. Using strategy rm

  7. Applied expm1-log1p-u10.3

    \[\leadsto \color{blue}{\mathsf{expm1}\left(\mathsf{log1p}\left(\sin^{-1} \left(\sqrt{\frac{{1}^{3} - {\left({\left(\frac{Om}{Omc}\right)}^{2}\right)}^{3}}{\left({\left(\frac{Om}{Omc}\right)}^{2} \cdot \left({\left(\frac{Om}{Omc}\right)}^{2} + 1\right) + 1 \cdot 1\right) \cdot \left(1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}\right)}}\right)\right)\right)}\]

  8. Final simplification10.3

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

Reproduce

herbie shell --seed 2020113 +o rules:numerics
(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)))))))