\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 -5.19934999049700251 \cdot 10^{36}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\left(2 \cdot \left(\frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}} - \frac{t}{\sqrt{2} \cdot {x}^{2}}\right) - \sqrt{2} \cdot t\right) - 2 \cdot \frac{t}{\sqrt{2} \cdot x}}\\
\mathbf{elif}\;t \le 2.85675210770604991 \cdot 10^{-244}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \ell \cdot \frac{\ell}{x}\right)}}\\
\mathbf{elif}\;t \le 3.03005972100276706 \cdot 10^{-173}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{2 \cdot \left(\frac{t}{\sqrt{2} \cdot {x}^{2}} + \frac{t}{\sqrt{2} \cdot x}\right) + \left(\sqrt{2} \cdot t - 2 \cdot \frac{t}{{\left(\sqrt{2}\right)}^{3} \cdot {x}^{2}}\right)}\\
\mathbf{elif}\;t \le 2.2589321978357942 \cdot 10^{91}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{\sqrt{4 \cdot \frac{{t}^{2}}{x} + 2 \cdot \left({t}^{2} + \frac{\ell}{\frac{x}{\ell}}\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2} \cdot t}{2 \cdot \left(\frac{t}{\sqrt{2} \cdot {x}^{2}} + \frac{t}{\sqrt{2} \cdot x}\right) + \left(\sqrt{2} \cdot t - 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 VAR;
if ((t <= -5.1993499904970025e+36)) {
VAR = ((sqrt(2.0) * t) / (((2.0 * ((t / (pow(sqrt(2.0), 3.0) * pow(x, 2.0))) - (t / (sqrt(2.0) * pow(x, 2.0))))) - (sqrt(2.0) * t)) - (2.0 * (t / (sqrt(2.0) * x)))));
} else {
double VAR_1;
if ((t <= 2.85675210770605e-244)) {
VAR_1 = ((sqrt(2.0) * t) / sqrt(((4.0 * (pow(t, 2.0) / x)) + (2.0 * (pow(t, 2.0) + (l * (l / x)))))));
} else {
double VAR_2;
if ((t <= 3.030059721002767e-173)) {
VAR_2 = ((sqrt(2.0) * t) / ((2.0 * ((t / (sqrt(2.0) * pow(x, 2.0))) + (t / (sqrt(2.0) * x)))) + ((sqrt(2.0) * t) - (2.0 * (t / (pow(sqrt(2.0), 3.0) * pow(x, 2.0)))))));
} else {
double VAR_3;
if ((t <= 2.258932197835794e+91)) {
VAR_3 = ((sqrt(2.0) * t) / sqrt(((4.0 * (pow(t, 2.0) / x)) + (2.0 * (pow(t, 2.0) + (l / (x / l)))))));
} else {
VAR_3 = ((sqrt(2.0) * t) / ((2.0 * ((t / (sqrt(2.0) * pow(x, 2.0))) + (t / (sqrt(2.0) * x)))) + ((sqrt(2.0) * t) - (2.0 * (t / (pow(sqrt(2.0), 3.0) * pow(x, 2.0)))))));
}
VAR_2 = VAR_3;
}
VAR_1 = VAR_2;
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -5.1993499904970025e+36Initial program 43.1
Taylor expanded around -inf 4.7
Simplified4.7
if -5.1993499904970025e+36 < t < 2.85675210770605e-244Initial program 45.1
Taylor expanded around inf 20.0
Simplified20.0
rmApplied *-un-lft-identity20.0
Applied add-sqr-sqrt43.3
Applied unpow-prod-down43.3
Applied times-frac41.7
Simplified41.7
Simplified16.6
if 2.85675210770605e-244 < t < 3.030059721002767e-173 or 2.258932197835794e+91 < t Initial program 52.7
Taylor expanded around inf 8.4
Simplified8.4
if 3.030059721002767e-173 < t < 2.258932197835794e+91Initial program 28.3
Taylor expanded around inf 11.3
Simplified11.3
rmApplied unpow211.3
Applied associate-/l*6.5
Final simplification9.4
herbie shell --seed 2020106
(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)))))