\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 r545088 = x;
double r545089 = y;
double r545090 = z;
double r545091 = r545089 - r545090;
double r545092 = 1.0;
double r545093 = r545091 + r545092;
double r545094 = r545088 * r545093;
double r545095 = r545094 / r545090;
return r545095;
}
double f(double x, double y, double z) {
double r545096 = x;
double r545097 = z;
double r545098 = r545096 / r545097;
double r545099 = y;
double r545100 = 1.0;
double r545101 = r545099 + r545100;
double r545102 = r545098 * r545101;
double r545103 = r545102 - r545096;
return r545103;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.1 |
|---|---|
| Target | 0.5 |
| Herbie | 1.6 |
Initial program 10.1
Simplified9.1
rmApplied div-inv9.2
Applied associate-*l*3.5
Simplified3.4
Taylor expanded around 0 3.5
Simplified1.8
Taylor expanded around 0 3.5
Simplified1.6
Final simplification1.6
herbie shell --seed 2019195
(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))