[cig-commits] r16479 - in short/3D/PyLith/trunk: . libsrc/topology

[brad at geodynamics.org]

Author: brad
Date: 2010-03-31 14:26:12 -0700 (Wed, 31 Mar 2010)
New Revision: 16479

Modified:
   short/3D/PyLith/trunk/TODO
   short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc
Log:
Added debugging output for mesh reordering.

Modified: short/3D/PyLith/trunk/TODO
===================================================================
--- short/3D/PyLith/trunk/TODO	2010-03-31 21:25:48 UTC (rev 16478)
+++ short/3D/PyLith/trunk/TODO	2010-03-31 21:26:12 UTC (rev 16479)
@@ -5,9 +5,8 @@
 MAIN PRIORITIES
 
 * Mesh reordering
+  debug
 
-  Make optional in MeshImporter.py.
-
 * Better preconditioning
   + Need settings for Schur complement
 
@@ -40,11 +39,6 @@
 
 * Uniform global refinement for tets with faults
 
-* Use rate fields.
-  + AbsorbingDampers
-  + ElasticityExplicit
-  + ElasticityExpligitLgDeform
-
 * Cleanup
   + memory model
   + full-scale testing

Modified: short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc	2010-03-31 21:25:48 UTC (rev 16478)
+++ short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc	2010-03-31 21:26:12 UTC (rev 16479)
@@ -17,6 +17,7 @@
 #include "pylith/topology/Mesh.hh" // USES Mesh
 
 #include <cassert> // USES assert()
+#include <stdexcept> // USES std::exception
 
 // ----------------------------------------------------------------------
 typedef pylith::topology::Mesh::SieveMesh SieveMesh;
@@ -36,6 +37,10 @@
     new ALE::Ordering<>::perm_type(sieveMesh->comm(), sieveMesh->debug());
 
   ALE::Ordering<>::calculateMeshReordering(sieveMesh, reordering);
+
+  reordering->view("REORDERING");
+  sieveMesh->view("MESH BEFORE RELABEL");
+
   sieveMesh->relabel(*reordering);
 
   //logger.stagePop();