\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;z \le -5.3424544320379512 \cdot 10^{-68}:\\
\;\;\;\;\frac{\frac{x}{z}}{\frac{1}{\frac{\sin y}{y}}}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{\frac{\sin y}{y}}{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 ((z <= -5.342454432037951e-68)) {
temp = ((x / z) / (1.0 / (sin(y) / y)));
} else {
temp = (x * ((sin(y) / y) / z));
}
return temp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.7 |
|---|---|
| Target | 0.3 |
| Herbie | 1.6 |
if z < -5.342454432037951e-68Initial program 0.3
rmApplied associate-/l*4.5
rmApplied div-inv4.5
Applied associate-/r*0.4
if -5.342454432037951e-68 < z Initial program 3.9
rmApplied *-un-lft-identity3.9
Applied times-frac2.2
Simplified2.2
Final simplification1.6
herbie shell --seed 2020058 +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))