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

Merge pull request #755 from blattms/istl-consistent-rhs

Browse Source 
Makes right hand side passed to linear solver consistent.
This commit is contained in:
Atgeirr Flø Rasmussen 2015-02-16 12:48:28 +01:00
commit 2f7a87aa4b
1 changed files with 2 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
opm/core/linalg/LinearSolverIstl.cpp
View File

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