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 e59db97128
commit 87f9cf8c77
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)
: max_it_(100) ,
rtol_ (1.0e-6),
is_spd_(false)
rtol_ (1.0e-8),
is_spd_(false),
linsolver_verbosity_(0)
{
max_it_ = param.getDefault("max_it", max_it_);
rtol_ = param.getDefault("rtol" , rtol_ );
is_spd_ = param.getDefault("is_spd", is_spd_);
linsolver_verbosity_ = param.getDefault("linsolver_verbosity",
linsolver_verbosity_);
}
LinearSolverAGMG::~LinearSolverAGMG() {}
@ -121,7 +124,7 @@ namespace Opm
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,
& 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.
*/
LinearSolverAGMG(const int max_it = 100 ,
const double rtol = 1.0e-6,
const double rtol = 1.0e-8,
const bool is_spd = false);
/**
@ -104,6 +104,7 @@ namespace Opm
int max_it_;
double rtol_ ;
bool is_spd_;
int linsolver_verbosity_;
};