dune-geometry  2.6-git
Public Types | Public Member Functions | Static Public Attributes | Protected Types | Protected Member Functions | Static Protected Member Functions | List of all members
Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits > Class Template Reference

Implement a MultiLinearGeometry with additional caching. More...

#include <dune/geometry/multilineargeometry.hh>

Inheritance diagram for Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >:
Inheritance graph

Public Types

typedef Base::ReferenceElement ReferenceElement
 
typedef Base::ctype ctype
 
typedef Base::LocalCoordinate LocalCoordinate
 
typedef Base::GlobalCoordinate GlobalCoordinate
 
typedef Base::JacobianTransposed JacobianTransposed
 
typedef Base::JacobianInverseTransposed JacobianInverseTransposed
 

Public Member Functions

template<class CornerStorage >
 CachedMultiLinearGeometry (const ReferenceElement &referenceElement, const CornerStorage &cornerStorage)
 
template<class CornerStorage >
 CachedMultiLinearGeometry (Dune::GeometryType gt, const CornerStorage &cornerStorage)
 
bool affine () const
 is this mapping affine? More...
 
GlobalCoordinate center () const
 obtain the centroid of the mapping's image More...
 
GlobalCoordinate global (const LocalCoordinate &local) const
 evaluate the mapping More...
 
 LocalCoordinate (const GlobalCoordinate &global) const
 evaluate the inverse mapping More...
 
ctype integrationElement (const LocalCoordinate &local) const
 obtain the integration element More...
 
ctype volume () const
 obtain the volume of the mapping's image More...
 
JacobianTransposed jacobianTransposed (const LocalCoordinate &local) const
 obtain the transposed of the Jacobian More...
 
JacobianInverseTransposed jacobianInverseTransposed (const LocalCoordinate &local) const
 obtain the transposed of the Jacobian's inverse More...
 
GlobalCoordinate corner (int i) const
 obtain coordinates of the i-th corner More...
 
bool affine () const
 is this mapping affine? More...
 
Dune::GeometryType type () const
 obtain the name of the reference element More...
 
int corners () const
 obtain number of corners of the corresponding reference element More...
 
GlobalCoordinate corner (int i) const
 obtain coordinates of the i-th corner More...
 
GlobalCoordinate center () const
 obtain the centroid of the mapping's image More...
 
ctype volume () const
 obtain the volume of the mapping's image More...
 

Static Public Attributes

static const int mydimension
 geometry dimension More...
 
static const int coorddimension
 coordinate dimension More...
 

Protected Types

typedef Base::MatrixHelper MatrixHelper
 
typedef Dune::ReferenceElements< ctype, mydimensionReferenceElements
 
typedef std::conditional< hasSingleGeometryType, std::integral_constant< unsigned int, MultiLinearGeometryTraits< ct > ::template hasSingleGeometryType< mydimension >::topologyId >, unsigned int >::type TopologyId
 

Protected Member Functions

ReferenceElement refElement () const
 
bool affine (JacobianTransposed &jacobianT) const
 
TopologyId topologyId () const
 

Static Protected Member Functions

static bool affine (TopologyId topologyId, std::integral_constant< int, dim >, CornerIterator &cit, JacobianTransposed &jt)
 
static bool affine (TopologyId topologyId, std::integral_constant< int, 0 >, CornerIterator &cit, JacobianTransposed &jt)
 
static void global (TopologyId topologyId, std::integral_constant< int, dim >, CornerIterator &cit, const ctype &df, const LocalCoordinate &x, const ctype &rf, GlobalCoordinate &y)
 
static void global (TopologyId topologyId, std::integral_constant< int, 0 >, CornerIterator &cit, const ctype &df, const LocalCoordinate &x, const ctype &rf, GlobalCoordinate &y)
 
static void jacobianTransposed (TopologyId topologyId, std::integral_constant< int, dim >, CornerIterator &cit, const ctype &df, const LocalCoordinate &x, const ctype &rf, FieldMatrix< ctype, rows, cdim > &jt)
 
static void jacobianTransposed (TopologyId topologyId, std::integral_constant< int, 0 >, CornerIterator &cit, const ctype &df, const LocalCoordinate &x, const ctype &rf, FieldMatrix< ctype, rows, cdim > &jt)
 

Detailed Description

template<class ct, int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
class Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >

Implement a MultiLinearGeometry with additional caching.

This class implements the same interface and functionality as MultiLinearGeometry. However, it additionally implements caching for various results.

Template Parameters
ctcoordinate type
mydimgeometry dimension
cdimcoordinate dimension
Traitstraits allowing to tweak some implementation details (optional)

Member Typedef Documentation

◆ ctype

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::ctype Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::ctype

◆ GlobalCoordinate

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::GlobalCoordinate Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::GlobalCoordinate

◆ JacobianInverseTransposed

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::JacobianInverseTransposed Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::JacobianInverseTransposed

◆ JacobianTransposed

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::JacobianTransposed Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::JacobianTransposed

◆ LocalCoordinate

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::LocalCoordinate Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::LocalCoordinate

◆ MatrixHelper

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::MatrixHelper Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::MatrixHelper
protected

◆ ReferenceElement

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
typedef Base::ReferenceElement Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::ReferenceElement

◆ ReferenceElements

◆ TopologyId

typedef std::conditional< hasSingleGeometryType, std::integral_constant< unsigned int, MultiLinearGeometryTraits< ct > ::template hasSingleGeometryType< mydimension >::topologyId >, unsigned int >::type Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::TopologyId
protectedinherited

Constructor & Destructor Documentation

◆ CachedMultiLinearGeometry() [1/2]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
template<class CornerStorage >
Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::CachedMultiLinearGeometry ( const ReferenceElement referenceElement,
const CornerStorage &  cornerStorage 
)
inline

◆ CachedMultiLinearGeometry() [2/2]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
template<class CornerStorage >
Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::CachedMultiLinearGeometry ( Dune::GeometryType  gt,
const CornerStorage &  cornerStorage 
)
inline

Member Function Documentation

◆ affine() [1/5]

bool Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::affine
inlineinherited

is this mapping affine?

◆ affine() [2/5]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
bool Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::affine ( ) const
inline

is this mapping affine?

◆ affine() [3/5]

bool Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::affine ( JacobianTransposed jacobianT) const
inlineprotectedinherited

◆ affine() [4/5]

bool Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::affine ( TopologyId  topologyId,
std::integral_constant< int, 0 >  ,
CornerIterator &  cit,
JacobianTransposed jt 
)
inlinestaticprotectedinherited

◆ affine() [5/5]

bool Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::affine ( TopologyId  topologyId,
std::integral_constant< int, dim >  ,
CornerIterator &  cit,
JacobianTransposed jt 
)
inlinestaticprotectedinherited

◆ center() [1/2]

GlobalCoordinate Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::center
inlineinherited

obtain the centroid of the mapping's image

◆ center() [2/2]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
GlobalCoordinate Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::center ( ) const
inline

obtain the centroid of the mapping's image

◆ corner() [1/2]

GlobalCoordinate Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::corner ( int  i) const
inlineinherited

obtain coordinates of the i-th corner

◆ corner() [2/2]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
GlobalCoordinate Dune::MultiLinearGeometry< ct, mydim, cdim, Traits >::corner
inline

obtain coordinates of the i-th corner

◆ corners()

int Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::corners
inlineinherited

obtain number of corners of the corresponding reference element

◆ global() [1/3]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
GlobalCoordinate Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::global ( const LocalCoordinate local) const
inline

evaluate the mapping

Parameters
[in]locallocal coordinate to map
Returns
corresponding global coordinate

◆ global() [2/3]

void Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::global ( TopologyId  topologyId,
std::integral_constant< int, 0 >  ,
CornerIterator &  cit,
const ctype df,
const LocalCoordinate x,
const ctype rf,
GlobalCoordinate y 
)
inlinestaticprotectedinherited

◆ global() [3/3]

void Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::global ( TopologyId  topologyId,
std::integral_constant< int, dim >  ,
CornerIterator &  cit,
const ctype df,
const LocalCoordinate x,
const ctype rf,
GlobalCoordinate y 
)
inlinestaticprotectedinherited

◆ integrationElement()

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
ctype Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::integrationElement ( const LocalCoordinate local) const
inline

obtain the integration element

If the Jacobian of the mapping is denoted by $J(x)$, the integration integration element $\mu(x)$ is given by

\[ \mu(x) = \sqrt{|\det (J^T(x) J(x))|}.\]

Parameters
[in]locallocal coordinate to evaluate the integration element in
Returns
the integration element $\mu(x)$.
Note
For affine mappings, it is more efficient to call jacobianInverseTransposed before integrationElement, if both are required.

◆ jacobianInverseTransposed()

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
JacobianInverseTransposed Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::jacobianInverseTransposed ( const LocalCoordinate local) const
inline

obtain the transposed of the Jacobian's inverse

The Jacobian's inverse is defined as a pseudo-inverse. If we denote the Jacobian by $J(x)$, the following condition holds:

\[J^{-1}(x) J(x) = I.\]

◆ jacobianTransposed() [1/3]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
JacobianTransposed Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::jacobianTransposed ( const LocalCoordinate local) const
inline

obtain the transposed of the Jacobian

Parameters
[in]locallocal coordinate to evaluate Jacobian in
Returns
a reference to the transposed of the Jacobian
Note
The returned reference is reused on the next call to JacobianTransposed, destroying the previous value.

◆ jacobianTransposed() [2/3]

void Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::jacobianTransposed ( TopologyId  topologyId,
std::integral_constant< int, 0 >  ,
CornerIterator &  cit,
const ctype df,
const LocalCoordinate x,
const ctype rf,
FieldMatrix< ctype, rows, cdim > &  jt 
)
inlinestaticprotectedinherited

◆ jacobianTransposed() [3/3]

void Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::jacobianTransposed ( TopologyId  topologyId,
std::integral_constant< int, dim >  ,
CornerIterator &  cit,
const ctype df,
const LocalCoordinate x,
const ctype rf,
FieldMatrix< ctype, rows, cdim > &  jt 
)
inlinestaticprotectedinherited

◆ LocalCoordinate()

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::LocalCoordinate ( const GlobalCoordinate global) const
inline

evaluate the inverse mapping

Parameters
[in]globalglobal coordinate to map
Returns
corresponding local coordinate
Note
For given global coordinate y the returned local coordinate x that minimizes the following function over the local coordinate space spanned by the reference element.
(global( x ) - y).two_norm()

◆ refElement()

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
ReferenceElement Dune::MultiLinearGeometry< ct, mydim, cdim, Traits >::refElement
inlineprotected

◆ topologyId()

TopologyId Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::topologyId
inlineprotectedinherited

◆ type()

Dune::GeometryType Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::type
inlineinherited

obtain the name of the reference element

◆ volume() [1/2]

ctype Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::volume
inlineinherited

obtain the volume of the mapping's image

Note
The current implementation just returns
integrationElement( refElement().position( 0, 0 ) ) * refElement().volume()
which is wrong for n-linear surface maps and other nonlinear maps.

◆ volume() [2/2]

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
ctype Dune::CachedMultiLinearGeometry< ct, mydim, cdim, Traits >::volume ( ) const
inline

obtain the volume of the mapping's image

Member Data Documentation

◆ coorddimension

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
const int Dune::MultiLinearGeometry< ct, mydim, cdim, Traits >::coorddimension
static

coordinate dimension

◆ mydimension

template<class ct , int mydim, int cdim, class Traits = MultiLinearGeometryTraits< ct >>
const int Dune::MultiLinearGeometry< ct, mydim, cdim, Traits >::mydimension
static

geometry dimension


The documentation for this class was generated from the following file:
Dune::CachedMultiLinearGeometry::integrationElement
ctype integrationElement(const LocalCoordinate &local) const
obtain the integration element
Definition: multilineargeometry.hh:588
Dune::CachedMultiLinearGeometry::global
GlobalCoordinate global(const LocalCoordinate &local) const
evaluate the mapping
Definition: multilineargeometry.hh:535
Dune::MultiLinearGeometry< ct, mydim, cdim, MultiLinearGeometryTraits< ct > >::refElement
ReferenceElement refElement() const
Definition: multilineargeometry.hh:383