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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.4 |
|---|---|
| Target | 0.3 |
| Herbie | 1.4 |
Initial program 2.4
rmApplied clear-num_binary64_126942.4
rmApplied un-div-inv_binary64_126932.4
rmApplied *-un-lft-identity_binary64_126952.4
Applied add-cube-cbrt_binary64_127302.7
Applied add-cube-cbrt_binary64_127303.3
Applied times-frac_binary64_127013.3
Applied times-frac_binary64_127011.4
Final simplification1.4
herbie shell --seed 2021015
(FPCore (x y z)
:name "Linear.Quaternion:$ctanh from linear-1.19.1.3"
:precision binary64
:herbie-target
(if (< z -4.2173720203427147e-29) (/ (* x (/ 1.0 (/ y (sin y)))) z) (if (< z 4.446702369113811e+64) (/ x (* z (/ y (sin y)))) (/ (* x (/ 1.0 (/ y (sin y)))) z)))
(/ (* x (/ (sin y) y)) z))