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

Add communication when multiplying with the matrix duneB_ in apply

Browse Source 
Here we go from cells to segments, and everything concerning segments is stored globally.
This commit is contained in:
Lisa Julia Nebel 2024-10-16 14:44:11 +02:00
parent 6784db2d1d
commit 81a5da8b63
1 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/simulators/wells/MultisegmentWellEquations.cpp
View File

@ -152,6 +152,11 @@ apply(const BVector& x, BVector& Ax) const
duneB_.mv(x, Bx); duneB_.mv(x, Bx);
if (this->pw_info_.communication().size() == 1) {
// We need to communicate here to get the contributions from all segments
this->pw_info_.communication().sum(Bx.data(), Bx.size());
}
// invDBx = duneD^-1 * Bx_ // invDBx = duneD^-1 * Bx_
const BVectorWell invDBx = mswellhelpers::applyUMFPack(*duneDSolver_, Bx); const BVectorWell invDBx = mswellhelpers::applyUMFPack(*duneDSolver_, Bx);