\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\frac{x}{z} \cdot \left(y + 1\right) - xdouble f(double x, double y, double z) {
double r418782 = x;
double r418783 = y;
double r418784 = z;
double r418785 = r418783 - r418784;
double r418786 = 1.0;
double r418787 = r418785 + r418786;
double r418788 = r418782 * r418787;
double r418789 = r418788 / r418784;
return r418789;
}
double f(double x, double y, double z) {
double r418790 = x;
double r418791 = z;
double r418792 = r418790 / r418791;
double r418793 = y;
double r418794 = 1.0;
double r418795 = r418793 + r418794;
double r418796 = r418792 * r418795;
double r418797 = r418796 - r418790;
return r418797;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.1 |
|---|---|
| Target | 0.4 |
| Herbie | 1.5 |
Initial program 10.1
Taylor expanded around 0 3.5
Simplified1.5
Final simplification1.5
herbie shell --seed 2019304 +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.7148310671343599e-162) (- (* (+ 1 y) (/ x z)) x) (if (< x 3.87410881643954616e-197) (* (* x (+ (- y z) 1)) (/ 1 z)) (- (* (+ 1 y) (/ x z)) x)))
(/ (* x (+ (- y z) 1)) z))