OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-04 13:36:57 -06:00
Code Issues Packages Projects Releases Wiki Activity

Assert refactored to include same-parent-leaf-cells

Browse Source
This commit is contained in:
Antonella Ritorto 2023-12-18 14:26:33 +01:00
parent acc08799a4
commit 1fb708a239
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
ebos/ecltransmissibility_impl.hh
View File

@ -311,8 +311,10 @@ update(bool global, const std::function<unsigned int(unsigned int)>& map, const
const auto& outsideElem = intersection.outside();
unsigned outsideElemIdx = elemMapper.index(outsideElem);
unsigned insideCartElemIdx = cartMapper_.cartesianIndex(elemIdx);
unsigned outsideCartElemIdx = cartMapper_.cartesianIndex(outsideElemIdx);
// Get the Cartesian indices of the origen cells (parent or equivalent cell on level zero), for CpGrid with LGRs.
// For genral grids and no LGRs, get the usual Cartesian Index.
unsigned insideCartElemIdx = this-> lookUpCartesianData_.template getFieldPropCartesianIdx<Grid>(elemIdx);
unsigned outsideCartElemIdx = this-> lookUpCartesianData_.template getFieldPropCartesianIdx<Grid>(outsideElemIdx);
// we only need to calculate a face's transmissibility
// once...
@ -401,7 +403,9 @@ update(bool global, const std::function<unsigned int(unsigned int)>& map, const
};
int kup = find_layer(insideCartElemIdx);
int kdown=find_layer(outsideCartElemIdx);
assert(kup != kdown);
// When a grid is a CpGrid with LGRs, insideCartElemIdx coincides with outsideCartElemIdx
// for cells on the leaf with the same parent cell on level zero.
assert((kup != kdown) || (insideCartElemIdx == outsideCartElemIdx));
if(std::abs(kup -kdown) > 1){
trans = 0.0;
}