[cig-commits] [commit] knepley/upgrade-petsc-interface: Faults: Added a check for hybrid edges in verifyConfiguration() (c7807de)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

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

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

commit c7807debaf1eb022eec4046e40c282c769a695ca
Author: Matthew G. Knepley <knepley at gmail.com>
Date:   Sun Nov 10 15:34:48 2013 -0600

    Faults: Added a check for hybrid edges in verifyConfiguration()


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

c7807debaf1eb022eec4046e40c282c769a695ca
 libsrc/pylith/faults/FaultCohesiveLagrange.cc | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/libsrc/pylith/faults/FaultCohesiveLagrange.cc b/libsrc/pylith/faults/FaultCohesiveLagrange.cc
index cb318f6..6043de8 100644
--- a/libsrc/pylith/faults/FaultCohesiveLagrange.cc
+++ b/libsrc/pylith/faults/FaultCohesiveLagrange.cc
@@ -904,6 +904,11 @@ pylith::faults::FaultCohesiveLagrange::verifyConfiguration(const topology::Mesh&
   PetscInt eMax;
 
   err = DMPlexGetHybridBounds(dmMesh, NULL, NULL, &eMax, NULL);PYLITH_CHECK_ERROR(err);
+  if (eMax < 0) {
+    std::ostringstream msg;
+    msg << "No hybrid edges found in mesh.";
+    throw std::runtime_error(msg.str());
+  } // if  
 
   // Check for fault groups
   PetscBool hasLabel;