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

PI/II Loop: Check That We Process All Active Connections

Browse Source 
Assert that the 'subsetPerfID' is equal to number_of_perforations_
at the end of the loop.

Suggested by: [at]atgeirr
This commit is contained in:
Bård Skaflestad 2020-11-06 21:15:25 +01:00
parent a46a732f9e
commit c8a43b982c
2 changed files with 10 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
opm/simulators/wells/MultisegmentWell_impl.hpp
View File

@ -1224,7 +1224,8 @@ namespace Opm
const auto& allConn = this->well_ecl_.getConnections();
const auto nPerf = allConn.size();
for (auto allPerfID = 0*nPerf, subsetPerfID = 0*nPerf; allPerfID < nPerf; ++allPerfID) {
auto subsetPerfID = 0*nPerf;
for (auto allPerfID = 0*nPerf; allPerfID < nPerf; ++allPerfID) {
if (allConn[allPerfID].state() == Connection::State::SHUT) {
continue;
}
@ -1259,6 +1260,9 @@ namespace Opm
++subsetPerfID;
connPI += np;
}
assert (static_cast<int>(subsetPerfID) == this->number_of_perforations_ &&
"Internal logic error in processing connections for PI/II");
}

6
opm/simulators/wells/StandardWell_impl.hpp
View File

@ -2330,7 +2330,8 @@ namespace Opm
const auto& allConn = this->well_ecl_.getConnections();
const auto nPerf = allConn.size();
for (auto allPerfID = 0*nPerf, subsetPerfID = 0*nPerf; allPerfID < nPerf; ++allPerfID) {
auto subsetPerfID = 0*nPerf;
for (auto allPerfID = 0*nPerf; allPerfID < nPerf; ++allPerfID) {
if (allConn[allPerfID].state() == Connection::State::SHUT) {
continue;
}
@ -2363,6 +2364,9 @@ namespace Opm
++subsetPerfID;
connPI += np;
}
assert (static_cast<int>(subsetPerfID) == this->number_of_perforations_ &&
"Internal logic error in processing connections for PI/II");
}