Makes right hand side passed to linear solver consistent.

Due to the size of the overlap layer and the discretization scheme
the rhs might not contain correct values for overlap cells. This
commit makes sure they are correct by an additional communication step.
This commit is contained in:
Markus Blatt 2015-02-12 21:00:32 +01:00
parent d6bf9d3a36
commit a31fae4507

View File

@ -195,6 +195,8 @@ namespace Opm
// System RHS
Vector b(opA.getmat().N());
std::copy(rhs, rhs+b.size(), b.begin());
// Make rhs consistent in the parallel case
comm.copyOwnerToOverlap(b,b);
// System solution
Vector x(opA.getmat().M());
x = 0.0;