\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;x \cdot \frac{\sin y}{y} \leq -3.058329899204841 \cdot 10^{-279}:\\
\;\;\;\;\frac{\frac{x}{\frac{y}{\sin y}}}{z}\\
\mathbf{elif}\;x \cdot \frac{\sin y}{y} \leq 5.031922861997968 \cdot 10^{-125}:\\
\;\;\;\;\frac{x}{\frac{y}{\sin y} \cdot 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 (/ (sin y) y)) -3.058329899204841e-279)
(/ (/ x (/ y (sin y))) z)
(if (<= (* x (/ (sin y) y)) 5.031922861997968e-125)
(/ x (* (/ y (sin y)) 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 * (sin(y) / y)) <= -3.058329899204841e-279) {
tmp = (x / (y / sin(y))) / z;
} else if ((x * (sin(y) / y)) <= 5.031922861997968e-125) {
tmp = x / ((y / sin(y)) * 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.5 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
| Alternative 1 | |
|---|---|
| Error | 0.3 |
| Cost | 7176 |
| Alternative 2 | |
|---|---|
| Error | 0.6 |
| Cost | 7176 |
| Alternative 3 | |
|---|---|
| Error | 3.1 |
| Cost | 6848 |
| Alternative 4 | |
|---|---|
| Error | 22.5 |
| Cost | 832 |
| Alternative 5 | |
|---|---|
| Error | 22.9 |
| Cost | 1346 |
| Alternative 6 | |
|---|---|
| Error | 22.7 |
| Cost | 834 |
| Alternative 7 | |
|---|---|
| Error | 47.9 |
| Cost | 64 |
| Alternative 8 | |
|---|---|
| Error | 61.7 |
| Cost | 64 |

if (*.f64 x (/.f64 (sin.f64 y) y)) < -3.0583298992048409e-279Initial program 0.2
rmApplied clear-num_binary64_143990.2
rmApplied un-div-inv_binary64_143980.2
Simplified0.2
if -3.0583298992048409e-279 < (*.f64 x (/.f64 (sin.f64 y) y)) < 5.03192286199796776e-125Initial program 7.4
rmApplied clear-num_binary64_143997.4
rmApplied un-div-inv_binary64_143987.4
Applied associate-/l/_binary64_143470.8
Simplified0.8
if 5.03192286199796776e-125 < (*.f64 x (/.f64 (sin.f64 y) y)) Initial program 0.2
Simplified0.2
Final simplification0.4
herbie shell --seed 2021044
(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))