diff --git a/opm/core/linalg/LinearSolverIstl.cpp b/opm/core/linalg/LinearSolverIstl.cpp index 1cf66f35d..eb3a84578 100644 --- a/opm/core/linalg/LinearSolverIstl.cpp +++ b/opm/core/linalg/LinearSolverIstl.cpp @@ -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.copyOwnerToAll(b,b); // System solution Vector x(opA.getmat().M()); x = 0.0;