\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \leq -0.0029089008421326404 \lor \neg \left(\varepsilon \leq 2.0310490321845372 \cdot 10^{-06}\right):\\
\;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\
\mathbf{else}:\\
\;\;\;\;-2 \cdot \left(\sin \left(\frac{\varepsilon}{2}\right) \cdot \sin \left(\frac{x + \left(\varepsilon + x\right)}{2}\right)\right)\\
\end{array}(FPCore (x eps) :precision binary64 (- (cos (+ x eps)) (cos x)))
(FPCore (x eps) :precision binary64 (if (or (<= eps -0.0029089008421326404) (not (<= eps 2.0310490321845372e-06))) (- (- (* (cos x) (cos eps)) (* (sin x) (sin eps))) (cos x)) (* -2.0 (* (sin (/ eps 2.0)) (sin (/ (+ x (+ eps x)) 2.0))))))
double code(double x, double eps) {
return ((double) (((double) cos(((double) (x + eps)))) - ((double) cos(x))));
}
double code(double x, double eps) {
double tmp;
if (((eps <= -0.0029089008421326404) || !(eps <= 2.0310490321845372e-06))) {
tmp = ((double) (((double) (((double) (((double) cos(x)) * ((double) cos(eps)))) - ((double) (((double) sin(x)) * ((double) sin(eps)))))) - ((double) cos(x))));
} else {
tmp = ((double) (-2.0 * ((double) (((double) sin((eps / 2.0))) * ((double) sin((((double) (x + ((double) (eps + x)))) / 2.0)))))));
}
return tmp;
}



Bits error versus x



Bits error versus eps
Results
if eps < -0.0029089008421326404 or 2.03104903218453725e-6 < eps Initial program 30.3
rmApplied cos-sum_binary641.0
if -0.0029089008421326404 < eps < 2.03104903218453725e-6Initial program 49.6
rmApplied diff-cos_binary6437.9
Simplified0.5
Final simplification0.7
herbie shell --seed 2020210
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))