\tan \left(x + \varepsilon\right) - \tan x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -5.246119333473003277498362100151355132534 \cdot 10^{-164} \lor \neg \left(\varepsilon \le 5.307417583108514168007488766475587468342 \cdot 10^{-118}\right):\\
\;\;\;\;\frac{\left(\left(\tan x + \tan \varepsilon\right) \cdot \cos x - \sin x\right) - \left(-\sin x \cdot \frac{\sin x \cdot \tan \varepsilon}{\cos x}\right)}{\cos x + \left(-\cos x \cdot \frac{\sin x \cdot \tan \varepsilon}{\cos x}\right)}\\
\mathbf{else}:\\
\;\;\;\;\varepsilon + \left(x \cdot \varepsilon\right) \cdot \left(\varepsilon + x\right)\\
\end{array}double f(double x, double eps) {
double r82935 = x;
double r82936 = eps;
double r82937 = r82935 + r82936;
double r82938 = tan(r82937);
double r82939 = tan(r82935);
double r82940 = r82938 - r82939;
return r82940;
}
double f(double x, double eps) {
double r82941 = eps;
double r82942 = -5.246119333473003e-164;
bool r82943 = r82941 <= r82942;
double r82944 = 5.307417583108514e-118;
bool r82945 = r82941 <= r82944;
double r82946 = !r82945;
bool r82947 = r82943 || r82946;
double r82948 = x;
double r82949 = tan(r82948);
double r82950 = tan(r82941);
double r82951 = r82949 + r82950;
double r82952 = cos(r82948);
double r82953 = r82951 * r82952;
double r82954 = sin(r82948);
double r82955 = r82953 - r82954;
double r82956 = r82954 * r82950;
double r82957 = r82956 / r82952;
double r82958 = r82954 * r82957;
double r82959 = -r82958;
double r82960 = r82955 - r82959;
double r82961 = r82952 * r82957;
double r82962 = -r82961;
double r82963 = r82952 + r82962;
double r82964 = r82960 / r82963;
double r82965 = r82948 * r82941;
double r82966 = r82941 + r82948;
double r82967 = r82965 * r82966;
double r82968 = r82941 + r82967;
double r82969 = r82947 ? r82964 : r82968;
return r82969;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.5 |
|---|---|
| Target | 15.4 |
| Herbie | 16.2 |
if eps < -5.246119333473003e-164 or 5.307417583108514e-118 < eps Initial program 32.4
rmApplied tan-sum11.4
rmApplied tan-quot11.4
Applied associate-*l/11.4
rmApplied add-cube-cbrt11.5
rmApplied tan-quot11.5
Applied frac-sub11.5
Simplified10.7
Simplified10.5
if -5.246119333473003e-164 < eps < 5.307417583108514e-118Initial program 51.2
rmApplied tan-sum51.2
Taylor expanded around 0 31.9
Simplified31.7
Final simplification16.2
herbie shell --seed 2019294
(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)))