[cig-commits] r16484 - short/3D/PyLith/trunk/unittests/libtests/faults

[brad at geodynamics.org]

Author: brad
Date: 2010-04-01 12:53:17 -0700 (Thu, 01 Apr 2010)
New Revision: 16484

Modified:
   short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc
Log:
More debugging output.

Modified: short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc	2010-04-01 19:52:46 UTC (rev 16483)
+++ short/3D/PyLith/trunk/unittests/libtests/faults/TestFaultCohesiveDyn.cc	2010-04-01 19:53:17 UTC (rev 16484)
@@ -358,8 +358,9 @@
         const int index = iVertex * spaceDim + i;
         const double valE = valsE[index];
 #if 1 // DEBUGGING
-	std::cout << "valE: " << valE
+	std::cout << "SOLUTION valE: " << valE
 		  << ", val: " << vals[i]
+		  << ", error: " << fabs(1.0-vals[i]/valE)
 		  << std::endl;
 #endif // DEBUGGING
 	if (fabs(valE) > tolerance)
@@ -390,6 +391,8 @@
       fault._fields->get("slip").section();
     CPPUNIT_ASSERT(!slipSection.isNull());
 
+    slipSection->view("SLIP");
+
     // Get expected values
     const double* valsE = _data->slipSlipE;
     int iVertex = 0; // variable to use as index into valsE array
@@ -409,8 +412,9 @@
         const int index = iVertex * spaceDim + i;
         const double valE = valsE[index];
 #if 1 // DEBUGGING
-	std::cout << "valE: " << valE
+	std::cout << "SLIP valE: " << valE
 		  << ", val: " << vals[i]
+		  << ", error: " << fabs(1.0-vals[i]/valE)
 		  << std::endl;
 #endif // DEBUGGING
         if (fabs(valE) > tolerance)