[cig-commits] [commit] knepley/fix-faults-parallel: Fault: Small fix for fault halo (d3610af)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/pylith

On branch  : knepley/fix-faults-parallel
Link       : https://github.com/geodynamics/pylith/compare/80029c389e447a372a70929ff189ad24b6fbf5b4...422d933ebf5e19a69537cf35cab1c4497c03fe90

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

commit d3610af3972f293388849911dfbd544eba3bd8b5
Author: Matthew G. Knepley <knepley at gmail.com>
Date:   Sat Apr 19 14:09:32 2014 -0500

    Fault: Small fix for fault halo


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

d3610af3972f293388849911dfbd544eba3bd8b5
 libsrc/pylith/faults/FaultCohesiveLagrange.cc | 1 +
 1 file changed, 1 insertion(+)

diff --git a/libsrc/pylith/faults/FaultCohesiveLagrange.cc b/libsrc/pylith/faults/FaultCohesiveLagrange.cc
index 2fb25f2..8b07bbf 100644
--- a/libsrc/pylith/faults/FaultCohesiveLagrange.cc
+++ b/libsrc/pylith/faults/FaultCohesiveLagrange.cc
@@ -1039,6 +1039,7 @@ pylith::faults::FaultCohesiveLagrange::checkConstraints(const topology::Field& s
   PetscInt fiberDim = 0, numConstraints = 0;
   for(int iVertex = 0; iVertex < numVertices; ++iVertex) {
 
+    if (_cohesiveVertices[iVertex].lagrange < 0) continue;
     const int v_negative = _cohesiveVertices[iVertex].negative;
     const int e_lagrange = _cohesiveVertices[iVertex].lagrange;