![]() |
MFEM v4.8.0
Finite element discretization library
|
#include <bilinearform.hpp>
Public Member Functions | |
MixedBilinearForm (FiniteElementSpace *tr_fes, FiniteElementSpace *te_fes) | |
Construct a MixedBilinearForm on the given trial, tr_fes, and test, te_fes, FiniteElementSpaces. | |
MixedBilinearForm (FiniteElementSpace *tr_fes, FiniteElementSpace *te_fes, MixedBilinearForm *mbf) | |
Create a MixedBilinearForm on the given trial, tr_fes, and test, te_fes, FiniteElementSpaces, using the same integrators as the MixedBilinearForm mbf. | |
real_t & | Elem (int i, int j) override |
Returns a reference to: \( M_{ij} \). | |
const real_t & | Elem (int i, int j) const override |
Returns a reference to: \( M_{ij} \). | |
void | Mult (const Vector &x, Vector &y) const override |
Matrix multiplication: \( y = M x \). | |
void | AddMult (const Vector &x, Vector &y, const real_t a=1.0) const override |
Add the matrix vector multiple to a vector: \( y += a M x \). | |
void | MultTranspose (const Vector &x, Vector &y) const override |
Matrix transpose vector multiplication: \( y = M^T x \). | |
void | AddMultTranspose (const Vector &x, Vector &y, const real_t a=1.0) const override |
Add the matrix transpose vector multiplication: \( y += a M^T x \). | |
MatrixInverse * | Inverse () const override |
Returns a pointer to (approximation) of the matrix inverse: \( M^{-1} \) (currently unimplemented and returns NULL) | |
void | Finalize (int skip_zeros=1) override |
Finalizes the matrix initialization if the AssemblyLevel is AssemblyLevel::LEGACY. | |
void | GetBlocks (Array2D< SparseMatrix * > &blocks) const |
Extract the associated matrix as SparseMatrix blocks. The number of block rows and columns is given by the vector dimensions (vdim) of the test and trial spaces, respectively. | |
const SparseMatrix & | SpMat () const |
Returns a const reference to the sparse matrix: \( M \). | |
SparseMatrix & | SpMat () |
Returns a reference to the sparse matrix: \( M \). | |
SparseMatrix * | LoseMat () |
Nullifies the internal matrix \( M \) and returns a pointer to it. Used for transferring ownership. | |
const SparseMatrix & | SpMatElim () const |
Returns a const reference to the sparse matrix of eliminated b.c.: \( M_e \). | |
SparseMatrix & | SpMatElim () |
Returns a reference to the sparse matrix of eliminated b.c.: \( M_e \). | |
void | AddDomainIntegrator (BilinearFormIntegrator *bfi) |
Adds a domain integrator. Assumes ownership of bfi. | |
void | AddDomainIntegrator (BilinearFormIntegrator *bfi, Array< int > &elem_marker) |
Adds a domain integrator. Assumes ownership of bfi. | |
void | AddBoundaryIntegrator (BilinearFormIntegrator *bfi) |
Adds a boundary integrator. Assumes ownership of bfi. | |
void | AddBoundaryIntegrator (BilinearFormIntegrator *bfi, Array< int > &bdr_marker) |
Adds a boundary integrator. Assumes ownership of bfi. | |
void | AddInteriorFaceIntegrator (BilinearFormIntegrator *bfi) |
Adds an interior face integrator. Assumes ownership of bfi. | |
void | AddBdrFaceIntegrator (BilinearFormIntegrator *bfi) |
Adds a boundary face integrator. Assumes ownership of bfi. | |
void | AddBdrFaceIntegrator (BilinearFormIntegrator *bfi, Array< int > &bdr_marker) |
Adds a boundary face integrator. Assumes ownership of bfi. | |
void | AddTraceFaceIntegrator (BilinearFormIntegrator *bfi) |
Add a trace face integrator. Assumes ownership of bfi. | |
void | AddBdrTraceFaceIntegrator (BilinearFormIntegrator *bfi) |
Adds a boundary trace face integrator. Assumes ownership of bfi. | |
void | AddBdrTraceFaceIntegrator (BilinearFormIntegrator *bfi, Array< int > &bdr_marker) |
Adds a boundary trace face integrator. Assumes ownership of bfi. | |
Array< BilinearFormIntegrator * > * | GetDBFI () |
Access all integrators added with AddDomainIntegrator(). | |
Array< Array< int > * > * | GetDBFI_Marker () |
Access all domain markers added with AddDomainIntegrator(). If no marker was specified when the integrator was added, the corresponding pointer (to Array<int>) will be NULL. | |
Array< BilinearFormIntegrator * > * | GetBBFI () |
Access all integrators added with AddBoundaryIntegrator(). | |
Array< Array< int > * > * | GetBBFI_Marker () |
Access all boundary markers added with AddBoundaryIntegrator(). | |
Array< BilinearFormIntegrator * > * | GetFBFI () |
Access all integrators added with AddInteriorFaceIntegrator(). | |
Array< BilinearFormIntegrator * > * | GetBFBFI () |
Access all integrators added with AddBdrFaceIntegrator(). | |
Array< Array< int > * > * | GetBFBFI_Marker () |
Access all boundary markers added with AddBdrFaceIntegrator(). If no marker was specified when the integrator was added, the corresponding pointer (to Array<int>) will be NULL. | |
Array< BilinearFormIntegrator * > * | GetTFBFI () |
Access all integrators added with AddTraceFaceIntegrator(). | |
Array< BilinearFormIntegrator * > * | GetBTFBFI () |
Access all integrators added with AddBdrTraceFaceIntegrator(). | |
Array< Array< int > * > * | GetBTFBFI_Marker () |
Access all boundary markers added with AddBdrTraceFaceIntegrator() | |
void | operator= (const real_t a) |
Sets all sparse values of \( M \) to a. | |
void | SetAssemblyLevel (AssemblyLevel assembly_level) |
Set the desired assembly level. The default is AssemblyLevel::LEGACY. | |
void | Assemble (int skip_zeros=1) |
void | AssembleDiagonal_ADAt (const Vector &D, Vector &diag) const |
Assemble the diagonal of ADA^T into diag, where A is this mixed bilinear form and D is a diagonal. | |
const Operator * | GetProlongation () const override |
Get the input finite element space prolongation matrix. | |
const Operator * | GetRestriction () const override |
Get the input finite element space restriction matrix. | |
const Operator * | GetOutputProlongation () const override |
Get the test finite element space prolongation matrix. | |
const Operator * | GetOutputRestriction () const override |
Get the test finite element space restriction matrix. | |
void | ConformingAssemble () |
For partially conforming trial and/or test FE spaces, complete the assembly process by performing \( P2^t A P1 \) where \( A \) is the internal sparse matrix; \( P1 \) and \( P2 \) are the conforming prolongation matrices of the trial and test FE spaces, respectively. After this call the MixedBilinearForm becomes an operator on the conforming FE spaces. | |
void | ComputeElementMatrix (int i, DenseMatrix &elmat) const |
Compute the element matrix of the given element. | |
void | ComputeBdrElementMatrix (int i, DenseMatrix &elmat) const |
Compute the boundary element matrix of the given boundary element. | |
void | ComputeTraceFaceMatrix (int i, DenseMatrix &elmat) const |
Compute the trace face matrix of the given face element. | |
void | ComputeBdrTraceFaceMatrix (int i, DenseMatrix &elmat) const |
Compute the boundary trace face matrix of the given boundary element. | |
void | ComputeFaceMatrix (int i, DenseMatrix &elmat) const |
Compute the face matrix of the given face element. | |
void | ComputeBdrFaceMatrix (int i, DenseMatrix &elmat) const |
Compute the boundary face matrix of the given boundary element. | |
void | AssembleElementMatrix (int i, const DenseMatrix &elmat, int skip_zeros=1) |
Assemble the given element matrix. | |
void | AssembleElementMatrix (int i, const DenseMatrix &elmat, Array< int > &trial_vdofs, Array< int > &test_vdofs, int skip_zeros=1) |
Assemble the given element matrix. | |
void | AssembleBdrElementMatrix (int i, const DenseMatrix &elmat, int skip_zeros=1) |
Assemble the given boundary element matrix. | |
void | AssembleBdrElementMatrix (int i, const DenseMatrix &elmat, Array< int > &trial_vdofs, Array< int > &test_vdofs, int skip_zeros=1) |
Assemble the given boundary element matrix. | |
void | EliminateTrialEssentialBC (const Array< int > &bdr_attr_is_ess, const Vector &sol, Vector &rhs) |
Eliminate essential boundary trial DOFs from the system. | |
void | EliminateTrialEssentialBC (const Array< int > &bdr_attr_is_ess) |
Eliminate essential boundary trial DOFs from the system matrix. | |
MFEM_DEPRECATED void | EliminateTrialDofs (const Array< int > &bdr_attr_is_ess, const Vector &sol, Vector &rhs) |
(DEPRECATED) Eliminate essential boundary trial DOFs from the system. | |
void | EliminateTrialVDofs (const Array< int > &vdofs, const Vector &sol, Vector &rhs) |
Eliminate the given trial vdofs. NOTE: here, vdofs is a list of DOFs. | |
void | EliminateTrialVDofs (const Array< int > &vdofs) |
Eliminate the given trial vdofs, storing the eliminated part internally in \( M_e \). | |
void | EliminateTrialVDofsInRHS (const Array< int > &vdofs, const Vector &x, Vector &b) |
Use the stored eliminated part of the matrix (see EliminateTrialVDofs(const Array<int> &)) to modify the r.h.s. b; vdofs is a list of DOFs (non-directional, i.e. >= 0). | |
void | EliminateEssentialBCFromTrialDofs (const Array< int > &marked_vdofs, const Vector &sol, Vector &rhs) |
Similar to EliminateTrialVDofs(const Array<int> &, const Vector &, Vector &) but here ess_dofs is a marker (boolean) array on all vector-dofs (ess_dofs[i] < 0 is true). | |
void | EliminateTestEssentialBC (const Array< int > &bdr_attr_is_ess) |
Eliminate essential boundary test DOFs from the system matrix. | |
virtual MFEM_DEPRECATED void | EliminateTestDofs (const Array< int > &bdr_attr_is_ess) |
(DEPRECATED) Eliminate essential boundary test DOFs from the system. | |
void | EliminateTestVDofs (const Array< int > &vdofs) |
Eliminate the given test vdofs. NOTE: here, vdofs is a list of DOFs. | |
virtual void | FormRectangularSystemMatrix (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, OperatorHandle &A) |
Return in A that is column-constrained. | |
template<typename OpType > | |
void | FormRectangularSystemMatrix (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, OpType &A) |
Form the column-constrained linear system matrix A. | |
virtual void | FormRectangularLinearSystem (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, Vector &x, Vector &b, OperatorHandle &A, Vector &X, Vector &B) |
Form the linear system A X = B, corresponding to this mixed bilinear form and the linear form b(.). | |
template<typename OpType > | |
void | FormRectangularLinearSystem (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, Vector &x, Vector &b, OpType &A, Vector &X, Vector &B) |
Form the linear system A X = B, corresponding to this bilinear form and the linear form b(.). | |
void | Update () |
Must be called after making changes to trial_fes or test_fes. | |
FiniteElementSpace * | TrialFESpace () |
Return the trial FE space associated with the BilinearForm. | |
const FiniteElementSpace * | TrialFESpace () const |
Read-only access to the associated trial FiniteElementSpace. | |
FiniteElementSpace * | TestFESpace () |
Return the test FE space associated with the BilinearForm. | |
const FiniteElementSpace * | TestFESpace () const |
Read-only access to the associated test FiniteElementSpace. | |
void | Print (std::ostream &os=mfem::out, int width_=4) const override |
Prints operator to stream os. | |
virtual | ~MixedBilinearForm () |
Deletes internal matrices, bilinear integrators, and the BilinearFormExtension. | |
![]() | |
Matrix (int s) | |
Creates a square matrix of size s. | |
Matrix (int h, int w) | |
Creates a matrix of the given height and width. | |
bool | IsSquare () const |
Returns whether the matrix is a square matrix. | |
virtual | ~Matrix () |
Destroys matrix. | |
![]() | |
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. | |
Operator (int s=0) | |
Construct a square Operator with given size s (default 0). | |
Operator (int h, int w) | |
Construct an Operator with the given height (output size) and width (input size). | |
int | Height () const |
Get the height (size of output) of the Operator. Synonym with NumRows(). | |
int | NumRows () const |
Get the number of rows (size of output) of the Operator. Synonym with Height(). | |
int | Width () const |
Get the width (size of input) of the Operator. Synonym with NumCols(). | |
int | NumCols () const |
Get the number of columns (size of input) of the Operator. Synonym with Width(). | |
virtual MemoryClass | GetMemoryClass () const |
Return the MemoryClass preferred by the Operator. | |
virtual void | ArrayMult (const Array< const Vector * > &X, Array< Vector * > &Y) const |
Operator application on a matrix: Y=A(X) . | |
virtual void | ArrayMultTranspose (const Array< const Vector * > &X, Array< Vector * > &Y) const |
Action of the transpose operator on a matrix: Y=A^t(X) . | |
virtual void | ArrayAddMult (const Array< const Vector * > &X, Array< Vector * > &Y, const real_t a=1.0) const |
Operator application on a matrix: Y+=A(X) (default) or Y+=a*A(X) . | |
virtual void | ArrayAddMultTranspose (const Array< const Vector * > &X, Array< Vector * > &Y, const real_t a=1.0) const |
Operator transpose application on a matrix: Y+=A^t(X) (default) or Y+=a*A^t(X) . | |
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. | |
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. | |
virtual const Operator * | GetOutputRestrictionTranspose () const |
Transpose of GetOutputRestriction, directly available in this form to facilitate matrix-free RAP-type operators. | |
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. | |
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. | |
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(). | |
void | FormSystemOperator (const Array< int > &ess_tdof_list, Operator *&A) |
Return in A a parallel (on truedofs) version of this square operator. | |
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). | |
void | FormDiscreteOperator (Operator *&A) |
Return in A a parallel (on truedofs) version of this rectangular operator. | |
void | PrintMatlab (std::ostream &out, int n, int m=0) const |
Prints operator with input size n and output size m in Matlab format. | |
virtual void | PrintMatlab (std::ostream &out) const |
Prints operator in Matlab format. | |
virtual | ~Operator () |
Virtual destructor. | |
Type | GetType () const |
Return the type ID of the Operator class. | |
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... | |
![]() | |
void | FormConstrainedSystemOperator (const Array< int > &ess_tdof_list, ConstrainedOperator *&Aout) |
see FormSystemOperator() | |
void | FormRectangularConstrainedSystemOperator (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, RectangularConstrainedOperator *&Aout) |
see FormRectangularSystemOperator() | |
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". | |
Class for assembling of bilinear forms a(u,v)
defined on different trial and test spaces. The assembled matrix M
is such that
a(u,v) = V^t M U
where U
and V
are the vectors representing the functions u
and v
, respectively. The first argument, u
, of a(,)
is in the trial space and the second argument, v
, is in the test space. Thus,
# of rows of M = dimension of the test space and # of cols of M = dimension of the trial space.
Both trial and test spaces should be defined on the same mesh.
Definition at line 755 of file bilinearform.hpp.
mfem::MixedBilinearForm::MixedBilinearForm | ( | FiniteElementSpace * | tr_fes, |
FiniteElementSpace * | te_fes ) |
Construct a MixedBilinearForm on the given trial, tr_fes, and test, te_fes, FiniteElementSpaces.
The pointers tr_fes and te_fes are not owned by the newly constructed object.
Definition at line 1308 of file bilinearform.cpp.
mfem::MixedBilinearForm::MixedBilinearForm | ( | FiniteElementSpace * | tr_fes, |
FiniteElementSpace * | te_fes, | ||
MixedBilinearForm * | mbf ) |
Create a MixedBilinearForm on the given trial, tr_fes, and test, te_fes, FiniteElementSpaces, using the same integrators as the MixedBilinearForm mbf.
The pointers tr_fes and te_fes are not owned by the newly constructed object.
The integrators in mbf are copied as pointers and they are not owned by the newly constructed MixedBilinearForm.
Definition at line 1321 of file bilinearform.cpp.
|
virtual |
Deletes internal matrices, bilinear integrators, and the BilinearFormExtension.
Definition at line 2353 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddBdrFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a boundary face integrator. Assumes ownership of bfi.
Definition at line 1492 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddBdrFaceIntegrator | ( | BilinearFormIntegrator * | bfi, |
Array< int > & | bdr_marker ) |
Adds a boundary face integrator. Assumes ownership of bfi.
Definition at line 1498 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddBdrTraceFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a boundary trace face integrator. Assumes ownership of bfi.
Definition at line 1510 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddBdrTraceFaceIntegrator | ( | BilinearFormIntegrator * | bfi, |
Array< int > & | bdr_marker ) |
Adds a boundary trace face integrator. Assumes ownership of bfi.
Definition at line 1517 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddBoundaryIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a boundary integrator. Assumes ownership of bfi.
Definition at line 1474 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddBoundaryIntegrator | ( | BilinearFormIntegrator * | bfi, |
Array< int > & | bdr_marker ) |
Adds a boundary integrator. Assumes ownership of bfi.
Definition at line 1480 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddDomainIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a domain integrator. Assumes ownership of bfi.
Definition at line 1461 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddDomainIntegrator | ( | BilinearFormIntegrator * | bfi, |
Array< int > & | elem_marker ) |
Adds a domain integrator. Assumes ownership of bfi.
Definition at line 1467 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddInteriorFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds an interior face integrator. Assumes ownership of bfi.
Definition at line 1487 of file bilinearform.cpp.
|
overridevirtual |
Add the matrix vector multiple to a vector: \( y += a M x \).
Reimplemented from mfem::Operator.
Definition at line 1395 of file bilinearform.cpp.
|
overridevirtual |
Add the matrix transpose vector multiplication: \( y += a M^T x \).
Reimplemented from mfem::Operator.
Definition at line 1414 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AddTraceFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Add a trace face integrator. Assumes ownership of bfi.
This type of integrator assembles terms over all faces of the mesh using the face FE from the trial space and the two adjacent volume FEs from the test space.
Definition at line 1505 of file bilinearform.cpp.
void mfem::MixedBilinearForm::Assemble | ( | int | skip_zeros = 1 | ) |
Definition at line 1524 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AssembleBdrElementMatrix | ( | int | i, |
const DenseMatrix & | elmat, | ||
Array< int > & | trial_vdofs, | ||
Array< int > & | test_vdofs, | ||
int | skip_zeros = 1 ) |
Assemble the given boundary element matrix.
The boundary element matrix elmat is assembled for the boundary element i, i.e. added to the system matrix. The vdofs of the element are returned in trial_vdofs and test_vdofs. The flag skip_zeros skips the zero elements of the matrix, unless they are breaking the symmetry of the system matrix.
Definition at line 2177 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AssembleBdrElementMatrix | ( | int | i, |
const DenseMatrix & | elmat, | ||
int | skip_zeros = 1 ) |
Assemble the given boundary element matrix.
The boundary element matrix elmat is assembled for the boundary element i, i.e. added to the system matrix. The flag skip_zeros skips the zero elements of the matrix, unless they are breaking the symmetry of the system matrix.
Definition at line 2171 of file bilinearform.cpp.
Assemble the diagonal of ADA^T into diag, where A is this mixed bilinear form and D is a diagonal.
Definition at line 1836 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AssembleElementMatrix | ( | int | i, |
const DenseMatrix & | elmat, | ||
Array< int > & | trial_vdofs, | ||
Array< int > & | test_vdofs, | ||
int | skip_zeros = 1 ) |
Assemble the given element matrix.
The element matrix elmat is assembled for the element i, i.e. added to the system matrix. The vdofs of the element are returned in trial_vdofs and test_vdofs. The flag skip_zeros skips the zero elements of the matrix, unless they are breaking the symmetry of the system matrix.
Definition at line 2158 of file bilinearform.cpp.
void mfem::MixedBilinearForm::AssembleElementMatrix | ( | int | i, |
const DenseMatrix & | elmat, | ||
int | skip_zeros = 1 ) |
Assemble the given element matrix.
The element matrix elmat is assembled for the element i, i.e. added to the system matrix. The flag skip_zeros skips the zero elements of the matrix, unless they are breaking the symmetry of the system matrix.
Definition at line 2152 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ComputeBdrElementMatrix | ( | int | i, |
DenseMatrix & | elmat ) const |
Compute the boundary element matrix of the given boundary element.
Definition at line 1944 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ComputeBdrFaceMatrix | ( | int | i, |
DenseMatrix & | elmat ) const |
Compute the boundary face matrix of the given boundary element.
Definition at line 2024 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ComputeBdrTraceFaceMatrix | ( | int | i, |
DenseMatrix & | elmat ) const |
Compute the boundary trace face matrix of the given boundary element.
Definition at line 2112 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ComputeElementMatrix | ( | int | i, |
DenseMatrix & | elmat ) const |
Compute the element matrix of the given element.
Definition at line 1917 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ComputeFaceMatrix | ( | int | i, |
DenseMatrix & | elmat ) const |
Compute the face matrix of the given face element.
Definition at line 1971 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ComputeTraceFaceMatrix | ( | int | i, |
DenseMatrix & | elmat ) const |
Compute the trace face matrix of the given face element.
Definition at line 2064 of file bilinearform.cpp.
void mfem::MixedBilinearForm::ConformingAssemble | ( | ) |
For partially conforming trial and/or test FE spaces, complete the assembly process by performing \( P2^t A P1 \) where \( A \) is the internal sparse matrix; \( P1 \) and \( P2 \) are the conforming prolongation matrices of the trial and test FE spaces, respectively. After this call the MixedBilinearForm becomes an operator on the conforming FE spaces.
Definition at line 1884 of file bilinearform.cpp.
|
overridevirtual |
Returns a reference to: \( M_{ij} \).
Implements mfem::Matrix.
Definition at line 1384 of file bilinearform.cpp.
|
overridevirtual |
Returns a reference to: \( M_{ij} \).
Implements mfem::Matrix.
Definition at line 1379 of file bilinearform.cpp.
void mfem::MixedBilinearForm::EliminateEssentialBCFromTrialDofs | ( | const Array< int > & | marked_vdofs, |
const Vector & | sol, | ||
Vector & | rhs ) |
Similar to EliminateTrialVDofs(const Array<int> &, const Vector &, Vector &) but here ess_dofs is a marker (boolean) array on all vector-dofs (ess_dofs[i] < 0 is true).
Definition at line 2235 of file bilinearform.cpp.
|
inlinevirtual |
(DEPRECATED) Eliminate essential boundary test DOFs from the system.
Definition at line 1128 of file bilinearform.hpp.
void mfem::MixedBilinearForm::EliminateTestEssentialBC | ( | const Array< int > & | bdr_attr_is_ess | ) |
Eliminate essential boundary test DOFs from the system matrix.
The array bdr_attr_is_ess marks boundary attributes that constitute the essential part of the boundary.
Definition at line 2241 of file bilinearform.cpp.
void mfem::MixedBilinearForm::EliminateTestVDofs | ( | const Array< int > & | vdofs | ) |
Eliminate the given test vdofs. NOTE: here, vdofs is a list of DOFs.
Definition at line 2262 of file bilinearform.cpp.
|
inline |
(DEPRECATED) Eliminate essential boundary trial DOFs from the system.
Definition at line 1092 of file bilinearform.hpp.
void mfem::MixedBilinearForm::EliminateTrialEssentialBC | ( | const Array< int > & | bdr_attr_is_ess | ) |
Eliminate essential boundary trial DOFs from the system matrix.
The array bdr_attr_is_ess marks boundary attributes that constitute the essential part of the boundary.
Definition at line 2198 of file bilinearform.cpp.
void mfem::MixedBilinearForm::EliminateTrialEssentialBC | ( | const Array< int > & | bdr_attr_is_ess, |
const Vector & | sol, | ||
Vector & | rhs ) |
Eliminate essential boundary trial DOFs from the system.
The array bdr_attr_is_ess marks boundary attributes that constitute the essential part of the boundary.
Definition at line 2190 of file bilinearform.cpp.
void mfem::MixedBilinearForm::EliminateTrialVDofs | ( | const Array< int > & | vdofs | ) |
Eliminate the given trial vdofs, storing the eliminated part internally in \( M_e \).
This method works in conjunction with EliminateTrialVDofsInRHS() and allows elimination of boundary conditions in multiple right-hand sides. In this method, vdofs is a list of DOFs.
Definition at line 2215 of file bilinearform.cpp.
void mfem::MixedBilinearForm::EliminateTrialVDofs | ( | const Array< int > & | vdofs, |
const Vector & | sol, | ||
Vector & | rhs ) |
Eliminate the given trial vdofs. NOTE: here, vdofs is a list of DOFs.
In this case the eliminations are applied to the internal \( M \) and rhs without storing the elimination matrix \( M_e \).
Definition at line 2206 of file bilinearform.cpp.
void mfem::MixedBilinearForm::EliminateTrialVDofsInRHS | ( | const Array< int > & | vdofs, |
const Vector & | x, | ||
Vector & | b ) |
Use the stored eliminated part of the matrix (see EliminateTrialVDofs(const Array<int> &)) to modify the r.h.s. b; vdofs is a list of DOFs (non-directional, i.e. >= 0).
Definition at line 2229 of file bilinearform.cpp.
|
overridevirtual |
Finalizes the matrix initialization if the AssemblyLevel is AssemblyLevel::LEGACY.
Reimplemented from mfem::Matrix.
Definition at line 1441 of file bilinearform.cpp.
|
virtual |
Form the linear system A X = B, corresponding to this mixed bilinear form and the linear form b(.).
Return in A a reference to the system matrix that is column-constrained. The reference will be invalidated when SetOperatorType(), Update(), or the destructor is called.
Reimplemented in mfem::ParMixedBilinearForm.
Definition at line 2312 of file bilinearform.cpp.
|
inline |
Form the linear system A X = B, corresponding to this bilinear form and the linear form b(.).
Version of the method FormRectangularLinearSystem() where the system matrix is returned in the variable A, of type OpType, holding a reference* to the system matrix (created with the method OpType::MakeRef()). The reference will be invalidated when SetOperatorType(), Update(), or the destructor is called.
Definition at line 1182 of file bilinearform.hpp.
|
virtual |
Return in A that is column-constrained.
This returns the same operator as FormRectangularLinearSystem(), but does without the transformations of the right-hand side.
Reimplemented in mfem::ParDiscreteLinearOperator, and mfem::ParMixedBilinearForm.
Definition at line 2270 of file bilinearform.cpp.
|
inline |
Form the column-constrained linear system matrix A.
Version of the method FormRectangularSystemMatrix() where the system matrix is returned in the variable A, of type OpType, holding a reference* to the system matrix (created with the method OpType::MakeRef()). The reference will be invalidated when SetOperatorType(), Update(), or the destructor is called.
Definition at line 1151 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddBoundaryIntegrator().
Definition at line 949 of file bilinearform.hpp.
Access all boundary markers added with AddBoundaryIntegrator().
If no marker was specified when the integrator was added, the corresponding pointer (to Array<int>) will be NULL.
Definition at line 955 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddBdrFaceIntegrator().
Definition at line 961 of file bilinearform.hpp.
Access all boundary markers added with AddBdrFaceIntegrator(). If no marker was specified when the integrator was added, the corresponding pointer (to Array<int>) will be NULL.
Definition at line 965 of file bilinearform.hpp.
void mfem::MixedBilinearForm::GetBlocks | ( | Array2D< SparseMatrix * > & | blocks | ) | const |
Extract the associated matrix as SparseMatrix blocks. The number of block rows and columns is given by the vector dimensions (vdim) of the test and trial spaces, respectively.
Definition at line 1449 of file bilinearform.cpp.
|
inline |
Access all integrators added with AddBdrTraceFaceIntegrator().
Definition at line 971 of file bilinearform.hpp.
Access all boundary markers added with AddBdrTraceFaceIntegrator()
If no marker was specified when the integrator was added, the corresponding pointer (to Array<int>) will be NULL.
Definition at line 978 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddDomainIntegrator().
Definition at line 942 of file bilinearform.hpp.
Access all domain markers added with AddDomainIntegrator(). If no marker was specified when the integrator was added, the corresponding pointer (to Array<int>) will be NULL.
Definition at line 946 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddInteriorFaceIntegrator().
Definition at line 958 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the test finite element space prolongation matrix.
Reimplemented from mfem::Operator.
Definition at line 1003 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the test finite element space restriction matrix.
Reimplemented from mfem::Operator.
Definition at line 1007 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the input finite element space prolongation matrix.
Reimplemented from mfem::Operator.
Definition at line 995 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the input finite element space restriction matrix.
Reimplemented from mfem::Operator.
Definition at line 999 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddTraceFaceIntegrator().
Definition at line 968 of file bilinearform.hpp.
|
overridevirtual |
Returns a pointer to (approximation) of the matrix inverse: \( M^{-1} \) (currently unimplemented and returns NULL)
Implements mfem::Matrix.
Definition at line 1427 of file bilinearform.cpp.
|
inline |
Nullifies the internal matrix \( M \) and returns a pointer to it. Used for transferring ownership.
Definition at line 887 of file bilinearform.hpp.
Matrix multiplication: \( y = M x \).
Implements mfem::Operator.
Definition at line 1389 of file bilinearform.cpp.
Matrix transpose vector multiplication: \( y = M^T x \).
Reimplemented from mfem::Operator.
Definition at line 1408 of file bilinearform.cpp.
|
inline |
Sets all sparse values of \( M \) to a.
Definition at line 982 of file bilinearform.hpp.
|
inlineoverridevirtual |
Prints operator to stream os.
Reimplemented from mfem::Matrix.
Definition at line 1211 of file bilinearform.hpp.
void mfem::MixedBilinearForm::SetAssemblyLevel | ( | AssemblyLevel | assembly_level | ) |
Set the desired assembly level. The default is AssemblyLevel::LEGACY.
This method must be called before assembly. See AssemblyLevel
Definition at line 1348 of file bilinearform.cpp.
|
inline |
Returns a reference to the sparse matrix: \( M \).
Definition at line 879 of file bilinearform.hpp.
|
inline |
Returns a const reference to the sparse matrix: \( M \).
This will segfault if the usual sparse mat is not defined like when static condensation is being used or AllocMat() has not yet been called.
Definition at line 872 of file bilinearform.hpp.
|
inline |
Returns a reference to the sparse matrix of eliminated b.c.: \( M_e \).
Definition at line 897 of file bilinearform.hpp.
|
inline |
Returns a const reference to the sparse matrix of eliminated b.c.: \( M_e \).
Definition at line 890 of file bilinearform.hpp.
|
inline |
Return the test FE space associated with the BilinearForm.
Definition at line 1205 of file bilinearform.hpp.
|
inline |
Read-only access to the associated test FiniteElementSpace.
Definition at line 1208 of file bilinearform.hpp.
|
inline |
Return the trial FE space associated with the BilinearForm.
Definition at line 1199 of file bilinearform.hpp.
|
inline |
Read-only access to the associated trial FiniteElementSpace.
Definition at line 1202 of file bilinearform.hpp.
void mfem::MixedBilinearForm::Update | ( | ) |
Must be called after making changes to trial_fes or test_fes.
Definition at line 2342 of file bilinearform.cpp.
|
protected |
The form assembly level (full, partial, etc.)
Definition at line 765 of file bilinearform.hpp.
|
protected |
Boundary face integrators.
Definition at line 792 of file bilinearform.hpp.
Entries are not owned.
Definition at line 794 of file bilinearform.hpp.
|
protected |
Boundary integrators.
Definition at line 784 of file bilinearform.hpp.
Entries are not owned.
Definition at line 786 of file bilinearform.hpp.
|
protected |
Boundary trace face (skeleton) integrators.
Definition at line 800 of file bilinearform.hpp.
Entries are not owned.
Definition at line 802 of file bilinearform.hpp.
|
protected |
Domain integrators.
Definition at line 779 of file bilinearform.hpp.
Entries are not owned.
Definition at line 781 of file bilinearform.hpp.
|
mutableprotected |
Definition at line 804 of file bilinearform.hpp.
|
protected |
Extension for supporting Full Assembly (FA), Element Assembly (EA), Partial Assembly (PA), or Matrix Free assembly (MF).
Definition at line 769 of file bilinearform.hpp.
|
protected |
Indicates the BilinearFormIntegrators stored in MixedBilinearForm::domain_integs, MixedBilinearForm::boundary_integs, MixedBilinearForm::trace_face_integs and MixedBilinearForm::boundary_trace_face_integs are owned by another MixedBilinearForm.
Definition at line 776 of file bilinearform.hpp.
|
protected |
Interior face integrators.
Definition at line 789 of file bilinearform.hpp.
|
protected |
Owned.
Definition at line 758 of file bilinearform.hpp.
|
protected |
Owned.
Definition at line 759 of file bilinearform.hpp.
|
protected |
Not owned.
Definition at line 762 of file bilinearform.hpp.
|
protected |
Definition at line 805 of file bilinearform.hpp.
|
protected |
Trace face (skeleton) integrators.
Definition at line 797 of file bilinearform.hpp.
|
protected |
Not owned.
Definition at line 761 of file bilinearform.hpp.
|
mutableprotected |
Definition at line 805 of file bilinearform.hpp.