[cig-commits] r20335 - short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/faults

[brad at geodynamics.org]

Author: brad
Date: 2012-06-08 12:45:53 -0700 (Fri, 08 Jun 2012)
New Revision: 20335

Modified:
   short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc
Log:
Updated unit tests to account for changes in friction formulation (updated solution with dispIncr adjust).

Modified: short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc	2012-06-08 19:35:16 UTC (rev 20334)
+++ short/3D/PyLith/branches/v1.7-trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc	2012-06-08 19:45:53 UTC (rev 20335)
@@ -233,6 +233,11 @@
   const PylithScalar dt = 0.01;
   fault.timeStep(dt);
   fault.constrainSolnSpace(&fields, t, jacobian);
+  
+  topology::Field<topology::Mesh>& solution = fields.solution();
+  const topology::Field<topology::Mesh>& dispIncrAdj = fields.get("dispIncr adjust");
+  solution += dispIncrAdj;
+
   fault.updateStateVars(t, &fields);
 
   { // Check solution values
@@ -352,6 +357,11 @@
   const PylithScalar dt = 0.01;
   fault.timeStep(dt);
   fault.constrainSolnSpace(&fields, t, jacobian);
+
+  topology::Field<topology::Mesh>& solution = fields.solution();
+  const topology::Field<topology::Mesh>& dispIncrAdj = fields.get("dispIncr adjust");
+  solution += dispIncrAdj;
+
   fault.updateStateVars(t, &fields);
 
   { // Check solution values
@@ -477,6 +487,11 @@
   const PylithScalar dt = 0.01;
   fault.timeStep(dt);
   fault.constrainSolnSpace(&fields, t, jacobian);
+
+  topology::Field<topology::Mesh>& solution = fields.solution();
+  const topology::Field<topology::Mesh>& dispIncrAdj = fields.get("dispIncr adjust");
+  solution += dispIncrAdj;
+
   fault.updateStateVars(t, &fields);
 
   //residual.view("RESIDUAL"); // DEBUGGING