\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\begin{array}{l}
\mathbf{if}\;z \le -2.1306842356386806 \cdot 10^{-13} \lor \neg \left(z \le 5.31693115780729608 \cdot 10^{-29}\right):\\
\;\;\;\;x \cdot \frac{\left(y - z\right) + 1}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{z} \cdot \left(\left(y - z\right) + 1\right)\\
\end{array}double code(double x, double y, double z) {
return ((x * ((y - z) + 1.0)) / z);
}
double code(double x, double y, double z) {
double VAR;
if (((z <= -2.1306842356386806e-13) || !(z <= 5.316931157807296e-29))) {
VAR = (x * (((y - z) + 1.0) / z));
} else {
VAR = ((x / z) * ((y - z) + 1.0));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 9.9 |
|---|---|
| Target | 0.5 |
| Herbie | 0.1 |
if z < -2.1306842356386806e-13 or 5.316931157807296e-29 < z Initial program 15.2
rmApplied *-un-lft-identity15.2
Applied times-frac0.1
Simplified0.1
if -2.1306842356386806e-13 < z < 5.316931157807296e-29Initial program 0.1
rmApplied associate-/l*8.1
rmApplied associate-/r/0.1
Final simplification0.1
herbie shell --seed 2020102 +o rules:numerics
(FPCore (x y z)
:name "Diagrams.TwoD.Segment.Bernstein:evaluateBernstein from diagrams-lib-1.3.0.3"
:precision binary64
:herbie-target
(if (< x -2.71483106713436e-162) (- (* (+ 1 y) (/ x z)) x) (if (< x 3.874108816439546e-197) (* (* x (+ (- y z) 1)) (/ 1 z)) (- (* (+ 1 y) (/ x z)) x)))
(/ (* x (+ (- y z) 1)) z))