\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \le -2.14002637841950426 \cdot 10^{-7}:\\
\;\;\;\;\frac{\sqrt[3]{{\left({\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right)}^{3}\right)}^{3}} - {\left(\cos x\right)}^{3}}{\left(\cos \varepsilon \cdot \cos x - \sin x \cdot \sin \varepsilon\right) \cdot \left(\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) + \cos x\right) + \cos x \cdot \cos x}\\
\mathbf{elif}\;\varepsilon \le 2.673419722608438 \cdot 10^{-8}:\\
\;\;\;\;\varepsilon \cdot \left(\left(\frac{1}{6} \cdot {x}^{3} - x\right) - \varepsilon \cdot \frac{1}{2}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\log \left(e^{\cos x \cdot \cos \varepsilon}\right) - \sin x \cdot \sin \varepsilon\right) - \cos x\\
\end{array}double f(double x, double eps) {
double r74388 = x;
double r74389 = eps;
double r74390 = r74388 + r74389;
double r74391 = cos(r74390);
double r74392 = cos(r74388);
double r74393 = r74391 - r74392;
return r74393;
}
double f(double x, double eps) {
double r74394 = eps;
double r74395 = -2.1400263784195043e-07;
bool r74396 = r74394 <= r74395;
double r74397 = x;
double r74398 = cos(r74397);
double r74399 = cos(r74394);
double r74400 = r74398 * r74399;
double r74401 = sin(r74397);
double r74402 = sin(r74394);
double r74403 = r74401 * r74402;
double r74404 = r74400 - r74403;
double r74405 = 3.0;
double r74406 = pow(r74404, r74405);
double r74407 = pow(r74406, r74405);
double r74408 = cbrt(r74407);
double r74409 = pow(r74398, r74405);
double r74410 = r74408 - r74409;
double r74411 = r74399 * r74398;
double r74412 = r74411 - r74403;
double r74413 = r74404 + r74398;
double r74414 = r74412 * r74413;
double r74415 = r74398 * r74398;
double r74416 = r74414 + r74415;
double r74417 = r74410 / r74416;
double r74418 = 2.6734197226084376e-08;
bool r74419 = r74394 <= r74418;
double r74420 = 0.16666666666666666;
double r74421 = pow(r74397, r74405);
double r74422 = r74420 * r74421;
double r74423 = r74422 - r74397;
double r74424 = 0.5;
double r74425 = r74394 * r74424;
double r74426 = r74423 - r74425;
double r74427 = r74394 * r74426;
double r74428 = exp(r74400);
double r74429 = log(r74428);
double r74430 = r74429 - r74403;
double r74431 = r74430 - r74398;
double r74432 = r74419 ? r74427 : r74431;
double r74433 = r74396 ? r74417 : r74432;
return r74433;
}



Bits error versus x



Bits error versus eps
Results
if eps < -2.1400263784195043e-07Initial program 31.0
rmApplied cos-sum1.1
rmApplied flip3--1.3
Simplified1.3
rmApplied add-cbrt-cube1.3
Simplified1.3
if -2.1400263784195043e-07 < eps < 2.6734197226084376e-08Initial program 48.7
Taylor expanded around 0 30.9
Simplified30.9
if 2.6734197226084376e-08 < eps Initial program 29.7
rmApplied cos-sum1.2
rmApplied add-log-exp1.4
Final simplification15.5
herbie shell --seed 2020056
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))