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 <= -6.583395804010651e-55)) {
VAR = ((double) (((double) (((double) (((double) tan(x)) + ((double) tan(eps)))) / ((double) (1.0 - ((double) (((double) (((double) tan(x)) * ((double) sin(eps)))) / ((double) cos(eps)))))))) - ((double) tan(x))));
} else {
double VAR_1;
if ((eps <= 2.6174394417726606e-73)) {
VAR_1 = ((double) (eps + ((double) (x * ((double) (eps * ((double) (eps + x))))))));
} else {
VAR_1 = ((double) (((double) (((double) (((double) (((double) (((double) cos(eps)) * ((double) sin(x)))) + ((double) (((double) sin(eps)) * ((double) cos(x)))))) / ((double) (((double) cos(eps)) * ((double) cos(x)))))) / ((double) (1.0 - ((double) (((double) tan(x)) * ((double) tan(eps)))))))) - ((double) tan(x))));
}
VAR = VAR_1;
}
return VAR;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.9 |
|---|---|
| Target | 15.1 |
| Herbie | 15.4 |
if eps < -6.5833958040106514e-55Initial program 30.1
rmApplied tan-sum4.1
rmApplied tan-quot4.2
Applied associate-*r/4.2
if -6.5833958040106514e-55 < eps < 2.61743944177266056e-73Initial program 47.8
Taylor expanded around 0 31.6
Simplified31.3
if 2.61743944177266056e-73 < eps Initial program 29.7
rmApplied tan-sum5.9
rmApplied tan-quot6.1
Applied tan-quot6.1
Applied frac-add6.1
Final simplification15.4
herbie shell --seed 2020191
(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)))