[cig-commits] [commit] knepley/upgrade-petsc-interface: Remove reordering after adjusting topology and refinement. (e2c0b27)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/upgrade-petsc-interface
Link       : https://github.com/geodynamics/pylith/compare/c7807debaf1eb022eec4046e40c282c769a695ca...258af47c24cade65c2145ea25750776577205221

>---------------------------------------------------------------

commit e2c0b274d60467ac89cb7634b7b839ad181fe240
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Sun Nov 10 20:34:41 2013 -0800

    Remove reordering after adjusting topology and refinement.
    
    Unmixing of hybrid cell stuff not implemented, so only reorder before adjusting topology.


>---------------------------------------------------------------

e2c0b274d60467ac89cb7634b7b839ad181fe240
 pylith/topology/MeshImporter.py | 12 +++---------
 1 file changed, 3 insertions(+), 9 deletions(-)

diff --git a/pylith/topology/MeshImporter.py b/pylith/topology/MeshImporter.py
index 60868d5..0498c80 100644
--- a/pylith/topology/MeshImporter.py
+++ b/pylith/topology/MeshImporter.py
@@ -101,7 +101,7 @@ class MeshImporter(MeshGenerator):
     if self.debug:
       mesh.view("Finite-element mesh.")
 
-    # Reorder mesh (serial) :TODO: Move this after adjusting topology?
+    # Reorder mesh
     if self.reorderMesh:
       logEvent2 = "%sreorder" % self._loggingPrefix
       self._eventLogger.eventBegin(logEvent2)
@@ -134,14 +134,8 @@ class MeshImporter(MeshGenerator):
       mesh.cleanup()
       newMesh.memLoggingStage = "RefinedMesh"
 
-    # Reorder mesh again (each processor independently)
-    if self.reorderMesh:
-      self._eventLogger.eventBegin(logEvent2)
-      self._debug.log(resourceUsageString())
-      if 0 == comm.rank:
-        self._info.log("Reordering cells and vertices.")
-      ordering.reorder(mesh)
-      self._eventLogger.eventEnd(logEvent2)
+    # Can't reorder mesh again, because we do not have routine to
+    # unmix normal and hybrid cells.
 
     # Nondimensionalize mesh (coordinates of vertices).
     from pylith.topology.topology import MeshOps_nondimensionalize