#include <bilininteg.hpp>

Inheritance diagram for mfem::DGDiffusionIntegrator:

Collaboration diagram for mfem::DGDiffusionIntegrator:

Public Member Functions
	DGDiffusionIntegrator (const real_t s, const real_t k)

	DGDiffusionIntegrator (Coefficient &q, const real_t s, const real_t k)

	DGDiffusionIntegrator (MatrixCoefficient &q, const real_t s, const real_t k)

void	AssembleFaceMatrix (const FiniteElement &el1, const FiniteElement &el2, FaceElementTransformations &Trans, DenseMatrix &elmat) override

bool	RequiresFaceNormalDerivatives () const override
	For bilinear forms on element faces, specifies if the normal derivatives are needed on the faces or just the face restriction.

void	AssemblePAInteriorFaces (const FiniteElementSpace &fes) override

void	AssemblePABoundaryFaces (const FiniteElementSpace &fes) override

void	AddMultPAFaceNormalDerivatives (const Vector &x, const Vector &dxdn, Vector &y, Vector &dydn) const override
	Method for partially assembled action.

const IntegrationRule &	GetRule (int order, FaceElementTransformations &T)

virtual void	AssembleFaceMatrix (const FiniteElement &trial_face_fe, const FiniteElement &test_fe1, const FiniteElement &test_fe2, FaceElementTransformations &Trans, DenseMatrix &elmat)

virtual void	AssemblePA (const FiniteElementSpace &fes)
	Method defining partial assembly.

virtual void	AssemblePA (const FiniteElementSpace &trial_fes, const FiniteElementSpace &test_fes)

Public Member Functions inherited from mfem::BilinearFormIntegrator
virtual void	AssembleNURBSPA (const FiniteElementSpace &fes)
	Method defining partial assembly on NURBS patches.

virtual void	AssemblePABoundary (const FiniteElementSpace &fes)

virtual void	AssembleDiagonalPA (Vector &diag)
	Assemble diagonal and add it to Vector diag.

virtual void	AssembleDiagonalPA_ADAt (const Vector &D, Vector &diag)
	Assemble diagonal of \(A D A^T\) ( \(A\) is this integrator) and add it to diag.

virtual void	AddMultPA (const Vector &x, Vector &y) const
	Method for partially assembled action.

virtual void	AddMultNURBSPA (const Vector &x, Vector &y) const
	Method for partially assembled action on NURBS patches.

virtual void	AddMultTransposePA (const Vector &x, Vector &y) const
	Method for partially assembled transposed action.

virtual void	AssembleEA (const FiniteElementSpace &fes, Vector &emat, const bool add=true)
	Method defining element assembly.

virtual void	AssembleMF (const FiniteElementSpace &fes)
	Method defining matrix-free assembly.

virtual void	AddMultMF (const Vector &x, Vector &y) const

virtual void	AddMultTransposeMF (const Vector &x, Vector &y) const

virtual void	AssembleDiagonalMF (Vector &diag)
	Assemble diagonal and add it to Vector diag.

virtual void	AssembleEAInteriorFaces (const FiniteElementSpace &fes, Vector &ea_data_int, Vector &ea_data_ext, const bool add=true)

virtual void	AssembleEABoundaryFaces (const FiniteElementSpace &fes, Vector &ea_data_bdr, const bool add=true)

virtual void	AssembleElementMatrix (const FiniteElement &el, ElementTransformation &Trans, DenseMatrix &elmat)
	Given a particular Finite Element computes the element matrix elmat.

virtual void	AssembleElementMatrix2 (const FiniteElement &trial_fe, const FiniteElement &test_fe, ElementTransformation &Trans, DenseMatrix &elmat)

virtual void	AssemblePatchMatrix (const int patch, const FiniteElementSpace &fes, SparseMatrix *&smat)

virtual void	AssembleTraceFaceMatrix (int elem, const FiniteElement &trial_face_fe, const FiniteElement &test_fe, FaceElementTransformations &Trans, DenseMatrix &elmat)

virtual void	AssembleElementVector (const FiniteElement &el, ElementTransformation &Tr, const Vector &elfun, Vector &elvect)
	Perform the local action of the BilinearFormIntegrator. Note that the default implementation in the base class is general but not efficient.

virtual void	AssembleFaceVector (const FiniteElement &el1, const FiniteElement &el2, FaceElementTransformations &Tr, const Vector &elfun, Vector &elvect)
	Perform the local action of the BilinearFormIntegrator resulting from a face integral term. Note that the default implementation in the base class is general but not efficient.

virtual void	AssembleElementGrad (const FiniteElement &el, ElementTransformation &Tr, const Vector &elfun, DenseMatrix &elmat)
	Assemble the local gradient matrix.

virtual void	AssembleFaceGrad (const FiniteElement &el1, const FiniteElement &el2, FaceElementTransformations &Tr, const Vector &elfun, DenseMatrix &elmat)
	Assemble the local action of the gradient of the NonlinearFormIntegrator resulting from a face integral term.

virtual void	ComputeElementFlux (const FiniteElement &el, ElementTransformation &Trans, Vector &u, const FiniteElement &fluxelem, Vector &flux, bool with_coef=true, const IntegrationRule *ir=NULL)
	Virtual method required for Zienkiewicz-Zhu type error estimators.

virtual real_t	ComputeFluxEnergy (const FiniteElement &fluxelem, ElementTransformation &Trans, Vector &flux, Vector *d_energy=NULL)
	Virtual method required for Zienkiewicz-Zhu type error estimators.

virtual	~BilinearFormIntegrator ()

Public Member Functions inherited from mfem::NonlinearFormIntegrator
virtual void	SetIntRule (const IntegrationRule *ir)
	Prescribe a fixed IntegrationRule to use (when ir != NULL) or let the integrator choose (when ir == NULL).

void	SetIntegrationMode (Mode m)

void	SetNURBSPatchIntRule (NURBSMeshRules *pr)
	For patchwise integration, SetNURBSPatchIntRule must be called.

bool	HasNURBSPatchIntRule () const

bool	Patchwise () const

void	SetIntegrationRule (const IntegrationRule &ir)
	Prescribe a fixed IntegrationRule to use.

void	SetPAMemoryType (MemoryType mt)

const IntegrationRule *	GetIntegrationRule () const
	Get the integration rule of the integrator (possibly NULL).

virtual real_t	GetElementEnergy (const FiniteElement &el, ElementTransformation &Tr, const Vector &elfun)
	Compute the local energy.

virtual void	AssembleGradPA (const Vector &x, const FiniteElementSpace &fes)
	Prepare the integrator for partial assembly (PA) gradient evaluations on the given FE space fes at the state x.

virtual real_t	GetLocalStateEnergyPA (const Vector &x) const
	Compute the local (to the MPI rank) energy with partial assembly.

virtual void	AddMultGradPA (const Vector &x, Vector &y) const
	Method for partially assembled gradient action.

virtual void	AssembleGradDiagonalPA (Vector &diag) const
	Method for computing the diagonal of the gradient with partial assembly.

virtual bool	SupportsCeed () const
	Indicates whether this integrator can use a Ceed backend.

ceed::Operator &	GetCeedOp ()

virtual	~NonlinearFormIntegrator ()

Protected Attributes
Coefficient *	Q

MatrixCoefficient *	MQ

real_t	sigma

real_t	kappa

Vector	shape1

Vector	shape2

Vector	dshape1dn

Vector	dshape2dn

Vector	nor

Vector	nh

Vector	ni

DenseMatrix	jmat

DenseMatrix	dshape1

DenseMatrix	dshape2

DenseMatrix	mq

DenseMatrix	adjJ

Vector	pa_data

const DofToQuad *	maps
	Not owned.

int	dim

int	nf

int	nq

int	dofs1D

int	quad1D

IntegrationRules	irs {0, Quadrature1D::GaussLobatto}

Protected Attributes inherited from mfem::NonlinearFormIntegrator
const IntegrationRule *	IntRule

Mode	integrationMode = Mode::ELEMENTWISE

NURBSMeshRules *	patchRules = nullptr

ceed::Operator *	ceedOp

MemoryType	pa_mt = MemoryType::DEFAULT

Additional Inherited Members
Public Types inherited from mfem::NonlinearFormIntegrator
enum	Mode { ELEMENTWISE = 0 , PATCHWISE = 1 , PATCHWISE_REDUCED = 2 }

Protected Member Functions inherited from mfem::BilinearFormIntegrator
	BilinearFormIntegrator (const IntegrationRule *ir=NULL)

Protected Member Functions inherited from mfem::NonlinearFormIntegrator
	NonlinearFormIntegrator (const IntegrationRule *ir=NULL)

Detailed Description

Integrator for the DG form:

\[ - \langle \{(Q \nabla u) \cdot n\}, [v] \rangle + \sigma \langle [u], \{(Q \nabla v) \cdot n \} \rangle + \kappa \langle \{h^{-1} Q\} [u], [v] \rangle \]

where \(Q\) is a scalar or matrix diffusion coefficient and \(u\), \(v\) are the trial and test spaces, respectively. The parameters \(\sigma\) and \(\kappa\) determine the DG method to be used (when this integrator is added to the "broken" DiffusionIntegrator):

\(\sigma = -1\), \(\kappa \geq \kappa_0\): symm. interior penalty (IP or SIPG) method,
\(\sigma = +1\), \(\kappa > 0\): non-symmetric interior penalty (NIPG) method,
\(\sigma = +1\), \(\kappa = 0\): the method of Baumann and Oden.

Definition at line 3254 of file bilininteg.hpp.

Constructor & Destructor Documentation

◆ DGDiffusionIntegrator() [1/3]

mfem::DGDiffusionIntegrator::DGDiffusionIntegrator	(	const real_t	s,
		const real_t	k )

inline

Definition at line 3273 of file bilininteg.hpp.

◆ DGDiffusionIntegrator() [2/3]

mfem::DGDiffusionIntegrator::DGDiffusionIntegrator	(	Coefficient &	q,
		const real_t	s,
		const real_t	k )

inline

Definition at line 3275 of file bilininteg.hpp.

◆ DGDiffusionIntegrator() [3/3]

mfem::DGDiffusionIntegrator::DGDiffusionIntegrator	(	MatrixCoefficient &	q,
		const real_t	s,
		const real_t	k )

inline

Definition at line 3277 of file bilininteg.hpp.

Member Function Documentation

◆ AddMultPAFaceNormalDerivatives()

void mfem::DGDiffusionIntegrator::AddMultPAFaceNormalDerivatives	(	const Vector &	x,
		const Vector &	dxdn,
		Vector &	y,
		Vector &	dydn ) const

overridevirtual

Method for partially assembled action.

For bilinear forms on element faces that depend on the normal derivative on the faces, computes the action of integrator to the face values x and reference-normal derivatives dxdn and adds the result to y and dydn.

This method can be called only after the method AssemblePA() has been called.

Parameters

[in]	x	E-vector of face values (provided by FaceRestriction::Mult)
[in]	dxdn	E-vector of face reference-normal derivatives (provided by FaceRestriction::NormalDerivativeMult)
[in,out]	y	E-vector of face values to add action to.
[in,out]	dydn	E-vector of face reference-normal derivative values to add action to.

Reimplemented from mfem::BilinearFormIntegrator.

◆ AssembleFaceMatrix() [1/2]

void mfem::DGDiffusionIntegrator::AssembleFaceMatrix	(	const FiniteElement &	el1,
		const FiniteElement &	el2,
		FaceElementTransformations &	Trans,
		DenseMatrix &	elmat )

overridevirtual

Reimplemented from mfem::BilinearFormIntegrator.

Definition at line 3428 of file bilininteg.cpp.

◆ AssembleFaceMatrix() [2/2]

void mfem::BilinearFormIntegrator::AssembleFaceMatrix	(	const FiniteElement &	trial_face_fe,
		const FiniteElement &	test_fe1,
		const FiniteElement &	test_fe2,
		FaceElementTransformations &	Trans,
		DenseMatrix &	elmat )

virtual

Abstract method used for assembling TraceFaceIntegrators in a MixedBilinearForm.

Reimplemented from mfem::BilinearFormIntegrator.

Definition at line 162 of file bilininteg.cpp.

◆ AssemblePA() [1/2]

void mfem::BilinearFormIntegrator::AssemblePA ( const FiniteElementSpace & fes )

virtual

Method defining partial assembly.

The result of the partial assembly is stored internally so that it can be used later in the methods AddMultPA() and AddMultTransposePA().

Reimplemented from mfem::BilinearFormIntegrator.

Definition at line 48 of file bilininteg.cpp.

◆ AssemblePA() [2/2]

void mfem::BilinearFormIntegrator::AssemblePA	(	const FiniteElementSpace &	trial_fes,
		const FiniteElementSpace &	test_fes )

virtual

Used with BilinearFormIntegrators that have different spaces.

Reimplemented from mfem::BilinearFormIntegrator.

Definition at line 50 of file bilininteg.cpp.

◆ AssemblePABoundaryFaces()

void mfem::DGDiffusionIntegrator::AssemblePABoundaryFaces ( const FiniteElementSpace & fes )

overridevirtual

Reimplemented from mfem::BilinearFormIntegrator.

◆ AssemblePAInteriorFaces()

void mfem::DGDiffusionIntegrator::AssemblePAInteriorFaces ( const FiniteElementSpace & fes )

overridevirtual

Reimplemented from mfem::BilinearFormIntegrator.

◆ GetRule()

const IntegrationRule & mfem::DGDiffusionIntegrator::GetRule	(	int	order,
		FaceElementTransformations &	T )

Definition at line 3655 of file bilininteg.cpp.

◆ RequiresFaceNormalDerivatives()

bool mfem::DGDiffusionIntegrator::RequiresFaceNormalDerivatives ( ) const

inlineoverridevirtual

For bilinear forms on element faces, specifies if the normal derivatives are needed on the faces or just the face restriction.

if RequiresFaceNormalDerivatives() == true, then AddMultPAFaceNormalDerivatives(...) should be invoked in place of AddMultPA(...) and L2NormalDerivativeFaceRestriction should be used to compute the normal derivatives. This is used for some DG integrators, for example DGDiffusionIntegrator.

Returns: whether normal derivatives appear in the bilinear form.

Reimplemented from mfem::BilinearFormIntegrator.

Definition at line 3285 of file bilininteg.hpp.

Member Data Documentation

◆ adjJ

DenseMatrix mfem::DGDiffusionIntegrator::adjJ

protected

Definition at line 3263 of file bilininteg.hpp.

◆ dim

int mfem::DGDiffusionIntegrator::dim

protected

Definition at line 3269 of file bilininteg.hpp.

◆ dofs1D

int mfem::DGDiffusionIntegrator::dofs1D

protected

Definition at line 3269 of file bilininteg.hpp.

◆ dshape1

DenseMatrix mfem::DGDiffusionIntegrator::dshape1

protected

Definition at line 3263 of file bilininteg.hpp.

◆ dshape1dn

Vector mfem::DGDiffusionIntegrator::dshape1dn

protected

Definition at line 3262 of file bilininteg.hpp.

◆ dshape2

DenseMatrix mfem::DGDiffusionIntegrator::dshape2

protected

Definition at line 3263 of file bilininteg.hpp.

◆ dshape2dn

Vector mfem::DGDiffusionIntegrator::dshape2dn

protected

Definition at line 3262 of file bilininteg.hpp.

◆ irs

IntegrationRules mfem::DGDiffusionIntegrator::irs {0, Quadrature1D::GaussLobatto}

protected

Definition at line 3270 of file bilininteg.hpp.

◆ jmat

DenseMatrix mfem::DGDiffusionIntegrator::jmat

protected

Definition at line 3263 of file bilininteg.hpp.

◆ kappa

real_t mfem::DGDiffusionIntegrator::kappa

protected

Definition at line 3259 of file bilininteg.hpp.

◆ maps

const DofToQuad* mfem::DGDiffusionIntegrator::maps

protected

Not owned.

Definition at line 3268 of file bilininteg.hpp.

◆ MQ

MatrixCoefficient* mfem::DGDiffusionIntegrator::MQ

protected

Definition at line 3258 of file bilininteg.hpp.

◆ mq

DenseMatrix mfem::DGDiffusionIntegrator::mq

protected

Definition at line 3263 of file bilininteg.hpp.

◆ nf

int mfem::DGDiffusionIntegrator::nf

protected

Definition at line 3269 of file bilininteg.hpp.

◆ nh

Vector mfem::DGDiffusionIntegrator::nh

protected

Definition at line 3262 of file bilininteg.hpp.

◆ ni

Vector mfem::DGDiffusionIntegrator::ni

protected

Definition at line 3262 of file bilininteg.hpp.

◆ nor

Vector mfem::DGDiffusionIntegrator::nor

protected

Definition at line 3262 of file bilininteg.hpp.

◆ nq

int mfem::DGDiffusionIntegrator::nq

protected

Definition at line 3269 of file bilininteg.hpp.

◆ pa_data

Vector mfem::DGDiffusionIntegrator::pa_data

protected

Definition at line 3267 of file bilininteg.hpp.

◆ Q

Coefficient* mfem::DGDiffusionIntegrator::Q

protected

Definition at line 3257 of file bilininteg.hpp.

◆ quad1D

int mfem::DGDiffusionIntegrator::quad1D

protected

Definition at line 3269 of file bilininteg.hpp.

◆ shape1

Vector mfem::DGDiffusionIntegrator::shape1

protected

Definition at line 3262 of file bilininteg.hpp.

◆ shape2

Vector mfem::DGDiffusionIntegrator::shape2

protected

Definition at line 3262 of file bilininteg.hpp.

◆ sigma

real_t mfem::DGDiffusionIntegrator::sigma

protected

Definition at line 3259 of file bilininteg.hpp.

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

fem/bilininteg.hpp
fem/bilininteg.cpp

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DGDiffusionIntegrator() [1/3]

◆ DGDiffusionIntegrator() [2/3]

◆ DGDiffusionIntegrator() [3/3]

Member Function Documentation

◆ AddMultPAFaceNormalDerivatives()

◆ AssembleFaceMatrix() [1/2]

◆ AssembleFaceMatrix() [2/2]

◆ AssemblePA() [1/2]

◆ AssemblePA() [2/2]

◆ AssemblePABoundaryFaces()

◆ AssemblePAInteriorFaces()

◆ GetRule()

◆ RequiresFaceNormalDerivatives()

Member Data Documentation

◆ adjJ

◆ dim

◆ dofs1D

◆ dshape1

◆ dshape1dn

◆ dshape2

◆ dshape2dn

◆ irs

◆ jmat

◆ kappa

◆ maps

◆ MQ

◆ mq

◆ nf

◆ nh

◆ ni

◆ nor

◆ nq

◆ pa_data

◆ Q

◆ quad1D

◆ shape1

◆ shape2

◆ sigma