![]() |
MFEM
v4.2.0
Finite element discretization library
|
Interface with UMFPack solver specialized for ComplexSparseMatrix This approach avoids forming a monolithic SparseMatrix which leads to increased memory and flops. More...
#include <complex_operator.hpp>
Public Member Functions | |
ComplexUMFPackSolver (bool _use_long_ints=false, bool transa_=false) | |
For larger matrices, if the solver fails, set the parameter _use_long_ints = true. More... | |
ComplexUMFPackSolver (ComplexSparseMatrix &A, bool _use_long_ints=false, bool transa_=false) | |
Factorize the given ComplexSparseMatrix using the defaults. For larger matrices, if the solver fails, set the parameter _use_long_ints = true. More... | |
virtual void | SetOperator (const Operator &op) |
Factorize the given Operator op which must be a ComplexSparseMatrix. More... | |
void | SetPrintLevel (int print_lvl) |
void | SetTransposeSolve (bool transa_) |
virtual void | Mult (const Vector &b, Vector &x) const |
This is solving the system A x = b. More... | |
virtual void | MultTranspose (const Vector &b, Vector &x) const |
This is solving the system: A^H x = b (when transa = false) This is equivalent to solving the transpose block system for the case of Convention = HERMITIAN A^T x = b (when transa = true) This is equivalent to solving the transpose block system for the case of Convention = BLOCK_SYMMETRIC. More... | |
virtual | ~ComplexUMFPackSolver () |
![]() | |
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 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 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 * | 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=0, int m=0) const |
Prints operator with input size n and output size m in Matlab format. More... | |
virtual | ~Operator () |
Virtual destructor. More... | |
Type | GetType () const |
Return the type ID of the Operator class. More... | |
Public Attributes | |
double | Control [UMFPACK_CONTROL] |
double | Info [UMFPACK_INFO] |
![]() | |
bool | iterative_mode |
If true, use the second argument of Mult() as an initial guess. More... | |
Protected Member Functions | |
void | Init () |
![]() | |
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, taking in input/output Prolongation matrices. More... | |
Protected Attributes | |
bool | use_long_ints |
bool | transa |
ComplexSparseMatrix * | mat |
void * | Numeric |
SuiteSparse_long * | AI |
SuiteSparse_long * | AJ |
![]() | |
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 } |
Enumeration defining IDs for some classes derived from Operator. More... | |
Interface with UMFPack solver specialized for ComplexSparseMatrix This approach avoids forming a monolithic SparseMatrix which leads to increased memory and flops.
Definition at line 182 of file complex_operator.hpp.
|
inline |
For larger matrices, if the solver fails, set the parameter _use_long_ints = true.
Definition at line 200 of file complex_operator.hpp.
|
inline |
Factorize the given ComplexSparseMatrix using the defaults. For larger matrices, if the solver fails, set the parameter _use_long_ints = true.
Definition at line 205 of file complex_operator.hpp.
|
virtual |
Definition at line 507 of file complex_operator.cpp.
|
protected |
Definition at line 273 of file complex_operator.cpp.
This is solving the system A x = b.
Implements mfem::Operator.
Definition at line 393 of file complex_operator.cpp.
This is solving the system: A^H x = b (when transa = false) This is equivalent to solving the transpose block system for the case of Convention = HERMITIAN A^T x = b (when transa = true) This is equivalent to solving the transpose block system for the case of Convention = BLOCK_SYMMETRIC.
Reimplemented from mfem::Operator.
Definition at line 446 of file complex_operator.cpp.
|
virtual |
Factorize the given Operator op which must be a ComplexSparseMatrix.
The factorization uses the parameters set in the Control data member.
Implements mfem::Solver.
Definition at line 288 of file complex_operator.cpp.
|
inline |
Definition at line 219 of file complex_operator.hpp.
|
inline |
Definition at line 222 of file complex_operator.hpp.
|
protected |
Definition at line 190 of file complex_operator.hpp.
|
protected |
Definition at line 190 of file complex_operator.hpp.
double mfem::ComplexUMFPackSolver::Control[UMFPACK_CONTROL] |
Definition at line 195 of file complex_operator.hpp.
|
mutable |
Definition at line 196 of file complex_operator.hpp.
|
protected |
Definition at line 187 of file complex_operator.hpp.
|
protected |
Definition at line 189 of file complex_operator.hpp.
|
protected |
Definition at line 186 of file complex_operator.hpp.
|
protected |
Definition at line 185 of file complex_operator.hpp.