![]() |
MFEM v4.10.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. | |
Public Member Functions inherited from mfem::Matrix | |
| 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. | |
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. | |
| 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 | AbsMult (const Vector &x, Vector &y) const |
Action of the absolute-value operator: y=|A|(x). The default behavior in class Operator is to generate an error. If the Operator is a composition of several operators, the composition unfold into a product of absolute-value operators too. | |
| virtual void | AbsMultTranspose (const Vector &x, Vector &y) const |
Action of the transpose absolute-value operator: y=|A|^t(x). The default behavior in class Operator is to generate an error. | |
| 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 void | MultMV (const MultiVector &x, MultiVector &y) const |
| Operator application, y = A(x), where the input x and the output y are MultiVector objects, i.e. they generally use non-contiguous memory representation. | |
| virtual void | MultTransposeMV (const MultiVector &x, MultiVector &y) const |
| Action of the transpose operator, y = A^t(x), where the input x and the output y are MultiVector objects, i.e. they generally use non-contiguous memory representation. | |
| 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 Operator & | GetGradientMV (const MultiVector &x) const |
| Evaluate the gradient operator at the point x. The input x is provided as a MultiVector, i.e. it generally uses non-contiguous memory representation. | |
| 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 | |
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::Operator | |
| 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 754 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 1359 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 1372 of file bilinearform.cpp.
|
virtual |
Deletes internal matrices, bilinear integrators, and the BilinearFormExtension.
Definition at line 2400 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::AddBdrFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a boundary face integrator. Assumes ownership of bfi.
Definition at line 1543 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 1549 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::AddBdrTraceFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a boundary trace face integrator. Assumes ownership of bfi.
Definition at line 1561 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 1568 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::AddBoundaryIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a boundary integrator. Assumes ownership of bfi.
Definition at line 1525 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 1531 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::AddDomainIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds a domain integrator. Assumes ownership of bfi.
Definition at line 1512 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 1518 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::AddInteriorFaceIntegrator | ( | BilinearFormIntegrator * | bfi | ) |
Adds an interior face integrator. Assumes ownership of bfi.
Definition at line 1538 of file bilinearform.cpp.
|
overridevirtual |
Add the matrix vector multiple to a vector: \( y += a M x \).
Reimplemented from mfem::Operator.
Definition at line 1446 of file bilinearform.cpp.
|
overridevirtual |
Add the matrix transpose vector multiplication: \( y += a M^T x \).
Reimplemented from mfem::Operator.
Definition at line 1465 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 1556 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::Assemble | ( | int | skip_zeros = 1 | ) |
Definition at line 1575 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 2224 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 2218 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 1883 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 2205 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 2199 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 1991 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 2071 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 2159 of file bilinearform.cpp.
| void mfem::MixedBilinearForm::ComputeElementMatrix | ( | int | i, |
| DenseMatrix & | elmat ) const |
Compute the element matrix of the given element.
Definition at line 1964 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 2018 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 2111 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 1931 of file bilinearform.cpp.
|
overridevirtual |
Returns a reference to: \( M_{ij} \).
Implements mfem::Matrix.
Definition at line 1435 of file bilinearform.cpp.
|
overridevirtual |
Returns a reference to: \( M_{ij} \).
Implements mfem::Matrix.
Definition at line 1430 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 2282 of file bilinearform.cpp.
|
inlinevirtual |
(DEPRECATED) Eliminate essential boundary test DOFs from the system.
Definition at line 1127 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 2288 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 2309 of file bilinearform.cpp.
|
inline |
(DEPRECATED) Eliminate essential boundary trial DOFs from the system.
Definition at line 1091 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 2245 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 2237 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 2262 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 2253 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 2276 of file bilinearform.cpp.
|
overridevirtual |
Finalizes the matrix initialization if the AssemblyLevel is AssemblyLevel::LEGACY.
Reimplemented from mfem::Matrix.
Definition at line 1492 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 2359 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 1181 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 2317 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 1150 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddBoundaryIntegrator().
Definition at line 948 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 954 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddBdrFaceIntegrator().
Definition at line 960 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 964 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 1500 of file bilinearform.cpp.
|
inline |
Access all integrators added with AddBdrTraceFaceIntegrator().
Definition at line 970 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 977 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddDomainIntegrator().
Definition at line 941 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 945 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddInteriorFaceIntegrator().
Definition at line 957 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the test finite element space prolongation matrix.
Reimplemented from mfem::Operator.
Definition at line 1002 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the test finite element space restriction matrix.
Reimplemented from mfem::Operator.
Definition at line 1006 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the input finite element space prolongation matrix.
Reimplemented from mfem::Operator.
Definition at line 994 of file bilinearform.hpp.
|
inlineoverridevirtual |
Get the input finite element space restriction matrix.
Reimplemented from mfem::Operator.
Definition at line 998 of file bilinearform.hpp.
|
inline |
Access all integrators added with AddTraceFaceIntegrator().
Definition at line 967 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 1478 of file bilinearform.cpp.
|
inline |
Nullifies the internal matrix \( M \) and returns a pointer to it. Used for transferring ownership.
Definition at line 886 of file bilinearform.hpp.
Matrix multiplication: \( y = M x \).
Implements mfem::Operator.
Definition at line 1440 of file bilinearform.cpp.
Matrix transpose vector multiplication: \( y = M^T x \).
Reimplemented from mfem::Operator.
Definition at line 1459 of file bilinearform.cpp.
|
inline |
Sets all sparse values of \( M \) to a.
Definition at line 981 of file bilinearform.hpp.
|
inlineoverridevirtual |
Prints operator to stream os.
Reimplemented from mfem::Matrix.
Definition at line 1210 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 1399 of file bilinearform.cpp.
|
inline |
Returns a reference to the sparse matrix: \( M \).
Definition at line 878 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 871 of file bilinearform.hpp.
|
inline |
Returns a reference to the sparse matrix of eliminated b.c.: \( M_e \).
Definition at line 896 of file bilinearform.hpp.
|
inline |
Returns a const reference to the sparse matrix of eliminated b.c.: \( M_e \).
Definition at line 889 of file bilinearform.hpp.
|
inline |
Return the test FE space associated with the BilinearForm.
Definition at line 1204 of file bilinearform.hpp.
|
inline |
Read-only access to the associated test FiniteElementSpace.
Definition at line 1207 of file bilinearform.hpp.
|
inline |
Return the trial FE space associated with the BilinearForm.
Definition at line 1198 of file bilinearform.hpp.
|
inline |
Read-only access to the associated trial FiniteElementSpace.
Definition at line 1201 of file bilinearform.hpp.
| void mfem::MixedBilinearForm::Update | ( | ) |
Must be called after making changes to trial_fes or test_fes.
Definition at line 2389 of file bilinearform.cpp.
|
protected |
The form assembly level (full, partial, etc.)
Definition at line 764 of file bilinearform.hpp.
|
protected |
Boundary face integrators.
Definition at line 791 of file bilinearform.hpp.
Entries are not owned.
Definition at line 793 of file bilinearform.hpp.
|
protected |
Boundary integrators.
Definition at line 783 of file bilinearform.hpp.
Entries are not owned.
Definition at line 785 of file bilinearform.hpp.
|
protected |
Boundary trace face (skeleton) integrators.
Definition at line 799 of file bilinearform.hpp.
Entries are not owned.
Definition at line 801 of file bilinearform.hpp.
|
protected |
Domain integrators.
Definition at line 778 of file bilinearform.hpp.
Entries are not owned.
Definition at line 780 of file bilinearform.hpp.
|
mutableprotected |
Definition at line 803 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 768 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 775 of file bilinearform.hpp.
|
protected |
Interior face integrators.
Definition at line 788 of file bilinearform.hpp.
|
protected |
Owned.
Definition at line 757 of file bilinearform.hpp.
|
protected |
Owned.
Definition at line 758 of file bilinearform.hpp.
|
protected |
Not owned.
Definition at line 761 of file bilinearform.hpp.
|
protected |
Definition at line 804 of file bilinearform.hpp.
|
protected |
Trace face (skeleton) integrators.
Definition at line 796 of file bilinearform.hpp.
|
protected |
Not owned.
Definition at line 760 of file bilinearform.hpp.
|
mutableprotected |
Definition at line 804 of file bilinearform.hpp.