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

#include <divdiv.hh>

Inheritance diagram for hp2D::DivDiv< Weight >:

Collaboration diagram for hp2D::DivDiv< Weight >:

Public Member Functions
virtual DivDiv *	clone () const
	Virtual constructor.
virtual MaxwellSharedData *	data () const
	Returns the data to be shared.
virtual void	data (concepts::Cloneable *d)
	Sets the data to be shared.
	DivDiv (SingularSet &singularities, const uint i=0, const uint j=0, const concepts::Formula< Real > *frm=0)
	Constructor.
	DivDiv (concepts::Array< Real > &jacobian, concepts::Array< concepts::MapReal2d > &jacobianInv, const Quad< Real > *oldElm, SingularSet &singularities, const uint i=0, const uint j=0, const concepts::Formula< Real > frm=0)
	Constructor.
void	operator() (const Quad< Real > &elmX, const Quad< Real > &elmY, concepts::ElementMatrix< Real > &em)
virtual void	operator() (const concepts::Element< Real > &elmX, const concepts::Element< Real > &elmY, concepts::ElementMatrix< Real > &em)
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)
	Evaluates the bilinear form for all shape functions on `elmX` and `elmY` and stores the result in the matrix `em`.
virtual void	operator() (const Element< typename Realtype< Real >::type > &elmX, const Element< typename Realtype< Real >::type > &elmY, ElementMatrix< Real > &em)=0
	Evaluates the bilinear form for all shape functions on `elmX` and `elmY` and stores the result in the matrix `em`.
virtual	~DivDiv ()
Static Public Member Functions
static void	setup (vectorial::BilinearForm< Real, Real > &bf, hp2D::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.
Protected Member Functions
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
Private Attributes
concepts::Array < concepts::MapReal2d >	coeffMatrix_
	Intermediate data for element matrix computation.
std::auto_ptr< const concepts::Formula< Real > >	frm_
	Formula for dielectricity $\epsilon$ .
uint	i_
	Row and column inside the vectorial bilinear form.
uint	j_
SingularSet &	singularities_
	Set of singular points.
Weight	weight_
	Weight of regularization.

concepts::Array< Real > *	jacobian_
	Intermediate data for element matrix computation.
concepts::Array < concepts::MapReal2d > *	jacobianInv_
	Intermediate data for element matrix computation.
const Quad< Real > **	oldElm_
	Intermediate data for element matrix computation.
MaxwellSharedData *	sharedData_
	Intermediate data for element matrix computation.

Detailed Description

template<typename Weight = TrivialWeight>
class hp2D::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, 2003

TODO: setupDivDiv()-function with singularities using SharedJacobianAdj instead of MaxwellSharedData.

Definition at line 57 of file divdiv.hh.

Constructor & Destructor Documentation

template<typename Weight = TrivialWeight>

hp2D::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$

template<typename Weight = TrivialWeight>

hp2D::DivDiv< Weight >::DivDiv	(	concepts::Array< Real > &	jacobian,
		concepts::Array< concepts::MapReal2d > &	jacobianInv,
		const Quad< Real > **	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 dielectic $\epsilon$

template<typename Weight = TrivialWeight>

virtual hp2D::DivDiv< Weight >::~DivDiv ( ) [virtual]

Member Function Documentation

template<typename Weight = TrivialWeight>

virtual DivDiv* hp2D::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.

template<typename Weight = TrivialWeight>

virtual void hp2D::DivDiv< Weight >::data ( concepts::Cloneable * ) [virtual]

Sets the data to be shared.

Implements vectorial::Vectorizable.

template<typename Weight = TrivialWeight>

virtual MaxwellSharedData* hp2D::DivDiv< Weight >::data ( ) const [virtual]

Returns the data to be shared.

Implements vectorial::Vectorizable.

template<typename Weight = TrivialWeight>

virtual std::ostream& hp2D::DivDiv< Weight >::info ( std::ostream & os ) const [protected, virtual]

Returns information in an output stream.

Reimplemented from concepts::BilinearForm< Real >.

virtual void concepts::BilinearForm< Real , typename Realtype<Real >::type >::operator()	(	const Element< typename Realtype<Real >::type > &	elmX,
		const Element< typename Realtype<Real >::type > &	elmY,
		ElementMatrix< Real > &	em
	)		`[pure virtual, inherited]`

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

virtual void concepts::BilinearForm< Real , typename Realtype<Real >::type >::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
	)		`[inline, virtual, inherited]`

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 53 of file bilinearForm.hh.

Here is the call graph for this function:

template<typename Weight = TrivialWeight>

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

template<typename Weight = TrivialWeight>

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

template<typename Weight = TrivialWeight>

static void hp2D::DivDiv< Weight >::setup	(	vectorial::BilinearForm< Real, Real > &	bf,
		hp2D::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 2 x 2.

Member Data Documentation

template<typename Weight = TrivialWeight>

concepts::Array<concepts::MapReal2d> hp2D::DivDiv< Weight >::coeffMatrix_ [private]

Intermediate data for element matrix computation.

Definition at line 124 of file divdiv.hh.

template<typename Weight = TrivialWeight>

std::auto_ptr<const concepts::Formula<Real> > hp2D::DivDiv< Weight >::frm_ [private]

Formula for dielectricity $\epsilon$ .

Definition at line 121 of file divdiv.hh.

template<typename Weight = TrivialWeight>

uint hp2D::DivDiv< Weight >::i_ [private]

Row and column inside the vectorial bilinear form.

Definition at line 113 of file divdiv.hh.

template<typename Weight = TrivialWeight>

uint hp2D::DivDiv< Weight >::j_ [private]

Definition at line 113 of file divdiv.hh.

template<typename Weight = TrivialWeight>

concepts::Array<Real>* hp2D::DivDiv< Weight >::jacobian_ [private]

Intermediate data for element matrix computation.

Definition at line 128 of file divdiv.hh.

template<typename Weight = TrivialWeight>

concepts::Array<concepts::MapReal2d>* hp2D::DivDiv< Weight >::jacobianInv_ [private]

Intermediate data for element matrix computation.

Definition at line 129 of file divdiv.hh.

template<typename Weight = TrivialWeight>

const Quad<Real>** hp2D::DivDiv< Weight >::oldElm_ [private]

Intermediate data for element matrix computation.

Definition at line 130 of file divdiv.hh.

template<typename Weight = TrivialWeight>

MaxwellSharedData* hp2D::DivDiv< Weight >::sharedData_ [private]

Intermediate data for element matrix computation.

Definition at line 131 of file divdiv.hh.

template<typename Weight = TrivialWeight>

SingularSet& hp2D::DivDiv< Weight >::singularities_ [private]

Set of singular points.

Definition at line 116 of file divdiv.hh.

template<typename Weight = TrivialWeight>

Weight hp2D::DivDiv< Weight >::weight_ [private]

Weight of regularization.

Definition at line 118 of file divdiv.hh.

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

hp2D/divdiv.hh

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

template<typename Weight = TrivialWeight> class hp2D::DivDiv< Weight >

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<typename Weight = TrivialWeight>
class hp2D::DivDiv< Weight >