![]() |
MFEM
v4.5.1
Finite element discretization library
|
FGMRES method. More...
#include <solvers.hpp>
Public Member Functions | |
FGMRESSolver () | |
FGMRESSolver (MPI_Comm comm_) | |
void | SetKDim (int dim) |
virtual void | Mult (const Vector &b, Vector &x) const |
Operator application: y=A(x) . More... | |
![]() | |
IterativeSolver () | |
IterativeSolver (MPI_Comm comm_) | |
virtual void | SetPreconditioner (Solver &pr) |
This should be called before SetOperator. More... | |
virtual void | SetOperator (const Operator &op) override |
Also calls SetOperator for the preconditioner if there is one. More... | |
void | SetMonitor (IterativeSolverMonitor &m) |
Set the iterative solver monitor. More... | |
MPI_Comm | GetComm () const |
Return the associated MPI communicator, or MPI_COMM_NULL if no communicator is set. More... | |
void | SetRelTol (double rtol) |
void | SetAbsTol (double atol) |
void | SetMaxIter (int max_it) |
virtual void | SetPrintLevel (int print_lvl) |
Legacy method to set the level of verbosity of the solver output. More... | |
virtual void | SetPrintLevel (PrintLevel) |
Set the level of verbosity of the solver output. More... | |
int | GetNumIterations () const |
bool | GetConverged () const |
double | GetFinalNorm () const |
![]() | |
Solver (int s=0, bool iter_mode=false) | |
Initialize a square Solver with size s. More... | |
Solver (int h, int w, bool iter_mode=false) | |
Initialize a Solver with height h and width w. More... | |
![]() | |
void | InitTVectors (const Operator *Po, const Operator *Ri, const Operator *Pi, Vector &x, Vector &b, Vector &X, Vector &B) const |
Initializes memory for true vectors of linear system. More... | |
Operator (int s=0) | |
Construct a square Operator with given size s (default 0). More... | |
Operator (int h, int w) | |
Construct an Operator with the given height (output size) and width (input size). More... | |
int | Height () const |
Get the height (size of output) of the Operator. Synonym with NumRows(). More... | |
int | NumRows () const |
Get the number of rows (size of output) of the Operator. Synonym with Height(). More... | |
int | Width () const |
Get the width (size of input) of the Operator. Synonym with NumCols(). More... | |
int | NumCols () const |
Get the number of columns (size of input) of the Operator. Synonym with Width(). More... | |
virtual MemoryClass | GetMemoryClass () const |
Return the MemoryClass preferred by the Operator. More... | |
virtual void | MultTranspose (const Vector &x, Vector &y) const |
Action of the transpose operator: y=A^t(x) . The default behavior in class Operator is to generate an error. More... | |
virtual Operator & | GetGradient (const Vector &x) const |
Evaluate the gradient operator at the point x. The default behavior in class Operator is to generate an error. More... | |
virtual void | AssembleDiagonal (Vector &diag) const |
Computes the diagonal entries into diag. Typically, this operation only makes sense for linear Operators. In some cases, only an approximation of the diagonal is computed. More... | |
virtual const Operator * | GetProlongation () const |
Prolongation operator from linear algebra (linear system) vectors, to input vectors for the operator. NULL means identity. More... | |
virtual const Operator * | GetRestriction () const |
Restriction operator from input vectors for the operator to linear algebra (linear system) vectors. NULL means identity. More... | |
virtual const Operator * | GetOutputProlongation () const |
Prolongation operator from linear algebra (linear system) vectors, to output vectors for the operator. NULL means identity. More... | |
virtual const Operator * | GetOutputRestrictionTranspose () const |
Transpose of GetOutputRestriction, directly available in this form to facilitate matrix-free RAP-type operators. More... | |
virtual const Operator * | GetOutputRestriction () const |
Restriction operator from output vectors for the operator to linear algebra (linear system) vectors. NULL means identity. More... | |
void | FormLinearSystem (const Array< int > &ess_tdof_list, Vector &x, Vector &b, Operator *&A, Vector &X, Vector &B, int copy_interior=0) |
Form a constrained linear system using a matrix-free approach. More... | |
void | FormRectangularLinearSystem (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, Vector &x, Vector &b, Operator *&A, Vector &X, Vector &B) |
Form a column-constrained linear system using a matrix-free approach. More... | |
virtual void | RecoverFEMSolution (const Vector &X, const Vector &b, Vector &x) |
Reconstruct a solution vector x (e.g. a GridFunction) from the solution X of a constrained linear system obtained from Operator::FormLinearSystem() or Operator::FormRectangularLinearSystem(). More... | |
void | FormSystemOperator (const Array< int > &ess_tdof_list, Operator *&A) |
Return in A a parallel (on truedofs) version of this square operator. More... | |
void | FormRectangularSystemOperator (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, Operator *&A) |
Return in A a parallel (on truedofs) version of this rectangular operator (including constraints). More... | |
void | FormDiscreteOperator (Operator *&A) |
Return in A a parallel (on truedofs) version of this rectangular operator. More... | |
void | PrintMatlab (std::ostream &out, int n, int m=0) const |
Prints operator with input size n and output size m in Matlab format. More... | |
virtual void | PrintMatlab (std::ostream &out) const |
Prints operator in Matlab format. More... | |
virtual | ~Operator () |
Virtual destructor. More... | |
Type | GetType () const |
Return the type ID of the Operator class. More... | |
Protected Attributes | |
int | m |
![]() | |
const Operator * | oper |
Solver * | prec |
IterativeSolverMonitor * | monitor = nullptr |
int | max_iter |
Limit for the number of iterations the solver is allowed to do. More... | |
double | rel_tol |
Relative tolerance. More... | |
double | abs_tol |
Absolute tolerance. More... | |
int | final_iter |
bool | converged |
double | final_norm |
int | print_level = -1 |
(DEPRECATED) Legacy print level definition, which is left for compatibility with custom iterative solvers. More... | |
PrintLevel | print_options |
Output behavior for the iterative solver. More... | |
![]() | |
int | height |
Dimension of the output / number of rows in the matrix. More... | |
int | width |
Dimension of the input / number of columns in the matrix. More... | |
Additional Inherited Members | |
![]() | |
enum | DiagonalPolicy { DIAG_ZERO, DIAG_ONE, DIAG_KEEP } |
Defines operator diagonal policy upon elimination of rows and/or columns. More... | |
enum | Type { ANY_TYPE, MFEM_SPARSEMAT, Hypre_ParCSR, PETSC_MATAIJ, PETSC_MATIS, PETSC_MATSHELL, PETSC_MATNEST, PETSC_MATHYPRE, PETSC_MATGENERIC, Complex_Operator, MFEM_ComplexSparseMat, Complex_Hypre_ParCSR, Complex_DenseMat, MFEM_Block_Matrix, MFEM_Block_Operator } |
Enumeration defining IDs for some classes derived from Operator. More... | |
![]() | |
bool | iterative_mode |
If true, use the second argument of Mult() as an initial guess. More... | |
![]() | |
double | Dot (const Vector &x, const Vector &y) const |
double | Norm (const Vector &x) const |
void | Monitor (int it, double norm, const Vector &r, const Vector &x, bool final=false) const |
PrintLevel | FromLegacyPrintLevel (int) |
Convert a legacy print level integer to a PrintLevel object. More... | |
![]() | |
void | FormConstrainedSystemOperator (const Array< int > &ess_tdof_list, ConstrainedOperator *&Aout) |
see FormSystemOperator() More... | |
void | FormRectangularConstrainedSystemOperator (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, RectangularConstrainedOperator *&Aout) |
see FormRectangularSystemOperator() More... | |
Operator * | SetupRAP (const Operator *Pi, const Operator *Po) |
Returns RAP Operator of this, using input/output Prolongation matrices Pi corresponds to "P", Po corresponds to "Rt". More... | |
![]() | |
static int | GuessLegacyPrintLevel (PrintLevel) |
Use some heuristics to guess a legacy print level corresponding to the given PrintLevel. More... | |
FGMRES method.
Definition at line 516 of file solvers.hpp.
|
inline |
Definition at line 522 of file solvers.hpp.
|
inline |
Definition at line 525 of file solvers.hpp.
Operator application: y=A(x)
.
Implements mfem::Operator.
Definition at line 1159 of file solvers.cpp.
|
inline |
Definition at line 528 of file solvers.hpp.
|
protected |
Definition at line 519 of file solvers.hpp.