[cig-commits] r22250 - short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/faults

[brad at geodynamics.org]

Author: brad
Date: 2013-06-13 09:54:36 -0700 (Thu, 13 Jun 2013)
New Revision: 22250

Modified:
   short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
Log:
Fixed bug (remove factor of 2 from Jacobian).

Modified: short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-06-13 16:53:51 UTC (rev 22249)
+++ short/3D/PyLith/branches/v1.7-trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-06-13 16:54:36 UTC (rev 22250)
@@ -1375,10 +1375,8 @@
       } // for
     } // for
     // Jacobian is diagonal and isotropic, so it is invariant with
-    // respect to rotation and contains one unique term.  Fault
-    // traction is equal and opposite, so the Jacobian for the change
-    // in traction with slip requires a factor of 0.5.
-    const PylithScalar jacobianShearVertex = -0.5 / (areaVertex * (1.0 / jacobianVertexN[0] + 1.0 / jacobianVertexP[0]));
+    // respect to rotation and contains one unique term.
+    const PylithScalar jacobianShearVertex = -1.0 / (areaVertex * (1.0 / jacobianVertexN[0] + 1.0 / jacobianVertexP[0]));
     
     // Get friction properties and state variables.
     _friction->retrievePropsStateVars(v_fault);