\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;x \leq 2.2974278247676085 \cdot 10^{+25}:\\
\;\;\;\;\frac{\sin y}{y} \cdot \frac{x}{z}\\
\mathbf{else}:\\
\;\;\;\;\frac{x \cdot \frac{\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 2.2974278247676085e+25) (* (/ (sin y) y) (/ x z)) (/ (* 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 <= 2.2974278247676085e+25) {
tmp = (sin(y) / y) * (x / z);
} else {
tmp = (x * (sin(y) / y)) / z;
}
return tmp;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.3 |
| Herbie | 1.5 |
if x < 2.29742782476760854e25Initial program 3.5
rmApplied add-cube-cbrt_binary64_130714.3
Applied associate-/r*_binary64_129804.3
Simplified2.6
rmApplied *-un-lft-identity_binary64_130362.6
Applied times-frac_binary64_130422.7
Simplified2.7
Simplified1.9
if 2.29742782476760854e25 < x Initial program 0.2
Final simplification1.5
herbie shell --seed 2021098
(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))