diff --git a/dumux/common/parameters.hh b/dumux/common/parameters.hh index 8e6e2637d7678dc9758630d635b2e720a57ba157..a3949028bc37a566498f0f45a89f545cf8e676a9 100644 --- a/dumux/common/parameters.hh +++ b/dumux/common/parameters.hh @@ -301,14 +301,6 @@ private: // parameters in the assembly group defaultParams["Assembly.NumericDifferenceMethod"] = "1"; - // parameters in the linear solver group - defaultParams["LinearSolver.GMResRestart"] = "10"; - defaultParams["LinearSolver.MaxIterations"] = "250"; - defaultParams["LinearSolver.PreconditionerIterations"] = "1"; - defaultParams["LinearSolver.PreconditionerRelaxation"] = "1.0"; - defaultParams["LinearSolver.ResidualReduction"] = "1e-13"; - defaultParams["LinearSolver.Verbosity"] = "0"; - // parameters in the problem group defaultParams["Problem.EnableGravity"] = "true"; defaultParams["Problem.EnableInertiaTerms"] = "true"; diff --git a/dumux/linear/seqsolverbackend.hh b/dumux/linear/seqsolverbackend.hh index a5c75b93513d087136afe23de54641f22228ddac..c52132bec4514986dd3bc5ed4527c680c606bf11 100644 --- a/dumux/linear/seqsolverbackend.hh +++ b/dumux/linear/seqsolverbackend.hh @@ -88,7 +88,7 @@ public: const std::string& modelParamGroup = "") { // get the restart threshold - const int restartGMRes = getParamFromGroup<double>(modelParamGroup, "LinearSolver.GMResRestart"); + const int restartGMRes = getParamFromGroup<int>(modelParamGroup, "LinearSolver.GMResRestart", 10); Preconditioner precond(A, s.precondIter(), s.relaxation()); @@ -131,7 +131,7 @@ public: const std::string& modelParamGroup = "") { // get the restart threshold - const int restartGMRes = getParamFromGroup<int>(modelParamGroup, "LinearSolver.GMResRestart"); + const int restartGMRes = getParamFromGroup<int>(modelParamGroup, "LinearSolver.GMResRestart", 10); Preconditioner precond(A, s.relaxation()); diff --git a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh index 2bec024454d46e4aafcc770d0567650c5f66c33b..4160138ba7ceb1a62381727359771ea4f90d4716 100644 --- a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh +++ b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressure.hh @@ -490,7 +490,7 @@ void MimeticPressure2P<TypeTag>::solve() { using Solver = GetPropType<TypeTag, Properties::LinearSolver>; - auto verboseLevelSolver = getParam<int>("LinearSolver.Verbosity"); + auto verboseLevelSolver = getParam<int>("LinearSolver.Verbosity", 0); if (verboseLevelSolver) std::cout << "MimeticPressure2P: solve for pressure" << std::endl; diff --git a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh index 3641886bd40b2891c594c9121bd7312973f89a07..a774ed7b1019b4cfc1a516a35195eca3b6a0f887 100644 --- a/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh +++ b/dumux/porousmediumflow/2p/sequential/diffusion/mimetic/pressureadaptive.hh @@ -504,7 +504,7 @@ void MimeticPressure2PAdaptive<TypeTag>::solve() { using Solver = GetPropType<TypeTag, Properties::LinearSolver>; - int verboseLevelSolver = getParam<int>("LinearSolver.Verbosity"); + int verboseLevelSolver = getParam<int>("LinearSolver.Verbosity", 0); if (verboseLevelSolver) std::cout << "MimeticPressure2PAdaptive: solve for pressure" << std::endl; diff --git a/dumux/porousmediumflow/sequential/cellcentered/pressure.hh b/dumux/porousmediumflow/sequential/cellcentered/pressure.hh index 12eacf2fea302f25272e4662c2a3e6af27645500..9d9df3d381132c344d87252d471237d536d6b57c 100644 --- a/dumux/porousmediumflow/sequential/cellcentered/pressure.hh +++ b/dumux/porousmediumflow/sequential/cellcentered/pressure.hh @@ -527,7 +527,7 @@ void FVPressure<TypeTag>::solve() { using Solver = GetPropType<TypeTag, Properties::LinearSolver>; - int verboseLevelSolver = getParam<int>("LinearSolver.Verbosity"); + int verboseLevelSolver = getParam<int>("LinearSolver.Verbosity", 0); if (verboseLevelSolver) std::cout << __FILE__ << ": solve for pressure" << std::endl;