[cig-commits] r6398 - in long/3D/Gale/trunk: . src/Underworld/plugins/EulerDeform

[walter at geodynamics.org]

Author: walter
Date: 2007-03-26 07:15:01 -0700 (Mon, 26 Mar 2007)
New Revision: 6398

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c
Log:
 r1747 at earth:  boo | 2007-03-26 07:12:32 -0700
 Make remeshing only happen once, not twice, for cell-centered grids



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1743
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:1747

Modified: long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c
===================================================================
--- long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c	2007-03-26 12:33:19 UTC (rev 6397)
+++ long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c	2007-03-26 14:15:01 UTC (rev 6398)
@@ -456,14 +456,14 @@
 	assert( edCtx );
 
 
-        /* We redo the the first system, because that is the cell
+        /* We do the second system first, because that is the vertex
            centered system.  We want the cell centered and vertex
-           centered systems to be compatible, but the vertex centered
-           coordinates are not available at first.  There is some
-           problem with doing the vertex centered first.  */
+           centered systems to be compatible, so the cell centered
+           system is based on the vertex centered one. */
 
-	for( sys_i = 0; sys_i < edCtx->nSystems+1; sys_i++ ) {
-		EulerDeform_System*	sys = edCtx->systems + sys_i%2;
+	for( sys_i = 1; sys_i < edCtx->nSystems+1; sys_i++ ) {
+		EulerDeform_System*	sys = edCtx->systems
+                  + sys_i%edCtx->nSystems;
 		double**		oldCrds;
 		double**		newCrds;
 		unsigned		nDomainNodes;