[cig-commits] r19533 - in short/3D/PyLith/trunk: libsrc/pylith/faults libsrc/pylith/friction unittests/libtests/friction unittests/libtests/friction/data

[brad at geodynamics.org]

Author: brad
Date: 2012-01-31 14:55:07 -0800 (Tue, 31 Jan 2012)
New Revision: 19533

Modified:
   short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
   short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakening.cc
   short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakeningTime.cc
   short/3D/PyLith/trunk/unittests/libtests/friction/TestFrictionModel.cc
   short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningData.cc
   short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningTimeData.cc
Log:
Manual patch (also pushed to stable) correcting updating of state variables for slip weakening friction.

Modified: short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2012-01-31 20:38:09 UTC (rev 19532)
+++ short/3D/PyLith/trunk/libsrc/pylith/faults/FaultCohesiveDyn.cc	2012-01-31 22:55:07 UTC (rev 19533)
@@ -981,7 +981,7 @@
       dDispTIncrVertexP[iDim] = +0.5*dDispRelVertex[iDim];
     } // for
 
-#if 1 // debugging
+#if 0 // debugging
     std::cout << "v_fault: " << v_fault;
     std::cout << ", tractionTpdtVertex: ";
     for (int iDim=0; iDim < spaceDim; ++iDim)

Modified: short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakening.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakening.cc	2012-01-31 20:38:09 UTC (rev 19532)
+++ short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakening.cc	2012-01-31 22:55:07 UTC (rev 19533)
@@ -323,7 +323,7 @@
   if (slipRate > tolerance && !_forceHealing) {
     const PylithScalar slipPrev = stateVars[s_slipPrev];
 
-    stateVars[s_slipPrev] = stateVars[s_slipCum];
+    stateVars[s_slipPrev] = slip;
     stateVars[s_slipCum] += fabs(slip - slipPrev);
   } else {
     // Sliding has stopped, so reset state variables.

Modified: short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakeningTime.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakeningTime.cc	2012-01-31 20:38:09 UTC (rev 19532)
+++ short/3D/PyLith/trunk/libsrc/pylith/friction/SlipWeakeningTime.cc	2012-01-31 22:55:07 UTC (rev 19533)
@@ -332,7 +332,7 @@
   if (slipRate > tolerance) {
     const PylithScalar slipPrev = stateVars[s_slipPrev];
 
-    stateVars[s_slipPrev] = stateVars[s_slipCum];
+    stateVars[s_slipPrev] = slip;
     stateVars[s_slipCum] += fabs(slip - slipPrev);
   } else {
     // Sliding has stopped, so reset state variables.

Modified: short/3D/PyLith/trunk/unittests/libtests/friction/TestFrictionModel.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/friction/TestFrictionModel.cc	2012-01-31 20:38:09 UTC (rev 19532)
+++ short/3D/PyLith/trunk/unittests/libtests/friction/TestFrictionModel.cc	2012-01-31 22:55:07 UTC (rev 19533)
@@ -360,7 +360,7 @@
 
     const int numStateVars = 2;
     const PylithScalar stateVars[2] = { 0.5, 0.1 };
-    const PylithScalar stateVarsUpdatedE[2] = { 0.65, 0.5 };
+    const PylithScalar stateVarsUpdatedE[2] = { 0.65, 0.25 };
     
 
     // Set state variables to given values

Modified: short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningData.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningData.cc	2012-01-31 20:38:09 UTC (rev 19532)
+++ short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningData.cc	2012-01-31 22:55:07 UTC (rev 19533)
@@ -142,9 +142,9 @@
 
 const PylithScalar pylith::friction::SlipWeakeningData::_stateVarsUpdated[] = {
   0.48,
-  0.4,
+  0.12,
   0.65,
-  0.5,
+  0.25,
 };
 
 pylith::friction::SlipWeakeningData::SlipWeakeningData(void)

Modified: short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningTimeData.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningTimeData.cc	2012-01-31 20:38:09 UTC (rev 19532)
+++ short/3D/PyLith/trunk/unittests/libtests/friction/data/SlipWeakeningTimeData.cc	2012-01-31 22:55:07 UTC (rev 19533)
@@ -150,9 +150,9 @@
 
 const PylithScalar pylith::friction::SlipWeakeningTimeData::_stateVarsUpdated[] = {
   0.48,
-  0.4,
+  0.12,
   0.65,
-  0.5,
+  0.25,
 };
 
 pylith::friction::SlipWeakeningTimeData::SlipWeakeningTimeData(void)