\tan \left(x + \varepsilon\right) - \tan x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -2.686281805292481 \cdot 10^{-23}:\\
\;\;\;\;\frac{\left(\frac{\tan \varepsilon + \tan x}{1 - \tan x \cdot \tan \varepsilon} + \tan x\right) \cdot \left(\frac{\tan \varepsilon + \tan x}{1 - \tan x \cdot \tan \varepsilon} - \tan x\right)}{\frac{\tan \varepsilon + \tan x}{1 - \tan x \cdot \tan \varepsilon} + \tan x}\\
\mathbf{elif}\;\varepsilon \le 7.203202361914666 \cdot 10^{-37}:\\
\;\;\;\;\varepsilon + \left(x + \varepsilon\right) \cdot \left(x \cdot \varepsilon\right)\\
\mathbf{else}:\\
\;\;\;\;\left(1 + \tan x \cdot \tan \varepsilon\right) \cdot \frac{\tan \varepsilon + \tan x}{1 - \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)} - \tan x\\
\end{array}double f(double x, double eps) {
double r2500008 = x;
double r2500009 = eps;
double r2500010 = r2500008 + r2500009;
double r2500011 = tan(r2500010);
double r2500012 = tan(r2500008);
double r2500013 = r2500011 - r2500012;
return r2500013;
}
double f(double x, double eps) {
double r2500014 = eps;
double r2500015 = -2.686281805292481e-23;
bool r2500016 = r2500014 <= r2500015;
double r2500017 = tan(r2500014);
double r2500018 = x;
double r2500019 = tan(r2500018);
double r2500020 = r2500017 + r2500019;
double r2500021 = 1.0;
double r2500022 = r2500019 * r2500017;
double r2500023 = r2500021 - r2500022;
double r2500024 = r2500020 / r2500023;
double r2500025 = r2500024 + r2500019;
double r2500026 = r2500024 - r2500019;
double r2500027 = r2500025 * r2500026;
double r2500028 = r2500027 / r2500025;
double r2500029 = 7.203202361914666e-37;
bool r2500030 = r2500014 <= r2500029;
double r2500031 = r2500018 + r2500014;
double r2500032 = r2500018 * r2500014;
double r2500033 = r2500031 * r2500032;
double r2500034 = r2500014 + r2500033;
double r2500035 = r2500021 + r2500022;
double r2500036 = r2500022 * r2500022;
double r2500037 = r2500021 - r2500036;
double r2500038 = r2500020 / r2500037;
double r2500039 = r2500035 * r2500038;
double r2500040 = r2500039 - r2500019;
double r2500041 = r2500030 ? r2500034 : r2500040;
double r2500042 = r2500016 ? r2500028 : r2500041;
return r2500042;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.0 |
|---|---|
| Target | 15.4 |
| Herbie | 15.0 |
if eps < -2.686281805292481e-23Initial program 30.1
rmApplied tan-sum1.5
rmApplied flip--1.6
Simplified1.5
if -2.686281805292481e-23 < eps < 7.203202361914666e-37Initial program 45.7
Taylor expanded around 0 31.0
Simplified31.0
if 7.203202361914666e-37 < eps Initial program 30.0
rmApplied tan-sum2.5
rmApplied flip--2.6
Applied associate-/r/2.6
Simplified2.6
Final simplification15.0
herbie shell --seed 2019152
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:herbie-target
(/ (sin eps) (* (cos x) (cos (+ x eps))))
(- (tan (+ x eps)) (tan x)))