\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 -3.794989872061021 \cdot 10^{143}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\frac{\frac{t}{x}}{\sqrt{2}} \cdot \left(\frac{1}{x} - \left(2 + \frac{2}{x}\right)\right) - t \cdot \sqrt{2}}\\
\mathbf{elif}\;t \le 3.3289760264210936 \cdot 10^{110}:\\
\;\;\;\;\frac{\left(\sqrt[3]{\sqrt{2}} \cdot \sqrt[3]{\sqrt{2}}\right) \cdot \left(\sqrt[3]{\sqrt{2}} \cdot t\right)}{\sqrt{{t}^{2} \cdot \left(2 + \frac{4}{x}\right) + 2 \cdot \frac{\ell}{\frac{x}{\ell}}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\frac{t}{{x}^{2} \cdot \sqrt{2}} \cdot \left(2 - 1\right) + t \cdot \left(\sqrt{2} + \frac{2}{x \cdot \sqrt{2}}\right)}\\
\end{array}double code(double x, double l, double t) {
return (((double) (((double) sqrt(2.0)) * t)) / ((double) sqrt(((double) (((double) ((((double) (x + 1.0)) / ((double) (x - 1.0))) * ((double) (((double) (l * l)) + ((double) (2.0 * ((double) (t * t)))))))) - ((double) (l * l)))))));
}
double code(double x, double l, double t) {
double VAR;
if ((t <= -3.794989872061021e+143)) {
VAR = (((double) (((double) sqrt(2.0)) * t)) / ((double) (((double) (((t / x) / ((double) sqrt(2.0))) * ((double) ((1.0 / x) - ((double) (2.0 + (2.0 / x))))))) - ((double) (t * ((double) sqrt(2.0)))))));
} else {
double VAR_1;
if ((t <= 3.3289760264210936e+110)) {
VAR_1 = (((double) (((double) (((double) cbrt(((double) sqrt(2.0)))) * ((double) cbrt(((double) sqrt(2.0)))))) * ((double) (((double) cbrt(((double) sqrt(2.0)))) * t)))) / ((double) sqrt(((double) (((double) (((double) pow(t, 2.0)) * ((double) (2.0 + (4.0 / x))))) + ((double) (2.0 * (l / (x / l)))))))));
} else {
VAR_1 = (((double) (((double) sqrt(2.0)) * t)) / ((double) (((double) ((t / ((double) (((double) pow(x, 2.0)) * ((double) sqrt(2.0))))) * ((double) (2.0 - 1.0)))) + ((double) (t * ((double) (((double) sqrt(2.0)) + (2.0 / ((double) (x * ((double) sqrt(2.0))))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -3.794989872061021e143Initial program 59.2
Taylor expanded around -inf 2.0
Simplified2.0
if -3.794989872061021e143 < t < 3.3289760264210936e110Initial program 35.8
Taylor expanded around inf 16.7
Simplified16.7
rmApplied unpow216.7
Applied associate-/l*12.1
rmApplied add-cube-cbrt12.1
Applied associate-*l*12.1
if 3.3289760264210936e110 < t Initial program 52.2
Taylor expanded around inf 2.9
Simplified2.9
Final simplification8.6
herbie shell --seed 2020182
(FPCore (x l t)
:name "Toniolo and Linder, Equation (7)"
:precision binary64
(/ (* (sqrt 2.0) t) (sqrt (- (* (/ (+ x 1.0) (- x 1.0)) (+ (* l l) (* 2.0 (* t t)))) (* l l)))))