x + \left(y - x\right) \cdot \frac{z}{t}\begin{array}{l}
\mathbf{if}\;\frac{z}{t} \le 3.2897350745799706 \cdot 10^{215}:\\
\;\;\;\;\mathsf{fma}\left(y - x, \frac{z}{t}, x\right)\\
\mathbf{else}:\\
\;\;\;\;x + \frac{\left(y - x\right) \cdot z}{t}\\
\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 (((z / t) <= 3.2897350745799706e+215)) {
temp = fma((y - x), (z / t), x);
} else {
temp = (x + (((y - x) * z) / t));
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t
Results
| Original | 2.2 |
|---|---|
| Target | 2.5 |
| Herbie | 1.5 |
if (/ z t) < 3.2897350745799706e+215Initial program 1.5
Simplified1.5
if 3.2897350745799706e+215 < (/ z t) Initial program 25.9
rmApplied associate-*r/0.5
Final simplification1.5
herbie shell --seed 2020060 +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))))