From b27e56e3864c5a480abdf067bf94743fc0540ee5 Mon Sep 17 00:00:00 2001 From: hnil Date: Mon, 18 Mar 2019 15:34:49 +0100 Subject: [PATCH] fixed bug in scaling of rhs for cpr --- opm/autodiff/ISTLSolverEbos.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/opm/autodiff/ISTLSolverEbos.hpp b/opm/autodiff/ISTLSolverEbos.hpp index 3440956d3..b92991835 100644 --- a/opm/autodiff/ISTLSolverEbos.hpp +++ b/opm/autodiff/ISTLSolverEbos.hpp @@ -769,11 +769,11 @@ protected: } block[pressureEqnIndex] = neweq; } - BlockVector newrhs(0.0); + Scalar newrhs(0.0); for (std::size_t ii = 0; ii < brhs.size(); ii++) { newrhs += bweights[ii]*brhs[ii]; } - brhs = newrhs; + brhs[pressureEqnIndex] = newrhs; } }