From 9eac2883f1e2552d3e16fa4baf60c967de4dceec Mon Sep 17 00:00:00 2001 From: Pavel Tomin Date: Sun, 1 Sep 2024 23:10:47 -0500 Subject: [PATCH] fix: Small fix (#3308) * Fix bug in SinglePhasePoromechanicsConformingFractures.cpp --------- --- .../SinglePhasePoromechanicsConformingFractures.cpp | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/coreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp b/src/coreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp index 97bfe909e3a..32f5980986e 100644 --- a/src/coreComponents/physicsSolvers/multiphysics/SinglePhasePoromechanicsConformingFractures.cpp +++ b/src/coreComponents/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,