\tan \left(x + \varepsilon\right) - \tan x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -6.095793693328762137298872934344524249228 \cdot 10^{-10}:\\
\;\;\;\;\frac{\cos x \cdot \left(\tan \varepsilon + \tan x\right) - \left(1 - \tan x \cdot \tan \varepsilon\right) \cdot \sin x}{\left(1 - \tan x \cdot \tan \varepsilon\right) \cdot \cos x}\\
\mathbf{elif}\;\varepsilon \le 1.066275681806479256320193537654211119112 \cdot 10^{-47}:\\
\;\;\;\;\varepsilon + \varepsilon \cdot \left(x \cdot \left(x + \varepsilon\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\left(\left(\tan \varepsilon + \tan x\right) \cdot \left(1 - \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) \cdot \cos x - \sin x \cdot \left(\left(1 - \tan x \cdot \tan \varepsilon\right) \cdot \left(1 - \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right)}{\left(\left(1 - \tan x \cdot \tan \varepsilon\right) \cdot \left(1 - \left(\tan x \cdot \tan \varepsilon\right) \cdot \left(\tan x \cdot \tan \varepsilon\right)\right)\right) \cdot \cos x}\\
\end{array}double f(double x, double eps) {
double r7032785 = x;
double r7032786 = eps;
double r7032787 = r7032785 + r7032786;
double r7032788 = tan(r7032787);
double r7032789 = tan(r7032785);
double r7032790 = r7032788 - r7032789;
return r7032790;
}
double f(double x, double eps) {
double r7032791 = eps;
double r7032792 = -6.095793693328762e-10;
bool r7032793 = r7032791 <= r7032792;
double r7032794 = x;
double r7032795 = cos(r7032794);
double r7032796 = tan(r7032791);
double r7032797 = tan(r7032794);
double r7032798 = r7032796 + r7032797;
double r7032799 = r7032795 * r7032798;
double r7032800 = 1.0;
double r7032801 = r7032797 * r7032796;
double r7032802 = r7032800 - r7032801;
double r7032803 = sin(r7032794);
double r7032804 = r7032802 * r7032803;
double r7032805 = r7032799 - r7032804;
double r7032806 = r7032802 * r7032795;
double r7032807 = r7032805 / r7032806;
double r7032808 = 1.0662756818064793e-47;
bool r7032809 = r7032791 <= r7032808;
double r7032810 = r7032794 + r7032791;
double r7032811 = r7032794 * r7032810;
double r7032812 = r7032791 * r7032811;
double r7032813 = r7032791 + r7032812;
double r7032814 = r7032801 * r7032801;
double r7032815 = r7032800 - r7032814;
double r7032816 = r7032798 * r7032815;
double r7032817 = r7032816 * r7032795;
double r7032818 = r7032802 * r7032815;
double r7032819 = r7032803 * r7032818;
double r7032820 = r7032817 - r7032819;
double r7032821 = r7032818 * r7032795;
double r7032822 = r7032820 / r7032821;
double r7032823 = r7032809 ? r7032813 : r7032822;
double r7032824 = r7032793 ? r7032807 : r7032823;
return r7032824;
}




Bits error versus x




Bits error versus eps
Results
| Original | 37.1 |
|---|---|
| Target | 15.1 |
| Herbie | 15.5 |
if eps < -6.095793693328762e-10Initial program 30.2
rmApplied tan-quot30.1
Applied tan-sum0.6
Applied frac-sub0.6
if -6.095793693328762e-10 < eps < 1.0662756818064793e-47Initial program 45.5
Taylor expanded around 0 31.3
Simplified31.3
if 1.0662756818064793e-47 < eps Initial program 30.0
rmApplied tan-sum3.7
rmApplied flip--3.7
Applied associate-/r/3.7
rmApplied tan-quot3.8
Applied flip-+3.8
Applied frac-times3.8
Applied frac-sub3.8
Final simplification15.5
herbie shell --seed 2019171
(FPCore (x eps)
:name "2tan (problem 3.3.2)"
:herbie-target
(/ (sin eps) (* (cos x) (cos (+ x eps))))
(- (tan (+ x eps)) (tan x)))