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




Bits error versus x




Bits error versus eps
Results
| Original | 37.0 |
|---|---|
| Target | 15.0 |
| Herbie | 0.4 |
if eps < -3.51660557575655631e-7Initial program 29.7
rmApplied tan-sum_binary64_19180.4
rmApplied tan-quot_binary64_19420.4
Applied associate-*r/_binary64_17250.4
Simplified0.4
rmApplied flip--_binary64_17580.5
Applied associate-/r/_binary64_17290.5
Simplified0.5
if -3.51660557575655631e-7 < eps < 4.45944387844819457e-10Initial program 44.6
Taylor expanded around 0 0.2
if 4.45944387844819457e-10 < eps Initial program 29.6
rmApplied tan-sum_binary64_19180.5
rmApplied tan-quot_binary64_19420.5
Applied associate-*r/_binary64_17250.5
Final simplification0.4
herbie shell --seed 2021028
(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)))