\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
\mathbf{if}\;\varepsilon \leq -0.0001675426203356274:\\
\;\;\;\;\cos x \cdot \cos \varepsilon - \left(\cos x + \sin x \cdot \sin \varepsilon\right)\\
\mathbf{elif}\;\varepsilon \leq 0.00015167359697961248:\\
\;\;\;\;\varepsilon \cdot \left(\varepsilon \cdot \left(\varepsilon \cdot \left(\sin x \cdot 0.16666666666666666\right) + \cos x \cdot -0.5\right) - \sin x\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\cos x \cdot \cos \varepsilon - \sin x \cdot \sin \varepsilon\right) - \cos x\\
\end{array}(FPCore (x eps) :precision binary64 (- (cos (+ x eps)) (cos x)))
(FPCore (x eps)
:precision binary64
(if (<= eps -0.0001675426203356274)
(- (* (cos x) (cos eps)) (+ (cos x) (* (sin x) (sin eps))))
(if (<= eps 0.00015167359697961248)
(*
eps
(-
(* eps (+ (* eps (* (sin x) 0.16666666666666666)) (* (cos x) -0.5)))
(sin x)))
(- (- (* (cos x) (cos eps)) (* (sin x) (sin eps))) (cos x)))))double code(double x, double eps) {
return cos(x + eps) - cos(x);
}
double code(double x, double eps) {
double tmp;
if (eps <= -0.0001675426203356274) {
tmp = (cos(x) * cos(eps)) - (cos(x) + (sin(x) * sin(eps)));
} else if (eps <= 0.00015167359697961248) {
tmp = eps * ((eps * ((eps * (sin(x) * 0.16666666666666666)) + (cos(x) * -0.5))) - sin(x));
} else {
tmp = ((cos(x) * cos(eps)) - (sin(x) * sin(eps))) - cos(x);
}
return tmp;
}














Bits error versus x














Bits error versus eps
Results
| Alternative 1 | |
|---|---|
| Error | 0.5 |
| Cost | 32904 |
| Alternative 2 | |
|---|---|
| Error | 0.5 |
| Cost | 26504 |
| Alternative 3 | |
|---|---|
| Error | 14.9 |
| Cost | 13632 |
| Alternative 4 | |
|---|---|
| Error | 14.6 |
| Cost | 14018 |
| Alternative 5 | |
|---|---|
| Error | 20.8 |
| Cost | 13512 |
| Alternative 6 | |
|---|---|
| Error | 20.7 |
| Cost | 13634 |
| Alternative 7 | |
|---|---|
| Error | 21.0 |
| Cost | 6984 |
| Alternative 8 | |
|---|---|
| Error | 39.5 |
| Cost | 6592 |
| Alternative 9 | |
|---|---|
| Error | 52.3 |
| Cost | 706 |
| Alternative 10 | |
|---|---|
| Error | 56.2 |
| Cost | 64 |
| Alternative 11 | |
|---|---|
| Error | 61.0 |
| Cost | 64 |

if eps < -1.6754262033562739e-4Initial program 29.7
rmApplied cos-sum_binary64_8940.9
Applied associate--l-_binary64_6980.9
Simplified0.9
Simplified0.9
if -1.6754262033562739e-4 < eps < 1.51673596979612482e-4Initial program 50.1
Taylor expanded around 0 0.2
Simplified0.2
Simplified0.2
if 1.51673596979612482e-4 < eps Initial program 29.5
rmApplied cos-sum_binary64_8940.8
Simplified0.8
Final simplification0.5
herbie shell --seed 2021044
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))