\frac{x \cdot \left(\left(y - z\right) + 1\right)}{z}\left(x \cdot \left(\sqrt[3]{y + 1} \cdot \frac{\sqrt[3]{y + 1}}{\sqrt[3]{z} \cdot \sqrt[3]{z}}\right)\right) \cdot \frac{\sqrt[3]{y + 1}}{\sqrt[3]{z}} - xdouble code(double x, double y, double z) {
return ((double) (((double) (x * ((double) (((double) (y - z)) + 1.0)))) / z));
}
double code(double x, double y, double z) {
return ((double) (((double) (((double) (x * ((double) (((double) cbrt(((double) (y + 1.0)))) * ((double) (((double) cbrt(((double) (y + 1.0)))) / ((double) (((double) cbrt(z)) * ((double) cbrt(z)))))))))) * ((double) (((double) cbrt(((double) (y + 1.0)))) / ((double) cbrt(z)))))) - x));
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 10.1 |
|---|---|
| Target | 0.4 |
| Herbie | 1.2 |
Initial program 10.1
Simplified3.2
rmApplied add-cube-cbrt3.8
Applied add-cube-cbrt3.8
Applied times-frac3.8
Applied associate-*r*1.2
Simplified1.2
Final simplification1.2
herbie shell --seed 2020184
(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.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))