\frac{\sqrt{2} \cdot t}{\sqrt{\frac{x + 1}{x - 1} \cdot \left(\ell \cdot \ell + 2 \cdot \left(t \cdot t\right)\right) - \ell \cdot \ell}}\begin{array}{l}
\mathbf{if}\;t \le -8.08541959572897258 \cdot 10^{99}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\mathsf{fma}\left(2, \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}}, -\mathsf{fma}\left(2, \frac{t}{\sqrt{2} \cdot {x}^{2}}, \mathsf{fma}\left(2, \frac{t}{\sqrt{2} \cdot x}, t \cdot \sqrt{2}\right)\right)\right)}\\
\mathbf{elif}\;t \le 5.31686358457741622 \cdot 10^{-279}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{\mathsf{fma}\left(2, {t}^{2}, \mathsf{fma}\left(2, \left(\sqrt[3]{{\left(\sqrt[3]{\ell}\right)}^{4}} \cdot \sqrt[3]{{\left(\sqrt[3]{\ell}\right)}^{4}}\right) \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\ell}\right)}^{4}} \cdot \frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{x}\right), 4 \cdot \frac{{t}^{2}}{x}\right)\right)}}\\
\mathbf{elif}\;t \le 8.3630275522259312 \cdot 10^{-198}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\mathsf{fma}\left(2, \frac{t}{\sqrt{2} \cdot {x}^{2}}, \mathsf{fma}\left(2, \frac{t}{\sqrt{2} \cdot x}, t \cdot \sqrt{2}\right) - 2 \cdot \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}}\right)}\\
\mathbf{elif}\;t \le 1.32955370740893664 \cdot 10^{133}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{\mathsf{fma}\left(2, {t}^{2}, \mathsf{fma}\left(2, \left(\sqrt[3]{{\left(\sqrt[3]{\ell}\right)}^{4}} \cdot \sqrt[3]{{\left(\sqrt[3]{\ell}\right)}^{4}}\right) \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\ell}\right)}^{4}} \cdot \frac{{\left(\sqrt[3]{\ell}\right)}^{2}}{x}\right), 4 \cdot \frac{{t}^{2}}{x}\right)\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\mathsf{fma}\left(2, \frac{t}{\sqrt{2} \cdot {x}^{2}}, \mathsf{fma}\left(2, \frac{t}{\sqrt{2} \cdot x}, t \cdot \sqrt{2}\right) - 2 \cdot \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}}\right)}\\
\end{array}double code(double x, double l, double t) {
return ((sqrt(2.0) * t) / sqrt(((((x + 1.0) / (x - 1.0)) * ((l * l) + (2.0 * (t * t)))) - (l * l))));
}
double code(double x, double l, double t) {
double temp;
if ((t <= -8.085419595728973e+99)) {
temp = ((sqrt(2.0) * t) / fma(2.0, (t / (pow(sqrt(2.0), 3.0) * pow(x, 2.0))), -fma(2.0, (t / (sqrt(2.0) * pow(x, 2.0))), fma(2.0, (t / (sqrt(2.0) * x)), (t * sqrt(2.0))))));
} else {
double temp_1;
if ((t <= 5.316863584577416e-279)) {
temp_1 = ((sqrt(2.0) * t) / sqrt(fma(2.0, pow(t, 2.0), fma(2.0, ((cbrt(pow(cbrt(l), 4.0)) * cbrt(pow(cbrt(l), 4.0))) * (cbrt(pow(cbrt(l), 4.0)) * (pow(cbrt(l), 2.0) / x))), (4.0 * (pow(t, 2.0) / x))))));
} else {
double temp_2;
if ((t <= 8.363027552225931e-198)) {
temp_2 = ((sqrt(2.0) * t) / fma(2.0, (t / (sqrt(2.0) * pow(x, 2.0))), (fma(2.0, (t / (sqrt(2.0) * x)), (t * sqrt(2.0))) - (2.0 * (t / (pow(sqrt(2.0), 3.0) * pow(x, 2.0)))))));
} else {
double temp_3;
if ((t <= 1.3295537074089366e+133)) {
temp_3 = ((sqrt(2.0) * t) / sqrt(fma(2.0, pow(t, 2.0), fma(2.0, ((cbrt(pow(cbrt(l), 4.0)) * cbrt(pow(cbrt(l), 4.0))) * (cbrt(pow(cbrt(l), 4.0)) * (pow(cbrt(l), 2.0) / x))), (4.0 * (pow(t, 2.0) / x))))));
} else {
temp_3 = ((sqrt(2.0) * t) / fma(2.0, (t / (sqrt(2.0) * pow(x, 2.0))), (fma(2.0, (t / (sqrt(2.0) * x)), (t * sqrt(2.0))) - (2.0 * (t / (pow(sqrt(2.0), 3.0) * pow(x, 2.0)))))));
}
temp_2 = temp_3;
}
temp_1 = temp_2;
}
temp = temp_1;
}
return temp;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -8.085419595728973e+99Initial program 50.0
Taylor expanded around -inf 3.3
Simplified3.3
if -8.085419595728973e+99 < t < 5.316863584577416e-279 or 8.363027552225931e-198 < t < 1.3295537074089366e+133Initial program 34.1
Taylor expanded around inf 15.2
Simplified15.2
rmApplied *-un-lft-identity15.2
Applied add-cube-cbrt15.3
Applied unpow-prod-down15.3
Applied times-frac12.0
Simplified12.0
rmApplied add-cube-cbrt12.0
Applied associate-*l*12.0
if 5.316863584577416e-279 < t < 8.363027552225931e-198 or 1.3295537074089366e+133 < t Initial program 58.2
Taylor expanded around inf 9.7
Simplified9.7
Final simplification9.8
herbie shell --seed 2020058 +o rules:numerics
(FPCore (x l t)
:name "Toniolo and Linder, Equation (7)"
:precision binary64
(/ (* (sqrt 2) t) (sqrt (- (* (/ (+ x 1) (- x 1)) (+ (* l l) (* 2 (* t t)))) (* l l)))))