\frac{x \cdot \frac{\sin y}{y}}{z}\begin{array}{l}
\mathbf{if}\;z \le -8.5463407620259075 \cdot 10^{83} \lor \neg \left(z \le 1.00828897081930849 \cdot 10^{-64}\right):\\
\;\;\;\;\frac{\frac{x}{z}}{\frac{1}{\frac{\sin y}{y}}}\\
\mathbf{else}:\\
\;\;\;\;\frac{x}{\frac{z}{\frac{\sin y}{y}}}\\
\end{array}double f(double x, double y, double z) {
double r1546 = x;
double r1547 = y;
double r1548 = sin(r1547);
double r1549 = r1548 / r1547;
double r1550 = r1546 * r1549;
double r1551 = z;
double r1552 = r1550 / r1551;
return r1552;
}
double f(double x, double y, double z) {
double r1553 = z;
double r1554 = -8.546340762025908e+83;
bool r1555 = r1553 <= r1554;
double r1556 = 1.0082889708193085e-64;
bool r1557 = r1553 <= r1556;
double r1558 = !r1557;
bool r1559 = r1555 || r1558;
double r1560 = x;
double r1561 = r1560 / r1553;
double r1562 = 1.0;
double r1563 = y;
double r1564 = sin(r1563);
double r1565 = r1564 / r1563;
double r1566 = r1562 / r1565;
double r1567 = r1561 / r1566;
double r1568 = r1553 / r1565;
double r1569 = r1560 / r1568;
double r1570 = r1559 ? r1567 : r1569;
return r1570;
}




Bits error versus x




Bits error versus y




Bits error versus z
Results
| Original | 2.7 |
|---|---|
| Target | 0.3 |
| Herbie | 0.4 |
if z < -8.546340762025908e+83 or 1.0082889708193085e-64 < z Initial program 0.2
rmApplied clear-num0.9
rmApplied div-inv1.0
Applied associate-/r*0.4
rmApplied *-un-lft-identity0.4
Applied times-frac0.5
Applied associate-/r*0.5
Simplified0.3
if -8.546340762025908e+83 < z < 1.0082889708193085e-64Initial program 5.4
rmApplied clear-num5.6
rmApplied div-inv6.1
Applied associate-/r*6.0
rmApplied *-un-lft-identity6.0
Applied *-un-lft-identity6.0
Applied times-frac6.0
Applied associate-/l*6.2
Simplified5.6
rmApplied *-un-lft-identity5.6
Applied times-frac1.0
Applied associate-/r*0.6
Simplified0.5
Final simplification0.4
herbie shell --seed 2020025 +o rules:numerics
(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))