OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-11 08:51:55 -06:00
Code Issues Packages Projects Releases Wiki Activity

Merge pull request #4310 from totto82/groupinfo_stage2

Browse Source 
BUGFIX we now use group_info in stage 2 and need to make sure its communicated
This commit is contained in:
Markus Blatt 2022-12-08 09:01:01 +01:00 committed by GitHub
commit b4acdedefd
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 0 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
opm/simulators/wells/BlackoilWellModel_impl.hpp
View File

@ -1053,10 +1053,6 @@ namespace Opm {
}
num_rates_to_sync = groups_to_sync.size();
}
// Since "group_info" is not used in stage2, there is no need to
// communicate rates if this is the last iteration...
if (i == (num_procs - 1))
break;
num_rates_to_sync = comm.sum(num_rates_to_sync);
if (num_rates_to_sync > 0) {
std::vector<int> group_indexes;
@ -1085,14 +1081,6 @@ namespace Opm {
group_water_rates.resize(num_rates_to_sync);
group_alq_rates.resize(num_rates_to_sync);
}
// TODO: We only need to broadcast to processors with index
// j > i since we do not use the "group_info" in stage 2. In
// this case we should use comm.send() and comm.receive()
// instead of comm.broadcast() to communicate with specific
// processes, and these processes only need to receive the
// data if they are going to check the group rates in stage1
// Another similar idea is to only communicate the rates to
// process j = i + 1
#if HAVE_MPI
EclMpiSerializer ser(comm);
ser.broadcast(i, group_indexes, group_oil_rates,