Toniolo and Linder, Equation (2)

?

Percentage Accurate: 84.1% → 98.8%
Time: 46.6s
Precision: binary64
Cost: 20872

?

\[\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right) \]
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+131}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{1 - \frac{Om}{Omc} \cdot \frac{Om}{Omc}} \cdot \frac{-\ell}{t \cdot \sqrt{2}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+97}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{\frac{1 - \frac{\frac{Om}{Omc}}{\frac{Omc}{Om}}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)\\ \end{array} \]
(FPCore (t l Om Omc)
 :precision binary64
 (asin
  (sqrt (/ (- 1.0 (pow (/ Om Omc) 2.0)) (+ 1.0 (* 2.0 (pow (/ t l) 2.0)))))))
(FPCore (t l Om Omc)
 :precision binary64
 (if (<= (/ t l) -5e+131)
   (asin
    (* (sqrt (- 1.0 (* (/ Om Omc) (/ Om Omc)))) (/ (- l) (* t (sqrt 2.0)))))
   (if (<= (/ t l) 5e+97)
     (asin
      (sqrt
       (/
        (- 1.0 (/ (/ Om Omc) (/ Omc Om)))
        (+ 1.0 (* 2.0 (pow (/ t l) 2.0))))))
     (asin (* (/ l t) (sqrt 0.5))))))
double code(double t, double l, double Om, double Omc) {
	return asin(sqrt(((1.0 - pow((Om / Omc), 2.0)) / (1.0 + (2.0 * pow((t / l), 2.0))))));
}
double code(double t, double l, double Om, double Omc) {
	double tmp;
	if ((t / l) <= -5e+131) {
		tmp = asin((sqrt((1.0 - ((Om / Omc) * (Om / Omc)))) * (-l / (t * sqrt(2.0)))));
	} else if ((t / l) <= 5e+97) {
		tmp = asin(sqrt(((1.0 - ((Om / Omc) / (Omc / Om))) / (1.0 + (2.0 * pow((t / l), 2.0))))));
	} else {
		tmp = asin(((l / t) * sqrt(0.5)));
	}
	return tmp;
}
real(8) function code(t, l, om, omc)
    real(8), intent (in) :: t
    real(8), intent (in) :: l
    real(8), intent (in) :: om
    real(8), intent (in) :: omc
    code = asin(sqrt(((1.0d0 - ((om / omc) ** 2.0d0)) / (1.0d0 + (2.0d0 * ((t / l) ** 2.0d0))))))
end function
real(8) function code(t, l, om, omc)
    real(8), intent (in) :: t
    real(8), intent (in) :: l
    real(8), intent (in) :: om
    real(8), intent (in) :: omc
    real(8) :: tmp
    if ((t / l) <= (-5d+131)) then
        tmp = asin((sqrt((1.0d0 - ((om / omc) * (om / omc)))) * (-l / (t * sqrt(2.0d0)))))
    else if ((t / l) <= 5d+97) then
        tmp = asin(sqrt(((1.0d0 - ((om / omc) / (omc / om))) / (1.0d0 + (2.0d0 * ((t / l) ** 2.0d0))))))
    else
        tmp = asin(((l / t) * sqrt(0.5d0)))
    end if
    code = tmp
end function
public static double code(double t, double l, double Om, double Omc) {
	return Math.asin(Math.sqrt(((1.0 - Math.pow((Om / Omc), 2.0)) / (1.0 + (2.0 * Math.pow((t / l), 2.0))))));
}
public static double code(double t, double l, double Om, double Omc) {
	double tmp;
	if ((t / l) <= -5e+131) {
		tmp = Math.asin((Math.sqrt((1.0 - ((Om / Omc) * (Om / Omc)))) * (-l / (t * Math.sqrt(2.0)))));
	} else if ((t / l) <= 5e+97) {
		tmp = Math.asin(Math.sqrt(((1.0 - ((Om / Omc) / (Omc / Om))) / (1.0 + (2.0 * Math.pow((t / l), 2.0))))));
	} else {
		tmp = Math.asin(((l / t) * Math.sqrt(0.5)));
	}
	return tmp;
}
def code(t, l, Om, Omc):
	return math.asin(math.sqrt(((1.0 - math.pow((Om / Omc), 2.0)) / (1.0 + (2.0 * math.pow((t / l), 2.0))))))
def code(t, l, Om, Omc):
	tmp = 0
	if (t / l) <= -5e+131:
		tmp = math.asin((math.sqrt((1.0 - ((Om / Omc) * (Om / Omc)))) * (-l / (t * math.sqrt(2.0)))))
	elif (t / l) <= 5e+97:
		tmp = math.asin(math.sqrt(((1.0 - ((Om / Omc) / (Omc / Om))) / (1.0 + (2.0 * math.pow((t / l), 2.0))))))
	else:
		tmp = math.asin(((l / t) * math.sqrt(0.5)))
	return tmp
function code(t, l, Om, Omc)
	return asin(sqrt(Float64(Float64(1.0 - (Float64(Om / Omc) ^ 2.0)) / Float64(1.0 + Float64(2.0 * (Float64(t / l) ^ 2.0))))))
end
function code(t, l, Om, Omc)
	tmp = 0.0
	if (Float64(t / l) <= -5e+131)
		tmp = asin(Float64(sqrt(Float64(1.0 - Float64(Float64(Om / Omc) * Float64(Om / Omc)))) * Float64(Float64(-l) / Float64(t * sqrt(2.0)))));
	elseif (Float64(t / l) <= 5e+97)
		tmp = asin(sqrt(Float64(Float64(1.0 - Float64(Float64(Om / Omc) / Float64(Omc / Om))) / Float64(1.0 + Float64(2.0 * (Float64(t / l) ^ 2.0))))));
	else
		tmp = asin(Float64(Float64(l / t) * sqrt(0.5)));
	end
	return tmp
end
function tmp = code(t, l, Om, Omc)
	tmp = asin(sqrt(((1.0 - ((Om / Omc) ^ 2.0)) / (1.0 + (2.0 * ((t / l) ^ 2.0))))));
end
function tmp_2 = code(t, l, Om, Omc)
	tmp = 0.0;
	if ((t / l) <= -5e+131)
		tmp = asin((sqrt((1.0 - ((Om / Omc) * (Om / Omc)))) * (-l / (t * sqrt(2.0)))));
	elseif ((t / l) <= 5e+97)
		tmp = asin(sqrt(((1.0 - ((Om / Omc) / (Omc / Om))) / (1.0 + (2.0 * ((t / l) ^ 2.0))))));
	else
		tmp = asin(((l / t) * sqrt(0.5)));
	end
	tmp_2 = tmp;
end
code[t_, l_, Om_, Omc_] := N[ArcSin[N[Sqrt[N[(N[(1.0 - N[Power[N[(Om / Omc), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(2.0 * N[Power[N[(t / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision]
code[t_, l_, Om_, Omc_] := If[LessEqual[N[(t / l), $MachinePrecision], -5e+131], N[ArcSin[N[(N[Sqrt[N[(1.0 - N[(N[(Om / Omc), $MachinePrecision] * N[(Om / Omc), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision] * N[((-l) / N[(t * N[Sqrt[2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision], If[LessEqual[N[(t / l), $MachinePrecision], 5e+97], N[ArcSin[N[Sqrt[N[(N[(1.0 - N[(N[(Om / Omc), $MachinePrecision] / N[(Omc / Om), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(1.0 + N[(2.0 * N[Power[N[(t / l), $MachinePrecision], 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]], $MachinePrecision]], $MachinePrecision], N[ArcSin[N[(N[(l / t), $MachinePrecision] * N[Sqrt[0.5], $MachinePrecision]), $MachinePrecision]], $MachinePrecision]]]
\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)
\begin{array}{l}
\mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+131}:\\
\;\;\;\;\sin^{-1} \left(\sqrt{1 - \frac{Om}{Omc} \cdot \frac{Om}{Omc}} \cdot \frac{-\ell}{t \cdot \sqrt{2}}\right)\\

\mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+97}:\\
\;\;\;\;\sin^{-1} \left(\sqrt{\frac{1 - \frac{\frac{Om}{Omc}}{\frac{Omc}{Om}}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\\

\mathbf{else}:\\
\;\;\;\;\sin^{-1} \left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)\\


\end{array}

Local Percentage Accuracy vs ?

The average percentage accuracy by input value. Horizontal axis shows value of an input variable; the variable is choosen in the title. Vertical axis is accuracy; higher is better. Red represent the original program, while blue represents Herbie's suggestion. These can be toggled with buttons below the plot. The line is an average while dots represent individual samples.

Herbie found 14 alternatives:

AlternativeAccuracySpeedup

Accuracy vs Speed

The accuracy (vertical axis) and speed (horizontal axis) of each alternatives. Up and to the right is better. The red square shows the initial program, and each blue circle shows an alternative.The line shows the best available speed-accuracy tradeoffs.

Bogosity?

Bogosity

Try it out?

Your Program's Arguments

Results

Enter valid numbers for all inputs

Derivation?

  1. Split input into 3 regimes
  2. if (/.f64 t l) < -4.99999999999999995e131

    1. Initial program 46.7%

      \[\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right) \]
    2. Applied egg-rr97.8%

      \[\leadsto \sin^{-1} \color{blue}{\left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right)} \]
      Step-by-step derivation

      [Start]46.7%

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

      sqrt-div [=>]46.6%

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

      div-inv [=>]46.6%

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

      add-sqr-sqrt [=>]46.6%

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

      hypot-1-def [=>]46.6%

      \[ \sin^{-1} \left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\color{blue}{\mathsf{hypot}\left(1, \sqrt{2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}\right)}}\right) \]

      *-commutative [=>]46.6%

      \[ \sin^{-1} \left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\mathsf{hypot}\left(1, \sqrt{\color{blue}{{\left(\frac{t}{\ell}\right)}^{2} \cdot 2}}\right)}\right) \]

      sqrt-prod [=>]46.6%

      \[ \sin^{-1} \left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\mathsf{hypot}\left(1, \color{blue}{\sqrt{{\left(\frac{t}{\ell}\right)}^{2}} \cdot \sqrt{2}}\right)}\right) \]

      unpow2 [=>]46.6%

      \[ \sin^{-1} \left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\mathsf{hypot}\left(1, \sqrt{\color{blue}{\frac{t}{\ell} \cdot \frac{t}{\ell}}} \cdot \sqrt{2}\right)}\right) \]

      sqrt-prod [=>]0.0%

      \[ \sin^{-1} \left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\mathsf{hypot}\left(1, \color{blue}{\left(\sqrt{\frac{t}{\ell}} \cdot \sqrt{\frac{t}{\ell}}\right)} \cdot \sqrt{2}\right)}\right) \]

      add-sqr-sqrt [<=]97.8%

      \[ \sin^{-1} \left(\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot \frac{1}{\mathsf{hypot}\left(1, \color{blue}{\frac{t}{\ell}} \cdot \sqrt{2}\right)}\right) \]
    3. Simplified97.9%

      \[\leadsto \sin^{-1} \color{blue}{\left(\frac{\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}}}{\mathsf{hypot}\left(1, \frac{t}{\frac{\ell}{\sqrt{2}}}\right)}\right)} \]
      Step-by-step derivation

      [Start]97.8%

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

      associate-*r/ [=>]97.8%

      \[ \sin^{-1} \color{blue}{\left(\frac{\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}} \cdot 1}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right)} \]

      unpow2 [=>]97.8%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \color{blue}{\frac{Om}{Omc} \cdot \frac{Om}{Omc}}} \cdot 1}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      times-frac [<=]76.0%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \color{blue}{\frac{Om \cdot Om}{Omc \cdot Omc}}} \cdot 1}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      unpow2 [<=]76.0%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \frac{\color{blue}{{Om}^{2}}}{Omc \cdot Omc}} \cdot 1}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      unpow2 [<=]76.0%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \frac{{Om}^{2}}{\color{blue}{{Omc}^{2}}}} \cdot 1}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      *-rgt-identity [=>]76.0%

      \[ \sin^{-1} \left(\frac{\color{blue}{\sqrt{1 - \frac{{Om}^{2}}{{Omc}^{2}}}}}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      unpow2 [=>]76.0%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \frac{\color{blue}{Om \cdot Om}}{{Omc}^{2}}}}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      unpow2 [=>]76.0%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \frac{Om \cdot Om}{\color{blue}{Omc \cdot Omc}}}}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      times-frac [=>]97.8%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \color{blue}{\frac{Om}{Omc} \cdot \frac{Om}{Omc}}}}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      unpow2 [<=]97.8%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - \color{blue}{{\left(\frac{Om}{Omc}\right)}^{2}}}}{\mathsf{hypot}\left(1, \frac{t}{\ell} \cdot \sqrt{2}\right)}\right) \]

      associate-*l/ [=>]97.8%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}}}{\mathsf{hypot}\left(1, \color{blue}{\frac{t \cdot \sqrt{2}}{\ell}}\right)}\right) \]

      associate-/l* [=>]97.9%

      \[ \sin^{-1} \left(\frac{\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}}}{\mathsf{hypot}\left(1, \color{blue}{\frac{t}{\frac{\ell}{\sqrt{2}}}}\right)}\right) \]
    4. Taylor expanded in t around -inf 77.0%

      \[\leadsto \sin^{-1} \color{blue}{\left(-1 \cdot \left(\frac{\ell}{\sqrt{2} \cdot t} \cdot \sqrt{1 - \frac{{Om}^{2}}{{Omc}^{2}}}\right)\right)} \]
    5. Simplified77.0%

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

      [Start]77.0%

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

      mul-1-neg [=>]77.0%

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

      *-commutative [=>]77.0%

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

      unpow2 [=>]77.0%

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

      unpow2 [=>]77.0%

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

      times-frac [=>]99.6%

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

      unpow2 [<=]99.6%

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

      unpow2 [=>]99.6%

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

      times-frac [<=]77.0%

      \[ \sin^{-1} \left(-\sqrt{1 - \color{blue}{\frac{Om \cdot Om}{Omc \cdot Omc}}} \cdot \frac{\ell}{\sqrt{2} \cdot t}\right) \]
    6. Applied egg-rr99.6%

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

      [Start]77.0%

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

      times-frac [=>]99.6%

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

    if -4.99999999999999995e131 < (/.f64 t l) < 4.99999999999999999e97

    1. Initial program 98.9%

      \[\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right) \]
    2. Applied egg-rr98.9%

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

      [Start]98.9%

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

      unpow2 [=>]98.9%

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

      clear-num [=>]98.9%

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

      un-div-inv [=>]98.9%

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

    if 4.99999999999999999e97 < (/.f64 t l)

    1. Initial program 41.5%

      \[\sin^{-1} \left(\sqrt{\frac{1 - {\left(\frac{Om}{Omc}\right)}^{2}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right) \]
    2. Taylor expanded in Om around 0 28.7%

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

      \[\leadsto \sin^{-1} \color{blue}{\left(\sqrt{\frac{1}{1 + 2 \cdot \frac{t \cdot t}{\ell \cdot \ell}}}\right)} \]
      Step-by-step derivation

      [Start]28.7%

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

      unpow2 [=>]28.7%

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

      unpow2 [=>]28.7%

      \[ \sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \frac{t \cdot t}{\color{blue}{\ell \cdot \ell}}}}\right) \]
    4. Applied egg-rr28.7%

      \[\leadsto \sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \color{blue}{{\left(\frac{\ell \cdot \ell}{t \cdot t}\right)}^{-1}}}}\right) \]
      Step-by-step derivation

      [Start]28.7%

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

      clear-num [=>]28.7%

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

      inv-pow [=>]28.7%

      \[ \sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \color{blue}{{\left(\frac{\ell \cdot \ell}{t \cdot t}\right)}^{-1}}}}\right) \]
    5. Simplified41.6%

      \[\leadsto \sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \color{blue}{\frac{1}{\frac{\ell}{t} \cdot \frac{\ell}{t}}}}}\right) \]
      Step-by-step derivation

      [Start]28.7%

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

      unpow-1 [=>]28.7%

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

      times-frac [=>]41.6%

      \[ \sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \frac{1}{\color{blue}{\frac{\ell}{t} \cdot \frac{\ell}{t}}}}}\right) \]
    6. Taylor expanded in l around 0 99.7%

      \[\leadsto \sin^{-1} \color{blue}{\left(\frac{\sqrt{0.5} \cdot \ell}{t}\right)} \]
    7. Simplified99.7%

      \[\leadsto \sin^{-1} \color{blue}{\left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)} \]
      Step-by-step derivation

      [Start]99.7%

      \[ \sin^{-1} \left(\frac{\sqrt{0.5} \cdot \ell}{t}\right) \]

      *-commutative [=>]99.7%

      \[ \sin^{-1} \left(\frac{\color{blue}{\ell \cdot \sqrt{0.5}}}{t}\right) \]

      associate-/l* [=>]99.6%

      \[ \sin^{-1} \color{blue}{\left(\frac{\ell}{\frac{t}{\sqrt{0.5}}}\right)} \]

      associate-/r/ [=>]99.7%

      \[ \sin^{-1} \color{blue}{\left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)} \]
  3. Recombined 3 regimes into one program.
  4. Final simplification99.2%

    \[\leadsto \begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+131}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{1 - \frac{Om}{Omc} \cdot \frac{Om}{Omc}} \cdot \frac{-\ell}{t \cdot \sqrt{2}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+97}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{\frac{1 - \frac{\frac{Om}{Omc}}{\frac{Omc}{Om}}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)\\ \end{array} \]

Alternatives

Alternative 1
Accuracy98.8%
Cost20872
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+131}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{1 - \frac{Om}{Omc} \cdot \frac{Om}{Omc}} \cdot \frac{-\ell}{t \cdot \sqrt{2}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+97}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{\frac{1 - \frac{\frac{Om}{Omc}}{\frac{Omc}{Om}}}{1 + 2 \cdot {\left(\frac{t}{\ell}\right)}^{2}}}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)\\ \end{array} \]
Alternative 2
Accuracy98.4%
Cost32832
\[\sin^{-1} \left(\frac{\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}}}{\mathsf{hypot}\left(1, \frac{t}{\frac{\ell}{\sqrt{2}}}\right)}\right) \]
Alternative 3
Accuracy98.4%
Cost32832
\[\sin^{-1} \left(\frac{\sqrt{1 - {\left(\frac{Om}{Omc}\right)}^{2}}}{\mathsf{hypot}\left(1, \sqrt{2} \cdot \frac{t}{\ell}\right)}\right) \]
Alternative 4
Accuracy98.0%
Cost20484
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+147}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{1 - \frac{Om}{Omc} \cdot \frac{Om}{Omc}} \cdot \frac{-\ell}{t \cdot \sqrt{2}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+80}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \frac{1}{\frac{\ell}{t} \cdot \frac{\ell}{t}}}}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell \cdot \sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 5
Accuracy97.8%
Cost14280
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+155}:\\ \;\;\;\;\sin^{-1} \left(\frac{\sqrt{0.5}}{t} \cdot \left(-\ell\right)\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+80}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \frac{1}{\frac{\ell}{t} \cdot \frac{\ell}{t}}}}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell \cdot \sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 6
Accuracy97.8%
Cost14152
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+155}:\\ \;\;\;\;\sin^{-1} \left(\frac{\sqrt{0.5}}{t} \cdot \left(-\ell\right)\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 5 \cdot 10^{+97}:\\ \;\;\;\;\sin^{-1} \left(\sqrt{\frac{1}{1 + 2 \cdot \left(\frac{t}{\ell} \cdot \frac{t}{\ell}\right)}}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)\\ \end{array} \]
Alternative 7
Accuracy97.0%
Cost13704
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -50:\\ \;\;\;\;\sin^{-1} \left(\frac{\sqrt{0.5}}{t} \cdot \left(-\ell\right)\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 0.02:\\ \;\;\;\;\sin^{-1} \left(1 - {\left(\frac{t}{\ell}\right)}^{2}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell \cdot \sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 8
Accuracy79.4%
Cost13640
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+213}:\\ \;\;\;\;\sin^{-1} \left(\frac{\sqrt{0.5}}{\frac{t}{\ell}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 0.02:\\ \;\;\;\;\sin^{-1} 1\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\ell \cdot \frac{\sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 9
Accuracy79.4%
Cost13640
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -5 \cdot 10^{+213}:\\ \;\;\;\;\sin^{-1} \left(\frac{\sqrt{0.5}}{\frac{t}{\ell}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 0.02:\\ \;\;\;\;\sin^{-1} 1\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell \cdot \sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 10
Accuracy96.6%
Cost13640
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -50:\\ \;\;\;\;\sin^{-1} \left(\frac{-\sqrt{0.5}}{\frac{t}{\ell}}\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 0.02:\\ \;\;\;\;\sin^{-1} 1\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell \cdot \sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 11
Accuracy96.8%
Cost13640
\[\begin{array}{l} \mathbf{if}\;\frac{t}{\ell} \leq -50:\\ \;\;\;\;\sin^{-1} \left(\frac{\sqrt{0.5}}{t} \cdot \left(-\ell\right)\right)\\ \mathbf{elif}\;\frac{t}{\ell} \leq 0.02:\\ \;\;\;\;\sin^{-1} 1\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} \left(\frac{\ell \cdot \sqrt{0.5}}{t}\right)\\ \end{array} \]
Alternative 12
Accuracy62.9%
Cost13385
\[\begin{array}{l} \mathbf{if}\;t \leq -1.5 \cdot 10^{+90} \lor \neg \left(t \leq 8.5 \cdot 10^{+103}\right):\\ \;\;\;\;\sin^{-1} \left(\ell \cdot \frac{\sqrt{0.5}}{t}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} 1\\ \end{array} \]
Alternative 13
Accuracy62.9%
Cost13385
\[\begin{array}{l} \mathbf{if}\;t \leq -2.5 \cdot 10^{+83} \lor \neg \left(t \leq 1.1 \cdot 10^{+104}\right):\\ \;\;\;\;\sin^{-1} \left(\frac{\ell}{t} \cdot \sqrt{0.5}\right)\\ \mathbf{else}:\\ \;\;\;\;\sin^{-1} 1\\ \end{array} \]
Alternative 14
Accuracy49.9%
Cost6464
\[\sin^{-1} 1 \]

Reproduce?

herbie shell --seed 2023263 
(FPCore (t l Om Omc)
  :name "Toniolo and Linder, Equation (2)"
  :precision binary64
  (asin (sqrt (/ (- 1.0 (pow (/ Om Omc) 2.0)) (+ 1.0 (* 2.0 (pow (/ t l) 2.0)))))))