MFEM
v4.6.0
Finite element discretization library
|
Data and methods for fully-assembled bilinear forms. More...
#include <bilinearform_ext.hpp>
Public Member Functions | |
FABilinearFormExtension (BilinearForm *form) | |
void | Assemble () |
Assemble at the level given for the BilinearFormExtension subclass. More... | |
void | RAP (OperatorHandle &A) |
void | EliminateBC (const Array< int > &ess_dofs, OperatorHandle &A) |
void | FormSystemMatrix (const Array< int > &ess_tdof_list, OperatorHandle &A) |
void | FormLinearSystem (const Array< int > &ess_tdof_list, Vector &x, Vector &b, OperatorHandle &A, Vector &X, Vector &B, int copy_interior=0) |
void | Mult (const Vector &x, Vector &y) const |
Operator application: y=A(x) . More... | |
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... | |
void | DGMult (const Vector &x, Vector &y) const |
void | DGMultTranspose (const Vector &x, Vector &y) const |
Public Member Functions inherited from mfem::EABilinearFormExtension | |
EABilinearFormExtension (BilinearForm *form) | |
Public Member Functions inherited from mfem::PABilinearFormExtension | |
PABilinearFormExtension (BilinearForm *) | |
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... | |
void | Update () |
Public Member Functions inherited from mfem::BilinearFormExtension | |
BilinearFormExtension (BilinearForm *form) | |
virtual MemoryClass | GetMemoryClass () const |
Return the MemoryClass preferred by the Operator. More... | |
virtual const Operator * | GetProlongation () const |
Get the finite element space prolongation matrix. More... | |
virtual const Operator * | GetRestriction () const |
Get the finite element space restriction matrix. More... | |
Public Member Functions inherited from mfem::Operator | |
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 void | AddMult (const Vector &x, Vector &y, const double a=1.0) const |
Operator application: y+=A(x) (default) or y+=a*A(x) . More... | |
virtual void | AddMultTranspose (const Vector &x, Vector &y, const double a=1.0) const |
Operator transpose application: y+=A^t(x) (default) or y+=a*A^t(x) . More... | |
virtual void | ArrayMult (const Array< const Vector *> &X, Array< Vector *> &Y) const |
Operator application on a matrix: Y=A(X) . More... | |
virtual void | ArrayMultTranspose (const Array< const Vector *> &X, Array< Vector *> &Y) const |
Action of the transpose operator on a matrix: Y=A^t(X) . More... | |
virtual void | ArrayAddMult (const Array< const Vector *> &X, Array< Vector *> &Y, const double a=1.0) const |
Operator application on a matrix: Y+=A(X) (default) or Y+=a*A(X) . More... | |
virtual void | ArrayAddMultTranspose (const Array< const Vector *> &X, Array< Vector *> &Y, const double a=1.0) const |
Operator transpose application on a matrix: Y+=A^t(X) (default) or Y+=a*A^t(X) . 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 * | 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... | |
Additional Inherited Members | |
Public Types inherited from mfem::Operator | |
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... | |
Protected Member Functions inherited from mfem::PABilinearFormExtension | |
void | SetupRestrictionOperators (const L2FaceValues m) |
void | AddMultWithMarkers (const BilinearFormIntegrator &integ, const Vector &x, const Array< int > *markers, const Array< int > &attributes, const bool transpose, Vector &y) const |
Accumulate the action (or transpose) of the integrator on x into y, taking into account the (possibly null) markers array. More... | |
Protected Member Functions inherited from mfem::Operator | |
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... | |
Protected Attributes inherited from mfem::EABilinearFormExtension | |
int | ne |
int | elemDofs |
Vector | ea_data |
int | nf_int |
int | nf_bdr |
int | faceDofs |
Vector | ea_data_int |
Vector | ea_data_ext |
Vector | ea_data_bdr |
bool | factorize_face_terms |
Protected Attributes inherited from mfem::PABilinearFormExtension | |
const FiniteElementSpace * | trial_fes |
const FiniteElementSpace * | test_fes |
Array< int > | elem_attributes |
Attributes of all mesh elements. More... | |
Array< int > | bdr_attributes |
Vector | tmp_evec |
Vector | localX |
Vector | localY |
Vector | int_face_X |
Vector | int_face_Y |
Vector | bdr_face_X |
Vector | bdr_face_Y |
const Operator * | elem_restrict |
const FaceRestriction * | int_face_restrict_lex |
const FaceRestriction * | bdr_face_restrict_lex |
Protected Attributes inherited from mfem::BilinearFormExtension | |
BilinearForm * | a |
Not owned. More... | |
Protected Attributes inherited from mfem::Operator | |
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... | |
Data and methods for fully-assembled bilinear forms.
Definition at line 140 of file bilinearform_ext.hpp.
mfem::FABilinearFormExtension::FABilinearFormExtension | ( | BilinearForm * | form | ) |
Definition at line 1027 of file bilinearform_ext.cpp.
|
virtual |
Assemble at the level given for the BilinearFormExtension subclass.
Reimplemented from mfem::EABilinearFormExtension.
Definition at line 1041 of file bilinearform_ext.cpp.
DGMult and DGMultTranspose use the extended L-vector to perform the computation.
Definition at line 1209 of file bilinearform_ext.cpp.
Definition at line 1270 of file bilinearform_ext.cpp.
void mfem::FABilinearFormExtension::EliminateBC | ( | const Array< int > & | ess_dofs, |
OperatorHandle & | A | ||
) |
DIAG_ONE
policy to be consistent with Operator::FormConstrainedSystemOperator
. Definition at line 1170 of file bilinearform_ext.cpp.
|
virtual |
Reimplemented from mfem::PABilinearFormExtension.
Definition at line 1197 of file bilinearform_ext.cpp.
|
virtual |
Reimplemented from mfem::PABilinearFormExtension.
Definition at line 1190 of file bilinearform_ext.cpp.
Operator application: y=A(x)
.
Reimplemented from mfem::EABilinearFormExtension.
Definition at line 1258 of file bilinearform_ext.cpp.
Action of the transpose operator: y=A^t(x)
. The default behavior in class Operator is to generate an error.
Reimplemented from mfem::EABilinearFormExtension.
Definition at line 1319 of file bilinearform_ext.cpp.
void mfem::FABilinearFormExtension::RAP | ( | OperatorHandle & | A | ) |
Definition at line 1156 of file bilinearform_ext.cpp.