OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-01-28 17:14:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add parameters controlling amg usage in cpr preconditioner.

Browse Source 
New parameters are:
 - cpr_use_amg      (default false) if true, use AMG preconditioner for elliptic part
 - cpr_use_bicgstab (default true)  if true, use BiCGStab (else use CG) for elliptic part
This commit is contained in:
Atgeirr Flø Rasmussen 2014-09-22 10:10:11 +02:00
parent 7874f177b6
commit 5d59a97972
2 changed files with 11 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
opm/autodiff/NewtonIterationBlackoilCPR.cpp
View File

@ -107,9 +107,10 @@ namespace Opm
/// Construct a system solver.
/// \param[in] linsolver linear solver to use
NewtonIterationBlackoilCPR::NewtonIterationBlackoilCPR(const parameter::ParameterGroup& /*param*/)
NewtonIterationBlackoilCPR::NewtonIterationBlackoilCPR(const parameter::ParameterGroup& param)
{
use_amg_ = param.getDefault("cpr_use_amg", false);
use_bicgstab_ = param.getDefault("cpr_use_bicgstab", true);
}
@ -185,7 +186,7 @@ namespace Opm
// typedef Dune::SeqILU0<Mat,Vector,Vector> Preconditioner;
typedef Opm::CPRPreconditioner<Mat,Vector,Vector> Preconditioner;
const double relax = 1.0;
Preconditioner precond(istlA, istlAe, relax);
Preconditioner precond(istlA, istlAe, relax, use_amg_, use_bicgstab_);
// Construct linear solver.
const double tolerance = 1e-3;

8
opm/autodiff/NewtonIterationBlackoilCPR.hpp
View File

@ -42,7 +42,9 @@ namespace Opm
/// \param[in] param parameters controlling the behaviour of
/// the preconditioning and choice of
/// linear solvers.
/// Note: parameters currently unused.
/// Parameters:
/// cpr_use_amg (default false) if true, use AMG preconditioner for elliptic part
/// cpr_use_bicgstab (default true) if true, use BiCGStab (else use CG) for elliptic part
NewtonIterationBlackoilCPR(const parameter::ParameterGroup& param);
/// Solve the system of linear equations Ax = b, with A being the
@ -51,6 +53,10 @@ namespace Opm
/// \param[in] residual residual object containing A and b.
/// \return the solution x
virtual SolutionVector computeNewtonIncrement(const LinearisedBlackoilResidual& residual) const;
private:
bool use_amg_;
bool use_bicgstab_;
};
} // namespace Opm