# HG changeset patch # User Carlo de Falco # Date 1383551363 -3600 # Mon Nov 04 08:49:23 2013 +0100 # Node ID da3a027cc9fbfc06f808c6d5656d0bc9cae38f83 # Parent 573f865aa8c8caf66f8b8032f86cca23f4ceaa8d add a test to mark regression in eigs * eigs.m : add a test to show failure when eigs falls back to eig diff --git a/scripts/sparse/eigs.m b/scripts/sparse/eigs.m --- a/scripts/sparse/eigs.m +++ b/scripts/sparse/eigs.m @@ -1104,7 +1104,15 @@ %! assert (max (abs ((A - d1(i)*eye (n))*v1(:,i))), 0, 1e-11); %! endfor +%!test +%! A = 2 * diag (ones (10, 1)) - diag (ones (9, 1), 1) - diag (ones (9, 1), -1); +%! B = diag (ones (10, 1)); +%! reseig = eig (A, B); +%! [~, idx] = sort (abs (reseig), 'ascend'); +%! assert (eigs (A, B, 10, 0), reseig (idx)) + %!assert (eigs (diag (1:5), 5, "sa"), [1;2;3;4;5]); %!assert (eigs (diag (1:5), 5, "la"), [5;4;3;2;1]); %!assert (eigs (diag (1:5), 3, "be"), [1;4;5]); +