[cig-commits] r14547 - long/3D/SNAC/trunk/Snac/plugins/plastic

[cstark at geodynamics.org]

Author: cstark
Date: 2009-03-30 08:39:22 -0700 (Mon, 30 Mar 2009)
New Revision: 14547

Modified:
   long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c
Log:
Minor reporting change to plastic plugin Constitutive.c

Flag tracks first instance of plastic failure (for local thread) and reports when it happens.

Helps to debug runs when plastic failure initiation is uncertain.




Modified: long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c
===================================================================
--- long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c	2009-03-30 11:22:23 UTC (rev 14546)
+++ long/3D/SNAC/trunk/Snac/plugins/plastic/Constitutive.c	2009-03-30 15:39:22 UTC (rev 14547)
@@ -88,6 +88,8 @@
 		const double		a2 = material->lambda ;
 		int                     ind=0;
 
+                static char             plasticStrainReportedFlag=0;
+
 		/* 
 		 *   Work out the plastic material properties of this element 
 		*/
@@ -179,6 +181,10 @@
 			ft = s[2] - st;
                         ind=0;
 			if( fs < 0.0f || ft > 0.0f ) {
+                                if(!plasticStrainReportedFlag) {
+                                        fprintf(stderr, "r=%d, ts=%d:  *** Plastic failure *** at (%d, %d, %d)\n",context->rank, context->timeStep, ijk[0],ijk[1],ijk[2]);
+                                        plasticStrainReportedFlag=1;
+                                }
 				/*! Failure: shear or tensile */
 				ind=1;
 				aP = sqrt( 1.0f + anphi * anphi ) + anphi;