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

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-06-04 13:15:03 -0700 (Mon, 04 Jun 2007)
New Revision: 7058

Modified:
   short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
Log:
I was using the wrong viscous multiplication factor for the Jacobian.


Modified: short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2007-06-04 20:13:18 UTC (rev 7057)
+++ short/3D/PyLith/trunk/libsrc/materials/MaxwellIsotropic3D.cc	2007-06-04 20:15:03 UTC (rev 7058)
@@ -327,8 +327,8 @@
     } else
       dq = maxwelltime*(1.0-exp(-_dt/maxwelltime))/_dt;
     const double visFac = mu*dq/3.0;
-    (*elasticConsts)[ 0] = bulkmodulus + 4.0*dq; // C1111
-    (*elasticConsts)[ 1] = bulkmodulus - 2.0*dq; // C1122
+    (*elasticConsts)[ 0] = bulkmodulus + 4.0*visFac; // C1111
+    (*elasticConsts)[ 1] = bulkmodulus - 2.0*visFac; // C1122
     (*elasticConsts)[ 2] = (*elasticConsts)[1]; // C1133
     (*elasticConsts)[ 3] = 0; // C1112
     (*elasticConsts)[ 4] = 0; // C1123