# HG changeset patch # User John W. Eaton # Date 1203123265 18000 # Node ID 6470f946a4250cbee4801ab9380162851497b395 # Parent 1e01db14700bb1d3a53592dc5daaf7fd14b20eda another small xGELSD workspace fix diff -r 1e01db14700b -r 6470f946a425 liboctave/CMatrix.cc --- 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; diff -r 1e01db14700b -r 6470f946a425 liboctave/dMatrix.cc --- 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);