\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 \leq -3.064774607274823 \cdot 10^{+89}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{-t \cdot \sqrt{\frac{2}{x - 1} + 2 \cdot \frac{x}{x - 1}}}\\
\mathbf{elif}\;t \leq -2.9474172478993817 \cdot 10^{-150}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{\sqrt{2 \cdot \frac{\ell \cdot \ell}{x} + \left(4 \cdot \left(\frac{t \cdot t}{x} + \frac{t \cdot t}{x \cdot x}\right) + 2 \cdot \left(t \cdot t + \frac{\ell}{\frac{x}{\frac{\ell}{x}}}\right)\right)}}\\
\mathbf{elif}\;t \leq -2.9041377634221667 \cdot 10^{-241}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{-t \cdot \sqrt{\frac{2}{x - 1} + 2 \cdot \frac{x}{x - 1}}}\\
\mathbf{elif}\;t \leq 2.1011173123428155 \cdot 10^{-259}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{2 \cdot \left(\frac{{t}^{2}}{\ell \cdot {x}^{2}} \cdot \sqrt{\frac{1}{2 \cdot \frac{1}{{x}^{2}} + 2 \cdot \frac{1}{x}}}\right) + \left(2 \cdot \left(\sqrt{\frac{1}{2 \cdot \frac{1}{{x}^{2}} + 2 \cdot \frac{1}{x}}} \cdot \frac{{t}^{2}}{x \cdot \ell}\right) + \left(\sqrt{\frac{1}{2 \cdot \frac{1}{{x}^{2}} + 2 \cdot \frac{1}{x}}} \cdot \frac{{t}^{2}}{\ell} + \ell \cdot \sqrt{2 \cdot \frac{1}{{x}^{2}} + 2 \cdot \frac{1}{x}}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{\sqrt{2}}{\sqrt{\frac{2}{x - 1} + 2 \cdot \frac{x}{x - 1}}}\\
\end{array}(FPCore (x l t) :precision binary64 (/ (* (sqrt 2.0) t) (sqrt (- (* (/ (+ x 1.0) (- x 1.0)) (+ (* l l) (* 2.0 (* t t)))) (* l l)))))
(FPCore (x l t)
:precision binary64
(if (<= t -3.064774607274823e+89)
(/
(* t (sqrt 2.0))
(- (* t (sqrt (+ (/ 2.0 (- x 1.0)) (* 2.0 (/ x (- x 1.0))))))))
(if (<= t -2.9474172478993817e-150)
(/
(* t (sqrt 2.0))
(sqrt
(+
(* 2.0 (/ (* l l) x))
(+
(* 4.0 (+ (/ (* t t) x) (/ (* t t) (* x x))))
(* 2.0 (+ (* t t) (/ l (/ x (/ l x)))))))))
(if (<= t -2.9041377634221667e-241)
(/
(* t (sqrt 2.0))
(- (* t (sqrt (+ (/ 2.0 (- x 1.0)) (* 2.0 (/ x (- x 1.0))))))))
(if (<= t 2.1011173123428155e-259)
(/
(* t (sqrt 2.0))
(+
(*
2.0
(*
(/ (pow t 2.0) (* l (pow x 2.0)))
(sqrt (/ 1.0 (+ (* 2.0 (/ 1.0 (pow x 2.0))) (* 2.0 (/ 1.0 x)))))))
(+
(*
2.0
(*
(sqrt (/ 1.0 (+ (* 2.0 (/ 1.0 (pow x 2.0))) (* 2.0 (/ 1.0 x)))))
(/ (pow t 2.0) (* x l))))
(+
(*
(sqrt (/ 1.0 (+ (* 2.0 (/ 1.0 (pow x 2.0))) (* 2.0 (/ 1.0 x)))))
(/ (pow t 2.0) l))
(* l (sqrt (+ (* 2.0 (/ 1.0 (pow x 2.0))) (* 2.0 (/ 1.0 x)))))))))
(/
(sqrt 2.0)
(sqrt (+ (/ 2.0 (- x 1.0)) (* 2.0 (/ x (- x 1.0)))))))))))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 tmp;
if (t <= -3.064774607274823e+89) {
tmp = (t * sqrt(2.0)) / -(t * sqrt((2.0 / (x - 1.0)) + (2.0 * (x / (x - 1.0)))));
} else if (t <= -2.9474172478993817e-150) {
tmp = (t * sqrt(2.0)) / sqrt((2.0 * ((l * l) / x)) + ((4.0 * (((t * t) / x) + ((t * t) / (x * x)))) + (2.0 * ((t * t) + (l / (x / (l / x)))))));
} else if (t <= -2.9041377634221667e-241) {
tmp = (t * sqrt(2.0)) / -(t * sqrt((2.0 / (x - 1.0)) + (2.0 * (x / (x - 1.0)))));
} else if (t <= 2.1011173123428155e-259) {
tmp = (t * sqrt(2.0)) / ((2.0 * ((pow(t, 2.0) / (l * pow(x, 2.0))) * sqrt(1.0 / ((2.0 * (1.0 / pow(x, 2.0))) + (2.0 * (1.0 / x)))))) + ((2.0 * (sqrt(1.0 / ((2.0 * (1.0 / pow(x, 2.0))) + (2.0 * (1.0 / x)))) * (pow(t, 2.0) / (x * l)))) + ((sqrt(1.0 / ((2.0 * (1.0 / pow(x, 2.0))) + (2.0 * (1.0 / x)))) * (pow(t, 2.0) / l)) + (l * sqrt((2.0 * (1.0 / pow(x, 2.0))) + (2.0 * (1.0 / x)))))));
} else {
tmp = sqrt(2.0) / sqrt((2.0 / (x - 1.0)) + (2.0 * (x / (x - 1.0))));
}
return tmp;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -3.064774607274823e89 or -2.9474172478993817e-150 < t < -2.9041377634221667e-241Initial program 52.0
Taylor expanded around -inf 9.5
Simplified9.5
if -3.064774607274823e89 < t < -2.9474172478993817e-150Initial program 26.4
Taylor expanded around inf 10.3
Simplified10.3
rmApplied associate-/l*_binary649.9
Simplified9.9
if -2.9041377634221667e-241 < t < 2.10111731234281555e-259Initial program 62.9
Taylor expanded around inf 36.1
Simplified36.1
Taylor expanded around inf 32.6
if 2.10111731234281555e-259 < t Initial program 41.6
Taylor expanded around inf 12.4
Simplified12.4
rmApplied associate-/l*_binary6412.4
Simplified12.3
Final simplification12.4
herbie shell --seed 2021173
(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)))))