Mercurial > octave-nkf

changeset 7488:6470f946a425

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 1e01db14700b
children 8e4592e49fa7
files liboctave/CMatrix.cc liboctave/dMatrix.cc
diffstat 2 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/liboctave/CMatrix.cc	Fri Feb 15 19:49:14 2008 -0500
+++ b/liboctave/CMatrix.cc	Fri Feb 15 19:54:25 2008 -0500
@@ -2464,10 +2464,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	Fri Feb 15 19:49:14 2008 -0500
+++ b/liboctave/dMatrix.cc	Fri Feb 15 19:54:25 2008 -0500
@@ -2067,10 +2067,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);