| Alternative 1 | |
|---|---|
| Error | 0.5 |
| Cost | 33288 |
(FPCore (x eps) :precision binary64 (- (cos (+ x eps)) (cos x)))
(FPCore (x eps)
:precision binary64
(let* ((t_0 (* (sin eps) (sin x))))
(if (<= eps -0.0056)
(- (* (cos x) (cos eps)) (+ (cos x) (* (sin x) (sin eps))))
(if (<= eps 0.0045)
(+
(* (sin x) (- (sin eps)))
(*
(cos x)
(+ (* 0.041666666666666664 (pow eps 4.0)) (* -0.5 (pow eps 2.0)))))
(+ (* t_0 -2.0) (+ t_0 (* (cos x) (+ (cos eps) -1.0))))))))double code(double x, double eps) {
return cos((x + eps)) - cos(x);
}
double code(double x, double eps) {
double t_0 = sin(eps) * sin(x);
double tmp;
if (eps <= -0.0056) {
tmp = (cos(x) * cos(eps)) - (cos(x) + (sin(x) * sin(eps)));
} else if (eps <= 0.0045) {
tmp = (sin(x) * -sin(eps)) + (cos(x) * ((0.041666666666666664 * pow(eps, 4.0)) + (-0.5 * pow(eps, 2.0))));
} else {
tmp = (t_0 * -2.0) + (t_0 + (cos(x) * (cos(eps) + -1.0)));
}
return tmp;
}
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
code = cos((x + eps)) - cos(x)
end function
real(8) function code(x, eps)
real(8), intent (in) :: x
real(8), intent (in) :: eps
real(8) :: t_0
real(8) :: tmp
t_0 = sin(eps) * sin(x)
if (eps <= (-0.0056d0)) then
tmp = (cos(x) * cos(eps)) - (cos(x) + (sin(x) * sin(eps)))
else if (eps <= 0.0045d0) then
tmp = (sin(x) * -sin(eps)) + (cos(x) * ((0.041666666666666664d0 * (eps ** 4.0d0)) + ((-0.5d0) * (eps ** 2.0d0))))
else
tmp = (t_0 * (-2.0d0)) + (t_0 + (cos(x) * (cos(eps) + (-1.0d0))))
end if
code = tmp
end function
public static double code(double x, double eps) {
return Math.cos((x + eps)) - Math.cos(x);
}
public static double code(double x, double eps) {
double t_0 = Math.sin(eps) * Math.sin(x);
double tmp;
if (eps <= -0.0056) {
tmp = (Math.cos(x) * Math.cos(eps)) - (Math.cos(x) + (Math.sin(x) * Math.sin(eps)));
} else if (eps <= 0.0045) {
tmp = (Math.sin(x) * -Math.sin(eps)) + (Math.cos(x) * ((0.041666666666666664 * Math.pow(eps, 4.0)) + (-0.5 * Math.pow(eps, 2.0))));
} else {
tmp = (t_0 * -2.0) + (t_0 + (Math.cos(x) * (Math.cos(eps) + -1.0)));
}
return tmp;
}
def code(x, eps): return math.cos((x + eps)) - math.cos(x)
def code(x, eps): t_0 = math.sin(eps) * math.sin(x) tmp = 0 if eps <= -0.0056: tmp = (math.cos(x) * math.cos(eps)) - (math.cos(x) + (math.sin(x) * math.sin(eps))) elif eps <= 0.0045: tmp = (math.sin(x) * -math.sin(eps)) + (math.cos(x) * ((0.041666666666666664 * math.pow(eps, 4.0)) + (-0.5 * math.pow(eps, 2.0)))) else: tmp = (t_0 * -2.0) + (t_0 + (math.cos(x) * (math.cos(eps) + -1.0))) return tmp
function code(x, eps) return Float64(cos(Float64(x + eps)) - cos(x)) end
function code(x, eps) t_0 = Float64(sin(eps) * sin(x)) tmp = 0.0 if (eps <= -0.0056) tmp = Float64(Float64(cos(x) * cos(eps)) - Float64(cos(x) + Float64(sin(x) * sin(eps)))); elseif (eps <= 0.0045) tmp = Float64(Float64(sin(x) * Float64(-sin(eps))) + Float64(cos(x) * Float64(Float64(0.041666666666666664 * (eps ^ 4.0)) + Float64(-0.5 * (eps ^ 2.0))))); else tmp = Float64(Float64(t_0 * -2.0) + Float64(t_0 + Float64(cos(x) * Float64(cos(eps) + -1.0)))); end return tmp end
function tmp = code(x, eps) tmp = cos((x + eps)) - cos(x); end
function tmp_2 = code(x, eps) t_0 = sin(eps) * sin(x); tmp = 0.0; if (eps <= -0.0056) tmp = (cos(x) * cos(eps)) - (cos(x) + (sin(x) * sin(eps))); elseif (eps <= 0.0045) tmp = (sin(x) * -sin(eps)) + (cos(x) * ((0.041666666666666664 * (eps ^ 4.0)) + (-0.5 * (eps ^ 2.0)))); else tmp = (t_0 * -2.0) + (t_0 + (cos(x) * (cos(eps) + -1.0))); end tmp_2 = tmp; end
code[x_, eps_] := N[(N[Cos[N[(x + eps), $MachinePrecision]], $MachinePrecision] - N[Cos[x], $MachinePrecision]), $MachinePrecision]
code[x_, eps_] := Block[{t$95$0 = N[(N[Sin[eps], $MachinePrecision] * N[Sin[x], $MachinePrecision]), $MachinePrecision]}, If[LessEqual[eps, -0.0056], N[(N[(N[Cos[x], $MachinePrecision] * N[Cos[eps], $MachinePrecision]), $MachinePrecision] - N[(N[Cos[x], $MachinePrecision] + N[(N[Sin[x], $MachinePrecision] * N[Sin[eps], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], If[LessEqual[eps, 0.0045], N[(N[(N[Sin[x], $MachinePrecision] * (-N[Sin[eps], $MachinePrecision])), $MachinePrecision] + N[(N[Cos[x], $MachinePrecision] * N[(N[(0.041666666666666664 * N[Power[eps, 4.0], $MachinePrecision]), $MachinePrecision] + N[(-0.5 * N[Power[eps, 2.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(t$95$0 * -2.0), $MachinePrecision] + N[(t$95$0 + N[(N[Cos[x], $MachinePrecision] * N[(N[Cos[eps], $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]
\cos \left(x + \varepsilon\right) - \cos x
\begin{array}{l}
t_0 := \sin \varepsilon \cdot \sin x\\
\mathbf{if}\;\varepsilon \leq -0.0056:\\
\;\;\;\;\cos x \cdot \cos \varepsilon - \left(\cos x + \sin x \cdot \sin \varepsilon\right)\\
\mathbf{elif}\;\varepsilon \leq 0.0045:\\
\;\;\;\;\sin x \cdot \left(-\sin \varepsilon\right) + \cos x \cdot \left(0.041666666666666664 \cdot {\varepsilon}^{4} + -0.5 \cdot {\varepsilon}^{2}\right)\\
\mathbf{else}:\\
\;\;\;\;t_0 \cdot -2 + \left(t_0 + \cos x \cdot \left(\cos \varepsilon + -1\right)\right)\\
\end{array}
Results
if eps < -0.00559999999999999994Initial program 30.9
Applied egg-rr0.8
Taylor expanded in x around -inf 0.8
if -0.00559999999999999994 < eps < 0.00449999999999999966Initial program 49.0
Applied egg-rr11.7
Simplified11.7
[Start]11.7 | \[ \left(\left(-\sin x \cdot \sin \varepsilon\right) + \left(-\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
|---|---|
rational_best-simplify-11 [=>]11.7 | \[ \left(\color{blue}{-1 \cdot \left(\sin x \cdot \sin \varepsilon\right)} + \left(-\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-2 [=>]11.7 | \[ \left(\color{blue}{\left(\sin x \cdot \sin \varepsilon\right) \cdot -1} + \left(-\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-11 [=>]11.7 | \[ \left(\left(\sin x \cdot \sin \varepsilon\right) \cdot -1 + \color{blue}{-1 \cdot \left(\sin x \cdot \sin \varepsilon\right)}\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-2 [=>]11.7 | \[ \left(\left(\sin x \cdot \sin \varepsilon\right) \cdot -1 + \color{blue}{\left(\sin x \cdot \sin \varepsilon\right) \cdot -1}\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-53 [=>]11.7 | \[ \color{blue}{\left(\sin x \cdot \sin \varepsilon\right) \cdot \left(-1 + -1\right)} + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-2 [=>]11.7 | \[ \color{blue}{\left(\sin \varepsilon \cdot \sin x\right)} \cdot \left(-1 + -1\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
metadata-eval [=>]11.7 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot \color{blue}{-2} + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-13 [=>]11.7 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \color{blue}{\left(0 - \sin x \cdot \sin \varepsilon\right)}\right)
\] |
rational_best-simplify-47 [=>]11.7 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \color{blue}{\left(\sin x \cdot \sin \varepsilon + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - 0\right)\right)}
\] |
rational_best-simplify-2 [=>]11.7 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \left(\color{blue}{\sin \varepsilon \cdot \sin x} + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - 0\right)\right)
\] |
rational_best-simplify-6 [=>]11.7 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \left(\sin \varepsilon \cdot \sin x + \color{blue}{\cos x \cdot \left(\cos \varepsilon + -1\right)}\right)
\] |
Taylor expanded in eps around inf 11.7
Simplified11.7
[Start]11.7 | \[ -2 \cdot \left(\sin x \cdot \sin \varepsilon\right) + \left(\cos x \cdot \left(\cos \varepsilon - 1\right) + \sin x \cdot \sin \varepsilon\right)
\] |
|---|---|
rational_best-simplify-1 [=>]11.7 | \[ -2 \cdot \left(\sin x \cdot \sin \varepsilon\right) + \color{blue}{\left(\sin x \cdot \sin \varepsilon + \cos x \cdot \left(\cos \varepsilon - 1\right)\right)}
\] |
rational_best-simplify-2 [=>]11.7 | \[ -2 \cdot \left(\sin x \cdot \sin \varepsilon\right) + \left(\color{blue}{\sin \varepsilon \cdot \sin x} + \cos x \cdot \left(\cos \varepsilon - 1\right)\right)
\] |
rational_best-simplify-2 [=>]11.7 | \[ -2 \cdot \left(\sin x \cdot \sin \varepsilon\right) + \left(\sin \varepsilon \cdot \sin x + \color{blue}{\left(\cos \varepsilon - 1\right) \cdot \cos x}\right)
\] |
rational_best-simplify-18 [<=]11.7 | \[ -2 \cdot \left(\sin x \cdot \sin \varepsilon\right) + \left(\sin \varepsilon \cdot \sin x + \color{blue}{\left(\cos \varepsilon + -1\right)} \cdot \cos x\right)
\] |
rational_best-simplify-43 [=>]11.7 | \[ \color{blue}{\left(\cos \varepsilon + -1\right) \cdot \cos x + \left(\sin \varepsilon \cdot \sin x + -2 \cdot \left(\sin x \cdot \sin \varepsilon\right)\right)}
\] |
rational_best-simplify-1 [=>]11.7 | \[ \color{blue}{\left(\sin \varepsilon \cdot \sin x + -2 \cdot \left(\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos \varepsilon + -1\right) \cdot \cos x}
\] |
Taylor expanded in eps around 0 0.1
Simplified0.1
[Start]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \left(0.041666666666666664 \cdot \left({\varepsilon}^{4} \cdot \cos x\right) + -0.5 \cdot \left({\varepsilon}^{2} \cdot \cos x\right)\right)
\] |
|---|---|
rational_best-simplify-1 [=>]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \color{blue}{\left(-0.5 \cdot \left({\varepsilon}^{2} \cdot \cos x\right) + 0.041666666666666664 \cdot \left({\varepsilon}^{4} \cdot \cos x\right)\right)}
\] |
rational_best-simplify-2 [=>]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \left(-0.5 \cdot \color{blue}{\left(\cos x \cdot {\varepsilon}^{2}\right)} + 0.041666666666666664 \cdot \left({\varepsilon}^{4} \cdot \cos x\right)\right)
\] |
rational_best-simplify-46 [=>]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \left(\color{blue}{\cos x \cdot \left(-0.5 \cdot {\varepsilon}^{2}\right)} + 0.041666666666666664 \cdot \left({\varepsilon}^{4} \cdot \cos x\right)\right)
\] |
rational_best-simplify-2 [=>]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \left(\cos x \cdot \left(-0.5 \cdot {\varepsilon}^{2}\right) + 0.041666666666666664 \cdot \color{blue}{\left(\cos x \cdot {\varepsilon}^{4}\right)}\right)
\] |
rational_best-simplify-46 [=>]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \left(\cos x \cdot \left(-0.5 \cdot {\varepsilon}^{2}\right) + \color{blue}{\cos x \cdot \left(0.041666666666666664 \cdot {\varepsilon}^{4}\right)}\right)
\] |
rational_best-simplify-53 [=>]0.1 | \[ \sin x \cdot \left(-\sin \varepsilon\right) + \color{blue}{\cos x \cdot \left(0.041666666666666664 \cdot {\varepsilon}^{4} + -0.5 \cdot {\varepsilon}^{2}\right)}
\] |
if 0.00449999999999999966 < eps Initial program 29.4
Applied egg-rr0.8
Simplified0.8
[Start]0.8 | \[ \left(\left(-\sin x \cdot \sin \varepsilon\right) + \left(-\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
|---|---|
rational_best-simplify-11 [=>]0.8 | \[ \left(\color{blue}{-1 \cdot \left(\sin x \cdot \sin \varepsilon\right)} + \left(-\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-2 [=>]0.8 | \[ \left(\color{blue}{\left(\sin x \cdot \sin \varepsilon\right) \cdot -1} + \left(-\sin x \cdot \sin \varepsilon\right)\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-11 [=>]0.8 | \[ \left(\left(\sin x \cdot \sin \varepsilon\right) \cdot -1 + \color{blue}{-1 \cdot \left(\sin x \cdot \sin \varepsilon\right)}\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-2 [=>]0.8 | \[ \left(\left(\sin x \cdot \sin \varepsilon\right) \cdot -1 + \color{blue}{\left(\sin x \cdot \sin \varepsilon\right) \cdot -1}\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-53 [=>]0.8 | \[ \color{blue}{\left(\sin x \cdot \sin \varepsilon\right) \cdot \left(-1 + -1\right)} + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-2 [=>]0.8 | \[ \color{blue}{\left(\sin \varepsilon \cdot \sin x\right)} \cdot \left(-1 + -1\right) + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
metadata-eval [=>]0.8 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot \color{blue}{-2} + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \left(-\sin x \cdot \sin \varepsilon\right)\right)
\] |
rational_best-simplify-13 [=>]0.8 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - \color{blue}{\left(0 - \sin x \cdot \sin \varepsilon\right)}\right)
\] |
rational_best-simplify-47 [=>]0.8 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \color{blue}{\left(\sin x \cdot \sin \varepsilon + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - 0\right)\right)}
\] |
rational_best-simplify-2 [=>]0.8 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \left(\color{blue}{\sin \varepsilon \cdot \sin x} + \left(\cos x \cdot \left(\cos \varepsilon + -1\right) - 0\right)\right)
\] |
rational_best-simplify-6 [=>]0.8 | \[ \left(\sin \varepsilon \cdot \sin x\right) \cdot -2 + \left(\sin \varepsilon \cdot \sin x + \color{blue}{\cos x \cdot \left(\cos \varepsilon + -1\right)}\right)
\] |
Final simplification0.5
| Alternative 1 | |
|---|---|
| Error | 0.5 |
| Cost | 33288 |
| Alternative 2 | |
|---|---|
| Error | 0.5 |
| Cost | 32708 |
| Alternative 3 | |
|---|---|
| Error | 0.5 |
| Cost | 26568 |
| Alternative 4 | |
|---|---|
| Error | 0.6 |
| Cost | 26440 |
| Alternative 5 | |
|---|---|
| Error | 14.3 |
| Cost | 20168 |
| Alternative 6 | |
|---|---|
| Error | 14.5 |
| Cost | 19976 |
| Alternative 7 | |
|---|---|
| Error | 15.1 |
| Cost | 13640 |
| Alternative 8 | |
|---|---|
| Error | 20.9 |
| Cost | 13256 |
| Alternative 9 | |
|---|---|
| Error | 32.5 |
| Cost | 7184 |
| Alternative 10 | |
|---|---|
| Error | 20.9 |
| Cost | 6920 |
| Alternative 11 | |
|---|---|
| Error | 35.8 |
| Cost | 6856 |
| Alternative 12 | |
|---|---|
| Error | 48.6 |
| Cost | 520 |
| Alternative 13 | |
|---|---|
| Error | 51.8 |
| Cost | 328 |
| Alternative 14 | |
|---|---|
| Error | 57.9 |
| Cost | 64 |
herbie shell --seed 2023097
(FPCore (x eps)
:name "2cos (problem 3.3.5)"
:precision binary64
(- (cos (+ x eps)) (cos x)))