![]() |
MFEM
v4.6.0
Finite element discretization library
|
Operator that extracts face degrees of freedom for H1, ND, or RT FiniteElementSpaces. More...
#include <restriction.hpp>
Public Member Functions | |
ConformingFaceRestriction (const FiniteElementSpace &fes, const ElementDofOrdering f_ordering, const FaceType type) | |
Construct a ConformingFaceRestriction. More... | |
void | Mult (const Vector &x, Vector &y) const override |
Scatter the degrees of freedom, i.e. goes from L-Vector to face E-Vector. More... | |
void | AddMultTranspose (const Vector &x, Vector &y, const double a=1.0) const override |
Gather the degrees of freedom, i.e. goes from face E-Vector to L-Vector. More... | |
void | AddMultTransposeUnsigned (const Vector &x, Vector &y, const double a=1.0) const override |
Gather the degrees of freedom, i.e. goes from face E-Vector to L-Vector not taking into account signs from DOF orientations. More... | |
virtual void | AddMultTransposeInPlace (Vector &x, Vector &y) const |
Add the face degrees of freedom x to the element degrees of freedom y. Perform the same computation as AddMultTranspose, but x is invalid after calling this method. More... | |
![]() | |
FaceRestriction () | |
FaceRestriction (int h, int w) | |
virtual | ~FaceRestriction () |
virtual void | AddMultTransposeInPlace (Vector &x, Vector &y) const |
Add the face degrees of freedom x to the element degrees of freedom y. Perform the same computation as AddMultTranspose, but x is invalid after calling this method. More... | |
void | MultTranspose (const Vector &x, Vector &y) const override |
Set the face degrees of freedom in the element degrees of freedom y to the values given in x. 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 | 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 | 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 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 Member Functions | |
ConformingFaceRestriction (const FiniteElementSpace &fes, const ElementDofOrdering f_ordering, const FaceType type, bool build) | |
Construct a ConformingFaceRestriction. More... | |
void | CheckFESpace (const ElementDofOrdering f_ordering) |
Verify that ConformingFaceRestriction is built from a supported finite element space. More... | |
void | SetFaceDofsScatterIndices (const Mesh::FaceInformation &face, const int face_index, const ElementDofOrdering f_ordering) |
Set the scattering indices of elem1, and increment the offsets for the face described by the face. More... | |
void | SetFaceDofsGatherIndices (const Mesh::FaceInformation &face, const int face_index, const ElementDofOrdering f_ordering) |
Set the gathering indices of elem1 for the interior face described by the face. 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... | |
Protected Attributes | |
const FiniteElementSpace & | fes |
const int | nf |
const int | vdim |
const bool | byvdim |
const int | face_dofs |
const int | elem_dofs |
const int | nfdofs |
const int | ndofs |
Array< int > | scatter_indices |
Array< int > | gather_offsets |
Array< int > | gather_indices |
Array< int > | vol_dof_map |
Array< int > | face_map |
![]() | |
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... | |
Operator that extracts face degrees of freedom for H1, ND, or RT FiniteElementSpaces.
Objects of this type are typically created and owned by FiniteElementSpace objects, see FiniteElementSpace::GetFaceRestriction().
Definition at line 230 of file restriction.hpp.
|
protected |
Construct a ConformingFaceRestriction.
[in] | fes | The FiniteElementSpace on which this operates |
[in] | f_ordering | Request a specific face dof ordering |
[in] | type | Request internal or boundary faces dofs |
[in] | build | Request the NCL2FaceRestriction to compute the scatter/gather indices. False should only be used when inheriting from ConformingFaceRestriction. |
Definition at line 601 of file restriction.cpp.
mfem::ConformingFaceRestriction::ConformingFaceRestriction | ( | const FiniteElementSpace & | fes, |
const ElementDofOrdering | f_ordering, | ||
const FaceType | type | ||
) |
Construct a ConformingFaceRestriction.
[in] | fes | The FiniteElementSpace on which this operates |
[in] | f_ordering | Request a specific face dof ordering |
[in] | type | Request internal or boundary faces dofs |
Definition at line 647 of file restriction.cpp.
|
overridevirtual |
Gather the degrees of freedom, i.e. goes from face E-Vector to L-Vector.
[in] | x | The face E-Vector degrees of freedom with the given format: face_dofs x vdim x nf where nf is the number of interior or boundary faces requested by type in the constructor. The face_dofs should be ordered according to the given ElementDofOrdering |
[in,out] | y | The L-vector degrees of freedom. |
[in] | a | Scalar coefficient for addition. |
Implements mfem::FaceRestriction.
Definition at line 717 of file restriction.cpp.
|
inline |
Add the face degrees of freedom x to the element degrees of freedom y. Perform the same computation as AddMultTranspose, but x is invalid after calling this method.
[in,out] | x | The face degrees of freedom on the face. |
[in,out] | y | The L-vector of degrees of freedom to which we add the face degrees of freedom. |
Definition at line 206 of file restriction.hpp.
|
overridevirtual |
Gather the degrees of freedom, i.e. goes from face E-Vector to L-Vector not taking into account signs from DOF orientations.
Reimplemented from mfem::FaceRestriction.
Definition at line 725 of file restriction.cpp.
|
protected |
Verify that ConformingFaceRestriction is built from a supported finite element space.
[in] | f_ordering | The requested face dof ordering. |
Definition at line 733 of file restriction.cpp.
Scatter the degrees of freedom, i.e. goes from L-Vector to face E-Vector.
[in] | x | The L-vector degrees of freedom. |
[out] | y | The face E-Vector degrees of freedom with the given format: face_dofs x vdim x nf where nf is the number of interior or boundary faces requested by type in the constructor. The face_dofs are ordered according to the given ElementDofOrdering. |
Implements mfem::FaceRestriction.
Definition at line 654 of file restriction.cpp.
|
protected |
Set the gathering indices of elem1 for the interior face described by the face.
[in] | face | The face information of the current face. |
[in] | face_index | The interior/boundary face index. |
[in] | f_ordering | Request a specific face dof ordering. |
Definition at line 876 of file restriction.cpp.
|
protected |
Set the scattering indices of elem1, and increment the offsets for the face described by the face.
[in] | face | The face information of the current face. |
[in] | face_index | The interior/boundary face index. |
[in] | f_ordering | Request a specific face dof ordering. |
Definition at line 846 of file restriction.cpp.
|
protected |
Definition at line 236 of file restriction.hpp.
|
protected |
Definition at line 238 of file restriction.hpp.
|
protected |
Definition at line 237 of file restriction.hpp.
|
mutableprotected |
Definition at line 325 of file restriction.hpp.
|
protected |
Definition at line 233 of file restriction.hpp.
|
protected |
Definition at line 243 of file restriction.hpp.
|
protected |
Definition at line 242 of file restriction.hpp.
|
protected |
Definition at line 240 of file restriction.hpp.
|
protected |
Definition at line 234 of file restriction.hpp.
|
protected |
Definition at line 239 of file restriction.hpp.
|
protected |
Definition at line 241 of file restriction.hpp.
|
protected |
Definition at line 235 of file restriction.hpp.
|
protected |
Definition at line 244 of file restriction.hpp.