DUNE-FEM (unstable)

Iterative Solvers Template Library (ISTL)

Iterative Solvers supporting block recursive matrix and vector classes at compile time. More...

Modules

 Iterative Solvers
 
 Communication Interface
 
 Sparse Matrix and Vector classes
 Matrix and Vector classes that support a block recursive structure capable of representing the natural structure from Finite Element discretisations.
 

Files

file  ldl.hh
 Class for using LDL with ISTL matrices.
 
file  spqr.hh
 Class for using SPQR with ISTL matrices.
 
file  superlu.hh
 Classes for using SuperLU with ISTL matrices.
 
file  umfpack.hh
 Classes for using UMFPack with ISTL matrices.
 

Classes

class  Dune::ILUSubdomainSolver< M, X, Y >
 base class encapsulating common algorithms of ILU0SubdomainSolver and ILUNSubdomainSolver. More...
 
class  Dune::ILU0SubdomainSolver< M, X, Y >
 Exact subdomain solver using ILU(p) with appropriate p. More...
 
class  Dune::ISTLError
 derive error class from the base class in common More...
 
class  Dune::BCRSMatrixError
 Error specific to BCRSMatrix. More...
 
class  Dune::ImplicitModeCompressionBufferExhausted
 Thrown when the compression buffer used by the implicit BCRSMatrix construction is exhausted. More...
 
class  Dune::SolverAbort
 Thrown when a solver aborts due to some problem. More...
 
class  Dune::MatrixBlockError
 Error when performing an operation on a matrix block. More...
 
class  Dune::SeqOverlappingSchwarz< M, X, TM, TD, TA >
 Sequential overlapping Schwarz preconditioner. More...
 
class  Dune::LDL< Matrix >
 Use the LDL package to directly solve linear systems – empty default class. More...
 
class  Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >
 The LDL direct sparse solver for matrices of type BCRSMatrix. More...
 
class  Dune::SPQR< Matrix >
 Use the SPQR package to directly solve linear systems – empty default class. More...
 
class  Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >
 The SPQR direct sparse solver for matrices of type BCRSMatrix. More...
 
class  Dune::SuperLU< M >
 SuperLu Solver. More...
 
class  Dune::UMFPack< M >
 The UMFPack direct sparse solver. More...
 

Typedefs

typedef Dune::BCRSMatrix< FieldMatrix< T, n, m >, A > Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::Matrix
 The matrix type.
 
typedef Dune::ISTL::Impl::BCCSMatrix< T, int > Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDLMatrix
 The corresponding SuperLU Matrix type.
 
typedef ISTL::Impl::BCCSMatrixInitializer< BCRSMatrix< FieldMatrix< T, n, m >, A >, int > Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::MatrixInitializer
 Type of an associated initializer class.
 
typedef Dune::BlockVector< FieldVector< T, m >, typename std::allocator_traits< A >::template rebind_alloc< FieldVector< T, m > > > Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::domain_type
 The type of the domain of the solver.
 
typedef Dune::BlockVector< FieldVector< T, n >, typename std::allocator_traits< A >::template rebind_alloc< FieldVector< T, n > > > Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::range_type
 The type of the range of the solver.
 
typedef Dune::BCRSMatrix< FieldMatrix< T, n, m >, A > Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::Matrix
 The matrix type.
 
typedef ISTL::Impl::BCCSMatrix< T, int > Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQRMatrix
 The corresponding SuperLU Matrix type.
 
typedef ISTL::Impl::BCCSMatrixInitializer< BCRSMatrix< FieldMatrix< T, n, m >, A >, int > Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::MatrixInitializer
 Type of an associated initializer class.
 
typedef Dune::BlockVector< FieldVector< T, m >, typename std::allocator_traits< A >::template rebind_alloc< FieldVector< T, m > > > Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::domain_type
 The type of the domain of the solver.
 
typedef Dune::BlockVector< FieldVector< T, n >, typename std::allocator_traits< A >::template rebind_alloc< FieldVector< T, n > > > Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::range_type
 The type of the range of the solver.
 
using Dune::UMFPack< M >::Matrix = M
 The matrix type.
 
using Dune::UMFPack< M >::UMFPackMatrix = ISTL::Impl::BCCSMatrix< typename Matrix::field_type, size_type >
 The corresponding (scalar) UMFPack matrix type.
 
using Dune::UMFPack< M >::MatrixInitializer = ISTL::Impl::BCCSMatrixInitializer< M, size_type >
 Type of an associated initializer class.
 
using Dune::UMFPack< M >::domain_type = Impl::UMFPackDomainType< M >
 The type of the domain of the solver.
 
using Dune::UMFPack< M >::range_type = Impl::UMFPackRangeType< M >
 The type of the range of the solver.
 

Functions

SolverCategory::Category Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::category () const override
 Category of the solver (see SolverCategory::Category)
 
 Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL (const Matrix &matrix, int verbose=0)
 Construct a solver object from a BCRSMatrix. More...
 
 Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL (const Matrix &matrix, int verbose, bool)
 Constructor for compatibility with SuperLU standard constructor. More...
 
 Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL (const Matrix &matrix, const ParameterTree &config)
 Constructs the LDL solver. More...
 
 Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL ()
 Default constructor.
 
virtual Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::~LDL ()
 Default constructor.
 
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply (domain_type &x, range_type &b, InverseOperatorResult &res) override
 Apply inverse operator,. More...
 
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply (domain_type &x, range_type &b, double reduction, InverseOperatorResult &res) override
 apply inverse operator, with given convergence criteria. More...
 
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply (T *x, T *b)
 Additional apply method with c-arrays in analogy to superlu. More...
 
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::setMatrix (const Matrix &matrix)
 Initialize data from given matrix.
 
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::setVerbosity (int v)
 Sets the verbosity level for the solver. More...
 
LDLMatrixDune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getInternalMatrix ()
 Return the column compress matrix. More...
 
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::free ()
 Free allocated space. More...
 
const char * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::name ()
 Get method name.
 
double * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getD ()
 Get factorization diagonal matrix D. More...
 
int * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getLp ()
 Get factorization Lp. More...
 
int * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getLi ()
 Get factorization Li. More...
 
double * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getLx ()
 Get factorization Lx. More...
 
SolverCategory::Category Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::category () const override
 Category of the solver (see SolverCategory::Category)
 
 Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR (const Matrix &matrix, int verbose=0)
 Construct a solver object from a BCRSMatrix. More...
 
 Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR (const Matrix &matrix, int verbose, bool)
 Constructor for compatibility with SuperLU standard constructor. More...
 
 Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR (const Matrix &matrix, const ParameterTree &config)
 Constructs the SPQR solver. More...
 
 Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR ()
 Default constructor.
 
virtual Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::~SPQR ()
 Destructor.
 
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply (domain_type &x, range_type &b, InverseOperatorResult &res) override
 Apply inverse operator,. More...
 
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply (domain_type &x, range_type &b, double reduction, InverseOperatorResult &res) override
 apply inverse operator, with given convergence criteria. More...
 
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::setMatrix (const Matrix &matrix)
 Initialize data from given matrix.
 
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::setVerbosity (int v)
 Sets the verbosity level for the solver. More...
 
SuiteSparseQR_factorization< T > * Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getFactorization ()
 Return the matrix factorization. More...
 
SPQRMatrixDune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getInternalMatrix ()
 Return the column coppressed matrix. More...
 
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::free ()
 Free allocated space. More...
 
const char * Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::name ()
 Get method name.
 
SolverCategory::Category Dune::UMFPack< M >::category () const override
 Category of the solver (see SolverCategory::Category)
 
 Dune::UMFPack< M >::UMFPack (const Matrix &matrix, int verbose=0)
 Construct a solver object from a matrix. More...
 
 Dune::UMFPack< M >::UMFPack (const Matrix &matrix, int verbose, bool)
 Constructor for compatibility with SuperLU standard constructor. More...
 
 Dune::UMFPack< M >::UMFPack (const Matrix &mat_, const ParameterTree &config)
 Construct a solver object from a matrix. More...
 
 Dune::UMFPack< M >::UMFPack ()
 default constructor
 
 Dune::UMFPack< M >::UMFPack (const Matrix &mat_, const char *file, int verbose=0)
 Try loading a decomposition from file and do a decomposition if unsuccessful. More...
 
 Dune::UMFPack< M >::UMFPack (const char *file, int verbose=0)
 try loading a decomposition from file More...
 
void Dune::UMFPack< M >::apply (domain_type &x, range_type &b, InverseOperatorResult &res) override
 Apply inverse operator,. More...
 
void Dune::UMFPack< M >::apply (domain_type &x, range_type &b, double reduction, InverseOperatorResult &res) override
 apply inverse operator, with given convergence criteria. More...
 
void Dune::UMFPack< M >::apply (T *x, T *b)
 additional apply method with c-arrays in analogy to superlu More...
 
void Dune::UMFPack< M >::setOption (unsigned int option, double value)
 Set UMFPack-specific options. More...
 
void Dune::UMFPack< M >::saveDecomposition (const char *file)
 saves a decomposition to a file More...
 
template<class BitVector = Impl::NoBitVector>
void Dune::UMFPack< M >::setMatrix (const Matrix &matrix, const BitVector &bitVector={})
 Initialize data from given matrix. More...
 
void Dune::UMFPack< M >::setVerbosity (int v)
 sets the verbosity level for the UMFPack solver More...
 
void * Dune::UMFPack< M >::getFactorization ()
 Return the matrix factorization. More...
 
UMFPackMatrixDune::UMFPack< M >::getInternalMatrix ()
 Return the column compress matrix from UMFPack. More...
 
void Dune::UMFPack< M >::free ()
 free allocated space. More...
 
template<class S >
std::size_t Dune::ILUSubdomainSolver< M, X, Y >::copyToLocalMatrix (const M &A, S &rowset)
 Copy the local part of the global matrix to ILU. More...
 
template<class S >
void Dune::ILU0SubdomainSolver< M, X, Y >::setSubMatrix (const M &A, S &rowset)
 Set the data of the local problem. More...
 
template<class S >
void Dune::ILUNSubdomainSolver< M, X, Y >::setSubMatrix (const M &A, S &rowset)
 Set the data of the local problem. More...
 
void Dune::SuperLU< M >::free ()
 free allocated space. More...
 
 Dune::SuperLU< M >::SuperLU (const Matrix &mat, bool verbose=false, bool reusevector=true)
 Constructs the SuperLU solver. More...
 
 Dune::SuperLU< M >::SuperLU ()
 Empty default constructor. More...
 
void Dune::SuperLU< M >::setMatrix (const Matrix &mat)
 Initialize data from given matrix.
 
void Dune::SuperLU< M >::apply (domain_type &x, range_type &b, InverseOperatorResult &res) override
 Apply inverse operator,. More...
 
void Dune::SuperLU< M >::apply (T *x, T *b)
 Apply SuperLu to C arrays.
 

Detailed Description

Iterative Solvers supporting block recursive matrix and vector classes at compile time.

The Iterative Solver Template Library applies generic programming in C++ to the domain of iterative solvers of linear systems stemming from finite element discretizations. Those discretizations exhibit a lot of structure, e.g:

  1. Certain discretizations for systems of PDEs or higher order methods result in matrices where individual entries are replaced by small blocks, say of size \(2\times 2\) or \(4\times 4\). Dense blocks of different sizes e.g. arise in \(hp\) Discontinuous Galerkin discretization methods. It is straightforward and efficient to treat these small dense blocks as fully coupled and solve them with direct methods within the iterative method.
  2. Equation-wise ordering for systems results in matrices having an \(n\times n\) block structure where \(n\) corresponds to the number of variables in the PDE and the blocks themselves are large and sparse. As an example we mention the Stokes system.
  3. Other discretisation, e.~g. those of reaction/diffusion systems, produce sparse matrices whose blocks are sparse matrices of small dense blocks,

Our matrix and vector interface supports a block recursive structure. Each sparse matrix entry can itself be either a sparse or a small dense matrix.

The solvers use this recursive block structure via template meta programming at compile time.

ISTL consists of the matrix and vector API and the solvers which use the Preconditioners preconditioners.

Function Documentation

◆ apply() [1/9]

template<typename T , typename A , int n, int m>
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply ( domain_type x,
range_type b,
double  reduction,
InverseOperatorResult res 
)
inlineoverride

apply inverse operator, with given convergence criteria.

Warning
Right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
reductionThe minimum defect reduction to achieve.
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

◆ apply() [2/9]

template<typename T , typename A , int n, int m>
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply ( domain_type x,
range_type b,
double  reduction,
InverseOperatorResult res 
)
inlineoverride

apply inverse operator, with given convergence criteria.

Warning
Right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
reductionThe minimum defect reduction to achieve.
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

◆ apply() [3/9]

template<typename M >
void Dune::UMFPack< M >::apply ( domain_type x,
range_type b,
double  reduction,
InverseOperatorResult res 
)
inlineoverride

apply inverse operator, with given convergence criteria.

Warning
Right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
reductionThe minimum defect reduction to achieve.
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

References Dune::UMFPack< M >::apply().

◆ apply() [4/9]

template<typename T , typename A , int n, int m>
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply ( domain_type x,
range_type b,
InverseOperatorResult res 
)
inlineoverride

Apply inverse operator,.

Warning
Note: right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

References Dune::InverseOperatorResult::converged, and Dune::InverseOperatorResult::iterations.

◆ apply() [5/9]

template<typename T , typename A , int n, int m>
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply ( domain_type x,
range_type b,
InverseOperatorResult res 
)
inlineoverride

Apply inverse operator,.

Warning
Note: right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

References Dune::InverseOperatorResult::converged, and Dune::InverseOperatorResult::iterations.

◆ apply() [6/9]

template<typename M >
void Dune::SuperLU< M >::apply ( domain_type x,
range_type b,
InverseOperatorResult res 
)
override

Apply inverse operator,.

Warning
Note: right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

Referenced by Dune::SuperLU< M >::apply().

◆ apply() [7/9]

template<typename M >
void Dune::UMFPack< M >::apply ( domain_type x,
range_type b,
InverseOperatorResult res 
)
inlineoverride

Apply inverse operator,.

Warning
Note: right hand side b may be overwritten!
Parameters
xThe left hand side to store the result in.
bThe right hand side
resObject to store the statistics about applying the operator.
Exceptions
SolverAbortWhen the solver detects a problem and cannot continue

References Dune::InverseOperatorResult::converged, DUNE_THROW, Dune::InverseOperatorResult::elapsed, Dune::flatVectorForEach(), and Dune::InverseOperatorResult::iterations.

Referenced by Dune::UMFPack< M >::apply().

◆ apply() [8/9]

template<typename T , typename A , int n, int m>
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::apply ( T *  x,
T *  b 
)
inline

Additional apply method with c-arrays in analogy to superlu.

Parameters
xsolution array
brhs array

◆ apply() [9/9]

template<typename M >
void Dune::UMFPack< M >::apply ( T *  x,
T *  b 
)
inline

additional apply method with c-arrays in analogy to superlu

Parameters
xsolution array
brhs array

NOTE If the user hands over a pure pointer, we assume that they know what they are doing, hence no copy to flat structures

◆ copyToLocalMatrix()

template<class M , class X , class Y >
template<class S >
std::size_t Dune::ILUSubdomainSolver< M, X, Y >::copyToLocalMatrix ( const M &  A,
S &  rowset 
)
protected

Copy the local part of the global matrix to ILU.

Parameters
AThe global matrix.
rowsetThe global indices of the local problem.

◆ free() [1/4]

template<typename T , typename A , int n, int m>
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::free ( )
inline

Free allocated space.

Warning
Later calling apply will result in an error.

◆ free() [2/4]

template<typename T , typename A , int n, int m>
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::free ( )
inline

Free allocated space.

Warning
Later calling apply will result in an error.

◆ free() [3/4]

template<typename M >
void Dune::SuperLU< M >::free

free allocated space.

Warning
later calling apply will result in an error.

◆ free() [4/4]

template<typename M >
void Dune::UMFPack< M >::free ( )
inline

free allocated space.

Warning
later calling apply will result in an error.

◆ getD()

template<typename T , typename A , int n, int m>
double * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getD ( )
inline

Get factorization diagonal matrix D.

Warning
It is up to the user to preserve consistency.

◆ getFactorization() [1/2]

template<typename T , typename A , int n, int m>
SuiteSparseQR_factorization< T > * Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getFactorization ( )
inline

Return the matrix factorization.

Warning
It is up to the user to keep consistency.

◆ getFactorization() [2/2]

template<typename M >
void * Dune::UMFPack< M >::getFactorization ( )
inline

Return the matrix factorization.

Warning
It is up to the user to keep consistency.

◆ getInternalMatrix() [1/3]

template<typename T , typename A , int n, int m>
LDLMatrix & Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getInternalMatrix ( )
inline

Return the column compress matrix.

Warning
It is up to the user to keep consistency.

◆ getInternalMatrix() [2/3]

template<typename T , typename A , int n, int m>
SPQRMatrix & Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getInternalMatrix ( )
inline

Return the column coppressed matrix.

Warning
It is up to the user to keep consistency.

◆ getInternalMatrix() [3/3]

template<typename M >
UMFPackMatrix & Dune::UMFPack< M >::getInternalMatrix ( )
inline

Return the column compress matrix from UMFPack.

Warning
It is up to the user to keep consistency.

◆ getLi()

template<typename T , typename A , int n, int m>
int * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getLi ( )
inline

Get factorization Li.

Warning
It is up to the user to preserve consistency.

◆ getLp()

template<typename T , typename A , int n, int m>
int * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getLp ( )
inline

Get factorization Lp.

Warning
It is up to the user to preserve consistency.

◆ getLx()

template<typename T , typename A , int n, int m>
double * Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::getLx ( )
inline

Get factorization Lx.

Warning
It is up to the user to preserve consistency.

◆ LDL() [1/3]

template<typename T , typename A , int n, int m>
Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL ( const Matrix matrix,
const ParameterTree config 
)
inline

Constructs the LDL solver.

Parameters
matrixThe matrix of the system to solve.
configParameterTree containing solver parameters.
ParameterTree Key Meaning
verbose The verbosity level. default=0

◆ LDL() [2/3]

template<typename T , typename A , int n, int m>
Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL ( const Matrix matrix,
int  verbose,
bool   
)
inline

Constructor for compatibility with SuperLU standard constructor.

This computes the matrix decomposition, and may take a long time (and use a lot of memory).

Parameters
matrixthe matrix to solve for
verbose0 or 1 set the verbosity level, defaults to 0

◆ LDL() [3/3]

template<typename T , typename A , int n, int m>
Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::LDL ( const Matrix matrix,
int  verbose = 0 
)
inline

Construct a solver object from a BCRSMatrix.

This computes the matrix decomposition, and may take a long time (and use a lot of memory).

Parameters
matrixthe matrix to solve for
verbose0 or 1 set the verbosity level, defaults to 0

◆ saveDecomposition()

template<typename M >
void Dune::UMFPack< M >::saveDecomposition ( const char *  file)
inline

saves a decomposition to a file

Parameters
filethe filename to save to

References DUNE_THROW.

Referenced by Dune::UMFPack< M >::UMFPack().

◆ setMatrix()

template<typename M >
template<class BitVector = Impl::NoBitVector>
void Dune::UMFPack< M >::setMatrix ( const Matrix matrix,
const BitVector &  bitVector = {} 
)
inline

Initialize data from given matrix.

Template Parameters
BitVectora compatible bitvector to the domain_type/range_type. Defaults to NoBitVector for backwards compatibility

A positive bit indices that the corresponding matrix row/column is excluded from the UMFPACK decomposition. WARNING This is an opposite behavior of the previous implementation in setSubMatrix.

Referenced by Dune::UMFPack< M >::UMFPack().

◆ setOption()

template<typename M >
void Dune::UMFPack< M >::setOption ( unsigned int  option,
double  value 
)
inline

Set UMFPack-specific options.

This method allows to set various options that control the UMFPack solver. More specifically, it allows to set values in the UMF_Control array. Please see the UMFPack documentation for a list of possible options and values.

Parameters
optionEntry in the UMF_Control array, e.g., UMFPACK_IRSTEP
valueCorresponding value
Exceptions
RangeErrorIf nonexisting option was requested

References DUNE_THROW.

◆ setSubMatrix() [1/2]

template<class M , class X , class Y >
template<class S >
void Dune::ILU0SubdomainSolver< M, X, Y >::setSubMatrix ( const M &  A,
S &  rowset 
)

Set the data of the local problem.

Parameters
AThe global matrix.
rowsetThe global indices of the local problem.
Template Parameters
SThe type of the set with the indices.

◆ setSubMatrix() [2/2]

template<class M , class X , class Y >
template<class S >
void Dune::ILUNSubdomainSolver< M, X, Y >::setSubMatrix ( const M &  A,
S &  rowset 
)

Set the data of the local problem.

Parameters
AThe global matrix.
rowsetThe global indices of the local problem.
Template Parameters
SThe type of the set with the indices.

◆ setVerbosity() [1/3]

template<typename T , typename A , int n, int m>
void Dune::LDL< BCRSMatrix< FieldMatrix< T, n, m >, A > >::setVerbosity ( int  v)
inline

Sets the verbosity level for the solver.

Parameters
vverbosity level: 0 only error messages, 1 a bit of statistics.

◆ setVerbosity() [2/3]

template<typename T , typename A , int n, int m>
void Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::setVerbosity ( int  v)
inline

Sets the verbosity level for the solver.

Parameters
vverbosity level: 0 only error messages, 1 a bit of statistics.

◆ setVerbosity() [3/3]

template<typename M >
void Dune::UMFPack< M >::setVerbosity ( int  v)
inline

sets the verbosity level for the UMFPack solver

Parameters
vverbosity level The following levels are implemented: 0 - only error messages 1 - a bit of statistics on decomposition and solution 2 - lots of statistics on decomposition and solution

Referenced by Dune::UMFPack< M >::UMFPack().

◆ SPQR() [1/3]

template<typename T , typename A , int n, int m>
Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR ( const Matrix matrix,
const ParameterTree config 
)
inline

Constructs the SPQR solver.

Parameters
matrixThe matrix of the system to solve.
configParameterTree containing solver parameters.
ParameterTree Key Meaning
verbose The verbosity level. default=0

◆ SPQR() [2/3]

template<typename T , typename A , int n, int m>
Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR ( const Matrix matrix,
int  verbose,
bool   
)
inline

Constructor for compatibility with SuperLU standard constructor.

This computes the matrix decomposition, and may take a long time (and use a lot of memory).

Parameters
matrixthe matrix to solve for
verbose0 or 1, set the verbosity level, defaults to 0

◆ SPQR() [3/3]

template<typename T , typename A , int n, int m>
Dune::SPQR< BCRSMatrix< FieldMatrix< T, n, m >, A > >::SPQR ( const Matrix matrix,
int  verbose = 0 
)
inline

Construct a solver object from a BCRSMatrix.

This computes the matrix decomposition, and may take a long time (and use a lot of memory).

Parameters
matrixthe matrix to solve for
verbose0 or 1, set the verbosity level, defaults to 0

◆ SuperLU() [1/2]

template<typename M >
Dune::SuperLU< M >::SuperLU

Empty default constructor.

Use setMatrix to tell SuperLU for what matrix it solves. Using this constructor no vectors will be reused.

◆ SuperLU() [2/2]

template<typename M >
Dune::SuperLU< M >::SuperLU ( const Matrix mat,
bool  verbose = false,
bool  reusevector = true 
)
explicit

Constructs the SuperLU solver.

During the construction the matrix will be decomposed. That means that in each apply call forward and backward substitutions take place (and no decomposition).

Parameters
matThe matrix of the system to solve.
verboseIf true some statistics are printed.
reusevectorDefault value is true. If true the two vectors are allocate in the first call to apply. These get resused in subsequent calls to apply and are deallocated in the destructor. If false these vectors are allocated at the beginning and deallocated at the end of each apply method. This allows using the same instance of superlu from different threads.

References Dune::SuperLU< M >::setMatrix().

◆ UMFPack() [1/5]

template<typename M >
Dune::UMFPack< M >::UMFPack ( const char *  file,
int  verbose = 0 
)
inline

try loading a decomposition from file

Parameters
filethe decomposition file
verbosethe verbosity level
Exceptions
Dune::ExceptionWhen not being able to load the file. Does not need knowledge of the actual matrix!

References DUNE_THROW, and Dune::UMFPack< M >::setVerbosity().

◆ UMFPack() [2/5]

template<typename M >
Dune::UMFPack< M >::UMFPack ( const Matrix mat_,
const char *  file,
int  verbose = 0 
)
inline

Try loading a decomposition from file and do a decomposition if unsuccessful.

Parameters
mat_the matrix to decompose when no decoposition file found
filethe decomposition file
verbosethe verbosity level

Use saveDecomposition(char* file) for manually storing a decomposition. This constructor will decompose mat_ and store the result to file if no file wasn't found in the first place. Thus, if you always use this you will only compute the decomposition once (and when you manually deleted the decomposition file).

References Dune::UMFPack< M >::saveDecomposition(), Dune::UMFPack< M >::setMatrix(), and Dune::UMFPack< M >::setVerbosity().

◆ UMFPack() [3/5]

template<typename M >
Dune::UMFPack< M >::UMFPack ( const Matrix mat_,
const ParameterTree config 
)
inline

Construct a solver object from a matrix.

Parameters
mat_the matrix to solve for
configParameterTree containing solver parameters.
ParameterTree Key Meaning
verbose The verbosity level. default=0

◆ UMFPack() [4/5]

template<typename M >
Dune::UMFPack< M >::UMFPack ( const Matrix matrix,
int  verbose,
bool   
)
inline

Constructor for compatibility with SuperLU standard constructor.

This computes the matrix decomposition, and may take a long time (and use a lot of memory).

Parameters
matrixthe matrix to solve for
verbose[0..2] set the verbosity level, defaults to 0

References Dune::UMFPack< M >::setMatrix(), and Dune::UMFPack< M >::setVerbosity().

◆ UMFPack() [5/5]

template<typename M >
Dune::UMFPack< M >::UMFPack ( const Matrix matrix,
int  verbose = 0 
)
inline

Construct a solver object from a matrix.

This computes the matrix decomposition, and may take a long time (and use a lot of memory).

Parameters
matrixthe matrix to solve for
verbose[0..2] set the verbosity level, defaults to 0

References Dune::UMFPack< M >::setMatrix(), and Dune::UMFPack< M >::setVerbosity().

Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.111.3 (Nov 20, 23:30, 2024)