[cig-commits] r16061 - mc/3D/CitcomCU/trunk/src

[tan2 at geodynamics.org]

Author: tan2
Date: 2009-12-02 18:18:49 -0800 (Wed, 02 Dec 2009)
New Revision: 16061

Modified:
   mc/3D/CitcomCU/trunk/src/Drive_solvers.c
Log:
Fixed a bug in non-newtonian solver, reported by Magali Billen.

Modified: mc/3D/CitcomCU/trunk/src/Drive_solvers.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/Drive_solvers.c	2009-12-02 18:44:58 UTC (rev 16060)
+++ mc/3D/CitcomCU/trunk/src/Drive_solvers.c	2009-12-03 02:18:49 UTC (rev 16061)
@@ -77,15 +77,15 @@
 	      damp = 0;
 	    }
 	  } /* end powerlaw */
-		oldU = (double *)malloc((neq + 2) * sizeof(double));
-		for(i = 1; i <= neq; i++)
+		oldU = (double *)malloc(neq * sizeof(double));
+		for(i = 0; i < neq; i++)
 			oldU[i] = 0.0;
 		visits++;
 	}
 
 	dUdot_mag = 0.0;
 
-	delta_U = (double *)malloc((neq + 2) * sizeof(double));
+	delta_U = (double *)malloc(neq * sizeof(double));
 
 	/* FIRST store the old velocity field */
 
@@ -124,10 +124,10 @@
 
 		  if(damp){
 		    /* add some of the old solution */
-		    for(i = 1; i <= neq; i++)
+		    for(i = 0; i < neq; i++)
 		      E->U[i] = alpha * E->U[i] + alpha1 * oldU[i];
 		  }
-			for(i = 1; i <= neq; i++)
+			for(i = 0; i < neq; i++)
 			{
 				delta_U[i] = E->U[i] - oldU[i];
 				oldU[i] = E->U[i];