fix: Small fix (#3308)

* Fix bug in SinglePhasePoromechanicsConformingFractures.cpp ---------
GEOS-DEV · Sep 17, 2024 · 9eac288 · 9eac288
1 parent dde1f49
commit 9eac288
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/...oreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp b/...oreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp
@@ -409,7 +409,7 @@ addTransmissibilityCouplingPattern( DomainPartition const & domain,
             // to the nodes of the adjacent cell
             localIndex const rowIndex = presDofNumber[sei[iconn][1-kf]] - rankOffset;
 
-            if( rowIndex > 0 && rowIndex < pattern.numRows() )
+            if( rowIndex >= 0 && rowIndex < pattern.numRows() )
             {
 
               // Get fracture, face and region/subregion/element indices (for elements on both sides)
@@ -655,7 +655,7 @@ assembleFluidMassResidualDerivativeWrtDisplacement( MeshLevel const & mesh,
 
         localIndex const localRow = LvArray::integerConversion< localIndex >( elemDOF[0] - rankOffset );
 
-        if( localRow > 0 && localRow < localMatrix.numRows() )
+        if( localRow >= 0 && localRow < localMatrix.numRows() )
         {
 
           localMatrix.addToRowBinarySearchUnsorted< serialAtomic >( localRow,
@@ -711,7 +711,7 @@ assembleFluidMassResidualDerivativeWrtDisplacement( MeshLevel const & mesh,
         {
           localIndex const localRow = LvArray::integerConversion< localIndex >( elemDOF[0] - rankOffset );
 
-          if( localRow > 0 && localRow < localMatrix.numRows() )
+          if( localRow >= 0 && localRow < localMatrix.numRows() )
           {
             localMatrix.addToRowBinarySearchUnsorted< serialAtomic >( localRow,
                                                                       nodeDOF,