Mercurial > hg > octave-nkf

diff liboctave/dMatrix.cc @ 10813:2c2d4a2f1047

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
fix workspace bug workaround for xGELSD (since 7486:6a6d2abe51ff)
author Jaroslav Hajek <highegg@gmail.com>
date Thu, 22 Jul 2010 09:54:03 +0200
parents 7c542263a92a
children 3622db30ff05
line wrap: on
line diff
--- a/liboctave/dMatrix.cc
+++ b/liboctave/dMatrix.cc
@@ -2294,7 +2294,7 @@
       // through 3.1.1 when n >= mnthr.  The obtuse formula below
       // should provide sufficient workspace for DGELSD to operate
       // efficiently.
-      if (n >= mnthr)
+      if (n > m && n >= mnthr)
         {
           const octave_idx_type wlalsd
             = 9*m + 2*m*smlsiz + 8*m*nlvl + m*nrhs + (smlsiz+1)*(smlsiz+1);