[cig-commits] r15226 - short/3D/PyLith/trunk/libsrc/materials

[brad at geodynamics.org]

Author: brad
Date: 2009-06-12 20:35:11 -0700 (Fri, 12 Jun 2009)
New Revision: 15226

Modified:
   short/3D/PyLith/trunk/libsrc/materials/PowerLaw3D.cc
Log:
Optimization.

Modified: short/3D/PyLith/trunk/libsrc/materials/PowerLaw3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/PowerLaw3D.cc	2009-06-13 02:09:26 UTC (rev 15225)
+++ short/3D/PyLith/trunk/libsrc/materials/PowerLaw3D.cc	2009-06-13 03:35:11 UTC (rev 15226)
@@ -1363,8 +1363,7 @@
     stateVars[s_stress+iComp] = devStressTpdt + diag[iComp] *
       (meanStressTpdt + meanStressInitial);
     devStressTau = (1.0 - alpha) * devStressT[iComp] + alpha * devStressTpdt;
-    deltaVisStrain = _dt * gammaTau * devStressTau;
-    stateVars[s_viscousStrain+iComp] = visStrainT[iComp] + deltaVisStrain;
+    stateVars[s_viscousStrain+iComp] += _dt * gammaTau * devStressTau;
   } // for
 
   _needNewJacobian = true;