From 77f330ee8aa20ed970f136d0a2e388e6b3b02210 Mon Sep 17 00:00:00 2001
From: Simon Emmert <simon.emmert@iws.uni-stuttgart.de>
Date: Wed, 29 Apr 2020 12:02:26 +0200
Subject: [PATCH] [2p][test] use ILU0 with GMRes as Linear Solver for 2p tests
Needed to ensure convergence with minimal setup (no UMFPack) for
test_2p_incompressible_tpfa_oilwet. The coarse grid solver is changed to
SSOR when UMFPack is not available, which did not converge as well.
---
test/porousmediumflow/2p/implicit/incompressible/main.cc | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/test/porousmediumflow/2p/implicit/incompressible/main.cc b/test/porousmediumflow/2p/implicit/incompressible/main.cc
index ce4d1b633c..59d9b0ab49 100644
--- a/test/porousmediumflow/2p/implicit/incompressible/main.cc
+++ b/test/porousmediumflow/2p/implicit/incompressible/main.cc
@@ -40,7 +40,7 @@
#include <dumux/common/valgrind.hh>
#include <dumux/common/dumuxmessage.hh>
-#include <dumux/linear/amgbackend.hh>
+#include <dumux/linear/seqsolverbackend.hh>
#include <dumux/linear/linearsolvertraits.hh>
#include <dumux/nonlinear/newtonsolver.hh>
@@ -181,8 +181,8 @@ int main(int argc, char** argv) try
auto assembler = std::make_shared<Assembler>(problem, gridGeometry, gridVariables, timeLoop, xOld);
// the linear solver
- using LinearSolver = AMGBiCGSTABBackend<LinearSolverTraits<GridGeometry>>;
- auto linearSolver = std::make_shared<LinearSolver>(leafGridView, gridGeometry->dofMapper());
+ using LinearSolver = ILU0RestartedGMResBackend;
+ auto linearSolver = std::make_shared<LinearSolver>();
// the non-linear solver
using NewtonSolver = Dumux::NewtonSolver<Assembler, LinearSolver>;
--
GitLab