\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;z \le -2.74989504026001211 \cdot 10^{-85} \lor \neg \left(z \le 5.2866687389386045 \cdot 10^{26}\right):\\
\;\;\;\;\left(x \cdot \frac{\sin y}{y}\right) \cdot \frac{1}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{1}{\frac{\frac{\sin y}{z}}{y}}}\\
\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 VAR;
if (((z <= -2.749895040260012e-85) || !(z <= 5.2866687389386045e+26))) {
VAR = ((x * (sin(y) / y)) * (1.0 / z));
} else {
VAR = (x / (1.0 / ((sin(y) / z) / y)));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
if z < -2.749895040260012e-85 or 5.2866687389386045e+26 < z Initial program 0.3
rmApplied div-inv0.4
if -2.749895040260012e-85 < z < 5.2866687389386045e+26Initial program 6.5
rmApplied associate-/l*0.2
rmApplied clear-num0.3
Simplified0.4
Final simplification0.4
herbie shell --seed 2020078 +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))