\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\left(\frac{x \cdot y}{z} + 1 \cdot \frac{x}{z}\right) - xdouble f(double x, double y, double z) {
double r481802 = x;
double r481803 = y;
double r481804 = z;
double r481805 = r481803 - r481804;
double r481806 = 1.0;
double r481807 = r481805 + r481806;
double r481808 = r481802 * r481807;
double r481809 = r481808 / r481804;
return r481809;
}
double f(double x, double y, double z) {
double r481810 = x;
double r481811 = y;
double r481812 = r481810 * r481811;
double r481813 = z;
double r481814 = r481812 / r481813;
double r481815 = 1.0;
double r481816 = r481810 / r481813;
double r481817 = r481815 * r481816;
double r481818 = r481814 + r481817;
double r481819 = r481818 - r481810;
return r481819;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.3 |
|---|---|
| Target | 0.4 |
| Herbie | 3.6 |
Initial program 10.3
Taylor expanded around 0 3.6
Simplified1.6
rmApplied distribute-lft-in1.6
Simplified3.6
Simplified3.6
Final simplification3.6
herbie shell --seed 2019306 +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))