\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}\;U \le -9.95931678236544397 \cdot 10^{54}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\ell}{\frac{Om}{\ell}}\right) - 0\right)}\\
\mathbf{elif}\;U \le 7.7008530070172207 \cdot 10^{77}:\\
\;\;\;\;\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{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \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)}} \cdot \sqrt{\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \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)}}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r215899 = 2.0;
double r215900 = n;
double r215901 = r215899 * r215900;
double r215902 = U;
double r215903 = r215901 * r215902;
double r215904 = t;
double r215905 = l;
double r215906 = r215905 * r215905;
double r215907 = Om;
double r215908 = r215906 / r215907;
double r215909 = r215899 * r215908;
double r215910 = r215904 - r215909;
double r215911 = r215905 / r215907;
double r215912 = pow(r215911, r215899);
double r215913 = r215900 * r215912;
double r215914 = U_;
double r215915 = r215902 - r215914;
double r215916 = r215913 * r215915;
double r215917 = r215910 - r215916;
double r215918 = r215903 * r215917;
double r215919 = sqrt(r215918);
return r215919;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r215920 = U;
double r215921 = -9.959316782365444e+54;
bool r215922 = r215920 <= r215921;
double r215923 = 2.0;
double r215924 = n;
double r215925 = r215923 * r215924;
double r215926 = r215925 * r215920;
double r215927 = t;
double r215928 = l;
double r215929 = Om;
double r215930 = r215929 / r215928;
double r215931 = r215928 / r215930;
double r215932 = r215923 * r215931;
double r215933 = r215927 - r215932;
double r215934 = 0.0;
double r215935 = r215933 - r215934;
double r215936 = r215926 * r215935;
double r215937 = sqrt(r215936);
double r215938 = 7.700853007017221e+77;
bool r215939 = r215920 <= r215938;
double r215940 = r215928 / r215929;
double r215941 = 2.0;
double r215942 = r215923 / r215941;
double r215943 = pow(r215940, r215942);
double r215944 = r215924 * r215943;
double r215945 = U_;
double r215946 = r215920 - r215945;
double r215947 = r215943 * r215946;
double r215948 = r215944 * r215947;
double r215949 = r215933 - r215948;
double r215950 = r215920 * r215949;
double r215951 = r215925 * r215950;
double r215952 = sqrt(r215951);
double r215953 = r215926 * r215949;
double r215954 = sqrt(r215953);
double r215955 = sqrt(r215954);
double r215956 = r215955 * r215955;
double r215957 = r215939 ? r215952 : r215956;
double r215958 = r215922 ? r215937 : r215957;
return r215958;
}



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 U < -9.959316782365444e+54Initial program 29.3
rmApplied associate-/l*26.8
Taylor expanded around 0 28.1
if -9.959316782365444e+54 < U < 7.700853007017221e+77Initial program 36.9
rmApplied associate-/l*34.0
rmApplied sqr-pow34.0
Applied associate-*r*33.0
rmApplied associate-*l*32.5
rmApplied associate-*l*28.9
if 7.700853007017221e+77 < U Initial program 29.8
rmApplied associate-/l*27.4
rmApplied sqr-pow27.4
Applied associate-*r*27.1
rmApplied associate-*l*27.7
rmApplied add-sqr-sqrt27.9
Final simplification28.7
herbie shell --seed 2020081
(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*))))))