\tan \left(x + \varepsilon\right) - \tan x
\frac{\sin \varepsilon}{\cos \varepsilon \cdot \left(1 - \frac{\sin \varepsilon \cdot \sin x}{\cos \varepsilon \cdot \cos x}\right)} + \frac{\sin x}{\cos x} \cdot \log \left(\frac{e^{\frac{1}{1 - \frac{\sin \varepsilon \cdot \sin x}{\cos \varepsilon \cdot \cos x}}}}{e}\right)(FPCore (x eps) :precision binary64 (- (tan (+ x eps)) (tan x)))
(FPCore (x eps)
:precision binary64
(+
(/
(sin eps)
(* (cos eps) (- 1.0 (/ (* (sin eps) (sin x)) (* (cos eps) (cos x))))))
(*
(/ (sin x) (cos x))
(log
(/
(exp (/ 1.0 (- 1.0 (/ (* (sin eps) (sin x)) (* (cos eps) (cos x))))))
E)))))double code(double x, double eps) {
return ((double) (((double) tan(((double) (x + eps)))) - ((double) tan(x))));
}
double code(double x, double eps) {
return ((double) ((((double) sin(eps)) / ((double) (((double) cos(eps)) * ((double) (1.0 - (((double) (((double) sin(eps)) * ((double) sin(x)))) / ((double) (((double) cos(eps)) * ((double) cos(x)))))))))) + ((double) ((((double) sin(x)) / ((double) cos(x))) * ((double) log((((double) exp((1.0 / ((double) (1.0 - (((double) (((double) sin(eps)) * ((double) sin(x)))) / ((double) (((double) cos(eps)) * ((double) cos(x)))))))))) / ((double) M_E))))))));
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.0 |
|---|---|
| Target | 15.3 |
| Herbie | 12.9 |
Initial program 37.0
rmApplied tan-sum_binary6421.6
rmApplied tan-quot_binary6421.6
Applied tan-quot_binary6421.6
Applied frac-times_binary6421.6
Taylor expanded around inf 21.7
Simplified12.8
rmApplied add-log-exp_binary6412.8
Applied add-log-exp_binary6412.9
Applied sum-log_binary6412.9
Simplified12.9
Final simplification12.9
herbie shell --seed 2020219
(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)))