\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;x \le -9.1854363141215809 \cdot 10^{55}:\\
\;\;\;\;\frac{x \cdot \frac{1}{\frac{y}{\sin y}}}{z}\\
\mathbf{else}:\\
\;\;\;\;x \cdot \frac{\frac{\sin y}{y}}{z}\\
\end{array}double code(double x, double y, double z) {
return ((double) (((double) (x * ((double) (((double) sin(y)) / y)))) / z));
}
double code(double x, double y, double z) {
double VAR;
if ((x <= -9.185436314121581e+55)) {
VAR = ((double) (((double) (x * ((double) (1.0 / ((double) (y / ((double) sin(y)))))))) / z));
} else {
VAR = ((double) (x * ((double) (((double) (((double) sin(y)) / y)) / z))));
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.3 |
| Herbie | 1.6 |
if x < -9.1854363141215809e55Initial program 0.2
rmApplied clear-num0.3
if -9.1854363141215809e55 < x Initial program 3.3
rmApplied *-un-lft-identity3.3
Applied times-frac1.9
Simplified1.9
Final simplification1.6
herbie shell --seed 2020150
(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))