Commit 92833d3c authored by MichelJuillard's avatar MichelJuillard
Browse files

In CHECK, use the eigenvalues as computed by the reordered real

generalized Schur decomposition, rather than the ones computed by
eig()
parent a22d1d41
......@@ -462,9 +462,6 @@ for i = 1:Size;
data(i).rank = rank(w(nd-nyf+1:end,nd-nyf+1:end));
dr.full_rank = dr.full_rank && (rcond(w(nd-nyf+1:end,nd- ...
nyf+1:end)) > 1e-9);
if ~exist('OCTAVE_VERSION','builtin')
data(i).eigval = eig(E,D);
end
dr.eigval = [dr.eigval ; data(i).eigval];
end
if (verbose)
......
......@@ -238,11 +238,6 @@ else
else
dr.full_rank = 1;
end
% Under Octave, eig(A,B) doesn't exist, and
% lambda = qz(A,B) won't return infinite eigenvalues
if ~exist('OCTAVE_VERSION')
dr.eigval = eig(E,D);
end
return
end
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment