[cig-commits] r17113 - short/3D/PyLith/trunk/libsrc/friction

[brad at geodynamics.org]

Author: brad
Date: 2010-08-23 16:53:05 -0700 (Mon, 23 Aug 2010)
New Revision: 17113

Modified:
   short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc
   short/3D/PyLith/trunk/libsrc/friction/TimeWeakening.cc
Log:
Switched from zero tolerance for healing to very small value to account for roundoff error in residual calculation.

Modified: short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc	2010-08-23 23:52:04 UTC (rev 17112)
+++ short/3D/PyLith/trunk/libsrc/friction/SlipWeakening.cc	2010-08-23 23:53:05 UTC (rev 17113)
@@ -304,7 +304,8 @@
   assert(0 != numStateVars);
   assert(0 != numProperties);
 
-  if (slipRate != 0.0) {
+  const double tolerance = 1.0e-12;
+  if (slipRate > tolerance) {
     const double slipPrev = stateVars[s_slipPrev];
 
     stateVars[s_slipPrev] = stateVars[s_slipCum];

Modified: short/3D/PyLith/trunk/libsrc/friction/TimeWeakening.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/friction/TimeWeakening.cc	2010-08-23 23:52:04 UTC (rev 17112)
+++ short/3D/PyLith/trunk/libsrc/friction/TimeWeakening.cc	2010-08-23 23:53:05 UTC (rev 17113)
@@ -294,7 +294,8 @@
   assert(0 != numStateVars);
   assert(0 != numProperties);
 
-  if (slipRate != 0.0) {
+  const double tolerance = 1.0e-12;
+  if (slipRate > tolerance) {
     const double dt = _dt;
 
     stateVars[s_time] += dt;