Template Numerical Library version\ main:ab3532f7
|
Iterative solver of linear systems based on the IDR(s) method. More...
#include <TNL/Solvers/Linear/IDRs.h>
Public Types | |
using | ConstVectorViewType = typename Base::ConstVectorViewType |
Type for constant vector view. | |
using | DeviceType = typename Base::DeviceType |
Device where the solver will run on and auxillary data will alloacted on. | |
using | IndexType = typename Base::IndexType |
Type for indexing. | |
using | RealType = typename Base::RealType |
Floating point type used for computations. | |
using | VectorViewType = typename Base::VectorViewType |
Type for vector view. | |
Public Types inherited from TNL::Solvers::Linear::LinearSolver< Matrix > | |
using | ConstVectorViewType = typename Traits< Matrix >::ConstVectorViewType |
Type for constant vector view. | |
using | DeviceType = typename Matrix::DeviceType |
Device where the solver will run on and auxillary data will alloacted on. | |
using | IndexType = typename Matrix::IndexType |
Type for indexing. | |
using | MatrixPointer = std::shared_ptr< std::add_const_t< MatrixType > > |
Type of shared pointer to the matrix. | |
using | MatrixType = Matrix |
Type of the matrix representing the linear system. | |
using | PreconditionerPointer = std::shared_ptr< std::add_const_t< PreconditionerType > > |
Type of shared pointer to the preconditioner. | |
using | PreconditionerType = Preconditioners::Preconditioner< MatrixType > |
Type of preconditioner. | |
using | RealType = typename Matrix::RealType |
Floating point type used for computations. | |
using | VectorViewType = typename Traits< Matrix >::VectorViewType |
Type for vector view. | |
Public Types inherited from TNL::Solvers::IterativeSolver< Matrix::RealType, Matrix::IndexType > | |
using | SolverMonitorType |
Type of an object used for monitoring of the convergence. | |
Public Member Functions | |
void | setResidualSmoothing (bool smoothing) |
Enables or disables the residual smoothing procedure. | |
void | setShadowSpaceDimension (int s) |
Sets the dimension of the shadow space. | |
bool | setup (const Config::ParameterContainer ¶meters, const String &prefix="") override |
Method for setup of the linear iterative solver based on configuration parameters. | |
bool | solve (ConstVectorViewType b, VectorViewType x) override |
Method for solving of a linear system. | |
Public Member Functions inherited from TNL::Solvers::Linear::LinearSolver< Matrix > | |
virtual | ~LinearSolver () |
Default destructor. | |
void | setMatrix (const MatrixPointer &matrix) |
Set the matrix of the linear system. | |
void | setPreconditioner (const PreconditionerPointer &preconditioner) |
Set the preconditioner. | |
Public Member Functions inherited from TNL::Solvers::IterativeSolver< Matrix::RealType, Matrix::IndexType > | |
IterativeSolver ()=default | |
Default constructor. | |
bool | checkConvergence () |
Checks whether the convergence occurred already. | |
bool | checkNextIteration () |
Checks if the solver is allowed to do the next iteration. | |
const Matrix::RealType & | getConvergenceResidue () const |
Gets the the convergence threshold. | |
const Matrix::RealType & | getDivergenceResidue () const |
Gets the limit for the divergence criterion. | |
const Matrix::IndexType & | getIterations () const |
Gets the number of iterations performed by the solver so far. | |
const Matrix::IndexType & | getMaxIterations () const |
Gets the maximal number of iterations the solver is allowed to perform. | |
const Matrix::IndexType & | getMinIterations () const |
Gets the minimal number of iterations the solver is supposed to do. | |
const Matrix::RealType & | getResidue () const |
Gets the residue reached at the current iteration. | |
bool | nextIteration () |
Proceeds to the next iteration. | |
void | resetIterations () |
Sets the the number of the current iterations to zero. | |
void | setConvergenceResidue (const Matrix::RealType &convergenceResidue) |
Sets the threshold for the convergence. | |
void | setDivergenceResidue (const Matrix::RealType &divergenceResidue) |
Sets the residue limit for the divergence criterion. | |
void | setMaxIterations (const Matrix::IndexType &maxIterations) |
Sets the maximal number of iterations the solver is allowed to perform. | |
void | setMinIterations (const Matrix::IndexType &minIterations) |
Sets the minimal number of iterations the solver is supposed to do. | |
void | setRefreshRate (const Matrix::IndexType &refreshRate) |
Sets the refresh rate (in milliseconds) for the solver monitor. | |
void | setResidue (const Matrix::RealType &residue) |
Sets the residue reached at the current iteration. | |
void | setSolverMonitor (SolverMonitorType &solverMonitor) |
Sets the solver monitor object. | |
bool | setup (const Config::ParameterContainer ¶meters, const std::string &prefix="") |
Method for setup of the iterative solver based on configuration parameters. | |
Static Public Member Functions | |
static void | configSetup (Config::ConfigDescription &config, const String &prefix="") |
This is method defines configuration entries for setup of the linear iterative solver. | |
Static Public Member Functions inherited from TNL::Solvers::Linear::LinearSolver< Matrix > | |
static void | configSetup (Config::ConfigDescription &config, const String &prefix="") |
This method defines configuration entries for setup of the linear iterative solver. | |
Static Public Member Functions inherited from TNL::Solvers::IterativeSolver< Matrix::RealType, Matrix::IndexType > | |
static void | configSetup (Config::ConfigDescription &config, const std::string &prefix="") |
This method defines configuration entries for setup of the iterative solver. | |
Protected Types | |
using | DeviceVector = Containers::Vector< RealType, DeviceType, IndexType > |
using | HostMatrix = Matrices::DenseMatrix< RealType, Devices::Sequential, IndexType > |
using | HostVector = Containers::Vector< RealType, Devices::Host, IndexType > |
using | VectorType = typename Traits::VectorType |
Protected Member Functions | |
void | matvec (ConstVectorViewType src, VectorViewType dst) |
void | psolve (ConstVectorViewType src, VectorViewType dst) |
void | setSize (const VectorViewType &x) |
Protected Attributes | |
HostVector | f |
DeviceVector | G |
HostMatrix | M |
DeviceVector | P |
VectorType | r |
VectorType | r_s |
int | s = 4 |
IndexType | size = 0 |
IndexType | sizeWithGhosts = 0 |
bool | smoothing = false |
VectorType | t |
DeviceVector | U |
VectorType | v |
VectorType | x_s |
Protected Attributes inherited from TNL::Solvers::Linear::LinearSolver< Matrix > | |
MatrixPointer | matrix = nullptr |
PreconditionerPointer | preconditioner = nullptr |
Protected Attributes inherited from TNL::Solvers::IterativeSolver< Matrix::RealType, Matrix::IndexType > | |
Matrix::RealType | convergenceResidue |
Matrix::IndexType | currentIteration |
Matrix::RealType | currentResidue |
Matrix::RealType | divergenceResidue |
Matrix::IndexType | maxIterations |
Matrix::IndexType | minIterations |
Matrix::IndexType | refreshRate |
std::ofstream | residualHistoryFile |
std::string | residualHistoryFileName |
IterativeSolverMonitor< Matrix::RealType, Matrix::IndexType > * | solverMonitor |
Iterative solver of linear systems based on the IDR(s) method.
IDRs implements an iterative solver for non-symmetric linear systems, using the Induced Dimension Reduction algorithm, denoted as IDR(s), according to the description in [1] and [2].
[1] Peter Sonneveld and Martin Van Gijzen. "IDR(s): A Family of Simple and Fast Algorithms for Solving Large Nonsymmetric Systems of Linear Equations", SIAM Journal on Scientific Computing 31.2 (2009): 1035-1062. [2] Martin Van Gijzen and Peter Sonneveld. "Algorithm 913: An elegant IDR(s) variant that efficiently exploits biorthogonality properties." ACM Transactions on Mathematical Software (TOMS) 38.1 (2011): 1-19.
See TNL::Solvers::Linear::LinearSolver for example of showing how to use the linear solvers.
Matrix | is type of matrix describing the linear system. |
|
static |
This is method defines configuration entries for setup of the linear iterative solver.
In addition to config entries defined by IterativeSolver::configSetup, this method defines the following:
bicgstab-ell - number of Bi-CG iterations before the MR part starts.
bicgstab-exact-residue - says whether the BiCGstab should compute the exact residue in each step (true) or to use a cheap approximation (false).
config | contains description of configuration parameters. |
prefix | is a prefix of particular configuration entries. |
|
overridevirtual |
Method for setup of the linear iterative solver based on configuration parameters.
parameters | contains values of the define configuration entries. |
prefix | is a prefix of particular configuration entries. |
Reimplemented from TNL::Solvers::Linear::LinearSolver< Matrix >.
|
overridevirtual |
Method for solving of a linear system.
See LinearSolver::solve for more details.
b | vector with the right-hand side of the linear system. |
x | vector for the solution of the linear system. |
Implements TNL::Solvers::Linear::LinearSolver< Matrix >.