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





Bits error versus x





Bits error versus y





Bits error versus z
Results
| Original | 10.3 |
|---|---|
| Target | 0.4 |
| Herbie | 0.6 |
| Alternative 1 | |
|---|---|
| Accuracy | 1.5 |
| Cost | 576 |
Initial program 10.3
Taylor expanded around 0 3.6
Simplified1.5
rmApplied add-cube-cbrt_binary64_168221.8
Applied add-cube-cbrt_binary64_168221.9
Applied times-frac_binary64_167931.9
Applied associate-*l*_binary64_167280.6
Simplified0.6
Final simplification0.6
herbie shell --seed 2020322
(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))