You are here: Concepts>Concepts Web>Class documentation

Class documentation of Concepts

Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | Protected Member Functions | List of all members
hp2D::DivDiv< Weight > Class Template Referenceabstract
Maxwell's Equations » Two Dimensions » Weighted Regularization

#include <divdiv.hh>

Inheritance diagram for hp2D::DivDiv< Weight >:
concepts::BilinearForm< Real > vectorial::Vectorizable concepts::Cloneable concepts::OutputOperator

Public Member Functions

 DivDiv (SingularSet &singularities, const uint i=0, const uint j=0, const concepts::Formula< Real > *frm=0)
 
 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)
 
virtual DivDiv * clone () const
 
virtual void operator() (const concepts::Element< Real > &elmX, const concepts::Element< Real > &elmY, concepts::ElementMatrix< Real > &em) const
 
void operator() (const Quad< Real > &elmX, const Quad< Real > &elmY, concepts::ElementMatrix< Real > &em) const
 
virtual void data (concepts::Cloneable *d)
 Sets the data to be shared.
 
virtual MaxwellSharedData * data () const
 Returns the data to be shared.
 
virtual void operator() (const Element< typename Realtype< Real >::type > &elmX, const Element< typename Realtype< Real >::type > &elmY, ElementMatrix< Real > &em) const=0
 
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
 

Static Public Member Functions

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

Protected Member Functions

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

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 59 of file divdiv.hh.

Constructor & Destructor Documentation

◆ DivDiv() [1/2]

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
singularitiesList of singular edges and vertices.
iIndex of the bilinear form in the vector valued bilinear form
jIndex of the bilinear form in the vector valued bilinear form
frmFormula for the dielectricity $\epsilon$

◆ DivDiv() [2/2]

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
jacobianReference to the storage space for the Jacobian. Here, also other quantities are stored (premultiplied).
jacobianInvReference to the storage space for the inverse of the Jacobian.
oldElmPointer to the element of the last invocation to detect if the information in jacobian and jacobianInv is still valid.
singularitiesList of singular edges and vertices.
iIndex of the bilinear form in the vector valued bilinear form
jIndex of the bilinear form in the vector valued bilinear form
frmFormula for the dielectic $\epsilon$

Member Function Documentation

◆ clone()

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::BilinearForm< Real >.

◆ data() [1/2]

template<typename Weight = TrivialWeight>
virtual MaxwellSharedData * hp2D::DivDiv< Weight >::data ( ) const
virtual

Returns the data to be shared.

Implements vectorial::Vectorizable.

◆ data() [2/2]

template<typename Weight = TrivialWeight>
virtual void hp2D::DivDiv< Weight >::data ( concepts::Cloneable *  )
virtual

Sets the data to be shared.

Implements vectorial::Vectorizable.

◆ info()

template<typename Weight = TrivialWeight>
virtual std::ostream & hp2D::DivDiv< Weight >::info ( std::ostream &  os) const
protectedvirtual

Returns information in an output stream.

Reimplemented from concepts::BilinearForm< Real >.

◆ operator()() [1/2]

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
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
elmXLeft element (test functions)
elmYRight element (trial functions)
emReturn element matrix

◆ operator()() [2/2]

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 
) 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
elmXLeft element
elmYRight element
emReturn element matrix
epElement pair holding more information on the pair elmX and elmY

Definition at line 57 of file bilinearForm.hh.

◆ setup()

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
bfResulting bilinear form (use this to compute the matrix)
singularitiesSingularities needed for computing the weight
frmFormula for the dielectricity (if none is given, 1 is assumed)
Precondition
The number of components of bf needs to be 2 x 2.
The documentation for this class was generated from the following file:
Generated on Wed Sep 13 2023 21:06:52 for Concepts by doxygen 1.9.8