[cig-commits] r14245 - long/3D/SNAC/trunk/Snac/snac2vtk

[cstark at geodynamics.org]

Author: cstark
Date: 2009-03-07 14:10:09 -0800 (Sat, 07 Mar 2009)
New Revision: 14245

Modified:
   long/3D/SNAC/trunk/Snac/snac2vtk/snac2vtk.c
Log:
Testing for NaNs in stress elements added in snac2vtk in DeriveStressMeasures().  

Error is reported if any element (of six) is NaN and then program bails.



Modified: long/3D/SNAC/trunk/Snac/snac2vtk/snac2vtk.c
===================================================================
--- long/3D/SNAC/trunk/Snac/snac2vtk/snac2vtk.c	2009-03-07 21:40:52 UTC (rev 14244)
+++ long/3D/SNAC/trunk/Snac/snac2vtk/snac2vtk.c	2009-03-07 22:10:09 UTC (rev 14245)
@@ -1015,6 +1015,13 @@
 	float	stressTensorArray[3][3];
 	fread( stressTensorArray, sizeof(float), 9, stressTensorIn );
 	/*
+	 *  Report error and bail if we pick up NaNs in any of the stress components
+	 */
+	if(isnan(stressTensorArray[tetra_I][0][0]) || isnan(stressTensorArray[tetra_I][1][1]) 
+	   || isnan(stressTensorArray[tetra_I][2][2]) || isnan(stressTensorArray[tetra_I][0][1]) 
+	   || isnan(stressTensorArray[tetra_I][0][2]) || isnan(stressTensorArray[tetra_I][1][2])) 
+	    fprintf(stderr,"NaN in stress tensor file\n");
+	/*
 	 *  Build average stress tensor for element by summing tetrahedral tensor components
 	 *   - even though it's symmetric, do for all 9 components in case we pick the wrong ones before diagonalization
 	 */