Mercurial > octave-nkf

changeset 17845:0f912af348e0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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	Mon Nov 04 16:51:35 2013 +0000
+++ b/scripts/sparse/eigs.m	Mon Nov 04 08:49:23 2013 +0100
@@ -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]);
 
+