[cig-commits] [commit] knepley/fix-faults-parallel: Faults: Fix the fault boundary label (01ad665)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/fix-faults-parallel
Link       : https://github.com/geodynamics/pylith/compare/cba3ba8659a67c89fc6b016a201455b52b4bdbd4...01ad6651de6adc20322c52870ee78cc04843d655

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

commit 01ad6651de6adc20322c52870ee78cc04843d655
Author: Matthew G. Knepley <knepley at gmail.com>
Date:   Tue Apr 29 21:00:21 2014 -0500

    Faults: Fix the fault boundary label


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

01ad6651de6adc20322c52870ee78cc04843d655
 libsrc/pylith/faults/CohesiveTopology.cc | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/libsrc/pylith/faults/CohesiveTopology.cc b/libsrc/pylith/faults/CohesiveTopology.cc
index 4930d0b..59ce5a1 100644
--- a/libsrc/pylith/faults/CohesiveTopology.cc
+++ b/libsrc/pylith/faults/CohesiveTopology.cc
@@ -27,6 +27,8 @@
 
 #include "pylith/utils/error.h" // USES PYLITH_CHECK_ERROR
 
+extern "C" PetscErrorCode DMPlexMarkBoundaryFaces_Internal(DM, PetscInt, DMLabel);
+
 // ----------------------------------------------------------------------
 void
 pylith::faults::CohesiveTopology::createFault(topology::Mesh* faultMesh,
@@ -59,7 +61,7 @@ pylith::faults::CohesiveTopology::createFault(topology::Mesh* faultMesh,
 
     err = DMPlexCreateLabel(subdm, labelName);PYLITH_CHECK_ERROR(err);
     err = DMPlexGetLabel(subdm, labelName, &label);PYLITH_CHECK_ERROR(err);
-    err = DMPlexMarkBoundaryFaces(subdm, label);PYLITH_CHECK_ERROR(err);
+    err = DMPlexMarkBoundaryFaces_Internal(subdm, 1, label);PYLITH_CHECK_ERROR(err);
     err = DMPlexLabelComplete(subdm, label);PYLITH_CHECK_ERROR(err);
     err = DMPlexCreateSubmesh(subdm, label, 1, &faultBoundary);PYLITH_CHECK_ERROR(err);
     std::string submeshLabel = "fault_" + std::string(groupName);