x + \frac{y \cdot \left(z - t\right)}{z - a}\begin{array}{l}
\mathbf{if}\;\frac{y \cdot \left(z - t\right)}{z - a} = -\infty \lor \neg \left(\frac{y \cdot \left(z - t\right)}{z - a} \le 2.7824488815679915 \cdot 10^{176}\right):\\
\;\;\;\;x + y \cdot \frac{z - t}{z - a}\\
\mathbf{else}:\\
\;\;\;\;x + \left(y \cdot \left(z - t\right)\right) \cdot \frac{1}{z - a}\\
\end{array}double code(double x, double y, double z, double t, double a) {
return (x + ((y * (z - t)) / (z - a)));
}
double code(double x, double y, double z, double t, double a) {
double temp;
if (((((y * (z - t)) / (z - a)) <= -inf.0) || !(((y * (z - t)) / (z - a)) <= 2.7824488815679915e+176))) {
temp = (x + (y * ((z - t) / (z - a))));
} else {
temp = (x + ((y * (z - t)) * (1.0 / (z - a))));
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z




Bits error versus t




Bits error versus a
Results
| Original | 11.3 |
|---|---|
| Target | 1.4 |
| Herbie | 0.6 |
if (/ (* y (- z t)) (- z a)) < -inf.0 or 2.7824488815679915e+176 < (/ (* y (- z t)) (- z a)) Initial program 53.4
rmApplied *-un-lft-identity53.4
Applied times-frac1.7
Simplified1.7
if -inf.0 < (/ (* y (- z t)) (- z a)) < 2.7824488815679915e+176Initial program 0.3
rmApplied div-inv0.3
Final simplification0.6
herbie shell --seed 2020060 +o rules:numerics
(FPCore (x y z t a)
:name "Graphics.Rendering.Plot.Render.Plot.Axis:renderAxisTicks from plot-0.2.3.4, A"
:precision binary64
:herbie-target
(+ x (/ y (/ (- z a) (- z t))))
(+ x (/ (* y (- z t)) (- z a))))