\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;x \leq -9.153932362511965 \cdot 10^{-41}:\\
\;\;\;\;\frac{x \cdot \frac{\sin y}{y}}{z}\\
\mathbf{elif}\;x \leq 4.5434271798398006 \cdot 10^{-18}:\\
\;\;\;\;\frac{x}{z \cdot \frac{y}{\sin y}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{x \cdot \sin y}{y}}{z}\\
\end{array}(FPCore (x y z) :precision binary64 (/ (* x (/ (sin y) y)) z))
(FPCore (x y z)
:precision binary64
(if (<= x -9.153932362511965e-41)
(/ (* x (/ (sin y) y)) z)
(if (<= x 4.5434271798398006e-18)
(/ x (* z (/ y (sin y))))
(/ (/ (* x (sin y)) y) z))))double code(double x, double y, double z) {
return (x * (sin(y) / y)) / z;
}
double code(double x, double y, double z) {
double tmp;
if (x <= -9.153932362511965e-41) {
tmp = (x * (sin(y) / y)) / z;
} else if (x <= 4.5434271798398006e-18) {
tmp = x / (z * (y / sin(y)));
} else {
tmp = ((x * sin(y)) / y) / z;
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.2 |
if x < -9.1539323625119649e-41Initial program 0.3
if -9.1539323625119649e-41 < x < 4.5434271798398006e-18Initial program 5.1
rmApplied clear-num_binary64_152145.1
rmApplied un-div-inv_binary64_152135.1
Applied associate-/l/_binary64_151620.1
if 4.5434271798398006e-18 < x Initial program 0.2
rmApplied associate-*r/_binary64_151570.3
Final simplification0.2
herbie shell --seed 2020292
(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))