\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}\;n \le -4.579846573878255 \cdot 10^{170}:\\
\;\;\;\;\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}\;n \le -9.1730180109330935 \cdot 10^{-133}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-\left(U - U*\right)\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right)\right)\right)\right)}\\
\mathbf{elif}\;n \le -3.37883772685991123 \cdot 10^{-162}:\\
\;\;\;\;\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 {\left(\sqrt[3]{\frac{\ell}{Om}} \cdot \sqrt[3]{\frac{\ell}{Om}}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot {\left(\sqrt[3]{\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}\;n \le 1.8079991333180533 \cdot 10^{-150}:\\
\;\;\;\;\sqrt{{\left(\left(\left(2 \cdot n\right) \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-\left(U - U*\right)\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right)\right)\right)\right) \cdot U\right)}^{1}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\left(2 \cdot n\right) \cdot \left(U \cdot \left(t - \left(2 \cdot \left(\ell \cdot \frac{\ell}{Om}\right) - \left(-\left(U - U*\right)\right) \cdot \left(n \cdot {\left(\frac{\ell}{Om}\right)}^{\left(2 \cdot \frac{2}{2}\right)}\right)\right)\right)\right)}\\
\end{array}double f(double n, double U, double t, double l, double Om, double U_) {
double r252469 = 2.0;
double r252470 = n;
double r252471 = r252469 * r252470;
double r252472 = U;
double r252473 = r252471 * r252472;
double r252474 = t;
double r252475 = l;
double r252476 = r252475 * r252475;
double r252477 = Om;
double r252478 = r252476 / r252477;
double r252479 = r252469 * r252478;
double r252480 = r252474 - r252479;
double r252481 = r252475 / r252477;
double r252482 = pow(r252481, r252469);
double r252483 = r252470 * r252482;
double r252484 = U_;
double r252485 = r252472 - r252484;
double r252486 = r252483 * r252485;
double r252487 = r252480 - r252486;
double r252488 = r252473 * r252487;
double r252489 = sqrt(r252488);
return r252489;
}
double f(double n, double U, double t, double l, double Om, double U_) {
double r252490 = n;
double r252491 = -4.5798465738782547e+170;
bool r252492 = r252490 <= r252491;
double r252493 = 2.0;
double r252494 = r252493 * r252490;
double r252495 = U;
double r252496 = r252494 * r252495;
double r252497 = sqrt(r252496);
double r252498 = t;
double r252499 = l;
double r252500 = Om;
double r252501 = r252499 / r252500;
double r252502 = r252499 * r252501;
double r252503 = r252493 * r252502;
double r252504 = r252498 - r252503;
double r252505 = 2.0;
double r252506 = r252493 / r252505;
double r252507 = pow(r252501, r252506);
double r252508 = r252490 * r252507;
double r252509 = r252508 * r252507;
double r252510 = U_;
double r252511 = r252495 - r252510;
double r252512 = r252509 * r252511;
double r252513 = r252504 - r252512;
double r252514 = sqrt(r252513);
double r252515 = r252497 * r252514;
double r252516 = -9.173018010933093e-133;
bool r252517 = r252490 <= r252516;
double r252518 = -r252511;
double r252519 = r252505 * r252506;
double r252520 = pow(r252501, r252519);
double r252521 = r252490 * r252520;
double r252522 = r252518 * r252521;
double r252523 = r252503 - r252522;
double r252524 = r252498 - r252523;
double r252525 = r252495 * r252524;
double r252526 = r252494 * r252525;
double r252527 = sqrt(r252526);
double r252528 = -3.378837726859911e-162;
bool r252529 = r252490 <= r252528;
double r252530 = cbrt(r252501);
double r252531 = r252530 * r252530;
double r252532 = pow(r252531, r252506);
double r252533 = r252490 * r252532;
double r252534 = pow(r252530, r252506);
double r252535 = r252533 * r252534;
double r252536 = r252535 * r252507;
double r252537 = r252536 * r252511;
double r252538 = r252504 - r252537;
double r252539 = r252496 * r252538;
double r252540 = sqrt(r252539);
double r252541 = 1.8079991333180533e-150;
bool r252542 = r252490 <= r252541;
double r252543 = r252494 * r252524;
double r252544 = r252543 * r252495;
double r252545 = 1.0;
double r252546 = pow(r252544, r252545);
double r252547 = sqrt(r252546);
double r252548 = r252542 ? r252547 : r252527;
double r252549 = r252529 ? r252540 : r252548;
double r252550 = r252517 ? r252527 : r252549;
double r252551 = r252492 ? r252515 : r252550;
return r252551;
}



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 n < -4.5798465738782547e+170Initial program 39.0
rmApplied *-un-lft-identity39.0
Applied times-frac38.0
Simplified38.0
rmApplied sqr-pow38.0
Applied associate-*r*35.8
rmApplied sqrt-prod46.0
if -4.5798465738782547e+170 < n < -9.173018010933093e-133 or 1.8079991333180533e-150 < n Initial program 31.2
rmApplied *-un-lft-identity31.2
Applied times-frac28.3
Simplified28.3
rmApplied sqr-pow28.3
Applied associate-*r*27.4
rmApplied associate-*l*27.7
Simplified28.7
if -9.173018010933093e-133 < n < -3.378837726859911e-162Initial program 33.4
rmApplied *-un-lft-identity33.4
Applied times-frac30.0
Simplified30.0
rmApplied sqr-pow30.0
Applied associate-*r*29.9
rmApplied add-cube-cbrt30.0
Applied unpow-prod-down30.0
Applied associate-*r*30.0
if -3.378837726859911e-162 < n < 1.8079991333180533e-150Initial program 39.7
rmApplied *-un-lft-identity39.7
Applied times-frac37.1
Simplified37.1
rmApplied sqr-pow37.1
Applied associate-*r*36.4
rmApplied pow136.4
Applied pow136.4
Applied pow136.4
Applied pow136.4
Applied pow-prod-down36.4
Applied pow-prod-down36.4
Applied pow-prod-down36.4
Simplified31.2
Final simplification30.7
herbie shell --seed 2020056 +o rules:numerics
(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*))))))