[cig-commits] r16799 - in short/3D/PyLith/trunk: libsrc/faults unittests/libtests/faults

[brad at geodynamics.org]

Author: brad
Date: 2010-05-26 20:22:15 -0700 (Wed, 26 May 2010)
New Revision: 16799

Modified:
   short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
Log:
Updated unit tests.

Modified: short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-05-27 02:53:11 UTC (rev 16798)
+++ short/3D/PyLith/trunk/libsrc/faults/FaultCohesiveLagrange.cc	2010-05-27 03:22:15 UTC (rev 16799)
@@ -1724,8 +1724,6 @@
     assemblyWtSection->updatePoint(v_fault, &wt);
   } // for
 
-  assemblyWtSection->view("ASSEMBLY WEIGHT");
-
 #if 0 // DEBUGGING
   area.view("AREA");
   //_faultMesh->getSendOverlap()->view("Send fault overlap");

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc	2010-05-27 02:53:11 UTC (rev 16798)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveKin.cc	2010-05-27 03:22:15 UTC (rev 16799)
@@ -364,7 +364,7 @@
     fault.useSolnIncr(false);
     fault.integrateResidual(residual, t, &fields);
 
-    residual.view("RESIDUAL"); // DEBUGGING
+    //residual.view("RESIDUAL"); // DEBUGGING
 
     // Check values
     const double* valsE = _data->residual;
@@ -395,7 +395,7 @@
     fault.useSolnIncr(true);
     fault.integrateResidual(residual, t, &fields);
 
-    residual.view("RESIDUAL"); // DEBUGGING
+    //residual.view("RESIDUAL"); // DEBUGGING
 
     // Check values
     const double* valsE = _data->residualIncr;
@@ -722,8 +722,11 @@
   } // setup disp
   jacobian.complete();
 
-  const topology::Field<topology::Mesh>& solution = fields.get("dispIncr(t->t+dt)");
+  topology::Field<topology::Mesh>& solution = fields.get("dispIncr(t->t+dt)");
   fault.adjustSolnLumped(&fields, jacobian);
+  const topology::Field<topology::Mesh>& dispIncrAdj = 
+    fields.get("dispIncr adjust");
+  solution += dispIncrAdj;
 
   //solution.view("SOLUTION AFTER ADJUSTMENT"); // DEBUGGING
 
@@ -1104,6 +1107,7 @@
   fields->add("residual", "residual");
   fields->add("disp(t)", "displacement");
   fields->add("dispIncr(t->t+dt)", "displacement_increment");
+  fields->add("dispIncr adjust", "displacement_adjust");
   fields->solutionName("dispIncr(t->t+dt)");
   
   const int spaceDim = _data->spaceDim;