Mercurial > hg > octave-nkf
changeset 17845:0f912af348e0
add a test to mark regression in eigs
* eigs.m : add a test to show failure when eigs falls back to eig
author | Carlo de Falco <cdf@users.sourceforge.net> |
---|---|
date | Mon, 04 Nov 2013 08:49:23 +0100 |
parents | e88a8fa25125 |
children | 42453dcfa95e |
files | scripts/sparse/eigs.m |
diffstat | 1 files changed, 8 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- 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]); +