\tan \left(x + \varepsilon\right) - \tan x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -1.58790806524755116 \cdot 10^{-31} \lor \neg \left(\varepsilon \le 2.27794334900708537 \cdot 10^{-160}\right):\\
\;\;\;\;\frac{\tan x + \tan \varepsilon}{1 - \frac{\sin x \cdot \sin \varepsilon}{\cos x \cdot \cos \varepsilon}} - \tan x\\
\mathbf{else}:\\
\;\;\;\;\mathsf{fma}\left({\varepsilon}^{2}, x, \mathsf{fma}\left(\varepsilon, {x}^{2}, \varepsilon\right)\right)\\
\end{array}double code(double x, double eps) {
return (tan((x + eps)) - tan(x));
}
double code(double x, double eps) {
double VAR;
if (((eps <= -1.5879080652475512e-31) || !(eps <= 2.2779433490070854e-160))) {
VAR = (((tan(x) + tan(eps)) / (1.0 - ((sin(x) * sin(eps)) / (cos(x) * cos(eps))))) - tan(x));
} else {
VAR = fma(pow(eps, 2.0), x, fma(eps, pow(x, 2.0), eps));
}
return VAR;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.4 |
|---|---|
| Target | 15.1 |
| Herbie | 15.8 |
if eps < -1.5879080652475512e-31 or 2.2779433490070854e-160 < eps Initial program 30.8
rmApplied tan-sum8.2
rmApplied tan-quot8.2
Applied tan-quot8.2
Applied frac-times8.2
if -1.5879080652475512e-31 < eps < 2.2779433490070854e-160Initial program 47.6
Taylor expanded around 0 30.8
Simplified30.8
Final simplification15.8
herbie shell --seed 2020092 +o rules:numerics
(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)))