\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 -4.458146329554135732988892774226101337842 \cdot 10^{62}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \frac{\frac{\ell}{\sqrt[3]{Om}} \cdot \frac{\ell}{\sqrt[3]{Om}}}{\sqrt[3]{Om}}\right) - \left(U - U*\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{2}\right)\right)}\\
\mathbf{elif}\;U \le 1.43075486013443336421549580762136466865 \cdot 10^{-67}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(\left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \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 n\right)\right) \cdot U\right)}\\
\mathbf{elif}\;U \le 3.638498443778453206704785881478597925308 \cdot 10^{-26}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \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 n\right)}\\
\mathbf{elif}\;U \le 1.219991701569961735518499105488225569099 \cdot 10^{291}:\\
\;\;\;\;\sqrt{\sqrt{\left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \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 n\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}} \cdot \sqrt{\sqrt{\left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \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 n\right)\right) \cdot \left(\left(2 \cdot n\right) \cdot U\right)}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(U - U*\right) \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 n\right)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r3785986 = 2.0;
double r3785987 = n;
double r3785988 = r3785986 * r3785987;
double r3785989 = U;
double r3785990 = r3785988 * r3785989;
double r3785991 = t;
double r3785992 = l;
double r3785993 = r3785992 * r3785992;
double r3785994 = Om;
double r3785995 = r3785993 / r3785994;
double r3785996 = r3785986 * r3785995;
double r3785997 = r3785991 - r3785996;
double r3785998 = r3785992 / r3785994;
double r3785999 = pow(r3785998, r3785986);
double r3786000 = r3785987 * r3785999;
double r3786001 = U_;
double r3786002 = r3785989 - r3786001;
double r3786003 = r3786000 * r3786002;
double r3786004 = r3785997 - r3786003;
double r3786005 = r3785990 * r3786004;
double r3786006 = sqrt(r3786005);
return r3786006;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r3786007 = U;
double r3786008 = -4.458146329554136e+62;
bool r3786009 = r3786007 <= r3786008;
double r3786010 = 2.0;
double r3786011 = n;
double r3786012 = r3786010 * r3786011;
double r3786013 = r3786012 * r3786007;
double r3786014 = t;
double r3786015 = l;
double r3786016 = Om;
double r3786017 = cbrt(r3786016);
double r3786018 = r3786015 / r3786017;
double r3786019 = r3786018 * r3786018;
double r3786020 = r3786019 / r3786017;
double r3786021 = r3786010 * r3786020;
double r3786022 = r3786014 - r3786021;
double r3786023 = U_;
double r3786024 = r3786007 - r3786023;
double r3786025 = r3786015 / r3786016;
double r3786026 = pow(r3786025, r3786010);
double r3786027 = r3786011 * r3786026;
double r3786028 = r3786024 * r3786027;
double r3786029 = r3786022 - r3786028;
double r3786030 = r3786013 * r3786029;
double r3786031 = sqrt(r3786030);
double r3786032 = 1.4307548601344334e-67;
bool r3786033 = r3786007 <= r3786032;
double r3786034 = r3786015 * r3786025;
double r3786035 = r3786010 * r3786034;
double r3786036 = r3786014 - r3786035;
double r3786037 = 2.0;
double r3786038 = r3786010 / r3786037;
double r3786039 = pow(r3786025, r3786038);
double r3786040 = r3786024 * r3786039;
double r3786041 = r3786039 * r3786011;
double r3786042 = r3786040 * r3786041;
double r3786043 = r3786036 - r3786042;
double r3786044 = r3786043 * r3786007;
double r3786045 = r3786012 * r3786044;
double r3786046 = sqrt(r3786045);
double r3786047 = 3.638498443778453e-26;
bool r3786048 = r3786007 <= r3786047;
double r3786049 = sqrt(r3786013);
double r3786050 = sqrt(r3786043);
double r3786051 = r3786049 * r3786050;
double r3786052 = 1.2199917015699617e+291;
bool r3786053 = r3786007 <= r3786052;
double r3786054 = r3786043 * r3786013;
double r3786055 = sqrt(r3786054);
double r3786056 = sqrt(r3786055);
double r3786057 = r3786056 * r3786056;
double r3786058 = r3786053 ? r3786057 : r3786051;
double r3786059 = r3786048 ? r3786051 : r3786058;
double r3786060 = r3786033 ? r3786046 : r3786059;
double r3786061 = r3786009 ? r3786031 : r3786060;
return r3786061;
}



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 < -4.458146329554136e+62Initial program 30.0
rmApplied add-cube-cbrt30.1
Applied associate-/r*30.1
Simplified27.8
if -4.458146329554136e+62 < U < 1.4307548601344334e-67Initial program 37.2
rmApplied *-un-lft-identity37.2
Applied times-frac34.4
Simplified34.4
rmApplied sqr-pow34.4
Applied associate-*r*33.3
rmApplied associate-*l*32.9
rmApplied associate-*l*29.3
if 1.4307548601344334e-67 < U < 3.638498443778453e-26 or 1.2199917015699617e+291 < U Initial program 30.9
rmApplied *-un-lft-identity30.9
Applied times-frac28.6
Simplified28.6
rmApplied sqr-pow28.6
Applied associate-*r*27.6
rmApplied associate-*l*28.2
rmApplied sqrt-prod41.8
if 3.638498443778453e-26 < U < 1.2199917015699617e+291Initial program 29.6
rmApplied *-un-lft-identity29.6
Applied times-frac26.1
Simplified26.1
rmApplied sqr-pow26.1
Applied associate-*r*25.4
rmApplied associate-*l*26.0
rmApplied add-sqr-sqrt26.2
Final simplification29.0
herbie shell --seed 2019168
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
(sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))