Do not filter connections on the schedule of the loadbalanced grid.

It would remove perforated cells from wells that cross the local
domain's border. That would make it impossible to figure out the
first connection. In addition we would not be able to check that
the connections exist (as rank 0 would have the complete information
-> inconsistency).
This commit is contained in:
Markus Blatt 2020-11-11 18:30:59 +01:00
parent 307315a03d
commit 46cb9012c0

View File

@ -234,15 +234,10 @@ public:
cartesianIndexMapper_.reset(); cartesianIndexMapper_.reset();
if ( ! equilGrid_ ) // Calling Schedule::filterConnectins would remove any perforated
{ // cells that exist on other ranks in the case of distributed wells
// for processes that do not hold the global grid we filter here using the local grid. // But we need them to figure out the first cell of a well (e.g. for
// If we would filter in filterConnection_ our partition would be empty and the connections of all // pressure). Hence this is now skipped. Rank 0 had everything even before.
// wells would be removed.
ActiveGridCells activeCells(grid().logicalCartesianSize(),
grid().globalCell().data(), grid().size(0));
this->schedule().filterConnections(activeCells);
}
} }
#endif #endif