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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 14.5 |
|---|---|
| Target | 4.2 |
| Herbie | 1.4 |
Initial program 14.5
rmApplied associate-/r*_binary6412.9
Simplified10.8
rmApplied add-cube-cbrt_binary6411.1
Applied times-frac_binary646.1
Applied associate-*r*_binary641.4
Final simplification1.4
herbie shell --seed 2020260
(FPCore (x y z)
:name "Statistics.Distribution.Beta:$cvariance from math-functions-0.1.5.2"
:precision binary64
:herbie-target
(if (< z 249.6182814532307) (/ (* y (/ x z)) (+ z (* z z))) (/ (* (/ (/ y z) (+ 1.0 z)) x) z))
(/ (* x y) (* (* z z) (+ z 1.0))))