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 \lor \neg \left(x + \frac{\left(y - x\right) \cdot z}{t} \leq 2.3361936796503275 \cdot 10^{+277}\right):\\
\;\;\;\;x + \left(y - x\right) \cdot \frac{z}{t}\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\left(y - x\right) \cdot z}{t}\\
\end{array}(FPCore (x y z t) :precision binary64 (+ x (/ (* (- y x) z) t)))
(FPCore (x y z t)
:precision binary64
(if (or (<= (+ x (/ (* (- y x) z) t)) (- INFINITY))
(not (<= (+ x (/ (* (- y x) z) t)) 2.3361936796503275e+277)))
(+ x (* (- y x) (/ z t)))
(+ x (/ (* (- y x) z) t))))double code(double x, double y, double z, double t) {
return ((double) (x + (((double) (((double) (y - x)) * z)) / t)));
}
double code(double x, double y, double z, double t) {
double VAR;
if (((((double) (x + (((double) (((double) (y - x)) * z)) / t))) <= ((double) -(((double) INFINITY)))) || !(((double) (x + (((double) (((double) (y - x)) * z)) / t))) <= 2.3361936796503275e+277))) {
VAR = ((double) (x + ((double) (((double) (y - x)) * (z / t)))));
} else {
VAR = ((double) (x + (((double) (((double) (y - x)) * z)) / t)));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 6.6 |
|---|---|
| Target | 1.9 |
| Herbie | 1.1 |
if (+ x (/ (* (- y x) z) t)) < -inf.0 or 2.33619367965032753e277 < (+ x (/ (* (- y x) z) t)) Initial program 49.4
Simplified2.0
if -inf.0 < (+ x (/ (* (- y x) z) t)) < 2.33619367965032753e277Initial program 0.9
Final simplification1.1
herbie shell --seed 2020198
(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)))