x + \left(y - x\right) \cdot \frac{z}{t}\begin{array}{l}
\mathbf{if}\;x \le -8.64219457899372803 \cdot 10^{-224} \lor \neg \left(x \le 1.6257691218099716 \cdot 10^{-308}\right):\\
\;\;\;\;\mathsf{fma}\left(y - x, \frac{z}{t}, x\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(y - x\right) \cdot z}{t} + x\\
\end{array}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 temp;
if (((x <= -8.642194578993728e-224) || !(x <= 1.6257691218099716e-308))) {
temp = fma((y - x), (z / t), x);
} else {
temp = ((((y - x) * z) / t) + x);
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.0 |
|---|---|
| Target | 2.2 |
| Herbie | 2.0 |
if x < -8.642194578993728e-224 or 1.6257691218099716e-308 < x Initial program 1.7
Simplified1.7
if -8.642194578993728e-224 < x < 1.6257691218099716e-308Initial program 6.0
Simplified6.0
rmApplied fma-udef6.0
rmApplied associate-*r/5.5
Final simplification2.0
herbie shell --seed 2020056 +o rules:numerics
(FPCore (x y z t)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:tickPosition from plot-0.2.3.4"
:precision binary64
:herbie-target
(if (< (* (- y x) (/ z t)) -1013646692435.887) (+ x (/ (- y x) (/ t z))) (if (< (* (- y x) (/ z t)) -0.0) (+ x (/ (* (- y x) z) t)) (+ x (/ (- y x) (/ t z)))))
(+ x (* (- y x) (/ z t))))