[cig-commits] r22253 - short/3D/PyLith/trunk/libsrc/pylith/faults

[brad at geodynamics.org]

Author: brad
Date: 2013-06-13 12:20:56 -0700 (Thu, 13 Jun 2013)
New Revision: 22253

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

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-06-13 19:20:12 UTC (rev 22252)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2013-06-13 19:20:56 UTC (rev 22253)
@@ -1200,10 +1200,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 / jacobianArray[jnoff+0] + 1.0 / jacobianArray[jpoff+0]));
+    // respect to rotation and contains one unique term.
+    const PylithScalar jacobianShearVertex = -1.0 / (areaVertex * (1.0 / jacobianArray[jnoff+0] + 1.0 / jacobianArray[jpoff+0]));
     
     // Get friction properties and state variables.
     _friction->retrievePropsStateVars(v_fault);