Class xlifepp::GeomMapData#

class GeomMapData#

Collaboration diagram for xlifepp::GeomMapData:

digraph { graph [bgcolor="#00000000"] node [shape=rectangle style=filled fillcolor="#FFFFFF" font=Helvetica padding=2] edge [color="#1414CE"] "40" [label="xlifepp::Collection< string_t >" tooltip="xlifepp::Collection< string_t >"] "37" [label="xlifepp::Matrix< real_t >" tooltip="xlifepp::Matrix< real_t >"] "24" [label="xlifepp::Vector< real_t >" tooltip="xlifepp::Vector< real_t >"] "46" [label="xlifepp::Vector< xlifepp::Vector< real_t > >" tooltip="xlifepp::Vector< xlifepp::Vector< real_t > >"] "55" [label="std::list< std::pair< number_t, real_t > >" tooltip="std::list< std::pair< number_t, real_t > >"] "19" [label="std::list< std::vector< PolynomialT< real_t > > >" tooltip="std::list< std::vector< PolynomialT< real_t > > >"] "34" [label="std::map< number_t, std::vector< number_t > >" tooltip="std::map< number_t, std::vector< number_t > >"] "33" [label="std::map< number_t, xlifepp::Geometry * >" tooltip="std::map< number_t, xlifepp::Geometry * >"] "54" [label="std::map< string_t, number_t >" tooltip="std::map< string_t, number_t >"] "62" [label="std::map< xlifepp::GeomElement *, std::set< number_t > >" tooltip="std::map< xlifepp::GeomElement *, std::set< number_t > >"] "56" [label="std::map< xlifepp::Quadrature *, std::vector< xlifepp::Point > >" tooltip="std::map< xlifepp::Quadrature *, std::vector< xlifepp::Point > >"] "22" [label="std::map< xlifepp::Quadrature *, std::vector< xlifepp::ShapeValues > >" tooltip="std::map< xlifepp::Quadrature *, std::vector< xlifepp::ShapeValues > >"] "20" [label="std::list< T >" tooltip="std::list< T >"] "23" [label="std::map< K, T >" tooltip="std::map< K, T >"] "4" [label="std::vector< T >" tooltip="std::vector< T >"] "43" [label="std::vector< T >" tooltip="std::vector< T >"] "26" [label="std::vector< K >" tooltip="std::vector< K >"] "30" [label="std::vector< RealPair >" tooltip="std::vector< RealPair >"] "10" [label="std::vector< ShapeType >" tooltip="std::vector< ShapeType >"] "61" [label="std::vector< const xlifepp::GeomDomain * >" tooltip="std::vector< const xlifepp::GeomDomain * >"] "39" [label="std::vector< const xlifepp::Transformation * >" tooltip="std::vector< const xlifepp::Transformation * >"] "5" [label="std::vector< number_t >" tooltip="std::vector< number_t >"] "6" [label="std::vector< real_t >" tooltip="std::vector< real_t >"] "51" [label="std::vector< std::vector< GeoNumPair > >" tooltip="std::vector< std::vector< GeoNumPair > >"] "12" [label="std::vector< std::vector< int_t > >" tooltip="std::vector< std::vector< int_t > >"] "11" [label="std::vector< std::vector< number_t > >" tooltip="std::vector< std::vector< number_t > >"] "63" [label="std::vector< std::vector< real_t > >" tooltip="std::vector< std::vector< real_t > >"] "41" [label="std::vector< string_t >" tooltip="std::vector< string_t >"] "50" [label="std::vector< xlifepp::GeomDomain * >" tooltip="std::vector< xlifepp::GeomDomain * >"] "49" [label="std::vector< xlifepp::GeomElement * >" tooltip="std::vector< xlifepp::GeomElement * >"] "13" [label="std::vector< xlifepp::GeomRefElement * >" tooltip="std::vector< xlifepp::GeomRefElement * >"] "15" [label="std::vector< xlifepp::Interpolation * >" tooltip="std::vector< xlifepp::Interpolation * >"] "53" [label="std::vector< xlifepp::Parameter * >" tooltip="std::vector< xlifepp::Parameter * >"] "3" [label="std::vector< xlifepp::Point * >" tooltip="std::vector< xlifepp::Point * >"] "32" [label="std::vector< xlifepp::Point >" tooltip="std::vector< xlifepp::Point >"] "21" [label="std::vector< xlifepp::PolynomialsBasisT >" tooltip="std::vector< xlifepp::PolynomialsBasisT >"] "16" [label="std::vector< xlifepp::RefDof * >" tooltip="std::vector< xlifepp::RefDof * >"] "17" [label="std::vector< xlifepp::RefElement * >" tooltip="std::vector< xlifepp::RefElement * >"] "44" [label="std::vector< xlifepp::Transformation * >" tooltip="std::vector< xlifepp::Transformation * >"] "47" [label="std::vector< xlifepp::Vector< real_t > >" tooltip="std::vector< xlifepp::Vector< real_t > >"] "29" [label="xlifepp::BoundingBox" tooltip="xlifepp::BoundingBox"] "42" [label="xlifepp::Collection< T >" tooltip="xlifepp::Collection< T >"] "60" [label="xlifepp::DomainInfo" tooltip="xlifepp::DomainInfo"] "58" [label="xlifepp::Extension" tooltip="xlifepp::Extension"] "57" [label="xlifepp::ExtensionData" tooltip="xlifepp::ExtensionData"] "35" [label="xlifepp::ExtrusionData" tooltip="xlifepp::ExtrusionData"] "45" [label="xlifepp::GeoNode" tooltip="xlifepp::GeoNode"] "59" [label="xlifepp::GeomDomain" tooltip="xlifepp::GeomDomain"] "1" [label="xlifepp::GeomMapData" tooltip="xlifepp::GeomMapData" fillcolor="#BFBFBF"] "9" [label="xlifepp::GeomRefElement" tooltip="xlifepp::GeomRefElement"] "28" [label="xlifepp::Geometry" tooltip="xlifepp::Geometry"] "14" [label="xlifepp::Interpolation" tooltip="xlifepp::Interpolation"] "38" [label="xlifepp::Matrix< K >" tooltip="xlifepp::Matrix< K >"] "48" [label="xlifepp::Mesh" tooltip="xlifepp::Mesh"] "2" [label="xlifepp::MeshElement" tooltip="xlifepp::MeshElement"] "31" [label="xlifepp::MinimalBox" tooltip="xlifepp::MinimalBox"] "52" [label="xlifepp::Parameters" tooltip="xlifepp::Parameters"] "27" [label="xlifepp::Parametrization" tooltip="xlifepp::Parametrization"] "7" [label="xlifepp::Point" tooltip="xlifepp::Point"] "18" [label="xlifepp::PolynomialsBasisT< K >" tooltip="xlifepp::PolynomialsBasisT< K >"] "8" [label="xlifepp::RefElement" tooltip="xlifepp::RefElement"] "36" [label="xlifepp::Transformation" tooltip="xlifepp::Transformation"] "25" [label="xlifepp::Vector< K >" tooltip="xlifepp::Vector< K >"] "40" -> "41" [dir=forward tooltip="public-inheritance"] "40" -> "42" [dir=forward tooltip="template-instance"] "37" -> "6" [dir=forward tooltip="public-inheritance"] "37" -> "38" [dir=forward tooltip="template-instance"] "24" -> "6" [dir=forward tooltip="public-inheritance"] "24" -> "25" [dir=forward tooltip="template-instance"] "46" -> "47" [dir=forward tooltip="public-inheritance"] "46" -> "25" [dir=forward tooltip="template-instance"] "55" -> "20" [dir=forward tooltip="template-instance"] "19" -> "20" [dir=forward tooltip="template-instance"] "34" -> "23" [dir=forward tooltip="template-instance"] "33" -> "23" [dir=forward tooltip="template-instance"] "54" -> "23" [dir=forward tooltip="template-instance"] "62" -> "23" [dir=forward tooltip="template-instance"] "56" -> "23" [dir=forward tooltip="template-instance"] "22" -> "23" [dir=forward tooltip="template-instance"] "26" -> "4" [dir=forward tooltip="template-instance"] "30" -> "4" [dir=forward tooltip="template-instance"] "10" -> "4" [dir=forward tooltip="template-instance"] "61" -> "4" [dir=forward tooltip="template-instance"] "39" -> "4" [dir=forward tooltip="template-instance"] "5" -> "4" [dir=forward tooltip="template-instance"] "6" -> "4" [dir=forward tooltip="template-instance"] "51" -> "4" [dir=forward tooltip="template-instance"] "12" -> "4" [dir=forward tooltip="template-instance"] "11" -> "4" [dir=forward tooltip="template-instance"] "63" -> "4" [dir=forward tooltip="template-instance"] "41" -> "4" [dir=forward tooltip="template-instance"] "50" -> "4" [dir=forward tooltip="template-instance"] "49" -> "4" [dir=forward tooltip="template-instance"] "13" -> "4" [dir=forward tooltip="template-instance"] "15" -> "4" [dir=forward tooltip="template-instance"] "53" -> "4" [dir=forward tooltip="template-instance"] "3" -> "4" [dir=forward tooltip="template-instance"] "32" -> "4" [dir=forward tooltip="template-instance"] "21" -> "4" [dir=forward tooltip="template-instance"] "16" -> "4" [dir=forward tooltip="template-instance"] "17" -> "4" [dir=forward tooltip="template-instance"] "44" -> "4" [dir=forward tooltip="template-instance"] "47" -> "4" [dir=forward tooltip="template-instance"] "29" -> "30" [dir=forward tooltip="usage"] "42" -> "43" [dir=forward tooltip="public-inheritance"] "60" -> "48" [dir=forward tooltip="usage"] "58" -> "59" [dir=forward tooltip="usage"] "58" -> "62" [dir=forward tooltip="usage"] "57" -> "58" [dir=forward tooltip="usage"] "57" -> "6" [dir=forward tooltip="usage"] "57" -> "63" [dir=forward tooltip="usage"] "57" -> "32" [dir=forward tooltip="usage"] "35" -> "36" [dir=forward tooltip="usage"] "35" -> "40" [dir=forward tooltip="usage"] "35" -> "5" [dir=forward tooltip="usage"] "35" -> "6" [dir=forward tooltip="usage"] "35" -> "24" [dir=forward tooltip="usage"] "35" -> "44" [dir=forward tooltip="usage"] "35" -> "7" [dir=forward tooltip="usage"] "45" -> "45" [dir=forward tooltip="usage"] "45" -> "28" [dir=forward tooltip="usage"] "59" -> "60" [dir=forward tooltip="usage"] "59" -> "59" [dir=forward tooltip="usage"] "59" -> "28" [dir=forward tooltip="usage"] "59" -> "61" [dir=forward tooltip="usage"] "1" -> "2" [dir=forward tooltip="usage"] "1" -> "7" [dir=forward tooltip="usage"] "1" -> "24" [dir=forward tooltip="usage"] "1" -> "27" [dir=forward tooltip="usage"] "1" -> "56" [dir=forward tooltip="usage"] "1" -> "57" [dir=forward tooltip="usage"] "1" -> "47" [dir=forward tooltip="usage"] "9" -> "6" [dir=forward tooltip="usage"] "9" -> "10" [dir=forward tooltip="usage"] "9" -> "11" [dir=forward tooltip="usage"] "9" -> "12" [dir=forward tooltip="usage"] "9" -> "13" [dir=forward tooltip="usage"] "28" -> "29" [dir=forward tooltip="usage"] "28" -> "31" [dir=forward tooltip="usage"] "28" -> "33" [dir=forward tooltip="usage"] "28" -> "34" [dir=forward tooltip="usage"] "28" -> "12" [dir=forward tooltip="usage"] "28" -> "35" [dir=forward tooltip="usage"] "28" -> "27" [dir=forward tooltip="usage"] "28" -> "28" [dir=forward tooltip="usage"] "28" -> "45" [dir=forward tooltip="usage"] "14" -> "15" [dir=forward tooltip="usage"] "38" -> "26" [dir=forward tooltip="public-inheritance"] "48" -> "28" [dir=forward tooltip="usage"] "48" -> "32" [dir=forward tooltip="usage"] "48" -> "49" [dir=forward tooltip="usage"] "48" -> "50" [dir=forward tooltip="usage"] "48" -> "5" [dir=forward tooltip="usage"] "48" -> "51" [dir=forward tooltip="usage"] "48" -> "48" [dir=forward tooltip="usage"] "2" -> "3" [dir=forward tooltip="usage"] "2" -> "5" [dir=forward tooltip="usage"] "2" -> "6" [dir=forward tooltip="usage"] "2" -> "7" [dir=forward tooltip="usage"] "2" -> "1" [dir=forward tooltip="usage"] "2" -> "2" [dir=forward tooltip="usage"] "2" -> "8" [dir=forward tooltip="usage"] "31" -> "32" [dir=forward tooltip="usage"] "52" -> "53" [dir=forward tooltip="usage"] "52" -> "54" [dir=forward tooltip="usage"] "27" -> "28" [dir=forward tooltip="usage"] "27" -> "24" [dir=forward tooltip="usage"] "27" -> "46" [dir=forward tooltip="usage"] "27" -> "48" [dir=forward tooltip="usage"] "27" -> "52" [dir=forward tooltip="usage"] "27" -> "55" [dir=forward tooltip="usage"] "7" -> "6" [dir=forward tooltip="public-inheritance"] "18" -> "19" [dir=forward tooltip="public-inheritance"] "8" -> "9" [dir=forward tooltip="usage"] "8" -> "14" [dir=forward tooltip="usage"] "8" -> "16" [dir=forward tooltip="usage"] "8" -> "17" [dir=forward tooltip="usage"] "8" -> "11" [dir=forward tooltip="usage"] "8" -> "18" [dir=forward tooltip="usage"] "8" -> "21" [dir=forward tooltip="usage"] "8" -> "22" [dir=forward tooltip="usage"] "36" -> "37" [dir=forward tooltip="usage"] "36" -> "24" [dir=forward tooltip="usage"] "36" -> "39" [dir=forward tooltip="usage"] "25" -> "26" [dir=forward tooltip="public-inheritance"] }

object handles useful computational data of a geometric element It stores jacobian, its determinant, … depending of a point

Public Functions

GeomMapData(const MeshElement*, bool withJ = false, bool withInvJ = false, bool withN = false)#

basic constructor (element centroid as point)

GeomMapData(const MeshElement*, const Point&, bool withJ = false, bool withInvJ = false, bool withN = false)#

basic constructor

GeomMapData(const MeshElement*, std::vector<real_t>::const_iterator, bool withJ = false, bool withInvJ = false, bool withN = false)#

basic constructor

void computeDifferentialElement()#

compute differential element assuming jacobian matrix is update

real_t computeJacobianDeterminant()#

compute jacobian determinant (jacobian already computed)

void computeJacobianMatrix(const ShapeValues &shv, number_t side = 0)#

compute jacobian matrix (shape functions given)

void computeJacobianMatrix(const std::vector<real_t>&, number_t side = 0)#

compute jacobian matrix at a point given as vector or Point

void computeJacobianMatrix(std::vector<real_t>::const_iterator, number_t side = 0)#

compute jacobian matrix at a point given by iterator

void computeMetricTensor()#

compute metric tensor for surface differential geometry

void computeNormalVector()#

computes normal vector at a point of element

void computeOrientedNormal()#

compute oriented unit normal vector, assuming jacobian matrix is update

void computeSurfaceGradient(real_t, real_t, std::vector<real_t>&)#

compute surface gradient gradS from 2D reference gradient (grad0, grad1)

void computeTangentVector()#

computes tangent (and bitangent) vector at a point of element

Matrix<real_t> contravariantPiolaMap(const Point& = Point())#

return the contravariant Piola map matrix at current point

return the contravariant Piola map matrix at current point: J/|J| when mapping from 1d->2d or 2d->3d the matrix J/|J| is replace by J/sqrt|Jt*J| note: in that case computeJacobianDeterminant() compute sqrt|Jt*J|

Matrix<real_t> covariantPiolaMap(const Point& = Point())#

return the covariant Piola map matrix at current point

return the covariant Piola map matrix at current point: J^{-T} when mapping from 1d->2d or 2d->3d the matrix J^{-T} is replace by J*(J{T}*J)^{-1}{T} note: in that case invertJacobian() compute J*(J{T}*J)^{-1}

inline const MeshElement &geomElement() const#

access to related MeshElement

Point geomMap(const ShapeValues &shv, number_t side)#

access to measure of element(s=0) or measure of side s>0

mapping from ref. elt onto geom. elt. Element (shape functions given, thread safe)

Point geomMap(const std::vector<real_t> &p, number_t side = 0)#

mapping from ref.

elt onto geom. elt. Element of point p

Point geomMap(std::vector<real_t>::const_iterator p, number_t side = 0)#

mapping from ref.

elt onto geom. elt. Element of point p

Point geomMapInverse(const std::vector<real_t>&, real_t eps = theTolerance, number_t maxIter = 50)#

mapping from elt onto ref.

elt. of point p

elt. (inverse of GeomMap) mapping may be non linear: order >1 or not simplicial element (quadrangle, hexahedron) so to inverse it we have to use an iterative method (Newton for instance) first step: initialize Newton’s method with first order map say P_0=(JF1)^-1(P-F1(0)) if element is a 1 order simplicial element return P_0 second step: Newton’s iteration P_k+1=P_k - JF(P_k)^-1(P-F(Pk)) stop if ||F(P_k+1)-X||< eps or ||P_k+1-P_k+1||< eps stop also if the number of iteration k>50 (error)

void invertJacobianMatrix()#

compute inverse of jacobian matrix (jacobian already computed)

void normalize()#

normalize normal vector to element

Point piolaMap(number_t side = 0)#

Piola mapping from ref.

elt onto geom. elt. Element (shape functions known)

void print(std::ostream&) const#

print GeoMapData

std::vector<Vector<real_t>> &sideNV() const#

get normals to the sides of the element

Public Members

Vector<real_t> bitangentVector#

second unit tangent vector at a given boundary point

Point currentPoint#

point in reference element used for last computation

real_t differentialElement#

differential element for elementary integrals (abs(jac.

det.))

dimen_t elementDim#

dim of reference space

ExtensionData *extdata#

pointer to additionnal data used by extension

Matrix<real_t> inverseJacobianMatrix#

inverse jacobian matrix (when needed)

const Parametrization *isoPar_p#

pointer to isoGeometric parametrization (default 0)

real_t jacobianDeterminant#

jacobian determinant

Matrix<real_t> jacobianMatrix#

jacobian matrix of mapping from reference element to current one

Matrix<real_t> metricTensor#

symetric metric tensor (t_i|t_j)

real_t metricTensorDeterminant#

metric_tensor determinant for differential geometry

Vector<real_t> normalVector#

unit outward normal vector at a given boundary point

std::map<Quadrature*, std::vector<Point>> phyPoints#

to store image of quadrature points in physical space

mutable std::vector<Vector<real_t>> sideNormalVectors#

normal to the sides of the element (build if required)

dimen_t spaceDim#

dim of physical space

Vector<real_t> tangentVector#

unit tangent vector at a given boundary point

bool useIsoNodes#

if true use MeshElement::isoNodes (default false)

bool useParametrization#

if true take into account isogeometric parametrization (default false)

Friends

friend std::ostream &operator<<(std::ostream&, const GeomMapData&)#

print operator