From 6f83ae4febe437efc971d711063ebaf65a791078 Mon Sep 17 00:00:00 2001
From: Max Kahnt <maxka@mi.fu-berlin.de>
Date: Wed, 23 Aug 2017 12:21:59 +0000
Subject: [PATCH] Disable regularization per default to avoid unintentional
 usage.

Thanks to Lasse for pointing that out.
---
 dune/solvers/iterationsteps/blockgssteps.hh | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/dune/solvers/iterationsteps/blockgssteps.hh b/dune/solvers/iterationsteps/blockgssteps.hh
index 76cae5f9..eec9dde1 100644
--- a/dune/solvers/iterationsteps/blockgssteps.hh
+++ b/dune/solvers/iterationsteps/blockgssteps.hh
@@ -208,7 +208,9 @@ auto truncateSymmetrically(LinearSolver&& linearSolver) {
 
 namespace LocalSolverRegularizer {
 
-static constexpr double defaultDiagRegularizeParameter = 1e-10;
+// a suitable value for defaultDiagRagularization could be 1e-10
+// nevertheless, the default is 0 to avoid unintentional regularization
+static constexpr double defaultDiagRegularizeParameter = 0;
 
 // Note: move-capture, auto-return and auto-lambda-arguments require C++14
 /**
-- 
GitLab