Mercurial > octave-nkf
changeset 11691:c6ce7274fb14 release-3-0-x
one more xGELSD workspace fix
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Wed, 12 Mar 2008 01:05:55 -0400 |
| parents | 42e936258b19 |
| children | 0cabd873264d |
| files | liboctave/ChangeLog liboctave/dMatrix.cc |
| diffstat | 2 files changed, 6 insertions(+), 2 deletions(-) [+] |
line wrap: on
line diff
--- a/liboctave/ChangeLog Tue Mar 11 23:38:00 2008 -0400 +++ b/liboctave/ChangeLog Wed Mar 12 01:05:55 2008 -0400 @@ -1,3 +1,7 @@ +2008-03-12 John W. Eaton <jwe@octave.org> + + * dMatrix.cc (Matrix::lssolve): One more xGELSD workspace fix. + 2008-03-07 John W. Eaton <jwe@octave.org> * idx-vector.cc (IDX_VEC_REP::idx_vector_rep (bool),
--- a/liboctave/dMatrix.cc Tue Mar 11 23:38:00 2008 -0400 +++ b/liboctave/dMatrix.cc Wed Mar 12 01:05:55 2008 -0400 @@ -2130,9 +2130,9 @@ double dminmn = static_cast<double> (minmn); double dsmlsizp1 = static_cast<double> (smlsiz+1); #if defined (HAVE_LOG2) - double tmp = log2 (dminmn) / dsmlsizp1; + double tmp = log2 (dminmn / dsmlsizp1); #else - double tmp = log (dminmn) / dsmlsizp1 / log (2.0); + double tmp = log (dminmn / dsmlsizp1) / log (2.0); #endif octave_idx_type nlvl = static_cast<octave_idx_type> (tmp) + 1; if (nlvl < 0)