\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 r687198 = x;
double r687199 = y;
double r687200 = z;
double r687201 = r687199 - r687200;
double r687202 = 1.0;
double r687203 = r687201 + r687202;
double r687204 = r687198 * r687203;
double r687205 = r687204 / r687200;
return r687205;
}
double f(double x, double y, double z) {
double r687206 = x;
double r687207 = z;
double r687208 = r687206 / r687207;
double r687209 = 1.0;
double r687210 = y;
double r687211 = r687209 + r687210;
double r687212 = r687208 * r687211;
double r687213 = r687212 - r687206;
return r687213;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.2 |
|---|---|
| Target | 0.5 |
| Herbie | 1.6 |
Initial program 10.2
rmApplied associate-/l*3.3
Taylor expanded around 0 3.6
Taylor expanded around 0 3.6
Simplified1.6
Final simplification1.6
herbie shell --seed 2019354
(FPCore (x y z)
:name "Diagrams.TwoD.Segment.Bernstein:evaluateBernstein from diagrams-lib-1.3.0.3"
:precision binary64
:herbie-target
(if (< x -2.71483106713436e-162) (- (* (+ 1 y) (/ x z)) x) (if (< x 3.874108816439546e-197) (* (* x (+ (- y z) 1)) (/ 1 z)) (- (* (+ 1 y) (/ x z)) x)))
(/ (* x (+ (- y z) 1)) z))