\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -3.658449891821199522487853550423153592419 \cdot 10^{-6} \lor \neg \left(\varepsilon \le 3.130998694967744134270873183900629754817 \cdot 10^{-8}\right):\\
\;\;\;\;\mathsf{fma}\left(\cos x, \cos \varepsilon, -\left(\sin x \cdot \sin \varepsilon + \cos x\right)\right)\\
\mathbf{else}:\\
\;\;\;\;\varepsilon \cdot \left(\left(\frac{1}{6} \cdot {x}^{3} - x\right) - \varepsilon \cdot \frac{1}{2}\right)\\
\end{array}double f(double x, double eps) {
double r58846 = x;
double r58847 = eps;
double r58848 = r58846 + r58847;
double r58849 = cos(r58848);
double r58850 = cos(r58846);
double r58851 = r58849 - r58850;
return r58851;
}
double f(double x, double eps) {
double r58852 = eps;
double r58853 = -3.6584498918211995e-06;
bool r58854 = r58852 <= r58853;
double r58855 = 3.130998694967744e-08;
bool r58856 = r58852 <= r58855;
double r58857 = !r58856;
bool r58858 = r58854 || r58857;
double r58859 = x;
double r58860 = cos(r58859);
double r58861 = cos(r58852);
double r58862 = sin(r58859);
double r58863 = sin(r58852);
double r58864 = r58862 * r58863;
double r58865 = r58864 + r58860;
double r58866 = -r58865;
double r58867 = fma(r58860, r58861, r58866);
double r58868 = 0.16666666666666666;
double r58869 = 3.0;
double r58870 = pow(r58859, r58869);
double r58871 = r58868 * r58870;
double r58872 = r58871 - r58859;
double r58873 = 0.5;
double r58874 = r58852 * r58873;
double r58875 = r58872 - r58874;
double r58876 = r58852 * r58875;
double r58877 = r58858 ? r58867 : r58876;
return r58877;
}



Bits error versus x



Bits error versus eps
if eps < -3.6584498918211995e-06 or 3.130998694967744e-08 < eps Initial program 30.1
rmApplied cos-sum1.1
Applied associate--l-1.1
Simplified1.1
rmApplied fma-neg1.1
rmApplied fma-udef1.1
if -3.6584498918211995e-06 < eps < 3.130998694967744e-08Initial program 49.7
Taylor expanded around 0 32.6
Simplified32.6
Final simplification16.5
herbie shell --seed 2019347 +o rules:numerics
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))