\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 -9.379719815826530589603530820932364731187 \cdot 10^{241}:\\
\;\;\;\;\sqrt{\left(\left(2 \cdot n\right) \cdot U\right) \cdot \left(t - \frac{\ell}{Om} \cdot \left(\ell \cdot 2\right)\right)}\\
\mathbf{elif}\;t \le -3.772245666195913577646904748734404208118 \cdot 10^{-180}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(\left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right)\right) - \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\right) \cdot \left({\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot \left(U - U*\right)\right)\right)\right)}\\
\mathbf{elif}\;t \le 1.39311700673337772156748504137624422583 \cdot 10^{-4}:\\
\;\;\;\;\sqrt{\left(\left(\left(-{\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\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) + \left(t - 2 \cdot \left(\ell \cdot \frac{\ell}{Om}\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(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)} \cdot n\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 r2791517 = 2.0;
double r2791518 = n;
double r2791519 = r2791517 * r2791518;
double r2791520 = U;
double r2791521 = r2791519 * r2791520;
double r2791522 = t;
double r2791523 = l;
double r2791524 = r2791523 * r2791523;
double r2791525 = Om;
double r2791526 = r2791524 / r2791525;
double r2791527 = r2791517 * r2791526;
double r2791528 = r2791522 - r2791527;
double r2791529 = r2791523 / r2791525;
double r2791530 = pow(r2791529, r2791517);
double r2791531 = r2791518 * r2791530;
double r2791532 = U_;
double r2791533 = r2791520 - r2791532;
double r2791534 = r2791531 * r2791533;
double r2791535 = r2791528 - r2791534;
double r2791536 = r2791521 * r2791535;
double r2791537 = sqrt(r2791536);
return r2791537;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r2791538 = t;
double r2791539 = -9.37971981582653e+241;
bool r2791540 = r2791538 <= r2791539;
double r2791541 = 2.0;
double r2791542 = n;
double r2791543 = r2791541 * r2791542;
double r2791544 = U;
double r2791545 = r2791543 * r2791544;
double r2791546 = l;
double r2791547 = Om;
double r2791548 = r2791546 / r2791547;
double r2791549 = r2791546 * r2791541;
double r2791550 = r2791548 * r2791549;
double r2791551 = r2791538 - r2791550;
double r2791552 = r2791545 * r2791551;
double r2791553 = sqrt(r2791552);
double r2791554 = -3.7722456661959136e-180;
bool r2791555 = r2791538 <= r2791554;
double r2791556 = r2791546 * r2791548;
double r2791557 = r2791541 * r2791556;
double r2791558 = r2791538 - r2791557;
double r2791559 = 2.0;
double r2791560 = r2791541 / r2791559;
double r2791561 = pow(r2791548, r2791560);
double r2791562 = r2791561 * r2791542;
double r2791563 = U_;
double r2791564 = r2791544 - r2791563;
double r2791565 = r2791561 * r2791564;
double r2791566 = r2791562 * r2791565;
double r2791567 = r2791558 - r2791566;
double r2791568 = r2791544 * r2791567;
double r2791569 = r2791543 * r2791568;
double r2791570 = sqrt(r2791569);
double r2791571 = 0.00013931170067333777;
bool r2791572 = r2791538 <= r2791571;
double r2791573 = -r2791561;
double r2791574 = r2791573 * r2791564;
double r2791575 = r2791574 * r2791562;
double r2791576 = r2791575 * r2791545;
double r2791577 = r2791558 * r2791545;
double r2791578 = r2791576 + r2791577;
double r2791579 = sqrt(r2791578);
double r2791580 = sqrt(r2791545);
double r2791581 = sqrt(r2791567);
double r2791582 = r2791580 * r2791581;
double r2791583 = r2791572 ? r2791579 : r2791582;
double r2791584 = r2791555 ? r2791570 : r2791583;
double r2791585 = r2791540 ? r2791553 : r2791584;
return r2791585;
}



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 < -9.37971981582653e+241Initial program 41.0
rmApplied sqr-pow41.0
Applied associate-*r*40.8
rmApplied *-un-lft-identity40.8
Applied times-frac38.5
Simplified38.5
Taylor expanded around 0 40.3
Simplified37.9
if -9.37971981582653e+241 < t < -3.7722456661959136e-180Initial program 32.6
rmApplied sqr-pow32.6
Applied associate-*r*31.8
rmApplied *-un-lft-identity31.8
Applied times-frac29.0
Simplified29.0
rmApplied associate-*l*29.1
rmApplied associate-*l*28.8
if -3.7722456661959136e-180 < t < 0.00013931170067333777Initial program 35.0
rmApplied sqr-pow35.0
Applied associate-*r*33.9
rmApplied *-un-lft-identity33.9
Applied times-frac30.6
Simplified30.6
rmApplied associate-*l*29.8
rmApplied sub-neg29.8
Applied distribute-lft-in29.8
if 0.00013931170067333777 < t Initial program 34.4
rmApplied sqr-pow34.4
Applied associate-*r*34.0
rmApplied *-un-lft-identity34.0
Applied times-frac31.7
Simplified31.7
rmApplied associate-*l*32.0
rmApplied sqrt-prod27.2
Final simplification29.2
herbie shell --seed 2019172
(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*))))))