\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 -1.8055396062121203 \cdot 10^{+120}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{\frac{t}{x \cdot x} \cdot \frac{1}{\sqrt{2}} - \left(t \cdot \sqrt{2} + 2 \cdot \left(\frac{t}{\sqrt{2} \cdot \left(x \cdot x\right)} + \frac{t}{\sqrt{2} \cdot x}\right)\right)}\\
\mathbf{elif}\;t \leq 7.744619413220778 \cdot 10^{+100}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{\sqrt{4 \cdot \frac{t \cdot t}{x} + 2 \cdot \left(t \cdot t + \ell \cdot \frac{\ell}{x}\right)}}\\
\mathbf{else}:\\
\;\;\;\;\frac{t \cdot \sqrt{2}}{\left(t \cdot \sqrt{2} + 2 \cdot \left(\frac{t}{\sqrt{2} \cdot \left(x \cdot x\right)} + \frac{t}{\sqrt{2} \cdot x}\right)\right) - \frac{t}{x \cdot x} \cdot \frac{1}{\sqrt{2}}}\\
\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 -1.8055396062121203e+120)
(/
(* t (sqrt 2.0))
(-
(* (/ t (* x x)) (/ 1.0 (sqrt 2.0)))
(+
(* t (sqrt 2.0))
(* 2.0 (+ (/ t (* (sqrt 2.0) (* x x))) (/ t (* (sqrt 2.0) x)))))))
(if (<= t 7.744619413220778e+100)
(/
(* t (sqrt 2.0))
(sqrt (+ (* 4.0 (/ (* t t) x)) (* 2.0 (+ (* t t) (* l (/ l x)))))))
(/
(* t (sqrt 2.0))
(-
(+
(* t (sqrt 2.0))
(* 2.0 (+ (/ t (* (sqrt 2.0) (* x x))) (/ t (* (sqrt 2.0) x)))))
(* (/ t (* x x)) (/ 1.0 (sqrt 2.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 <= -1.8055396062121203e+120) {
tmp = (t * sqrt(2.0)) / (((t / (x * x)) * (1.0 / sqrt(2.0))) - ((t * sqrt(2.0)) + (2.0 * ((t / (sqrt(2.0) * (x * x))) + (t / (sqrt(2.0) * x))))));
} else if (t <= 7.744619413220778e+100) {
tmp = (t * sqrt(2.0)) / sqrt((4.0 * ((t * t) / x)) + (2.0 * ((t * t) + (l * (l / x)))));
} else {
tmp = (t * sqrt(2.0)) / (((t * sqrt(2.0)) + (2.0 * ((t / (sqrt(2.0) * (x * x))) + (t / (sqrt(2.0) * x))))) - ((t / (x * x)) * (1.0 / sqrt(2.0))));
}
return tmp;
}



Bits error versus x



Bits error versus l



Bits error versus t
Results
if t < -1.80556e120Initial program 53.1
Taylor expanded around -inf 2.4
Simplified2.4
if -1.80556e120 < t < 7.74466e100Initial program 38.0
Taylor expanded around inf 18.0
Simplified18.0
rmApplied *-un-lft-identity_binary6418.0
Applied times-frac_binary6413.7
Simplified13.7
if 7.74466e100 < t Initial program 50.7
Taylor expanded around inf 2.8
Simplified2.8
Final simplification9.5
herbie shell --seed 2020231
(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)))))