Lgrs neighbors considered in removeSmallNonCart

This commit is contained in:
Antonella Ritorto 2023-11-13 15:51:03 +01:00
parent e33d27ca90
commit 33f8d3e7bc

View File

@ -628,7 +628,9 @@ removeSmallNonCartesianTransmissibilities_()
int gc2 = std::max(cartMapper_.cartesianIndex(elements.first), cartMapper_.cartesianIndex(elements.second));
// only adjust the NNCs
if (gc2 - gc1 == 1 || gc2 - gc1 == cartDims[0] || gc2 - gc1 == cartDims[0]*cartDims[1])
// When LGRs, all neighbors in the LGR are cartesian neighbours on the level grid representing the LGR.
// When elements on the leaf grid view have the same parent cell, gc1 and gc2 coincide.
if (gc2 - gc1 == 1 || gc2 - gc1 == cartDims[0] || gc2 - gc1 == cartDims[0]*cartDims[1] || gc2 - gc1 == 0)
continue;
//remove transmissibilities less than the threshold (by default 1e-6 in the deck's unit system)