x + \frac{\left(y - x\right) \cdot z}{t}\begin{array}{l}
\mathbf{if}\;x + \frac{\left(y - x\right) \cdot z}{t} \leq -\infty:\\
\;\;\;\;x + \frac{y - x}{\frac{t}{z}}\\
\mathbf{elif}\;x + \frac{\left(y - x\right) \cdot z}{t} \leq 3.2613229998891658 \cdot 10^{+299}:\\
\;\;\;\;x + \frac{\left(y - x\right) \cdot z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + z \cdot \frac{y - x}{t}\\
\end{array}(FPCore (x y z t) :precision binary64 (+ x (/ (* (- y x) z) t)))
(FPCore (x y z t)
:precision binary64
(if (<= (+ x (/ (* (- y x) z) t)) (- INFINITY))
(+ x (/ (- y x) (/ t z)))
(if (<= (+ x (/ (* (- y x) z) t)) 3.2613229998891658e+299)
(+ x (/ (* (- y x) z) t))
(+ x (* z (/ (- y x) t))))))double code(double x, double y, double z, double t) {
return x + (((y - x) * z) / t);
}
double code(double x, double y, double z, double t) {
double tmp;
if ((x + (((y - x) * z) / t)) <= -((double) INFINITY)) {
tmp = x + ((y - x) / (t / z));
} else if ((x + (((y - x) * z) / t)) <= 3.2613229998891658e+299) {
tmp = x + (((y - x) * z) / t);
} else {
tmp = x + (z * ((y - x) / t));
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.8 |
|---|---|
| Target | 2.1 |
| Herbie | 0.9 |
if (+.f64 x (/.f64 (*.f64 (-.f64 y x) z) t)) < -inf.0Initial program 64.0
rmApplied associate-/l*_binary640.2
if -inf.0 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) z) t)) < 3.26132299988916578e299Initial program 0.7
if 3.26132299988916578e299 < (+.f64 x (/.f64 (*.f64 (-.f64 y x) z) t)) Initial program 54.3
rmApplied associate-/l*_binary641.5
rmApplied associate-/r/_binary643.9
Final simplification0.9
herbie shell --seed 2020232
(FPCore (x y z t)
:name "Numeric.Histogram:binBounds from Chart-1.5.3"
:precision binary64
:herbie-target
(if (< x -9.025511195533005e-135) (- x (* (/ z t) (- x y))) (if (< x 4.275032163700715e-250) (+ x (* (/ (- y x) t) z)) (+ x (/ (- y x) (/ t z)))))
(+ x (/ (* (- y x) z) t)))