\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 -2.739921110099047161246117288663608441381 \cdot 10^{-223}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\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{elif}\;U* \le 6.078150703052965766771052983090706385029 \cdot 10^{-288}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - 0\right)}\\
\mathbf{elif}\;U* \le 1.246666962673744573461126104160207795017 \cdot 10^{-74}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\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{elif}\;U* \le 1.266602938111869964038515141872421636651 \cdot 10^{-33}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\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)}\\
\mathbf{elif}\;U* \le 6.570671311497278466304286938757794883486 \cdot 10^{73}:\\
\;\;\;\;\sqrt{{\left(\left(\left(2 \cdot n\right) \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)\right)\right) \cdot U\right)}^{1}}\\
\mathbf{elif}\;U* \le 7.601947622374171152660894727648559871279 \cdot 10^{174}:\\
\;\;\;\;\sqrt{\left(\sqrt[3]{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)} \cdot \sqrt[3]{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)}\right) \cdot \sqrt[3]{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left(\left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)}}\\
\mathbf{elif}\;U* \le 2.543190609382229994294004240158197999196 \cdot 10^{261}:\\
\;\;\;\;\sqrt{{\left(\left(\left(2 \cdot n\right) \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)\right)\right) \cdot U\right)}^{1}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\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)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r210708 = 2.0;
double r210709 = n;
double r210710 = r210708 * r210709;
double r210711 = U;
double r210712 = r210710 * r210711;
double r210713 = t;
double r210714 = l;
double r210715 = r210714 * r210714;
double r210716 = Om;
double r210717 = r210715 / r210716;
double r210718 = r210708 * r210717;
double r210719 = r210713 - r210718;
double r210720 = r210714 / r210716;
double r210721 = pow(r210720, r210708);
double r210722 = r210709 * r210721;
double r210723 = U_;
double r210724 = r210711 - r210723;
double r210725 = r210722 * r210724;
double r210726 = r210719 - r210725;
double r210727 = r210712 * r210726;
double r210728 = sqrt(r210727);
return r210728;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r210729 = U_;
double r210730 = -2.739921110099047e-223;
bool r210731 = r210729 <= r210730;
double r210732 = 2.0;
double r210733 = n;
double r210734 = r210732 * r210733;
double r210735 = U;
double r210736 = t;
double r210737 = l;
double r210738 = Om;
double r210739 = r210737 / r210738;
double r210740 = r210737 * r210739;
double r210741 = r210732 * r210740;
double r210742 = r210736 - r210741;
double r210743 = 2.0;
double r210744 = r210732 / r210743;
double r210745 = pow(r210739, r210744);
double r210746 = r210733 * r210745;
double r210747 = r210735 - r210729;
double r210748 = r210745 * r210747;
double r210749 = r210746 * r210748;
double r210750 = r210742 - r210749;
double r210751 = r210735 * r210750;
double r210752 = r210734 * r210751;
double r210753 = sqrt(r210752);
double r210754 = 6.078150703052966e-288;
bool r210755 = r210729 <= r210754;
double r210756 = r210734 * r210735;
double r210757 = 0.0;
double r210758 = r210742 - r210757;
double r210759 = r210756 * r210758;
double r210760 = sqrt(r210759);
double r210761 = 1.2466669626737446e-74;
bool r210762 = r210729 <= r210761;
double r210763 = 1.26660293811187e-33;
bool r210764 = r210729 <= r210763;
double r210765 = sqrt(r210756);
double r210766 = sqrt(r210750);
double r210767 = r210765 * r210766;
double r210768 = 6.570671311497278e+73;
bool r210769 = r210729 <= r210768;
double r210770 = r210743 * r210744;
double r210771 = pow(r210739, r210770);
double r210772 = r210733 * r210771;
double r210773 = -r210772;
double r210774 = r210773 * r210747;
double r210775 = r210741 - r210774;
double r210776 = r210736 - r210775;
double r210777 = r210734 * r210776;
double r210778 = r210777 * r210735;
double r210779 = 1.0;
double r210780 = pow(r210778, r210779);
double r210781 = sqrt(r210780);
double r210782 = 7.601947622374171e+174;
bool r210783 = r210729 <= r210782;
double r210784 = r210746 * r210745;
double r210785 = r210784 * r210747;
double r210786 = r210742 - r210785;
double r210787 = r210756 * r210786;
double r210788 = cbrt(r210787);
double r210789 = r210788 * r210788;
double r210790 = r210789 * r210788;
double r210791 = sqrt(r210790);
double r210792 = 2.54319060938223e+261;
bool r210793 = r210729 <= r210792;
double r210794 = r210793 ? r210781 : r210753;
double r210795 = r210783 ? r210791 : r210794;
double r210796 = r210769 ? r210781 : r210795;
double r210797 = r210764 ? r210767 : r210796;
double r210798 = r210762 ? r210753 : r210797;
double r210799 = r210755 ? r210760 : r210798;
double r210800 = r210731 ? r210753 : r210799;
return r210800;
}



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* < -2.739921110099047e-223 or 6.078150703052966e-288 < U* < 1.2466669626737446e-74 or 2.54319060938223e+261 < U* Initial program 34.3
rmApplied *-un-lft-identity34.3
Applied times-frac31.7
Simplified31.7
rmApplied sqr-pow31.7
Applied associate-*r*30.8
rmApplied associate-*l*30.6
rmApplied associate-*l*30.6
if -2.739921110099047e-223 < U* < 6.078150703052966e-288Initial program 32.8
rmApplied *-un-lft-identity32.8
Applied times-frac28.8
Simplified28.8
Taylor expanded around 0 27.5
if 1.2466669626737446e-74 < U* < 1.26660293811187e-33Initial program 35.7
rmApplied *-un-lft-identity35.7
Applied times-frac33.5
Simplified33.5
rmApplied sqr-pow33.5
Applied associate-*r*33.1
rmApplied associate-*l*32.6
rmApplied sqrt-prod44.3
if 1.26660293811187e-33 < U* < 6.570671311497278e+73 or 7.601947622374171e+174 < U* < 2.54319060938223e+261Initial program 35.4
rmApplied *-un-lft-identity35.4
Applied times-frac32.8
Simplified32.8
rmApplied sqr-pow32.8
Applied associate-*r*32.1
rmApplied pow132.1
Applied pow132.1
Applied pow132.1
Applied pow132.1
Applied pow-prod-down32.1
Applied pow-prod-down32.1
Applied pow-prod-down32.1
Simplified32.4
if 6.570671311497278e+73 < U* < 7.601947622374171e+174Initial program 34.1
rmApplied *-un-lft-identity34.1
Applied times-frac32.4
Simplified32.4
rmApplied sqr-pow32.4
Applied associate-*r*31.6
rmApplied add-cube-cbrt32.0
Final simplification31.2
herbie shell --seed 2019352
(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*))))))