A function class to calculate element matrices for the bilinear form. More...

#include <bf_graduv.hh>

Inheritance diagram for hp2Dedge::GraduvMatrix< F >:

Public Types
typedef concepts::ElementFormulaContainer< concepts::Mapping< F, 2 > >	FrmE_Matrix

typedef concepts::Combtype< Real, Real >::type	value_type

Public Member Functions
virtual GraduvMatrix< F > *	clone () const
	Virtual constructor. More...

virtual BilinearForm *	clone () const=0
	Virtual constructor. More...

concepts::RCP< concepts::SharedJacobianAdj< 2 > >	data () const
	Gets the pointer to the shared data. More...

void	data (const concepts::RCP< concepts::SharedJacobianAdj< 2 > > d)
	Set the pointer to the shared data. More...

	GraduvMatrix (const FrmE_Matrix frm=FrmE_Matrix())

void	operator() (const concepts::Element< Real > &elmX, const concepts::Element< Real > &elmY, concepts::ElementMatrix< F > &em) const

virtual void	operator() (const Element< typename Realtype< Real >::type > &elmX, const Element< typename Realtype< Real >::type > &elmY, ElementMatrix< Real > &em) const=0
	Evaluates the bilinear form for all shape functions on `elmX` and `elmY` and stores the result in the matrix `em`. More...

virtual void	operator() (const Element< typename Realtype< Real >::type > &elmX, const Element< typename Realtype< Real >::type > &elmY, ElementMatrix< Real > &em, const ElementPair< typename Realtype< Real >::type > &ep) const
	Evaluates the bilinear form for all shape functions on `elmX` and `elmY` and stores the result in the matrix `em`. More...

void	operator() (const hp2D::Quad< Real > &elmX, const Quad< Real > &elmY, concepts::ElementMatrix< F > &em) const

virtual	~GraduvMatrix ()

Protected Member Functions
bool	assemble_ (const hp2D::Quad< Real > elmX, const Quad< Real > elmY, concepts::ElementMatrix< value_type > &em) const

void	computeIntermediate_ (const BaseQuad< Real > &elm, const int i=-1, const int j=-1) const
	Compute the intermediate data for element matrix computation. More...

virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream. More...

Protected Attributes
bool	all_
	Parameter for the sum factorisation. More...

concepts::ElementFormulaContainer< Real >	frm_
	Element formula. More...

concepts::ElementFormulaContainer< concepts::Mapping< Real, 2 > >	frmM_
	Matrix element formula. More...

concepts::Array< concepts::Mapping< Real, 2 > >	intermediateMatrix_
	Intermediate matrix. More...

concepts::Array< Real >	intermediateValue_
	Intermediate value. More...

Private Member Functions
void	computeadjJ_adjJT_rank1_ (concepts::Array< concepts::Mapping< Real, 2 > > &intermediateMatrix, const int i, const int j) const
	Computes the either adj(J)*adj(J)^T or in the case of partial derivatives (i > NONE, j > NONE) the rank-1 product of j-th column of adj(J) with i-th row of adj(J)^T. More...

void	computeJacobianMatrix_ (const BaseQuad< Real > &elm, concepts::Array< concepts::Mapping< Real, 2 > > &J, concepts::Array< Real > &detJ_inv) const
	Compute the Jacobian matrix and the inverse of its determinant on each quadrature point. More...

Private Attributes
concepts::Array< Real >	factor_
	Intermediate data for element matrix computation. More...

concepts::Array< concepts::Mapping< Real, 2 > >	jacobianInv_

concepts::RCP< concepts::SharedJacobianAdj< 2 > >	sharedData_
	Shared data for vectorial bilinear forms. More...

Detailed Description

template<class F = Real>
class hp2Dedge::GraduvMatrix< F >

A function class to calculate element matrices for the bilinear form.

$\int\limits_{K}(A(\underline{\xi})\underline{u})^\top\;\nabla{v}\;d\underline{\xi} = \int\limits_{\hat{K}}\underline{\hat{u}}^\top J^{-1}A(\underline{\xi})^{\top}J^{-\top}\;\nabla{\hat{v}}|\det J|\;d\underline{\hat{\xi}}$

Definition at line 134 of file bf_graduv.hh.

Member Typedef Documentation

◆ FrmE_Matrix

template<class F = Real>

typedef concepts::ElementFormulaContainer<concepts::Mapping<F, 2> > hp2Dedge::GraduvMatrix< F >::FrmE_Matrix

Definition at line 139 of file bf_graduv.hh.

◆ value_type

typedef concepts::Combtype<Real , Real >::type hp2Dedge::GraduvBase< Real , Real >::value_type

inherited

Definition at line 66 of file bf_graduv.hh.

Constructor & Destructor Documentation

◆ GraduvMatrix()

template<class F = Real>

hp2Dedge::GraduvMatrix< F >::GraduvMatrix ( const FrmE_Matrix frm = FrmE_Matrix() )

◆ ~GraduvMatrix()

template<class F = Real>

virtual hp2Dedge::GraduvMatrix< F >::~GraduvMatrix ( )

inlinevirtual

Definition at line 142 of file bf_graduv.hh.

Member Function Documentation

◆ assemble_()

bool hp2Dedge::GraduvBase< Real , Real >::assemble_	(	const hp2D::Quad< Real > *	elmX,
		const Quad< Real > *	elmY,
		concepts::ElementMatrix< value_type > &	em
	)		const

protectedinherited

◆ clone() [1/2]

template<class F = Real>

virtual GraduvMatrix<F>* hp2Dedge::GraduvMatrix< F >::clone ( ) const

inlinevirtual

Virtual constructor.

Returns a pointer to a copy of itself. The caller is responsible to destroy this copy.

Implements concepts::Cloneable.

Definition at line 153 of file bf_graduv.hh.

◆ clone() [2/2]

virtual BilinearForm* concepts::BilinearForm< Real , typename Realtype<Real >::type >::clone ( ) const

pure virtualinherited

Virtual constructor.

Returns a pointer to a copy of itself. The caller is responsible to destroy this copy.

◆ computeadjJ_adjJT_rank1_()

void hp2D::BilinearFormHelper_1_1< Real , Real >::computeadjJ_adjJT_rank1_	(	concepts::Array< concepts::Mapping< Real, 2 > > &	intermediateMatrix,
		const int	i,
		const int	j
	)		const

privateinherited

Computes the either adj(J)*adj(J)^T or in the case of partial derivatives (i > NONE, j > NONE) the rank-1 product of j-th column of adj(J) with i-th row of adj(J)^T.

◆ computeIntermediate_()

void hp2D::BilinearFormHelper_1_1< Real , Real >::computeIntermediate_	(	const BaseQuad< Real > &	elm,
		const int	i = `-1`,
		const int	j = `-1`
	)		const

protectedinherited

Compute the intermediate data for element matrix computation.

Parameters

i	if i=0 or 1, then take only i-th column of Jacobian matrix (for test function)
j	if j=0 or 1, then take only j-th column of Jacobian matrix (for trial function)

The Jacobian matrices have to been taken both full (i,j = -1) or both partial (i,j = 0 or 1).

Matrix formulas and complex valued scalar formulas are only implemented for full Jacobians.

◆ computeJacobianMatrix_()

void hp2D::BilinearFormHelper_1_1< Real , Real >::computeJacobianMatrix_	(	const BaseQuad< Real > &	elm,
		concepts::Array< concepts::Mapping< Real, 2 > > &	J,
		concepts::Array< Real > &	detJ_inv
	)		const

privateinherited

Compute the Jacobian matrix and the inverse of its determinant on each quadrature point.

◆ data() [1/2]

concepts::RCP<concepts::SharedJacobianAdj<2> > hp2D::BilinearFormHelper_1_1< Real , Real >::data ( ) const

inherited

Gets the pointer to the shared data.

◆ data() [2/2]

void hp2D::BilinearFormHelper_1_1< Real , Real >::data ( const concepts::RCP< concepts::SharedJacobianAdj< 2 > > d )

inherited

Set the pointer to the shared data.

◆ info()

template<class F = Real>

virtual std::ostream& hp2Dedge::GraduvMatrix< F >::info ( std::ostream & os ) const

protectedvirtual

Returns information in an output stream.

Reimplemented from concepts::BilinearForm< Real >.

◆ operator()() [1/4]

template<class F = Real>

void hp2Dedge::GraduvMatrix< F >::operator()	(	const concepts::Element< Real > &	elmX,
		const concepts::Element< Real > &	elmY,
		concepts::ElementMatrix< F > &	em
	)		const

◆ operator()() [2/4]

virtual void concepts::BilinearForm< Real , typename Realtype<Real >::type >::operator()	(	const Element< G > &	elmX,
		const Element< G > &	elmY,
		ElementMatrix< F > &	em
	)		const

pure virtualinherited

Evaluates the bilinear form for all shape functions on elmX and elmY and stores the result in the matrix em.

Postcondition: The returned matrix em has the correct size.

Parameters

elmX	Left element (test functions)
elmY	Right element (trial functions)
em	Return element matrix

◆ operator()() [3/4]

virtual void concepts::BilinearForm< Real , typename Realtype<Real >::type >::operator()	(	const Element< G > &	elmX,
		const Element< G > &	elmY,
		ElementMatrix< F > &	em,
		const ElementPair< G > &	ep
	)		const

inlinevirtualinherited

Evaluates the bilinear form for all shape functions on elmX and elmY and stores the result in the matrix em.

If this method is not reimplemented in a derived class, the default behaviour is to call the application operator without ep.

Postcondition: The returned matrix em has the correct size.

Parameters

elmX	Left element
elmY	Right element
em	Return element matrix
ep	Element pair holding more information on the pair `elmX` and `elmY`

Definition at line 57 of file bilinearForm.hh.

◆ operator()() [4/4]

template<class F = Real>

void hp2Dedge::GraduvMatrix< F >::operator()	(	const hp2D::Quad< Real > &	elmX,
		const Quad< Real > &	elmY,
		concepts::ElementMatrix< F > &	em
	)		const

Member Data Documentation

◆ all_

bool hp2Dedge::GraduvBase< Real , Real >::all_

protectedinherited

Parameter for the sum factorisation.

Definition at line 79 of file bf_graduv.hh.

◆ factor_

concepts::Array<Real > hp2Dedge::GraduvBase< Real , Real >::factor_

mutableprivateinherited

Intermediate data for element matrix computation.

Definition at line 87 of file bf_graduv.hh.

◆ frm_

concepts::ElementFormulaContainer<Real > hp2D::BilinearFormHelper_1_1< Real , Real >::frm_

protectedinherited

Element formula.

Definition at line 193 of file bilinearFormHelper.hh.

◆ frmM_

concepts::ElementFormulaContainer<concepts::Mapping<Real ,2> > hp2D::BilinearFormHelper_1_1< Real , Real >::frmM_

protectedinherited

Matrix element formula.

Definition at line 195 of file bilinearFormHelper.hh.

◆ intermediateMatrix_

concepts::Array<concepts::Mapping<Real ,2> > hp2D::BilinearFormHelper_1_1< Real , Real >::intermediateMatrix_

mutableprotectedinherited

Intermediate matrix.

In case of a scalar formula:

$\mbox{adj}(J) \mbox{adj}(J)^\top$

In case of a matrix formula :

$\mbox{adj}(J) M \mbox{adj}(J)^\top$

In case of partial Jacobian:

$\mbox{adj}(J)_{\cdot,j} (\mbox{adj}(J)_{\cdot,i})^\top$

Definition at line 191 of file bilinearFormHelper.hh.

◆ intermediateValue_

concepts::Array<Real > hp2D::BilinearFormHelper_1_1< Real , Real >::intermediateValue_

mutableprotectedinherited

Intermediate value.

In case of a scalar formula:

$\frac{f(F_K(\xi))}{\det J}$

In case of a matrix formula:

$\frac{1}{\det J}$

Definition at line 179 of file bilinearFormHelper.hh.

◆ jacobianInv_

concepts::Array<concepts::Mapping<Real , 2> > hp2Dedge::GraduvBase< Real , Real >::jacobianInv_

mutableprivateinherited

Definition at line 88 of file bf_graduv.hh.

◆ sharedData_

concepts::RCP<concepts::SharedJacobianAdj<2> > hp2D::BilinearFormHelper_1_1< Real , Real >::sharedData_

privateinherited

Shared data for vectorial bilinear forms.

Definition at line 212 of file bilinearFormHelper.hh.

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

hp2Dedge/bf_graduv.hh

hp2Dedge::GraduvMatrix< F > Class Template Referenceabstract

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

template<class F = Real> class hp2Dedge::GraduvMatrix< F >

Member Typedef Documentation

◆ FrmE_Matrix

◆ value_type

Constructor & Destructor Documentation

◆ GraduvMatrix()

◆ ~GraduvMatrix()

Member Function Documentation

◆ assemble_()

◆ clone() [1/2]

◆ clone() [2/2]

◆ computeadjJ_adjJT_rank1_()

◆ computeIntermediate_()

◆ computeJacobianMatrix_()

◆ data() [1/2]

◆ data() [2/2]

◆ info()

◆ operator()() [1/4]

◆ operator()() [2/4]

◆ operator()() [3/4]

◆ operator()() [4/4]

Member Data Documentation

◆ all_

◆ factor_

◆ frm_

◆ frmM_

◆ intermediateMatrix_

◆ intermediateValue_

◆ jacobianInv_

◆ sharedData_

template<class F = Real>
class hp2Dedge::GraduvMatrix< F >