Mercurial > hg > octave-nkf
changeset 11652:ef95e842ba81 release-3-0-x
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);