[cig-commits] r16603 - short/3D/PyLith/trunk/libsrc/topology

[knepley at geodynamics.org]

Author: knepley
Date: 2010-04-30 08:56:13 -0700 (Fri, 30 Apr 2010)
New Revision: 16603

Modified:
   short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc
Log:
Fixed reordering


Modified: short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc	2010-04-30 00:28:50 UTC (rev 16602)
+++ short/3D/PyLith/trunk/libsrc/topology/ReverseCuthillMcKee.cc	2010-04-30 15:56:13 UTC (rev 16603)
@@ -33,15 +33,19 @@
 
   const ALE::Obj<SieveMesh>& sieveMesh = mesh->sieveMesh();
   assert(!sieveMesh.isNull());
+  ALE::Obj<ALE::Ordering<>::perm_type> perm = 
+    new ALE::Ordering<>::perm_type(sieveMesh->comm(), sieveMesh->debug());
   ALE::Obj<ALE::Ordering<>::perm_type> reordering = 
     new ALE::Ordering<>::perm_type(sieveMesh->comm(), sieveMesh->debug());
 
-  ALE::Ordering<>::calculateMeshReordering(sieveMesh, reordering);
+  ALE::Ordering<>::calculateMeshReordering(sieveMesh, perm, reordering);
 
-  reordering->view("REORDERING");
-  sieveMesh->view("MESH BEFORE RELABEL");
+  //perm->view("PERMUTATION");
+  //reordering->view("REORDERING");
+  //sieveMesh->view("MESH BEFORE RELABEL");
 
   sieveMesh->relabel(*reordering);
+  //sieveMesh->view("MESH AFTER RELABEL");
 
   //logger.stagePop();
 } // reorder