[cig-commits] r15303 - short/3D/PyLith/trunk/templates/materials

[brad at geodynamics.org]

Author: brad
Date: 2009-06-16 13:36:37 -0700 (Tue, 16 Jun 2009)
New Revision: 15303

Modified:
   short/3D/PyLith/trunk/templates/materials/PlaneStrainState.cc
Log:
Fixed small bug. Template for bulk material has state variables.

Modified: short/3D/PyLith/trunk/templates/materials/PlaneStrainState.cc
===================================================================
--- short/3D/PyLith/trunk/templates/materials/PlaneStrainState.cc	2009-06-16 19:10:37 UTC (rev 15302)
+++ short/3D/PyLith/trunk/templates/materials/PlaneStrainState.cc	2009-06-16 20:36:37 UTC (rev 15303)
@@ -258,7 +258,8 @@
   assert(0 != density);
   assert(0 != properties);
   assert(_numPropsQuadPt == numProperties);
-  assert(0 == numStateVars);
+  assert(0 != stateVars);
+  assert(_numVarsQuadPt == numStateVars);
 
   // Set density using physical properties (trivial since one our
   // physical properties is density).
@@ -287,7 +288,8 @@
   assert(_PlaneStrainState::tensorSize == stressSize);
   assert(0 != properties);
   assert(_numPropsQuadPt == numProperties);
-  assert(0 == numStateVars);
+  assert(0 != stateVars);
+  assert(_numVarsQuadPt == numStateVars);
   assert(0 != totalStrain);
   assert(_PlaneStrainState::tensorSize == strainSize);
   assert(0 != initialStress);
@@ -341,7 +343,8 @@
   assert(_PlaneStrainState::numElasticConsts == numElasticConsts);
   assert(0 != properties);
   assert(_numPropsQuadPt == numProperties);
-  assert(0 == numStateVars);
+  assert(0 != stateVars);
+  assert(_numVarsQuadPt == numStateVars);
   assert(0 != totalStrain);
   assert(_PlaneStrainState::tensorSize == strainSize);
   assert(0 != initialStress);