OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Transmissibility::update: add faceIdxDir lambda function

Browse Source 
short loop increases readability
This commit is contained in:
Arne Morten Kvarving 2025-01-17 12:30:50 +01:00
parent 2197a1ea7c
commit e10531513a
1 changed files with 18 additions and 22 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

40
opm/simulators/flow/Transmissibility_impl.hpp
View File

@ -250,6 +250,23 @@ update(bool global, const TransUpdateQuantities update_quantities,
: 1.0 / (1.0 / half1 + 1.0 / half2);
};
auto faceIdToDir = [](int insideFaceIdx) {
switch (insideFaceIdx) {
case 0:
case 1:
return FaceDir::XPlus;
case 2:
case 3:
return FaceDir::YPlus;
break;
case 4:
case 5:
return FaceDir::ZPlus;
default:
throw std::logic_error("Could not determine a face direction");
}
};
// compute the transmissibilities for all intersections
for (const auto& elem : elements(gridView_)) {
unsigned elemIdx = elemMapper.index(elem);
@ -410,30 +427,9 @@ update(bool global, const TransUpdateQuantities update_quantities,
}
// apply the region multipliers (cf. the MULTREGT keyword)
FaceDir::DirEnum faceDir;
switch (insideFaceIdx) {
case 0:
case 1:
faceDir = FaceDir::XPlus;
break;
case 2:
case 3:
faceDir = FaceDir::YPlus;
break;
case 4:
case 5:
faceDir = FaceDir::ZPlus;
break;
default:
throw std::logic_error("Could not determine a face direction");
}
trans *= transMult.getRegionMultiplier(insideCartElemIdx,
outsideCartElemIdx,
faceDir);
faceIdToDir(insideFaceIdx));
trans_.insert_or_assign(details::isId(elemIdx, outsideElemIdx), trans);