Mercurial > hg > octave-max
changeset 7488:6470f946a425
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 +++ b/liboctave/CMatrix.cc @@ -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 +++ b/liboctave/dMatrix.cc @@ -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);