\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;x \le -1.309325788103135 \cdot 10^{135}:\\
\;\;\;\;\frac{1}{\frac{z}{x \cdot \frac{\sin y}{y}}}\\
\mathbf{elif}\;x \le 1.15923258913015539 \cdot 10^{113}:\\
\;\;\;\;\frac{x}{\frac{z}{\frac{\sin y}{y}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{z}}{\frac{y}{x \cdot \sin 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 ((x <= -1.309325788103135e+135)) {
VAR = (1.0 / (z / (x * (sin(y) / y))));
} else {
double VAR_1;
if ((x <= 1.1592325891301554e+113)) {
VAR_1 = (x / (z / (sin(y) / y)));
} else {
VAR_1 = ((1.0 / z) / (y / (x * sin(y))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.8 |
|---|---|
| Target | 0.3 |
| Herbie | 0.8 |
if x < -1.309325788103135e+135Initial program 0.2
rmApplied clear-num0.5
if -1.309325788103135e+135 < x < 1.1592325891301554e+113Initial program 3.6
rmApplied associate-/l*1.0
if 1.1592325891301554e+113 < x Initial program 0.2
rmApplied clear-num0.6
rmApplied div-inv0.6
Applied associate-/r*0.4
Taylor expanded around inf 0.4
Final simplification0.8
herbie shell --seed 2020091
(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))