\frac{1 - \cos x}{\sin x}\begin{array}{l}
\mathbf{if}\;x \le -0.014503396267940102:\\
\;\;\;\;\frac{\frac{1}{\sin x} \cdot \frac{1}{\sin x} - \frac{\cos x}{\sin x} \cdot \frac{\cos x}{\sin x}}{\frac{1}{\sin x} + \frac{\cos x}{\sin x}}\\
\mathbf{elif}\;x \le 0.024683183633800634:\\
\;\;\;\;\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}:\\
\;\;\;\;1 \cdot \frac{1 - \cos x}{\sin x}\\
\end{array}double code(double x) {
return ((double) (((double) (1.0 - ((double) cos(x)))) / ((double) sin(x))));
}
double code(double x) {
double VAR;
if ((x <= -0.014503396267940102)) {
VAR = ((double) (((double) (((double) (((double) (1.0 / ((double) sin(x)))) * ((double) (1.0 / ((double) sin(x)))))) - ((double) (((double) (((double) cos(x)) / ((double) sin(x)))) * ((double) (((double) cos(x)) / ((double) sin(x)))))))) / ((double) (((double) (1.0 / ((double) sin(x)))) + ((double) (((double) cos(x)) / ((double) sin(x))))))));
} else {
double VAR_1;
if ((x <= 0.024683183633800634)) {
VAR_1 = ((double) fma(0.041666666666666664, ((double) pow(x, 3.0)), ((double) fma(0.004166666666666667, ((double) pow(x, 5.0)), ((double) (0.5 * x))))));
} else {
VAR_1 = ((double) (1.0 * ((double) (((double) (1.0 - ((double) cos(x)))) / ((double) sin(x))))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x
Results
| Original | 31.0 |
|---|---|
| Target | 0.0 |
| Herbie | 0.6 |
if x < -0.014503396267940102Initial program 0.9
rmApplied div-sub1.2
rmApplied flip--1.6
if -0.014503396267940102 < x < 0.024683183633800634Initial program 59.9
Taylor expanded around 0 0.0
Simplified0.0
if 0.024683183633800634 < x Initial program 0.9
rmApplied add-log-exp1.0
rmApplied pow11.0
Applied log-pow1.0
Simplified0.9
Final simplification0.6
herbie shell --seed 2020113 +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)))