Mercurial > hg > octave-max

changeset 11652:ef95e842ba81 release-3-0-x

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
another small xGELSD workspace fix
author John W. Eaton <jwe@octave.org>
date Fri, 15 Feb 2008 19:54:25 -0500
parents 74de76325d12
children 5d50ebf98273
files liboctave/CMatrix.cc liboctave/dMatrix.cc
diffstat 2 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/liboctave/CMatrix.cc
+++ b/liboctave/CMatrix.cc
@@ -2546,10 +2546,10 @@
 				 lwork, prwork, piwork, info));
 
       // The workspace query is broken in at least LAPACK 3.0.0
-      // through 3.1.1 when n > mnthr.  The obtuse formula below
+      // through 3.1.1 when n >= mnthr.  The obtuse formula below
       // should provide sufficient workspace for ZGELSD to operate
       // efficiently.
-      if (n > mnthr)
+      if (n >= mnthr)
 	{
 	  octave_idx_type addend = m;
 
--- a/liboctave/dMatrix.cc
+++ b/liboctave/dMatrix.cc
@@ -2149,10 +2149,10 @@
 				 lwork, piwork, info));
 
       // The workspace query is broken in at least LAPACK 3.0.0
-      // through 3.1.1 when n > mnthr.  The obtuse formula below
+      // 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 >= mnthr)
 	{
 	  const octave_idx_type wlalsd
 	    = 9*m + 2*m*smlsiz + 8*m*nlvl + m*nrhs + (smlsiz+1)*(smlsiz+1);