[cig-commits] r21947 - short/3D/PyLith/trunk/unittests/libtests/meshio

[brad at geodynamics.org]

Author: brad
Date: 2013-04-25 12:28:21 -0700 (Thu, 25 Apr 2013)
New Revision: 21947

Modified:
   short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
Log:
Debuggin output.

Modified: short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc
===================================================================
--- short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc	2013-04-25 19:16:47 UTC (rev 21946)
+++ short/3D/PyLith/trunk/unittests/libtests/meshio/TestDataWriterFaultMesh.cc	2013-04-25 19:28:21 UTC (rev 21947)
@@ -106,6 +106,7 @@
   PYLITH_METHOD_END;
 } // _initialize
 
+#include "pylith/topology/CoordsVisitor.hh" // DEBUGGING
 // ----------------------------------------------------------------------
 // Create vertex fields.
 void
@@ -124,6 +125,11 @@
   const PetscInt vStart = verticesStratum.begin();
   const PetscInt vEnd = verticesStratum.end();
 
+#if 1 // DEBUGGING
+  topology::CoordsVisitor coordsVisitor(dmMesh);
+  const PetscScalar* coordsArray = coordsVisitor.localArray();
+#endif
+
   // Set vertex fields
   for (int i=0; i < nfields; ++i) {
     const char* name = _data->vertexFieldsInfo[i].name;
@@ -143,6 +149,16 @@
       for(PetscInt d = 0; d < fiberDim; ++d, ++index) {
 	fieldArray[off+d] = _data->vertexFields[i][index];
       } // for
+
+#if 1 // DEBUGGING
+      const PetscInt coff = coordsVisitor.sectionOffset(v);
+      std::cout << "VERTEX " << v;
+      for(PetscInt d = 0; d < coordsVisitor.sectionDof(v); ++d) {
+	std::cout << " " << coordsArray[off+d];
+      } // for
+      std::cout << std::endl;
+#endif
+
     } // for
     CPPUNIT_ASSERT_EQUAL(_data->numVertices, vEnd-vStart);
   } // for