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

using C++17 style for to reduce nesting level of for loop

Browse Source
This commit is contained in:
Kai Bao 2022-06-29 14:20:35 +02:00
parent ed78d8c99a
commit 1b3cee8088
1 changed files with 8 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
opm/simulators/wells/MultisegmentWellEval.cpp
View File

@ -1554,16 +1554,15 @@ assembleICDPressureEq(const int seg,
// top segment can not be a spiral ICD device
assert(seg != 0);
if(this->segmentSet()[seg].segmentType() == Segment::SegmentType::VALVE) {
const Valve& valve = this->segmentSet()[seg].valve();
if (valve.status() == Opm::ICDStatus::SHUT) { // we switch here to be a zero rate equation
resWell_[seg][SPres] = this->primary_variables_evaluation_[seg][GTotal].value();
duneD_[seg][seg][SPres][GTotal] = 1.;
if (const auto& segment = this->segmentSet()[seg];
(segment.segmentType() == Segment::SegmentType::VALVE) &&
(segment.valve().status() == Opm::ICDStatus::SHUT) ) { // we use a zero rate equation to handle SHUT valve
resWell_[seg][SPres] = this->primary_variables_evaluation_[seg][GTotal].value();
duneD_[seg][seg][SPres][GTotal] = 1.;
auto& ws = well_state.well(baseif_.indexOfWell());
ws.segments.pressure_drop_friction[seg] = 0.;
return;
}
auto& ws = well_state.well(baseif_.indexOfWell());
ws.segments.pressure_drop_friction[seg] = 0.;
return;
}
// the pressure equation is something like