\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;x \cdot \frac{\sin y}{y} \le -1.4187362230985416 \cdot 10^{-273}:\\
\;\;\;\;\frac{x \cdot \frac{\sin y}{y}}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x}{\sqrt[3]{z} \cdot \sqrt[3]{z}}}{\frac{y}{\sin y} \cdot \sqrt[3]{z}}\\
\end{array}double code(double x, double y, double z) {
return ((x * (sin(y) / y)) / z);
}
double code(double x, double y, double z) {
double temp;
if (((x * (sin(y) / y)) <= -1.4187362230985416e-273)) {
temp = ((x * (sin(y) / y)) / z);
} else {
temp = ((x / (cbrt(z) * cbrt(z))) / ((y / sin(y)) * cbrt(z)));
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.3 |
| Herbie | 1.7 |
if (* x (/ (sin y) y)) < -1.4187362230985416e-273Initial program 0.2
if -1.4187362230985416e-273 < (* x (/ (sin y) y)) Initial program 4.7
rmApplied *-un-lft-identity4.7
Applied *-commutative4.7
Applied times-frac2.6
Simplified2.6
rmApplied add-cube-cbrt3.4
Applied associate-/r*3.4
Applied clear-num3.4
Applied frac-times2.8
Simplified2.8
Final simplification1.7
herbie shell --seed 2020066 +o rules:numerics
(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 (/ y (sin y)))) z) (if (< z 4.446702369113811e+64) (/ x (* z (/ y (sin y)))) (/ (* x (/ 1 (/ y (sin y)))) z)))
(/ (* x (/ (sin y) y)) z))