# HG changeset patch
# User jwe
# Date 869870400 0
# Node ID a6a00badcc12aba0bacc4ade5b83795fc60009f8
# Parent  42975c59d2a08f376b586e9530803a3fae62ba1e
[project @ 1997-07-25 22:39:56 by jwe]

diff --git a/liboctave/CMatrix.cc b/liboctave/CMatrix.cc
--- a/liboctave/CMatrix.cc
+++ b/liboctave/CMatrix.cc
@@ -1454,6 +1454,8 @@
       else
 	lwork = 2*n + (nrhs > m ? nrhs : m);
 
+      lwork *= 16;
+
       Array<Complex> work (lwork);
       Complex *pwork = work.fortran_vec ();
 
@@ -1534,6 +1536,8 @@
       else
 	lwork = 2*n + (nrhs > m ? nrhs : m);
 
+      lwork *= 16;
+
       Array<Complex> work (lwork);
       Complex *pwork = work.fortran_vec ();
 
diff --git a/liboctave/ChangeLog b/liboctave/ChangeLog
--- a/liboctave/ChangeLog
+++ b/liboctave/ChangeLog
@@ -1,3 +1,8 @@
+Fri Jul 25 17:31:26 1997  John W. Eaton  <jwe@bevo.che.wisc.edu>
+
+	* dMatrix.cc (Matrix::lssolve): Increase lwork by factor of 16.
+	* CMatrix.cc (ComplexMatrix::lssolve): Ditto.
+
 Thu Jul 24 14:32:48 1997  John W. Eaton  <jwe@bevo.che.wisc.edu>
 
 	* file-ops.cc (tilde_expand_word): Fix off-by-one error.
diff --git a/liboctave/dMatrix.cc b/liboctave/dMatrix.cc
--- a/liboctave/dMatrix.cc
+++ b/liboctave/dMatrix.cc
@@ -1166,6 +1166,8 @@
 		       ? (2*n > m ? 2*n : m)
 		       : (nrhs > m ? nrhs : m));
 
+      lwork *= 16;
+
       Array<double> work (lwork);
       double *pwork = work.fortran_vec ();
 
@@ -1267,6 +1269,8 @@
 		       ? (2*n > m ? 2*n : m)
 		       : (nrhs > m ? nrhs : m));
 
+      lwork *= 16;
+
       Array<double> work (lwork);
       double *pwork = work.fortran_vec ();