OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-22 15:33:29 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #5704 from bska/guaranteed-actionx-welpi-well-order

Browse Source 
Simplify Logic for Exchanging PI/II Values in ACTIONX/WELPI
This commit is contained in:
Bård Skaflestad 2024-11-01 11:53:55 +01:00 committed by GitHub
commit cdc96ac19e
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 7 additions and 28 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

35
opm/simulators/flow/ActionHandler.cpp
View File

@ -151,38 +151,17 @@ namespace {
return wellpi;
}
const auto num_wells = schedule[reportStep].well_order().size();
std::vector<Scalar> wellpi_vector(num_wells);
for (const auto& wname : wellpi_wells) {
if (wellModel.hasWell(wname)) {
const auto& well = schedule.getWell(wname, reportStep);
wellpi_vector[well.seqIndex()] = wellModel.wellPI(wname);
auto wellPI = std::vector<Scalar>(wellpi_wells.size());
for (auto i = 0*wellpi_wells.size(); i < wellpi_wells.size(); ++i) {
if (wellModel.hasWell(wellpi_wells[i])) {
wellPI[i] = wellModel.wellPI(wellpi_wells[i]);
}
}
if (comm.size() > 1) {
std::vector<Scalar> wellpi_buffer(num_wells * comm.size());
comm.gather(wellpi_vector.data(), wellpi_buffer.data(), num_wells, 0);
comm.max(wellPI.data(), wellPI.size());
if (comm.rank() == 0) {
for (int rank = 1; rank < comm.size(); ++rank) {
for (std::size_t well_index = 0; well_index < num_wells; ++well_index) {
const auto global_index = rank*num_wells + well_index;
const auto value = wellpi_buffer[global_index];
if (value != Scalar{0}) {
wellpi_vector[well_index] = value;
}
}
}
}
comm.broadcast(wellpi_vector.data(), wellpi_vector.size(), 0);
}
for (const auto& wname : wellpi_wells) {
const auto& well = schedule.getWell(wname, reportStep);
wellpi.insert_or_assign(wname, wellpi_vector[well.seqIndex()]);
for (auto i = 0*wellpi_wells.size(); i < wellpi_wells.size(); ++i) {
wellpi.emplace(wellpi_wells[i], wellPI[i]);
}
return wellpi;