\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}
t_1 := \sqrt{n \cdot \left(U \cdot \left(\frac{n \cdot U*}{Om \cdot Om} - \left(\frac{2}{Om} + \frac{n \cdot U}{Om \cdot Om}\right)\right)\right)} \cdot \left(\ell \cdot \sqrt{2}\right)\\
\mathbf{if}\;\ell \leq -1.0774931511781235 \cdot 10^{+189}:\\
\;\;\;\;-t_1\\
\mathbf{elif}\;\ell \leq -8.228016994005398 \cdot 10^{-40}:\\
\;\;\;\;{\left(\left(n \cdot 2\right) \cdot \left(U \cdot \mathsf{fma}\left(\frac{\ell}{Om}, \mathsf{fma}\left(\ell, -2, \left(U* - U\right) \cdot \left(n \cdot \frac{\ell}{Om}\right)\right), t\right)\right)\right)}^{0.5}\\
\mathbf{elif}\;\ell \leq 7.2207560999646215 \cdot 10^{-224}:\\
\;\;\;\;\sqrt{\left(U \cdot \left(n \cdot 2\right)\right) \cdot \left(t + \frac{\ell}{Om} \cdot \mathsf{fma}\left(\ell, -2, \frac{n \cdot \left(\ell \cdot U* - \ell \cdot U\right)}{Om}\right)\right)}\\
\mathbf{elif}\;\ell \leq 1.3083543023179113 \cdot 10^{+195}:\\
\;\;\;\;\sqrt{\left(n \cdot 2\right) \cdot \left(U \cdot \left(t + \frac{\ell}{Om} \cdot \mathsf{fma}\left(\ell, -2, \left(U* - U\right) \cdot \left(\left(n \cdot \sqrt{\ell}\right) \cdot \frac{\sqrt{\ell}}{Om}\right)\right)\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;t_1\\
\end{array}
(FPCore (n U t l Om U*) :precision binary64 (sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))
(FPCore (n U t l Om U*)
:precision binary64
(let* ((t_1
(*
(sqrt
(*
n
(*
U
(- (/ (* n U*) (* Om Om)) (+ (/ 2.0 Om) (/ (* n U) (* Om Om)))))))
(* l (sqrt 2.0)))))
(if (<= l -1.0774931511781235e+189)
(- t_1)
(if (<= l -8.228016994005398e-40)
(pow
(*
(* n 2.0)
(* U (fma (/ l Om) (fma l -2.0 (* (- U* U) (* n (/ l Om)))) t)))
0.5)
(if (<= l 7.2207560999646215e-224)
(sqrt
(*
(* U (* n 2.0))
(+ t (* (/ l Om) (fma l -2.0 (/ (* n (- (* l U*) (* l U))) Om))))))
(if (<= l 1.3083543023179113e+195)
(sqrt
(*
(* n 2.0)
(*
U
(+
t
(*
(/ l Om)
(fma
l
-2.0
(* (- U* U) (* (* n (sqrt l)) (/ (sqrt l) Om)))))))))
t_1))))))double code(double n, double U, double t, double l, double Om, double U_42_) {
return sqrt(((2.0 * n) * U) * ((t - (2.0 * ((l * l) / Om))) - ((n * pow((l / Om), 2.0)) * (U - U_42_))));
}
double code(double n, double U, double t, double l, double Om, double U_42_) {
double t_1 = sqrt(n * (U * (((n * U_42_) / (Om * Om)) - ((2.0 / Om) + ((n * U) / (Om * Om)))))) * (l * sqrt(2.0));
double tmp;
if (l <= -1.0774931511781235e+189) {
tmp = -t_1;
} else if (l <= -8.228016994005398e-40) {
tmp = pow(((n * 2.0) * (U * fma((l / Om), fma(l, -2.0, ((U_42_ - U) * (n * (l / Om)))), t))), 0.5);
} else if (l <= 7.2207560999646215e-224) {
tmp = sqrt((U * (n * 2.0)) * (t + ((l / Om) * fma(l, -2.0, ((n * ((l * U_42_) - (l * U))) / Om)))));
} else if (l <= 1.3083543023179113e+195) {
tmp = sqrt((n * 2.0) * (U * (t + ((l / Om) * fma(l, -2.0, ((U_42_ - U) * ((n * sqrt(l)) * (sqrt(l) / Om))))))));
} else {
tmp = t_1;
}
return tmp;
}



Bits error versus n



Bits error versus U



Bits error versus t



Bits error versus l



Bits error versus Om



Bits error versus U*
if l < -1.0774931511781235e189Initial program 64.0
Simplified52.7
Taylor expanded in l around -inf 32.7
Simplified32.7
if -1.0774931511781235e189 < l < -8.22801699400539806e-40Initial program 37.3
Simplified31.7
Applied associate-*l*_binary6430.7
Applied add-cube-cbrt_binary6431.1
Simplified31.1
Simplified31.1
Applied pow1/3_binary6432.6
Applied pow1/3_binary6433.0
Applied pow1/3_binary6433.3
Applied pow-prod-up_binary6433.3
Applied pow-prod-up_binary6430.7
Applied sqrt-pow1_binary6430.7
if -8.22801699400539806e-40 < l < 7.2207560999646215e-224Initial program 24.7
Simplified23.9
Taylor expanded in Om around -inf 24.4
Simplified24.4
if 7.2207560999646215e-224 < l < 1.3083543023179113e195Initial program 32.0
Simplified28.4
Applied associate-*l*_binary6428.1
Applied *-un-lft-identity_binary6428.1
Applied add-sqr-sqrt_binary6428.1
Applied times-frac_binary6428.1
Applied associate-*r*_binary6428.4
Simplified28.4
if 1.3083543023179113e195 < l Initial program 64.0
Simplified52.6
Taylor expanded in l around inf 32.6
Simplified32.6
Final simplification28.0
herbie shell --seed 2022081
(FPCore (n U t l Om U*)
:name "Toniolo and Linder, Equation (13)"
:precision binary64
(sqrt (* (* (* 2.0 n) U) (- (- t (* 2.0 (/ (* l l) Om))) (* (* n (pow (/ l Om) 2.0)) (- U U*))))))