Mercurial > hg > octave-lyh

diff scripts/linear-algebra/krylov.m @ 9872:72d6e0de76c7

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix qp, condest and krylov
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 26 Nov 2009 08:28:07 +0100
parents 71483d19204f
children 95c3e38098bf
line wrap: on
line diff
--- a/scripts/linear-algebra/krylov.m
+++ b/scripts/linear-algebra/krylov.m
@@ -81,10 +81,10 @@
     eps1 = defeps;
   endif
 
-  na = issquare (A);
-  if (! na)
-    error ("A(%d x %d) must be square", rows (A), columns (A));
+  if (! issquare (A) || isempty (A))
+    error ("A(%d x %d) must be non-empty square matrix", rows (A), columns (A));
   endif
+  na = rows (A);
 
   [m, kb] = size (V);
   if (m != na)