\tan \left(x + \varepsilon\right) - \tan x
\begin{array}{l}
\mathbf{if}\;\varepsilon \leq -5.846200722042183 \cdot 10^{-05}:\\
\;\;\;\;\frac{\tan x + \tan \varepsilon}{1 - \frac{\tan \varepsilon \cdot \sin x}{\cos x}} - \tan x\\
\mathbf{elif}\;\varepsilon \leq 2.2886471627852177 \cdot 10^{-05}:\\
\;\;\;\;\left(\frac{\varepsilon \cdot {\sin x}^{2}}{{\cos x}^{2}} + 1.3333333333333333 \cdot \frac{{\sin x}^{2} \cdot {\varepsilon}^{3}}{{\cos x}^{2}}\right) + \left(\varepsilon + \left({\varepsilon}^{3} \cdot \left(\frac{{\sin x}^{4}}{{\cos x}^{4}} + 0.3333333333333333\right) + \left(\varepsilon \cdot \varepsilon\right) \cdot \left(\frac{{\sin x}^{3}}{{\cos x}^{3}} + \frac{\sin x}{\cos x}\right)\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\tan x + \tan \varepsilon\right) \cdot \cos x - \sin x \cdot \left(1 - \tan x \cdot \tan \varepsilon\right)}{\cos x \cdot \left(1 - \tan x \cdot \tan \varepsilon\right)}\\
\end{array}(FPCore (x eps) :precision binary64 (- (tan (+ x eps)) (tan x)))
(FPCore (x eps)
:precision binary64
(if (<= eps -5.846200722042183e-05)
(-
(/ (+ (tan x) (tan eps)) (- 1.0 (/ (* (tan eps) (sin x)) (cos x))))
(tan x))
(if (<= eps 2.2886471627852177e-05)
(+
(+
(/ (* eps (pow (sin x) 2.0)) (pow (cos x) 2.0))
(*
1.3333333333333333
(/ (* (pow (sin x) 2.0) (pow eps 3.0)) (pow (cos x) 2.0))))
(+
eps
(+
(*
(pow eps 3.0)
(+ (/ (pow (sin x) 4.0) (pow (cos x) 4.0)) 0.3333333333333333))
(*
(* eps eps)
(+ (/ (pow (sin x) 3.0) (pow (cos x) 3.0)) (/ (sin x) (cos x)))))))
(/
(-
(* (+ (tan x) (tan eps)) (cos x))
(* (sin x) (- 1.0 (* (tan x) (tan eps)))))
(* (cos x) (- 1.0 (* (tan x) (tan eps))))))))double code(double x, double eps) {
return tan(x + eps) - tan(x);
}
double code(double x, double eps) {
double tmp;
if (eps <= -5.846200722042183e-05) {
tmp = ((tan(x) + tan(eps)) / (1.0 - ((tan(eps) * sin(x)) / cos(x)))) - tan(x);
} else if (eps <= 2.2886471627852177e-05) {
tmp = (((eps * pow(sin(x), 2.0)) / pow(cos(x), 2.0)) + (1.3333333333333333 * ((pow(sin(x), 2.0) * pow(eps, 3.0)) / pow(cos(x), 2.0)))) + (eps + ((pow(eps, 3.0) * ((pow(sin(x), 4.0) / pow(cos(x), 4.0)) + 0.3333333333333333)) + ((eps * eps) * ((pow(sin(x), 3.0) / pow(cos(x), 3.0)) + (sin(x) / cos(x))))));
} else {
tmp = (((tan(x) + tan(eps)) * cos(x)) - (sin(x) * (1.0 - (tan(x) * tan(eps))))) / (cos(x) * (1.0 - (tan(x) * tan(eps))));
}
return tmp;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.5 |
|---|---|
| Target | 15.2 |
| Herbie | 0.3 |
if eps < -5.84620072204218282e-5Initial program 29.9
rmApplied tan-sum_binary64_15770.4
rmApplied tan-quot_binary64_16010.4
Applied associate-*l/_binary64_13850.4
if -5.84620072204218282e-5 < eps < 2.28864716278521768e-5Initial program 44.2
Taylor expanded around 0 0.2
Simplified0.2
if 2.28864716278521768e-5 < eps Initial program 28.8
rmApplied tan-quot_binary64_160128.7
Applied tan-sum_binary64_15770.4
Applied frac-sub_binary64_14510.5
Final simplification0.3
herbie shell --seed 2021032
(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)))