\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 r177962 = 2.0;
double r177963 = n;
double r177964 = r177962 * r177963;
double r177965 = U;
double r177966 = r177964 * r177965;
double r177967 = t;
double r177968 = l;
double r177969 = r177968 * r177968;
double r177970 = Om;
double r177971 = r177969 / r177970;
double r177972 = r177962 * r177971;
double r177973 = r177967 - r177972;
double r177974 = r177968 / r177970;
double r177975 = pow(r177974, r177962);
double r177976 = r177963 * r177975;
double r177977 = U_;
double r177978 = r177965 - r177977;
double r177979 = r177976 * r177978;
double r177980 = r177973 - r177979;
double r177981 = r177966 * r177980;
double r177982 = sqrt(r177981);
return r177982;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r177983 = t;
double r177984 = -5.661273572231306e-263;
bool r177985 = r177983 <= r177984;
double r177986 = 2.0;
double r177987 = n;
double r177988 = r177986 * r177987;
double r177989 = U;
double r177990 = l;
double r177991 = Om;
double r177992 = r177991 / r177990;
double r177993 = r177990 / r177992;
double r177994 = r177986 * r177993;
double r177995 = r177983 - r177994;
double r177996 = r177990 / r177991;
double r177997 = 2.0;
double r177998 = r177986 / r177997;
double r177999 = pow(r177996, r177998);
double r178000 = r177987 * r177999;
double r178001 = U_;
double r178002 = r177989 - r178001;
double r178003 = r177999 * r178002;
double r178004 = r178000 * r178003;
double r178005 = r177995 - r178004;
double r178006 = r177989 * r178005;
double r178007 = r177988 * r178006;
double r178008 = sqrt(r178007);
double r178009 = r177988 * r177989;
double r178010 = sqrt(r178009);
double r178011 = sqrt(r178005);
double r178012 = r178010 * r178011;
double r178013 = r177985 ? r178008 : r178012;
return r178013;
}



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*))))))