\tan \left(x + \varepsilon\right) - \tan x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -5.17742371337321258 \cdot 10^{-51}:\\
\;\;\;\;\frac{\tan x + \tan \varepsilon}{1 - \sqrt[3]{{\left(\tan x \cdot \tan \varepsilon\right)}^{3}}} - \tan x\\
\mathbf{elif}\;\varepsilon \le 8.3497052348174857 \cdot 10^{-111}:\\
\;\;\;\;\left(\varepsilon \cdot x\right) \cdot \left(x + \varepsilon\right) + \varepsilon\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\left(\tan x \cdot \tan x - \tan \varepsilon \cdot \tan \varepsilon\right) \cdot \frac{1}{1 - \tan x \cdot \tan \varepsilon}\right) \cdot \cos x - \left(\tan x - \tan \varepsilon\right) \cdot \sin x}{\left(\tan x - \tan \varepsilon\right) \cdot \cos x}\\
\end{array}double code(double x, double eps) {
return ((double) (((double) tan(((double) (x + eps)))) - ((double) tan(x))));
}
double code(double x, double eps) {
double VAR;
if ((eps <= -5.1774237133732126e-51)) {
VAR = ((double) (((double) (((double) (((double) tan(x)) + ((double) tan(eps)))) / ((double) (1.0 - ((double) cbrt(((double) pow(((double) (((double) tan(x)) * ((double) tan(eps)))), 3.0)))))))) - ((double) tan(x))));
} else {
double VAR_1;
if ((eps <= 8.349705234817486e-111)) {
VAR_1 = ((double) (((double) (((double) (eps * x)) * ((double) (x + eps)))) + eps));
} else {
VAR_1 = ((double) (((double) (((double) (((double) (((double) (((double) (((double) tan(x)) * ((double) tan(x)))) - ((double) (((double) tan(eps)) * ((double) tan(eps)))))) * ((double) (1.0 / ((double) (1.0 - ((double) (((double) tan(x)) * ((double) tan(eps)))))))))) * ((double) cos(x)))) - ((double) (((double) (((double) tan(x)) - ((double) tan(eps)))) * ((double) sin(x)))))) / ((double) (((double) (((double) tan(x)) - ((double) tan(eps)))) * ((double) cos(x))))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.6 |
|---|---|
| Target | 15.2 |
| Herbie | 15.3 |
if eps < -5.1774237133732126e-51Initial program 29.1
rmApplied tan-sum3.4
rmApplied add-cbrt-cube3.4
Applied add-cbrt-cube3.5
Applied cbrt-unprod3.4
Simplified3.4
if -5.1774237133732126e-51 < eps < 8.349705234817486e-111Initial program 47.5
Taylor expanded around 0 31.0
Simplified30.8
if 8.349705234817486e-111 < eps Initial program 31.7
rmApplied tan-sum9.5
rmApplied div-inv9.5
rmApplied tan-quot9.5
Applied flip-+9.6
Applied associate-*l/9.6
Applied frac-sub9.6
Final simplification15.3
herbie shell --seed 2020122
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:precision binary64
:herbie-target
(/ (sin eps) (* (cos x) (cos (+ x eps))))
(- (tan (+ x eps)) (tan x)))