Index: mttroot/mtt/bin/trans/make_ode2odes
==================================================================
--- mttroot/mtt/bin/trans/make_ode2odes
+++ mttroot/mtt/bin/trans/make_ode2odes
@@ -7,10 +7,13 @@
 ###############################################################
 ## Version control history
 ###############################################################
 ## $Id$
 ## $Log$
+## Revision 1.50  2001/02/18 09:18:49  geraint
+## Removed temporary Matrices from mtt_implicit.cc
+##
 ## Revision 1.49  2001/02/14 06:06:34  geraint
 ## Removed octave_value_list wrappers from standalone.exe - speed improvements
 ##
 ## Revision 1.48  2001/02/11 07:08:59  geraint
 ## Static declarations of octave_value_lists: small .exe speed improvement
@@ -300,11 +303,11 @@
  
 extern ColumnVector Fmtt_implicit (
 	ColumnVector &x,
 	ColumnVector &dx,
 	Matrix &AA,
-	const ColumnVector &AAx,
+	ColumnVector &AAx,
 	const double &ddt,
 	const int &nx,
 	const ColumnVector &open_switches);
 
 extern ColumnVector F${sys}_${ode} (
@@ -394,11 +397,11 @@
 
 inline ColumnVector
 mtt_implicit (ColumnVector &x,
 	      ColumnVector &dx,
 	      Matrix &AA,
-	      const ColumnVector &AAx,
+	      ColumnVector &AAx,
 	      const double &ddt,
 	      const int &nx,
 	      const ColumnVector &open_switches)
 {
 #ifdef STANDALONE

Index: mttroot/mtt/lib/cc/mtt_implicit.cc
==================================================================
--- mttroot/mtt/lib/cc/mtt_implicit.cc
+++ mttroot/mtt/lib/cc/mtt_implicit.cc
@@ -3,11 +3,11 @@
 
 #ifdef STANDALONE
 ColumnVector Fmtt_implicit (      ColumnVector	&x,
 				  ColumnVector	&dx,
 			          Matrix	&AA,
-			    const ColumnVector	&AAx,
+				  ColumnVector	&AAx,
 			    const double	&t,
 			    const int		&Nx,
 			    const ColumnVector	&openx)
 {
 #else // !STANDALONE
@@ -37,14 +37,15 @@
 
   for (row = 0; row < Nx; row++)
     {
       if (0 != openx (row))
 	{
+	  AAx (row) = 0.0;
           dx (row) = 0.0;
 	  for (col = 0; col < Nx; col++)
 	    {
-	      AA (row,col) = AA (col,row) = 0.0;
+	      AA (row,col) = 0.0;
 	    }
 	}
     }
 
   x = static_cast<ColumnVector> (xleftdiv (AA, static_cast<Matrix>(AAx + dx * t)));