[cig-commits] [commit] knepley/upgrade-petsc-interface: Mesh: Fix DMPlexCheckSkeleton() for submeshes (fd50ef8)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/upgrade-petsc-interface
Link       : https://github.com/geodynamics/pylith/compare/056b66054f792da728de46f7d51d922eb18f9188...fd50ef895510b6352f372ba3b28e6b36059d957c

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

commit fd50ef895510b6352f372ba3b28e6b36059d957c
Author: Matthew G. Knepley <knepley at gmail.com>
Date:   Tue Nov 12 14:55:35 2013 -0600

    Mesh: Fix DMPlexCheckSkeleton() for submeshes


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

fd50ef895510b6352f372ba3b28e6b36059d957c
 libsrc/pylith/topology/MeshOps.cc | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/libsrc/pylith/topology/MeshOps.cc b/libsrc/pylith/topology/MeshOps.cc
index ae786b9..cf529fb 100644
--- a/libsrc/pylith/topology/MeshOps.cc
+++ b/libsrc/pylith/topology/MeshOps.cc
@@ -90,11 +90,14 @@ pylith::topology::MeshOps::checkTopology(const Mesh& mesh)
       (cellDim == 3 && numCorners == 8)) {
     isSimplexMesh = PETSC_FALSE;
   } // if
+  DMLabel subpointMap;
+  PetscErrorCode ierr = DMPlexGetSubpointMap(dmMesh, &subpointMap);PYLITH_CHECK_ERROR(ierr);
+  PetscInt cellHeight = subpointMap ? 1 : 0;
 
   PetscErrorCode err;
   err = DMPlexCheckSymmetry(dmMesh);PYLITH_CHECK_ERROR_MSG(err, "Error in topology of mesh associated with symmetry of adjacency information.");
 
-  err = DMPlexCheckSkeleton(dmMesh, isSimplexMesh);PYLITH_CHECK_ERROR_MSG(err, "Error in topology of mesh cells.");
+  err = DMPlexCheckSkeleton(dmMesh, isSimplexMesh, cellHeight);PYLITH_CHECK_ERROR_MSG(err, "Error in topology of mesh cells.");
 } // checkTopology