# HG changeset patch
# User John W. Eaton <jwe@octave.org>
# Date 1205298355 14400
# Node ID c6ce7274fb141682dfc394c6d6219316350a1dd7
# Parent  42e936258b19b70792b494e80b3a8335875a0ce5
one more xGELSD workspace fix

diff -r 42e936258b19 -r c6ce7274fb14 liboctave/ChangeLog
--- 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),
diff -r 42e936258b19 -r c6ce7274fb14 liboctave/dMatrix.cc
--- 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)