From 11f9eb9d88b62036fff1d44d3e04d61375ae8afe Mon Sep 17 00:00:00 2001 From: Markus Blatt Date: Mon, 7 Sep 2020 16:44:17 +0200 Subject: [PATCH] Base face processing decision of faces on cartesian indices This will process the same faces in serial and parallel. Hence it make the silent assumption that we only process faces from low cartesian index to high cartesian index hold again. This should fix PINCH MULTZ ALL in parallel. --- ebos/ecltransmissibility.hh | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/ebos/ecltransmissibility.hh b/ebos/ecltransmissibility.hh index 9af88dfb1..374434c4f 100644 --- a/ebos/ecltransmissibility.hh +++ b/ebos/ecltransmissibility.hh @@ -273,14 +273,14 @@ public: A * (n*d)/(d*d); } - // we only need to calculate a face's transmissibility - // once... - if (elemIdx > outsideElemIdx) - continue; - unsigned insideCartElemIdx = cartMapper.cartesianIndex(elemIdx); unsigned outsideCartElemIdx = cartMapper.cartesianIndex(outsideElemIdx); + // we only need to calculate a face's transmissibility + // once... + if (insideCartElemIdx > outsideCartElemIdx) + continue; + // local indices of the faces of the inside and // outside elements which contain the intersection int insideFaceIdx = intersection.indexInInside();