[cig-commits] r12471 - short/3D/PyLith/trunk/libsrc/feassemble

[willic3 at geodynamics.org]

Author: willic3
Date: 2008-07-24 10:48:22 -0700 (Thu, 24 Jul 2008)
New Revision: 12471

Modified:
   short/3D/PyLith/trunk/libsrc/feassemble/IntegratorElasticity.cc
Log:
Added missing calls to getPropertiesCell.
This fixed the stress output problem when using initial state variables,
and should also fix potential problems with the stress output when
properties vary within a material.



Modified: short/3D/PyLith/trunk/libsrc/feassemble/IntegratorElasticity.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/feassemble/IntegratorElasticity.cc	2008-07-24 14:29:01 UTC (rev 12470)
+++ short/3D/PyLith/trunk/libsrc/feassemble/IntegratorElasticity.cc	2008-07-24 17:48:22 UTC (rev 12471)
@@ -357,6 +357,7 @@
     if (!calcStress) {
       (*field)->updatePoint(*c_iter, &totalStrain[0]);
     } else {
+      _material->getPropertiesCell(*c_iter, numQuadPts);
       const double_array& stress = _material->calcStress(totalStrain);
       (*field)->updatePoint(*c_iter, &stress[0]);	
     } // else
@@ -415,6 +416,7 @@
       (*field)->restrictPoint(*c_iter);
     memcpy(&totalStrain[0], strainVals, sizeof(double)*totalStrain.size());
 
+    _material->getPropertiesCell(*c_iter, numQuadPts);
     const double_array& stress = _material->calcStress(totalStrain);
     (*field)->updatePoint(*c_iter, &stress[0]);	
   } // for