\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\frac{x \cdot \left(1 + y\right)}{z} - xdouble f(double x, double y, double z) {
double r485162 = x;
double r485163 = y;
double r485164 = z;
double r485165 = r485163 - r485164;
double r485166 = 1.0;
double r485167 = r485165 + r485166;
double r485168 = r485162 * r485167;
double r485169 = r485168 / r485164;
return r485169;
}
double f(double x, double y, double z) {
double r485170 = x;
double r485171 = 1.0;
double r485172 = y;
double r485173 = r485171 + r485172;
double r485174 = r485170 * r485173;
double r485175 = z;
double r485176 = r485174 / r485175;
double r485177 = r485176 - r485170;
return r485177;
}




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 pow11.6
Applied pow11.6
Applied pow-prod-down1.6
Simplified3.6
Final simplification3.6
herbie shell --seed 2019306
(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))