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




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.6 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
if (/ (* x (/ (sin y) y)) z) < -9.58817383040855e-299 or 0.0 < (/ (* x (/ (sin y) y)) z) Initial program 0.4
rmApplied clear-num0.4
if -9.58817383040855e-299 < (/ (* x (/ (sin y) y)) z) < 0.0Initial program 8.7
rmApplied *-un-lft-identity8.7
Applied times-frac0.3
Simplified0.3
rmApplied div-inv0.3
Applied associate-/l*0.3
Simplified0.3
Final simplification0.4
herbie shell --seed 2020105 +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))