OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-02 05:49:09 -06:00
Code Issues Packages Projects Releases Wiki Activity

Update ecltransmissibility_impl.hh

Browse Source 
Update the comparison。
This commit is contained in:
locture 2024-01-06 00:39:46 -06:00 committed by GitHub
parent 204ebacbf2
commit da2ce029ff
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 5 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
ebos/ecltransmissibility_impl.hh
View File

@ -318,13 +318,11 @@ update(bool global, const std::function<unsigned int(unsigned int)>& map, const
// we only need to calculate a face's transmissibility
// once...
// This comparison can be problematic for corner cases in mpi
// Test case observed elemIdx <= outsideElemIdx AND insideFaceIdx==4
// This means an assert failure in applyAllZMultipliers_() while useSmallestMultiplier==true
// Converting to comparison between insideCartElemIdx and outsideCartElemIdx fixes the problem
// if (elemIdx > outsideElemIdx)
if (insideCartElemIdx > outsideCartElemIdx)
// In a parallel run insideCartElemIdx>outsideCartElemIdx does not imply elemIdx>outsideElemIdx for
// ghost cells and we need to use the cartesian index as this will be used when applying Z multipliers
// We still need to cover the case where both cells are part of an LGR and as a consequence might have
// the same cartesian index
if (std::tie(insideCartElemIdx, elemIdx) > std::tie(outsideCartElemIdx, outsideElemIdx))
continue;
// local indices of the faces of the inside and