From a980993289ebd690b5e494367309e4d959f5f338 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Atgeirr=20Fl=C3=B8=20Rasmussen?= Date: Wed, 3 May 2017 14:29:09 +0200 Subject: [PATCH] Make reordering solver work with iterated implicit. Also tighten tolerances. They were too slack, giving significant differences from non-reordering implementation towards the end of SPE1CASE2 for example. --- opm/autodiff/BlackoilReorderingTransportModel.hpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/opm/autodiff/BlackoilReorderingTransportModel.hpp b/opm/autodiff/BlackoilReorderingTransportModel.hpp index 339c9b34f..700822536 100644 --- a/opm/autodiff/BlackoilReorderingTransportModel.hpp +++ b/opm/autodiff/BlackoilReorderingTransportModel.hpp @@ -324,7 +324,7 @@ namespace Opm { { auto rs = reservoir_state; auto ws = well_state; - tr_model_.nonlinearIteration(iteration, timer, nonlinear_solver, rs, ws); + tr_model_.nonlinearIteration(/*iteration*/ 0, timer, nonlinear_solver, rs, ws); } // Create report and exit. @@ -830,7 +830,7 @@ namespace Opm { bool getConvergence(const int cell, const Vec2& res) { - const double tol = 1e-5; + const double tol = 1e-7; // Compute scaled residuals (scaled like saturations). double sres[] = { res[0] / (cstate_[cell].b[Oil] * Base::pvdt_[cell]), res[1] / (cstate_[cell].b[Gas] * Base::pvdt_[cell]) };