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

[walter at geodynamics.org]

Author: walter
Date: 2007-03-26 05:33:19 -0700 (Mon, 26 Mar 2007)
New Revision: 6397

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c
Log:
 r1743 at earth:  boo | 2007-03-26 05:32:06 -0700
 Make the cell centered mesh agree with the vertex centered mesh



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

Modified: long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c
===================================================================
--- long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c	2007-03-26 06:01:18 UTC (rev 6396)
+++ long/3D/Gale/trunk/src/Underworld/plugins/EulerDeform/EulerDeform.c	2007-03-26 12:33:19 UTC (rev 6397)
@@ -455,8 +455,15 @@
 
 	assert( edCtx );
 
-	for( sys_i = 0; sys_i < edCtx->nSystems; sys_i++ ) {
-		EulerDeform_System*	sys = edCtx->systems + sys_i;
+
+        /* We redo the the first system, because that is the cell
+           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.  */
+
+	for( sys_i = 0; sys_i < edCtx->nSystems+1; sys_i++ ) {
+		EulerDeform_System*	sys = edCtx->systems + sys_i%2;
 		double**		oldCrds;
 		double**		newCrds;
 		unsigned		nDomainNodes;