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

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-06-05 11:34:29 -0700 (Tue, 05 Jun 2007)
New Revision: 7071

Modified:
   short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
Log:
Forgot to divide by 3 for meanStrainT in calcStress.
All existing unit tests now pass.
Still need a real unit test for updateState.


Modified: short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2007-06-05 18:20:20 UTC (rev 7070)
+++ short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2007-06-05 18:34:29 UTC (rev 7071)
@@ -224,9 +224,9 @@
   } else {
     const double diag[] = { 1.0, 1.0, 1.0, 0.0, 0.0, 0.0 };
     const double meanStressTpdt = bulkmodulus * traceStrainTpdt;
-    const double meanStrainT = parameters[_MaxwellIsotropic3D::pidStrainT][0] +
-      parameters[_MaxwellIsotropic3D::pidStrainT][1] +
-      parameters[_MaxwellIsotropic3D::pidStrainT][2];
+    const double meanStrainT = (parameters[_MaxwellIsotropic3D::pidStrainT][0] +
+				parameters[_MaxwellIsotropic3D::pidStrainT][1] +
+				parameters[_MaxwellIsotropic3D::pidStrainT][2])/3.0;
 
     // The code below should probably be in a separate function since it
     // is used more than once.  I should also probably cover the possibility