From 33f8d3e7bcff30d700b6acd8409e198ac9b5fc88 Mon Sep 17 00:00:00 2001 From: Antonella Ritorto Date: Mon, 13 Nov 2023 15:51:03 +0100 Subject: [PATCH] Lgrs neighbors considered in removeSmallNonCart --- ebos/ecltransmissibility_impl.hh | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/ebos/ecltransmissibility_impl.hh b/ebos/ecltransmissibility_impl.hh index 2d9580631..21fd6f5b8 100644 --- a/ebos/ecltransmissibility_impl.hh +++ b/ebos/ecltransmissibility_impl.hh @@ -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)