[cig-commits] r19269 - short/3D/PyLith/branches/v1.6-stable/unittests/libtests/faults

[brad at geodynamics.org]

Author: brad
Date: 2011-12-05 13:36:04 -0800 (Mon, 05 Dec 2011)
New Revision: 19269

Modified:
   short/3D/PyLith/branches/v1.6-stable/unittests/libtests/faults/TestFaultCohesiveDyn.cc
Log:
Added updating state vars to constrain solution space test. Updates slip values even if they aren't changed in constraining the solution space.

Modified: short/3D/PyLith/branches/v1.6-stable/unittests/libtests/faults/TestFaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/branches/v1.6-stable/unittests/libtests/faults/TestFaultCohesiveDyn.cc	2011-12-05 21:14:19 UTC (rev 19268)
+++ short/3D/PyLith/branches/v1.6-stable/unittests/libtests/faults/TestFaultCohesiveDyn.cc	2011-12-05 21:36:04 UTC (rev 19269)
@@ -216,6 +216,7 @@
   const double dt = 0.01;
   fault.timeStep(dt);
   fault.constrainSolnSpace(&fields, t, jacobian);
+  fault.updateStateVars(t, &fields);
 
   { // Check solution values
     // No change to Lagrange multipliers for stick case.
@@ -334,6 +335,7 @@
   const double dt = 0.01;
   fault.timeStep(dt);
   fault.constrainSolnSpace(&fields, t, jacobian);
+  fault.updateStateVars(t, &fields);
 
   { // Check solution values
     // Lagrange multipliers should be adjusted according to friction
@@ -457,6 +459,7 @@
   const double dt = 0.01;
   fault.timeStep(dt);
   fault.constrainSolnSpace(&fields, t, jacobian);
+  fault.updateStateVars(t, &fields);
 
   //residual.view("RESIDUAL"); // DEBUGGING