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

[walter at geodynamics.org]

Author: walter
Date: 2006-08-01 01:54:38 -0700 (Tue, 01 Aug 2006)
New Revision: 4159

Modified:
   long/3D/Gale/trunk/src/StgFEM/
   long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c
Log:
 r723 at earth:  boo | 2006-08-01 01:52:54 -0700
  r713 at earth (orig r617):  PatrickSunter | 2006-07-24 19:48:50 -0700
  Instructed the FeVariable ASCII file format save to use up
  to 15 significant figures to checkpoint field values of unknowns
  such as pressure and velocity.
  
  This is close to the limit of double precision on a 32-bit machine
  according to what I've read, so will minimise error when restarting
  (compared to the default 6 significant figures used by %g).
  
  Figure the file size increase is OK for now, since eventually we will probably
  use a binary dump format by default for FeVariables, to reduce size.
  
 



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

Modified: long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c
===================================================================
--- long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c	2006-08-01 08:54:35 UTC (rev 4158)
+++ long/3D/Gale/trunk/src/StgFEM/Discretisation/src/FeVariable.c	2006-08-01 08:54:38 UTC (rev 4159)
@@ -1998,10 +1998,10 @@
 		gNode_I = self->feMesh->nodeL2G[lNode_I];
 		fprintf( outputFile, "%u ", gNode_I );
 		coord = self->feMesh->nodeCoord[lNode_I];
-		fprintf( outputFile, "%g %g %g ", coord[0], coord[1], coord[2] );
+		fprintf( outputFile, "%.15g %.15g %.15g ", coord[0], coord[1], coord[2] );
 		FeVariable_GetValueAtNode( self, lNode_I, variableValues );
 		for ( dof_I = 0; dof_I < dofAtEachNodeCount; dof_I++ ) {
-			fprintf( outputFile, "%g ", variableValues[dof_I] );
+			fprintf( outputFile, "%.15g ", variableValues[dof_I] );
 		}	
 		fprintf( outputFile, "\n" );
 	}