DUNE PDELab (git)

SuperLu Solver. More...

#include <dune/istl/superlu.hh>

Public Types

using Matrix = M
 The matrix type.
 
typedef Dune::SuperLUMatrix< MatrixSuperLUMatrix
 The corresponding SuperLU Matrix type.
 
typedef SuperMatrixInitializer< MatrixMatrixInitializer
 Type of an associated initializer class.
 
using domain_type = typename Impl::SuperLUVectorChooser< M >::domain_type
 The type of the domain of the solver.
 
using range_type = typename Impl::SuperLUVectorChooser< M >::range_type
 The type of the range of the solver.
 

Public Member Functions

SolverCategory::Category category () const override
 Category of the solver (see SolverCategory::Category)
 
 SuperLU (const Matrix &mat, bool verbose=false, bool reusevector=true)
 Constructs the SuperLU solver. More...
 
 SuperLU (const Matrix &mat, const ParameterTree &config)
 Constructs the SuperLU solver. More...
 
 SuperLU ()
 Empty default constructor. More...
 
void apply (domain_type &x, range_type &b, InverseOperatorResult &res) override
 Apply inverse operator,. More...
 
void apply (domain_type &x, range_type &b, double reduction, InverseOperatorResult &res) override
 apply inverse operator, with given convergence criteria. More...
 
void apply (T *x, T *b)
 Apply SuperLu to C arrays.
 
void setMatrix (const Matrix &mat)
 Initialize data from given matrix.
 
void free ()
 free allocated space. More...
 

Detailed Description

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

SuperLu Solver.

Uses the well known SuperLU package to solve the system.

SuperLU supports single and double precision floating point and complex numbers. Unfortunately these cannot be used at the same time. Therefore users must set SUPERLU_NTYPE (0: float, 1: double, 2: std::complex<float>, 3: std::complex<double>) if the numeric type should be different from double.

Constructor & Destructor Documentation

◆ SuperLU()

template<typename M >
Dune::SuperLU< M >::SuperLU ( const Matrix mat,
const ParameterTree config 
)
inline

Constructs the SuperLU solver.

Parameters
matThe matrix of the system to solve.
configParameterTree containing solver parameters.
ParameterTree Key Meaning
verbose The verbosity level. default=false
reuseVector Reuse initially allocated vectors in apply. default=true

Member Function Documentation

◆ apply()

template<typename M >
void Dune::SuperLU< 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::SuperLU< M >::apply().


The documentation for this class was generated from the following file:
Creative Commons License   |  Legal Statements / Impressum  |  Hosted by TU Dresden  |  generated with Hugo v0.111.3 (Nov 23, 23:29, 2024)