A function class to calculate element matrices for the Div u*Div v Bilinearform. More...

#include <divdiv.hh>

Inheritance diagram for hp3D::DivDiv< Weight >:

Public Member Functions
virtual DivDiv *	clone () const
	Virtual constructor. More...

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

virtual MaxwellSharedData *	data () const
	Returns the data to be shared. More...

virtual void	data (concepts::Cloneable *d)
	Sets the data to be shared. More...

	DivDiv (concepts::Array< Real > &jacobian, concepts::Array< concepts::MapReal3d > &jacobianInv, const Hexahedron *oldElm, SingularSet &singularities, const uint i=0, const uint j=0, const concepts::Formula< Real > frm=0)
	Constructor. More...

	DivDiv (SingularSet &singularities, const uint i=0, const uint j=0, const concepts::Formula< Real > *frm=0)
	Constructor. More...

virtual void	operator() (const concepts::Element< Real > &elmX, const concepts::Element< Real > &elmY, concepts::ElementMatrix< Real > &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 Hexahedron &elmX, const Hexahedron &elmY, concepts::ElementMatrix< Real > &em) const

virtual	~DivDiv ()

Static Public Member Functions
static void	setup (vectorial::BilinearForm< Real, Real > &bf, hp3D::SingularSet &singularities, const concepts::Formula< Real > *frm=0)
	Sets up the bilinear form for the div-div part for Maxwell's equation using weighted regularization. More...

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

Private Attributes
concepts::Array< concepts::MapReal3d >	coeffMatrix_
	Intermediate data for element matrix computation. More...

std::unique_ptr< const concepts::Formula< Real > >	frm_
	Formula for dielectricity $\epsilon$ . More...

uint	i_
	row and column inside the vectorial Bilinearform More...

uint	j_

concepts::Array< Real > *	jacobian_
	Data shared between components of the div-div bilinear forms,. More...

concepts::Array< concepts::MapReal3d > *	jacobianInv_

const Hexahedron **	oldElm_

MaxwellSharedData *	sharedData_

SingularSet &	singularities_
	two lists, one with edges, one with vertices (which are marked as singular) More...

Weight	weight_
	weight of regularization More...

Detailed Description

template<class Weight = TrivialWeight>
class hp3D::DivDiv< Weight >

A function class to calculate element matrices for the Div u*Div v Bilinearform.

The Jacobian, its inverse and additional quantities only depending on the physical coordinates are only computed once for all components of the bilinear form. To this end, several additional parameters can be given to the constructor (see there). However, it is much simpler to use the DivDiv::setup method which sets up everything automatically.

Author: Kersten Schmidt, 2002

Definition at line 46 of file divdiv.hh.

Constructor & Destructor Documentation

◆ DivDiv() [1/2]

template<class Weight = TrivialWeight>

hp3D::DivDiv< Weight >::DivDiv	(	SingularSet &	singularities,
		const uint	i = `0`,
		const uint	j = `0`,
		const concepts::Formula< Real > *	frm = `0`
	)

Constructor.

Parameters

singularities	List of singular edges and vertices.
i	Index of the bilinear form in the vector valued bilinear form
j	Index of the bilinear form in the vector valued bilinear form
frm	Formula for the dielectricity $\epsilon$

◆ DivDiv() [2/2]

template<class Weight = TrivialWeight>

hp3D::DivDiv< Weight >::DivDiv	(	concepts::Array< Real > &	jacobian,
		concepts::Array< concepts::MapReal3d > &	jacobianInv,
		const Hexahedron **	oldElm,
		SingularSet &	singularities,
		const uint	i = `0`,
		const uint	j = `0`,
		const concepts::Formula< Real > *	frm = `0`
	)

Constructor.

Parameters

jacobian	Reference to the storage space for the Jacobian. Here, also other quantities are stored (premultiplied).
jacobianInv	Reference to the storage space for the inverse of the Jacobian.
oldElm	Pointer to the element of the last invocation to detect if the information in `jacobian` and `jacobianInv` is still valid.
singularities	List of singular edges and vertices.
i	Index of the bilinear form in the vector valued bilinear form
j	Index of the bilinear form in the vector valued bilinear form
frm	Formula for the dielectricity $\epsilon$

◆ ~DivDiv()

template<class Weight = TrivialWeight>

virtual hp3D::DivDiv< Weight >::~DivDiv ( )

virtual

Member Function Documentation

◆ clone() [1/2]

template<class Weight = TrivialWeight>

virtual DivDiv* hp3D::DivDiv< Weight >::clone ( ) const

virtual

Virtual constructor.

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

Implements concepts::Cloneable.

◆ 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.

◆ data() [1/2]

template<class Weight = TrivialWeight>

virtual MaxwellSharedData* hp3D::DivDiv< Weight >::data ( ) const

virtual

Returns the data to be shared.

Implements vectorial::Vectorizable.

◆ data() [2/2]

template<class Weight = TrivialWeight>

virtual void hp3D::DivDiv< Weight >::data ( concepts::Cloneable * )

virtual

Sets the data to be shared.

Implements vectorial::Vectorizable.

◆ info()

template<class Weight = TrivialWeight>

virtual std::ostream& hp3D::DivDiv< Weight >::info ( std::ostream & os ) const

protectedvirtual

Returns information in an output stream.

Reimplemented from concepts::BilinearForm< Real >.

◆ operator()() [1/4]

template<class Weight = TrivialWeight>

virtual void hp3D::DivDiv< Weight >::operator()	(	const concepts::Element< Real > &	elmX,
		const concepts::Element< Real > &	elmY,
		concepts::ElementMatrix< Real > &	em
	)		const

virtual

◆ 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 Weight = TrivialWeight>

void hp3D::DivDiv< Weight >::operator()	(	const Hexahedron &	elmX,
		const Hexahedron &	elmY,
		concepts::ElementMatrix< Real > &	em
	)		const

◆ setup()

template<class Weight = TrivialWeight>

static void hp3D::DivDiv< Weight >::setup	(	vectorial::BilinearForm< Real, Real > &	bf,
		hp3D::SingularSet &	singularities,
		const concepts::Formula< Real > *	frm = `0`
	)

static

Sets up the bilinear form for the div-div part for Maxwell's equation using weighted regularization.

Parameters

bf	Resulting bilinear form (use this to compute the matrix)
singularities	Singularities needed for computing the weight
frm	Formula for the dielectricity (if none is given, 1 is assumed)

Precondition: The number of components of bf needs to be 3 x 3.

Examples: hpFEM3d-EV.cc.

Member Data Documentation

◆ coeffMatrix_

template<class Weight = TrivialWeight>

concepts::Array<concepts::MapReal3d> hp3D::DivDiv< Weight >::coeffMatrix_

mutableprivate

Intermediate data for element matrix computation.

Definition at line 113 of file divdiv.hh.

◆ frm_

template<class Weight = TrivialWeight>

std::unique_ptr<const concepts::Formula<Real> > hp3D::DivDiv< Weight >::frm_

private

Formula for dielectricity $\epsilon$ .

Definition at line 110 of file divdiv.hh.

◆ i_

template<class Weight = TrivialWeight>

uint hp3D::DivDiv< Weight >::i_

private

row and column inside the vectorial Bilinearform

Definition at line 102 of file divdiv.hh.

◆ j_

template<class Weight = TrivialWeight>

uint hp3D::DivDiv< Weight >::j_

private

Definition at line 102 of file divdiv.hh.

◆ jacobian_

template<class Weight = TrivialWeight>

concepts::Array<Real>* hp3D::DivDiv< Weight >::jacobian_

private

Data shared between components of the div-div bilinear forms,.

Definition at line 117 of file divdiv.hh.

◆ jacobianInv_

template<class Weight = TrivialWeight>

concepts::Array<concepts::MapReal3d>* hp3D::DivDiv< Weight >::jacobianInv_

private

Definition at line 118 of file divdiv.hh.

◆ oldElm_

template<class Weight = TrivialWeight>

const Hexahedron** hp3D::DivDiv< Weight >::oldElm_

private

Definition at line 119 of file divdiv.hh.

◆ sharedData_

template<class Weight = TrivialWeight>

MaxwellSharedData* hp3D::DivDiv< Weight >::sharedData_

private

Definition at line 120 of file divdiv.hh.

◆ singularities_

template<class Weight = TrivialWeight>

SingularSet& hp3D::DivDiv< Weight >::singularities_

private

two lists, one with edges, one with vertices (which are marked as singular)

Definition at line 105 of file divdiv.hh.

◆ weight_

template<class Weight = TrivialWeight>

Weight hp3D::DivDiv< Weight >::weight_

private

weight of regularization

Definition at line 107 of file divdiv.hh.

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

hp3D/divdiv.hh

hp3D::DivDiv< Weight > Class Template ReferenceabstractMaxwell's Equations » Three Dimensions » Weighted Regularization

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

template<class Weight = TrivialWeight> class hp3D::DivDiv< Weight >

Constructor & Destructor Documentation

◆ DivDiv() [1/2]

◆ DivDiv() [2/2]

◆ ~DivDiv()

Member Function Documentation

◆ clone() [1/2]

◆ clone() [2/2]

◆ data() [1/2]

◆ data() [2/2]

◆ info()

◆ operator()() [1/4]

◆ operator()() [2/4]

◆ operator()() [3/4]

◆ operator()() [4/4]

◆ setup()

Member Data Documentation

◆ coeffMatrix_

◆ frm_

◆ i_

◆ j_

◆ jacobian_

◆ jacobianInv_

◆ oldElm_

◆ sharedData_

◆ singularities_

◆ weight_

hp3D::DivDiv< Weight > Class Template Referenceabstract
Maxwell's Equations » Three Dimensions » Weighted Regularization

template<class Weight = TrivialWeight>
class hp3D::DivDiv< Weight >