\frac{1 - \cos x}{\sin x}\begin{array}{l}
\mathbf{if}\;x \le -0.0231416446907326372:\\
\;\;\;\;\frac{1}{\sin x} - \frac{\cos x}{\sin x}\\
\mathbf{elif}\;x \le 0.0251863311991388956:\\
\;\;\;\;\mathsf{fma}\left(\frac{1}{24}, {x}^{3}, \mathsf{fma}\left(\frac{1}{240}, {x}^{5}, \frac{1}{2} \cdot x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{e^{\log \left({1}^{3} - \sqrt[3]{{\left({\left(\cos x\right)}^{3}\right)}^{3}}\right)}}{\mathsf{fma}\left(\cos x, \sqrt[3]{{\left(1 + \cos x\right)}^{3}}, 1 \cdot 1\right)}}{\sin x}\\
\end{array}double code(double x) {
return ((1.0 - cos(x)) / sin(x));
}
double code(double x) {
double VAR;
if ((x <= -0.023141644690732637)) {
VAR = ((1.0 / sin(x)) - (cos(x) / sin(x)));
} else {
double VAR_1;
if ((x <= 0.025186331199138896)) {
VAR_1 = fma(0.041666666666666664, pow(x, 3.0), fma(0.004166666666666667, pow(x, 5.0), (0.5 * x)));
} else {
VAR_1 = ((exp(log((pow(1.0, 3.0) - cbrt(pow(pow(cos(x), 3.0), 3.0))))) / fma(cos(x), cbrt(pow((1.0 + cos(x)), 3.0)), (1.0 * 1.0))) / sin(x));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x
Results
| Original | 30.3 |
|---|---|
| Target | 0 |
| Herbie | 0.6 |
if x < -0.023141644690732637Initial program 0.9
rmApplied div-sub1.2
if -0.023141644690732637 < x < 0.025186331199138896Initial program 60.0
Taylor expanded around 0 0.0
Simplified0.0
if 0.025186331199138896 < x Initial program 1.0
rmApplied flip3--1.1
Simplified1.1
rmApplied add-cbrt-cube1.2
Simplified1.2
rmApplied add-exp-log1.2
rmApplied add-cbrt-cube1.2
Simplified1.2
Final simplification0.6
herbie shell --seed 2020106 +o rules:numerics
(FPCore (x)
:name "tanhf (example 3.4)"
:precision binary64
:herbie-expected 2
:herbie-target
(tan (/ x 2))
(/ (- 1 (cos x)) (sin x)))