Mercurial > hg > octave-avbm
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 +++ b/liboctave/ChangeLog @@ -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 +++ b/liboctave/dMatrix.cc @@ -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)