[cig-commits] [commit] knepley/fix-faults-parallel: Fix bug in checking topology for submeshes. (ce91c33)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/fix-faults-parallel
Link       : https://github.com/geodynamics/pylith/compare/96a97b54c2e3426d62617703be018cb08f5f4dbb...ce4ce38e9b8191d09bdccb7e13f6d3f74dd68a69

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

commit ce91c339a2799829235c01b95076dda150ac94cf
Author: Brad Aagaard <baagaard at usgs.gov>
Date:   Mon Apr 28 15:04:52 2014 -0700

    Fix bug in checking topology for submeshes.


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

ce91c339a2799829235c01b95076dda150ac94cf
 libsrc/pylith/topology/Mesh.icc | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/libsrc/pylith/topology/Mesh.icc b/libsrc/pylith/topology/Mesh.icc
index 5a89202..8af2816 100644
--- a/libsrc/pylith/topology/Mesh.icc
+++ b/libsrc/pylith/topology/Mesh.icc
@@ -128,7 +128,8 @@ pylith::topology::Mesh::numCorners(void) const {
   if (_dmMesh) {
     PetscInt cStart, cEnd, vStart, vEnd, closureSize, *closure = NULL;
     PetscErrorCode err;
-    err = DMPlexGetHeightStratum(_dmMesh, 0, &cStart, &cEnd);PYLITH_CHECK_ERROR(err);
+    const int cellHeight = _isSubMesh ? 1 : 0;
+    err = DMPlexGetHeightStratum(_dmMesh, cellHeight, &cStart, &cEnd);PYLITH_CHECK_ERROR(err);
     err = DMPlexGetDepthStratum(_dmMesh, 0, &vStart, &vEnd);PYLITH_CHECK_ERROR(err);
     if (cEnd > cStart) {
       err = DMPlexGetTransitiveClosure(_dmMesh, cStart, PETSC_TRUE, &closureSize, &closure);PYLITH_CHECK_ERROR(err);