OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Use 1e-8 as residual tolerance similar to istl. Added verbosity mode.

Browse Source
This commit is contained in:
Markus Blatt 2013-03-18 09:57:13 +01:00
parent d86b4c09a4
commit bf04993820
2 changed files with 8 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
opm/core/linalg/LinearSolverAGMG.cpp
View File

@ -69,12 +69,15 @@ namespace Opm
LinearSolverAGMG::LinearSolverAGMG(const parameter::ParameterGroup& param) LinearSolverAGMG::LinearSolverAGMG(const parameter::ParameterGroup& param)
: max_it_(100) , : max_it_(100) ,
rtol_ (1.0e-6), rtol_ (1.0e-8),
is_spd_(false) is_spd_(false),
linsolver_verbosity_(0)
{ {
max_it_ = param.getDefault("max_it", max_it_); max_it_ = param.getDefault("max_it", max_it_);
rtol_ = param.getDefault("rtol" , rtol_ ); rtol_ = param.getDefault("rtol" , rtol_ );
is_spd_ = param.getDefault("is_spd", is_spd_); is_spd_ = param.getDefault("is_spd", is_spd_);
linsolver_verbosity_ = param.getDefault("linsolver_verbosity",
linsolver_verbosity_);
} }
LinearSolverAGMG::~LinearSolverAGMG() {} LinearSolverAGMG::~LinearSolverAGMG() {}
@ -121,7 +124,7 @@ namespace Opm
nrest = 10; // Suggested default number of GCR restarts. nrest = 10; // Suggested default number of GCR restarts.
} }
int iprint = 0; // Suppress most output int iprint = linsolver_verbosity_; // Suppress most output
DAGMG_(& size, & a[0], & j[0], & i[0], rhs, solution, DAGMG_(& size, & a[0], & j[0], & i[0], rhs, solution,
& ijob, & iprint, & nrest, & rpt.iterations, & rtol_); & ijob, & iprint, & nrest, & rpt.iterations, & rtol_);

3
opm/core/linalg/LinearSolverAGMG.hpp
View File

@ -60,7 +60,7 @@ namespace Opm
* systems are solved using GCR. * systems are solved using GCR.
*/ */
LinearSolverAGMG(const int max_it = 100 , LinearSolverAGMG(const int max_it = 100 ,
const double rtol = 1.0e-6, const double rtol = 1.0e-8,
const bool is_spd = false); const bool is_spd = false);
/** /**
@ -104,6 +104,7 @@ namespace Opm
int max_it_; int max_it_;
double rtol_ ; double rtol_ ;
bool is_spd_; bool is_spd_;
int linsolver_verbosity_;
}; };