\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\left(1 + y\right) \cdot \frac{x}{z} - xdouble f(double x, double y, double z) {
double r28847645 = x;
double r28847646 = y;
double r28847647 = z;
double r28847648 = r28847646 - r28847647;
double r28847649 = 1.0;
double r28847650 = r28847648 + r28847649;
double r28847651 = r28847645 * r28847650;
double r28847652 = r28847651 / r28847647;
return r28847652;
}
double f(double x, double y, double z) {
double r28847653 = 1.0;
double r28847654 = y;
double r28847655 = r28847653 + r28847654;
double r28847656 = x;
double r28847657 = z;
double r28847658 = r28847656 / r28847657;
double r28847659 = r28847655 * r28847658;
double r28847660 = r28847659 - r28847656;
return r28847660;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.4 |
|---|---|
| Target | 0.4 |
| Herbie | 1.8 |
Initial program 10.4
rmApplied *-un-lft-identity10.4
Applied times-frac3.3
Simplified3.3
Taylor expanded around 0 3.7
Simplified1.8
Final simplification1.8
herbie shell --seed 2019200
(FPCore (x y z)
:name "Diagrams.TwoD.Segment.Bernstein:evaluateBernstein from diagrams-lib-1.3.0.3"
:herbie-target
(if (< x -2.71483106713436e-162) (- (* (+ 1.0 y) (/ x z)) x) (if (< x 3.874108816439546e-197) (* (* x (+ (- y z) 1.0)) (/ 1.0 z)) (- (* (+ 1.0 y) (/ x z)) x)))
(/ (* x (+ (- y z) 1.0)) z))