\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 -20628.158740343424:\\
\;\;\;\;\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 -3.38886052338754529 \cdot 10^{-282}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \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)}\\
\mathbf{elif}\;U* \le 2.66582174134472983 \cdot 10^{-120}:\\
\;\;\;\;\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 6.37844133114075821 \cdot 10^{-28}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot U} \cdot \sqrt{\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)}\\
\mathbf{elif}\;U* \le 1.3559215670195436 \cdot 10^{141}:\\
\;\;\;\;\sqrt{\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(\left(n \cdot \ell\right) \cdot {\left(\frac{1}{{Om}^{1}}\right)}^{1}\right) \cdot {\left(\frac{\ell}{Om}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \left(U - U*\right)\right)}\\
\mathbf{elif}\;U* \le 7.30429266906966537 \cdot 10^{231}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \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)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\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)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r214543 = 2.0;
double r214544 = n;
double r214545 = r214543 * r214544;
double r214546 = U;
double r214547 = r214545 * r214546;
double r214548 = t;
double r214549 = l;
double r214550 = r214549 * r214549;
double r214551 = Om;
double r214552 = r214550 / r214551;
double r214553 = r214543 * r214552;
double r214554 = r214548 - r214553;
double r214555 = r214549 / r214551;
double r214556 = pow(r214555, r214543);
double r214557 = r214544 * r214556;
double r214558 = U_;
double r214559 = r214546 - r214558;
double r214560 = r214557 * r214559;
double r214561 = r214554 - r214560;
double r214562 = r214547 * r214561;
double r214563 = sqrt(r214562);
return r214563;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r214564 = U_;
double r214565 = -20628.158740343424;
bool r214566 = r214564 <= r214565;
double r214567 = 2.0;
double r214568 = n;
double r214569 = r214567 * r214568;
double r214570 = t;
double r214571 = l;
double r214572 = Om;
double r214573 = r214571 / r214572;
double r214574 = r214571 * r214573;
double r214575 = r214567 * r214574;
double r214576 = 2.0;
double r214577 = r214567 / r214576;
double r214578 = r214576 * r214577;
double r214579 = pow(r214573, r214578);
double r214580 = r214568 * r214579;
double r214581 = -r214580;
double r214582 = U;
double r214583 = r214582 - r214564;
double r214584 = r214581 * r214583;
double r214585 = r214575 - r214584;
double r214586 = r214570 - r214585;
double r214587 = r214569 * r214586;
double r214588 = r214587 * r214582;
double r214589 = 1.0;
double r214590 = pow(r214588, r214589);
double r214591 = sqrt(r214590);
double r214592 = -3.3888605233875453e-282;
bool r214593 = r214564 <= r214592;
double r214594 = r214582 * r214586;
double r214595 = r214569 * r214594;
double r214596 = sqrt(r214595);
double r214597 = 2.66582174134473e-120;
bool r214598 = r214564 <= r214597;
double r214599 = 6.378441331140758e-28;
bool r214600 = r214564 <= r214599;
double r214601 = r214569 * r214582;
double r214602 = sqrt(r214601);
double r214603 = r214570 - r214575;
double r214604 = pow(r214573, r214577);
double r214605 = r214568 * r214604;
double r214606 = r214605 * r214604;
double r214607 = r214606 * r214583;
double r214608 = r214603 - r214607;
double r214609 = sqrt(r214608);
double r214610 = r214602 * r214609;
double r214611 = 1.3559215670195436e+141;
bool r214612 = r214564 <= r214611;
double r214613 = r214568 * r214571;
double r214614 = 1.0;
double r214615 = pow(r214572, r214614);
double r214616 = r214589 / r214615;
double r214617 = pow(r214616, r214614);
double r214618 = r214613 * r214617;
double r214619 = r214618 * r214604;
double r214620 = r214619 * r214583;
double r214621 = r214603 - r214620;
double r214622 = r214601 * r214621;
double r214623 = sqrt(r214622);
double r214624 = 7.304292669069665e+231;
bool r214625 = r214564 <= r214624;
double r214626 = r214601 * r214608;
double r214627 = sqrt(r214626);
double r214628 = r214625 ? r214596 : r214627;
double r214629 = r214612 ? r214623 : r214628;
double r214630 = r214600 ? r214610 : r214629;
double r214631 = r214598 ? r214591 : r214630;
double r214632 = r214593 ? r214596 : r214631;
double r214633 = r214566 ? r214591 : r214632;
return r214633;
}



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* < -20628.158740343424 or -3.3888605233875453e-282 < U* < 2.66582174134473e-120Initial program 35.7
rmApplied *-un-lft-identity35.7
Applied times-frac33.4
Simplified33.4
rmApplied sqr-pow33.4
Applied associate-*r*32.4
rmApplied pow132.4
Applied pow132.4
Applied pow132.4
Applied pow132.4
Applied pow-prod-down32.4
Applied pow-prod-down32.4
Applied pow-prod-down32.4
Simplified33.2
if -20628.158740343424 < U* < -3.3888605233875453e-282 or 1.3559215670195436e+141 < U* < 7.304292669069665e+231Initial program 33.3
rmApplied *-un-lft-identity33.3
Applied times-frac30.7
Simplified30.7
rmApplied sqr-pow30.7
Applied associate-*r*29.8
rmApplied associate-*l*30.6
Simplified31.8
if 2.66582174134473e-120 < U* < 6.378441331140758e-28Initial program 34.9
rmApplied *-un-lft-identity34.9
Applied times-frac31.7
Simplified31.7
rmApplied sqr-pow31.7
Applied associate-*r*30.8
rmApplied sqrt-prod42.6
if 6.378441331140758e-28 < U* < 1.3559215670195436e+141Initial program 34.1
rmApplied *-un-lft-identity34.1
Applied times-frac31.2
Simplified31.2
rmApplied sqr-pow31.2
Applied associate-*r*30.0
Taylor expanded around 0 30.7
if 7.304292669069665e+231 < U* Initial program 38.5
rmApplied *-un-lft-identity38.5
Applied times-frac36.5
Simplified36.5
rmApplied sqr-pow36.5
Applied associate-*r*35.3
Final simplification33.2
herbie shell --seed 2020049
(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*))))))