[cig-commits] r3952 - in long/3D/Gale/trunk/src/StgFEM: . Discretisation/src

[walter at geodynamics.org]

Author: walter
Date: 2006-07-06 02:08:28 -0700 (Thu, 06 Jul 2006)
New Revision: 3952

Modified:
   long/3D/Gale/trunk/src/StgFEM/
   long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c
Log:
 r668 at earth:  boo | 2006-07-06 02:04:21 -0700
  r661 at earth (orig r602):  PatrickSunter | 2006-07-04 21:58:10 -0700
  Fixed a bug with loading a geometry from checkpoint:
   problems that get compressed in a particular dimension
   weren't being reloaded correctly. Fixed.
  
 



Property changes on: long/3D/Gale/trunk/src/StgFEM
___________________________________________________________________
Name: svk:merge
   - 38867592-cf10-0410-9e16-a142ea72ac34:/cig:667
db209038-57f2-0310-97fa-b160e0ae9d04:/trunk:601
   + 38867592-cf10-0410-9e16-a142ea72ac34:/cig:668
db209038-57f2-0310-97fa-b160e0ae9d04:/trunk:602

Modified: long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c
===================================================================
--- long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c	2006-07-06 09:08:25 UTC (rev 3951)
+++ long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c	2006-07-06 09:08:28 UTC (rev 3952)
@@ -1968,6 +1968,13 @@
 
 	dofAtEachNodeCount = self->fieldComponentCount;
 
+	/* Need to re-set the geometry here, in case we're loading from a checkpoint that had compression/squashing BCs,
+		and hence ended up with a smaller mesh than the original */
+	for ( dim_I = 0; dim_I < 3; dim_I++ ) {
+		geometry->min[dim_I] = HUGE_VAL;
+		geometry->max[dim_I] = -HUGE_VAL;
+	}
+
 	while ( !feof(inputFile) ) {
 		fscanf( inputFile, "%u ", &gNode_I );
 		lNode_I = Mesh_NodeMapGlobalToLocal( self->feMesh, gNode_I );