\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell \cdot \ell}{Om}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right) \cdot \left(U - U*\right)\right)}\begin{array}{l}
\mathbf{if}\;t \le -5.661273572231306187747529113901313439198 \cdot 10^{-263}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r215790 = 2.0;
double r215791 = n;
double r215792 = r215790 * r215791;
double r215793 = U;
double r215794 = r215792 * r215793;
double r215795 = t;
double r215796 = l;
double r215797 = r215796 * r215796;
double r215798 = Om;
double r215799 = r215797 / r215798;
double r215800 = r215790 * r215799;
double r215801 = r215795 - r215800;
double r215802 = r215796 / r215798;
double r215803 = pow(r215802, r215790);
double r215804 = r215791 * r215803;
double r215805 = U_;
double r215806 = r215793 - r215805;
double r215807 = r215804 * r215806;
double r215808 = r215801 - r215807;
double r215809 = r215794 * r215808;
double r215810 = sqrt(r215809);
return r215810;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r215811 = t;
double r215812 = -5.661273572231306e-263;
bool r215813 = r215811 <= r215812;
double r215814 = 2.0;
double r215815 = n;
double r215816 = r215814 * r215815;
double r215817 = U;
double r215818 = l;
double r215819 = Om;
double r215820 = r215819 / r215818;
double r215821 = r215818 / r215820;
double r215822 = r215814 * r215821;
double r215823 = r215811 - r215822;
double r215824 = r215818 / r215819;
double r215825 = 2.0;
double r215826 = r215814 / r215825;
double r215827 = pow(r215824, r215826);
double r215828 = r215815 * r215827;
double r215829 = U_;
double r215830 = r215817 - r215829;
double r215831 = r215827 * r215830;
double r215832 = r215828 * r215831;
double r215833 = r215823 - r215832;
double r215834 = r215817 * r215833;
double r215835 = r215816 * r215834;
double r215836 = sqrt(r215835);
double r215837 = r215816 * r215817;
double r215838 = sqrt(r215837);
double r215839 = sqrt(r215833);
double r215840 = r215838 * r215839;
double r215841 = r215813 ? r215836 : r215840;
return r215841;
}



Bits error versus n



Bits error versus U



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus U*
Results
if t < -5.661273572231306e-263Initial program 35.0
rmApplied associate-/l*32.4
rmApplied sqr-pow32.4
Applied associate-*r*31.5
rmApplied associate-*l*31.5
rmApplied associate-*l*30.7
if -5.661273572231306e-263 < t Initial program 35.0
rmApplied associate-/l*32.4
rmApplied sqr-pow32.4
Applied associate-*r*31.4
rmApplied associate-*l*31.2
rmApplied sqrt-prod29.9
Final simplification30.3
herbie shell --seed 2019291
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
:precision binary64
(sqrt (* (* (* 2 n) U) (- (- t (* 2 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2)) (- U U*))))))