\tan \left(x + \varepsilon\right) - \tan x
\left(\frac{\frac{\sin x}{\cos x}}{1 - \frac{\frac{\sin x \cdot \sin \varepsilon}{\cos x}}{\cos \varepsilon}} - \frac{\sin x}{\cos x}\right) + \frac{1}{\frac{\cos \varepsilon \cdot \left(1 - \frac{\frac{\sin \varepsilon}{\cos x} \cdot \sin x}{\cos \varepsilon}\right)}{\sin \varepsilon}}double f(double x, double eps) {
double r2976647 = x;
double r2976648 = eps;
double r2976649 = r2976647 + r2976648;
double r2976650 = tan(r2976649);
double r2976651 = tan(r2976647);
double r2976652 = r2976650 - r2976651;
return r2976652;
}
double f(double x, double eps) {
double r2976653 = x;
double r2976654 = sin(r2976653);
double r2976655 = cos(r2976653);
double r2976656 = r2976654 / r2976655;
double r2976657 = 1.0;
double r2976658 = eps;
double r2976659 = sin(r2976658);
double r2976660 = r2976654 * r2976659;
double r2976661 = r2976660 / r2976655;
double r2976662 = cos(r2976658);
double r2976663 = r2976661 / r2976662;
double r2976664 = r2976657 - r2976663;
double r2976665 = r2976656 / r2976664;
double r2976666 = r2976665 - r2976656;
double r2976667 = r2976659 / r2976655;
double r2976668 = r2976667 * r2976654;
double r2976669 = r2976668 / r2976662;
double r2976670 = r2976657 - r2976669;
double r2976671 = r2976662 * r2976670;
double r2976672 = r2976671 / r2976659;
double r2976673 = r2976657 / r2976672;
double r2976674 = r2976666 + r2976673;
return r2976674;
}




Bits error versus x




Bits error versus eps
Results
| Original | 36.1 |
|---|---|
| Target | 14.3 |
| Herbie | 12.8 |
Initial program 36.1
rmApplied tan-sum21.7
Taylor expanded around inf 21.8
Simplified12.7
rmApplied clear-num12.8
Taylor expanded around inf 12.8
Final simplification12.8
herbie shell --seed 2019158 +o rules:numerics
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:herbie-target
(/ (sin eps) (* (cos x) (cos (+ x eps))))
(- (tan (+ x eps)) (tan x)))