[cig-commits] [commit] knepley/fix-faults-parallel: Fix bug in checking topology for submeshes. (ce91c33)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Mon May 5 03:05:42 PDT 2014
Repository : ssh://geoshell/pylith
On branch : knepley/fix-faults-parallel
Link : https://github.com/geodynamics/pylith/compare/80029c389e447a372a70929ff189ad24b6fbf5b4...422d933ebf5e19a69537cf35cab1c4497c03fe90
>---------------------------------------------------------------
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);
More information about the CIG-COMMITS
mailing list