Class representing the parallel weak formulation. (Convenient for DPG Equations) More...

#include <pweakform.hpp>

Inheritance diagram for mfem::ParDPGWeakForm:

Collaboration diagram for mfem::ParDPGWeakForm:

Public Member Functions
	ParDPGWeakForm ()

	ParDPGWeakForm (Array< ParFiniteElementSpace * > &trial_pfes_, Array< FiniteElementCollection * > &fecol_)
	Creates bilinear form associated with FE spaces trial_pfes_.

void	SetParSpaces (Array< ParFiniteElementSpace * > &trial_pfes_, Array< FiniteElementCollection * > &fecol_)

void	Assemble (int skip_zeros=1)
	Assemble the local matrix.

void	ParallelAssemble (BlockMatrix *mat)
	Returns the matrix assembled on the true dofs, i.e. P^t A P.

void	FormLinearSystem (const Array< int > &ess_tdof_list, Vector &x, OperatorHandle &A, Vector &X, Vector &B, int copy_interior=0)
	Form the linear system A X = B, corresponding to this DPG weak form.

void	FormSystemMatrix (const Array< int > &ess_tdof_list, OperatorHandle &A)
	Form the linear system matrix A, see FormLinearSystem() for details.

virtual void	RecoverFEMSolution (const Vector &X, Vector &x)

virtual void	Update ()
	Update the DPGWeakForm after mesh modifications (AMR)

virtual	~ParDPGWeakForm ()
	Destroys bilinear form.

Public Member Functions inherited from mfem::DPGWeakForm
	DPGWeakForm ()
	Default constructor. User must call SetSpaces to setup the FE spaces.

	DPGWeakForm (Array< FiniteElementSpace * > &fes_, Array< FiniteElementCollection * > &fecol_)
	Creates bilinear form associated with FE spaces fes_.

void	SetTestFECollVdim (int test_fec, int vdim)

void	SetSpaces (Array< FiniteElementSpace * > &fes_, Array< FiniteElementCollection * > &fecol_)

int	Size () const
	Get the size of the bilinear form of the DPGWeakForm.

void	AllocateMatrix ()
	Pre-allocate the internal BlockMatrix before assembly.

void	Finalize (int skip_zeros=1)
	Finalizes the matrix initialization.

BlockMatrix &	BlockMat ()
	Returns a reference to the BlockMatrix: \( M \).

BlockMatrix &	BlockMatElim ()
	Returns a reference to the sparse matrix of eliminated b.c.: \( M_e \).

void	AddTrialIntegrator (BilinearFormIntegrator *bfi, int n, int m)
	Adds new Trial Integrator. Assumes ownership of bfi.

void	AddTestIntegrator (BilinearFormIntegrator *bfi, int n, int m)
	Adds new Test Integrator. Assumes ownership of bfi.

void	AddDomainLFIntegrator (LinearFormIntegrator *lfi, int n)
	Adds new Domain LF Integrator. Assumes ownership of bfi.

void	Assemble (int skip_zeros=1)
	Assembles the form i.e. sums over all integrators.

template<typename OpType >
void	FormLinearSystem (const Array< int > &ess_tdof_list, Vector &x, OpType &A, Vector &X, Vector &B, int copy_interior=0)
	Form the linear system A X = B, corresponding to this DPG weak form Version of the method FormLinearSystem() 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()).

template<typename OpType >
void	FormSystemMatrix (const Array< int > &ess_tdof_list, OpType &A)
	Form the linear system matrix A, see FormLinearSystem() for details.

void	EliminateVDofs (const Array< int > &vdofs, Operator::DiagonalPolicy dpolicy=Operator::DIAG_ONE)
	Eliminate the given vdofs, storing the eliminated part internally in \( M_e \).

void	EliminateVDofsInRHS (const Array< int > &vdofs, const Vector &x, Vector &b)
	Use the stored eliminated part of the matrix (see EliminateVDofs(const Array<int> &, DiagonalPolicy)) to modify the r.h.s. b; vdofs is a list of DOFs.

void	SetDiagonalPolicy (Operator::DiagonalPolicy policy)
	Sets diagonal policy used upon construction of the linear system.

void	StoreMatrices (bool store_matrices_=true)
	Store internal element matrices used for computation of residual after solve.

void	EnableStaticCondensation ()

Vector &	ComputeResidual (const BlockVector &x)
	Compute DPG residual based error estimator.

virtual	~DPGWeakForm ()

Protected Member Functions
void	FillEssTdofLists (const Array< int > &ess_tdof_list)

void	BuildProlongation ()

Protected Member Functions inherited from mfem::DPGWeakForm
void	Init ()

void	ReleaseInitMemory ()

void	AllocMat ()
	Allocate appropriate BlockMatrix and assign it to mat.

void	ConformingAssemble ()

void	ComputeOffsets ()

Protected Attributes
Array< ParFiniteElementSpace * >	trial_pfes

Array< Array< int > * >	ess_tdofs

BlockOperator *	P = nullptr

BlockMatrix *	R = nullptr

BlockOperator *	p_mat = nullptr

BlockOperator *	p_mat_e = nullptr

Protected Attributes inherited from mfem::DPGWeakForm
BlockStaticCondensation *	static_cond
	Owned.

bool	initialized = false

Mesh *	mesh = nullptr

int	height

int	width

int	nblocks

Array< int >	dof_offsets

Array< int >	tdof_offsets

BlockMatrix *	mat = nullptr
	Block matrix \( M \) to be associated with the Block bilinear form. Owned.

BlockVector *	y = nullptr
	Block vector \( y \) to be associated with the Block linear form.

BlockMatrix *	mat_e = nullptr
	Block Matrix \( M_e \) used to store the eliminations from the b.c. Owned. \( M + M_e = M_{original} \).

Array< FiniteElementSpace * >	trial_fes
	Trial FE spaces.

Array< int >	IsTraceFes
	Flags to determine if a FiniteElementSpace is Trace.

Array< FiniteElementCollection * >	test_fecols
	Test FE Collections (Broken)

Array< int >	test_fecols_vdims

Array2D< Array< BilinearFormIntegrator * > * >	trial_integs
	Set of Trial Integrators to be applied for matrix B.

Array2D< Array< BilinearFormIntegrator * > * >	test_integs
	Set of Test Space (broken) Integrators to be applied for matrix G.

Array< Array< LinearFormIntegrator * > * >	lfis
	Set of Linear Form Integrators to be applied.

BlockMatrix *	P = nullptr
	Block Prolongation.

BlockMatrix *	R = nullptr
	Block Restriction.

mfem::Operator::DiagonalPolicy	diag_policy

bool	store_matrices = false

Array< DenseMatrix * >	Bmat

Array< Vector * >	fvec

Vector	residuals

Detailed Description

Class representing the parallel weak formulation. (Convenient for DPG Equations)

Definition at line 25 of file pweakform.hpp.

Constructor & Destructor Documentation

◆ ParDPGWeakForm() [1/2]

mfem::ParDPGWeakForm::ParDPGWeakForm ( )

inline

Definition at line 53 of file pweakform.hpp.

◆ ParDPGWeakForm() [2/2]

mfem::ParDPGWeakForm::ParDPGWeakForm	(	Array< ParFiniteElementSpace * > &	trial_pfes_,
		Array< FiniteElementCollection * > &	fecol_ )

inline

Creates bilinear form associated with FE spaces trial_pfes_.

Definition at line 56 of file pweakform.hpp.

◆ ~ParDPGWeakForm()

mfem::ParDPGWeakForm::~ParDPGWeakForm ( )

virtual

Destroys bilinear form.

Definition at line 206 of file pweakform.cpp.

Member Function Documentation

◆ Assemble()

void mfem::ParDPGWeakForm::Assemble ( int skip_zeros = 1 )

Assemble the local matrix.

Definition at line 33 of file pweakform.cpp.

◆ BuildProlongation()

void mfem::ParDPGWeakForm::BuildProlongation ( )

protectedvirtual

Reimplemented from mfem::DPGWeakForm.

Definition at line 84 of file pweakform.cpp.

◆ FillEssTdofLists()

void mfem::ParDPGWeakForm::FillEssTdofLists ( const Array< int > & ess_tdof_list )

protected

Definition at line 19 of file pweakform.cpp.

◆ FormLinearSystem()

void mfem::ParDPGWeakForm::FormLinearSystem	(	const Array< int > &	ess_tdof_list,
		Vector &	x,
		OperatorHandle &	A,
		Vector &	X,
		Vector &	B,
		int	copy_interior = 0 )

virtual

Form the linear system A X = B, corresponding to this DPG weak form.

This method applies any necessary transformations to the linear system such as: eliminating boundary conditions; applying conforming constraints for non-conforming AMR; static condensation;

The GridFunction-size vector x must contain the essential b.c. The DPGWeakForm must be assembled.

The vector X is initialized with a suitable initial guess: the essential entries of X are set to the corresponding b.c. and all other entries are set to zero (copy_interior == 0) or copied from x (copy_interior != 0).

After solving the linear system, the finite element solution x can be recovered by calling RecoverFEMSolution() (with the same vectors X, and x).

NOTE: If there are no transformations, X simply reuses the data of x.

Reimplemented from mfem::DPGWeakForm.

Definition at line 100 of file pweakform.cpp.

◆ FormSystemMatrix()

void mfem::ParDPGWeakForm::FormSystemMatrix	(	const Array< int > &	ess_tdof_list,
		OperatorHandle &	A )

virtual

Form the linear system matrix A, see FormLinearSystem() for details.

Reimplemented from mfem::DPGWeakForm.

Definition at line 140 of file pweakform.cpp.

◆ ParallelAssemble()

void mfem::ParDPGWeakForm::ParallelAssemble ( BlockMatrix * mat )

Returns the matrix assembled on the true dofs, i.e. P^t A P.

The returned matrix has to be deleted by the caller.

Definition at line 38 of file pweakform.cpp.

◆ RecoverFEMSolution()

void mfem::ParDPGWeakForm::RecoverFEMSolution	(	const Vector &	X,
		Vector &	x )

virtual

Call this method after solving a linear system constructed using the FormLinearSystem method to recover the solution as a ParGridFunction-size vector in x. Use the same arguments as in the FormLinearSystem call.

Reimplemented from mfem::DPGWeakForm.

Definition at line 173 of file pweakform.cpp.

◆ SetParSpaces()

void mfem::ParDPGWeakForm::SetParSpaces	(	Array< ParFiniteElementSpace * > &	trial_pfes_,
		Array< FiniteElementCollection * > &	fecol_ )

inline

Definition at line 63 of file pweakform.hpp.

◆ Update()

void mfem::ParDPGWeakForm::Update ( )

virtual

Update the DPGWeakForm after mesh modifications (AMR)

Reimplemented from mfem::DPGWeakForm.

Definition at line 188 of file pweakform.cpp.

Member Data Documentation

◆ ess_tdofs

Array<Array<int> *> mfem::ParDPGWeakForm::ess_tdofs

protected

Definition at line 33 of file pweakform.hpp.

◆ P

BlockOperator* mfem::ParDPGWeakForm::P = nullptr

protected

Definition at line 40 of file pweakform.hpp.

◆ p_mat

BlockOperator* mfem::ParDPGWeakForm::p_mat = nullptr

protected

Definition at line 44 of file pweakform.hpp.

◆ p_mat_e

BlockOperator* mfem::ParDPGWeakForm::p_mat_e = nullptr

protected

Definition at line 45 of file pweakform.hpp.

◆ R

BlockMatrix* mfem::ParDPGWeakForm::R = nullptr

protected

Definition at line 41 of file pweakform.hpp.

◆ trial_pfes

Array<ParFiniteElementSpace * > mfem::ParDPGWeakForm::trial_pfes

protected

Definition at line 30 of file pweakform.hpp.

The documentation for this class was generated from the following files:

miniapps/dpg/util/pweakform.hpp
miniapps/dpg/util/pweakform.cpp

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ParDPGWeakForm() [1/2]

◆ ParDPGWeakForm() [2/2]

◆ ~ParDPGWeakForm()

Member Function Documentation

◆ Assemble()

◆ BuildProlongation()

◆ FillEssTdofLists()

◆ FormLinearSystem()

◆ FormSystemMatrix()

◆ ParallelAssemble()

◆ RecoverFEMSolution()

◆ SetParSpaces()

◆ Update()

Member Data Documentation

◆ ess_tdofs

◆ P

◆ p_mat

◆ p_mat_e

◆ R

◆ trial_pfes