\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\frac{x}{z} \cdot \left(1 + y\right) - xdouble f(double x, double y, double z) {
double r611251 = x;
double r611252 = y;
double r611253 = z;
double r611254 = r611252 - r611253;
double r611255 = 1.0;
double r611256 = r611254 + r611255;
double r611257 = r611251 * r611256;
double r611258 = r611257 / r611253;
return r611258;
}
double f(double x, double y, double z) {
double r611259 = x;
double r611260 = z;
double r611261 = r611259 / r611260;
double r611262 = 1.0;
double r611263 = y;
double r611264 = r611262 + r611263;
double r611265 = r611261 * r611264;
double r611266 = r611265 - r611259;
return r611266;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.3 |
|---|---|
| Target | 0.4 |
| Herbie | 1.6 |
Initial program 10.3
Taylor expanded around 0 3.7
Simplified1.6
Final simplification1.6
herbie shell --seed 2019303
(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))