3D geometry information More...

#include <Geometry3D.h>

Inheritance diagram for Nektar::SpatialDomains::Geometry3D:

Collaboration diagram for Nektar::SpatialDomains::Geometry3D:

Public Member Functions
	Geometry3D ()

	Geometry3D (const int coordim)

virtual	~Geometry3D ()

int	GetEid (int i) const
	Return the ID of edge i in this element. More...

const Geometry1DSharedPtr	GetEdge (int i) const

Geometry2DSharedPtr	GetFace (int i)
	Return face i in this element. More...

int	GetDir (const int faceidx, const int facedir) const
	Returns the element coordinate direction corresponding to a given face coordinate direction. More...

Public Member Functions inherited from Nektar::SpatialDomains::Geometry
	Geometry ()

	Geometry (int coordim)

virtual	~Geometry ()

bool	IsElmtConnected (int gvo_id, int locid) const

void	AddElmtConnected (int gvo_id, int locid)

int	NumElmtConnected () const

int	GetCoordim () const

void	SetCoordim (int coordim)

GeomFactorsSharedPtr	GetGeomFactors ()

GeomFactorsSharedPtr	GetRefGeomFactors (const Array< OneD, const LibUtilities::BasisSharedPtr > &tbasis)

GeomFactorsSharedPtr	GetMetricInfo ()

LibUtilities::ShapeType	GetShapeType (void)

int	GetGlobalID (void)

void	SetGlobalID (int globalid)

int	GetVid (int i) const

int	GetEid (int i) const

int	GetFid (int i) const

int	GetTid (int i) const

int	GetNumVerts () const

PointGeomSharedPtr	GetVertex (int i) const

StdRegions::Orientation	GetEorient (const int i) const

StdRegions::Orientation	GetPorient (const int i) const

StdRegions::Orientation	GetForient (const int i) const

int	GetNumEdges () const

int	GetNumFaces () const

int	GetShapeDim () const

StdRegions::StdExpansionSharedPtr	GetXmap () const

const Array< OneD, const NekDouble > &	GetCoeffs (const int i) const

bool	ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, NekDouble tol=0.0)

bool	ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, Array< OneD, NekDouble > &locCoord, NekDouble tol)

bool	ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, Array< OneD, NekDouble > &locCoord, NekDouble tol, NekDouble &resid)

int	GetVertexEdgeMap (int i, int j) const

int	GetVertexFaceMap (int i, int j) const
	return the id of the $j^{th}$ face attached to the $i^{th}$ vertex More...

int	GetEdgeFaceMap (int i, int j) const

void	FillGeom ()
	Put all quadrature information into face/edge structure and backward transform. More...

NekDouble	GetLocCoords (const Array< OneD, const NekDouble > &coords, Array< OneD, NekDouble > &Lcoords)

NekDouble	GetCoord (const int i, const Array< OneD, const NekDouble > &Lcoord)
	Given local collapsed coordinate Lcoord return the value of physical coordinate in direction i. More...

void	SetOwnData ()

const LibUtilities::BasisSharedPtr	GetBasis (const int i)
	Return the j-th basis of the i-th co-ordinate dimension. More...

const LibUtilities::PointsKeyVector	GetPointsKeys ()

void	Reset (CurveMap &curvedEdges, CurveMap &curvedFaces)

Protected Member Functions
void	NewtonIterationForLocCoord (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &ptsx, const Array< OneD, const NekDouble > &ptsy, const Array< OneD, const NekDouble > &ptsz, Array< OneD, NekDouble > &Lcoords, NekDouble &resid)

virtual void	v_FillGeom ()
	Put all quadrature information into face/edge structure and backward transform. More...

virtual NekDouble	v_GetCoord (const int i, const Array< OneD, const NekDouble > &Lcoord)
	Given local collapsed coordinate Lcoord return the value of physical coordinate in direction i. More...

virtual void	v_GenGeomFactors ()

virtual int	v_GetShapeDim () const
	Return the dimension of this element. More...

virtual int	v_GetVid (int i) const
	Return the vertex ID of vertex i. More...

virtual PointGeomSharedPtr	v_GetVertex (int i) const
	Return vertex i in this element. More...

virtual const SegGeomSharedPtr	v_GetEdge (int i) const
	Return edge i of this element. More...

virtual StdRegions::Orientation	v_GetEorient (const int i) const
	Return the orientation of edge i in this element. More...

virtual int	v_GetEid (int i) const
	Return the ID of edge i in this element. More...

virtual const Geometry2DSharedPtr	v_GetFace (int i) const
	Return face i in this element. More...

virtual StdRegions::Orientation	v_GetForient (const int i) const
	Return the orientation of face i in this element. More...

virtual int	v_GetFid (int i) const
	Return the ID of face i in this element. More...

virtual int	v_GetEid () const
	Return the ID of this element. More...

virtual int	v_WhichEdge (SegGeomSharedPtr edge)
	Return the local ID of a given edge. More...

virtual int	v_WhichFace (Geometry2DSharedPtr face)
	Return the local ID of a given face. More...

virtual const LibUtilities::BasisSharedPtr	v_GetBasis (const int i)
	Return the j-th basis of the i-th co-ordinate dimension. More...

virtual int	v_GetDir (const int faceidx, const int facedir) const =0

virtual void	v_AddElmtConnected (int gvo_id, int locid)

virtual bool	v_IsElmtConnected (int gvo_id, int locid) const

virtual int	v_NumElmtConnected () const

virtual void	v_SetOwnData ()

Protected Member Functions inherited from Nektar::SpatialDomains::Geometry
void	GenGeomFactors ()

virtual int	v_GetNumVerts () const

virtual StdRegions::Orientation	v_GetPorient (const int i) const

virtual int	v_GetNumEdges () const

virtual int	v_GetNumFaces () const

virtual StdRegions::StdExpansionSharedPtr	v_GetXmap () const

virtual int	v_GetCoordim () const

virtual bool	v_ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, NekDouble tol=0.0)

virtual bool	v_ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, Array< OneD, NekDouble > &locCoord, NekDouble tol)

virtual bool	v_ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, Array< OneD, NekDouble > &locCoord, NekDouble tol, NekDouble &resid)

virtual int	v_GetVertexEdgeMap (int i, int j) const

virtual int	v_GetVertexFaceMap (int i, int j) const

virtual int	v_GetEdgeFaceMap (int i, int j) const

virtual NekDouble	v_GetLocCoords (const Array< OneD, const NekDouble > &coords, Array< OneD, NekDouble > &Lcoords)

virtual void	v_Reset (CurveMap &curvedEdges, CurveMap &curvedFaces)
	Reset this geometry object: unset the current state and remove allocated GeomFactors. More...

void	SetUpCoeffs (const int nCoeffs)
	Initialise the m_coeffs array. More...

Protected Attributes
PointGeomVector	m_verts

SegGeomVector	m_edges

Geometry2DVector	m_faces

std::vector < StdRegions::Orientation >	m_eorient

std::vector < StdRegions::Orientation >	m_forient

std::list< CompToElmt >	m_elmtmap

bool	m_owndata

int	m_eid

bool	m_ownverts

Protected Attributes inherited from Nektar::SpatialDomains::Geometry
int	m_coordim
	coordinate dimension More...

GeomFactorsSharedPtr	m_geomFactors

GeomState	m_geomFactorsState

StdRegions::StdExpansionSharedPtr	m_xmap

GeomState	m_state
	enum identifier to determine if quad points are filled More...

GeomType	m_geomType

LibUtilities::ShapeType	m_shapeType

int	m_globalID

Array< OneD, Array< OneD, NekDouble > >	m_coeffs

Additional Inherited Members
Static Protected Member Functions inherited from Nektar::SpatialDomains::Geometry
static GeomFactorsSharedPtr	ValidateRegGeomFactor (GeomFactorsSharedPtr geomFactor)

Static Protected Attributes inherited from Nektar::SpatialDomains::Geometry
static GeomFactorsVector	m_regGeomFactorsManager

Detailed Description

3D geometry information

Definition at line 70 of file Geometry3D.h.

Constructor & Destructor Documentation

Nektar::SpatialDomains::Geometry3D::Geometry3D ( )

Definition at line 51 of file Geometry3D.cpp.

52 {

53 }

Nektar::SpatialDomains::Geometry3D::Geometry3D ( const int coordim )

Definition at line 55 of file Geometry3D.cpp.

References ASSERTL0, and Nektar::SpatialDomains::Geometry::m_coordim.

                                            :
       Geometry(coordim)
     {
       ASSERTL0(m_coordim > 2,
                "Coordinate dimension should be at least 3 for a 3D geometry.");
     }

Nektar::SpatialDomains::Geometry3D::~Geometry3D ( )

virtual

Definition at line 62 of file Geometry3D.cpp.

63 {

64 }

Member Function Documentation

int Nektar::SpatialDomains::Geometry3D::GetDir	(	const int	faceidx,
		const int	facedir
	)		const

Returns the element coordinate direction corresponding to a given face coordinate direction.

Definition at line 91 of file Geometry3D.cpp.

References v_GetDir().

Referenced by Nektar::SpatialDomains::MeshGraph3D::GetFaceBasisKey().

     {
       return v_GetDir(faceidx, facedir);
     }

const Geometry1DSharedPtr Nektar::SpatialDomains::Geometry3D::GetEdge ( int i ) const

Referenced by Nektar::SpatialDomains::PyrGeom::SetUpLocalEdges(), Nektar::SpatialDomains::TetGeom::SetUpLocalEdges(), Nektar::SpatialDomains::PrismGeom::SetUpLocalEdges(), and Nektar::SpatialDomains::HexGeom::SetUpLocalEdges().

int Nektar::SpatialDomains::Geometry3D::GetEid ( int i ) const

Return the ID of edge i in this element.

Definition at line 74 of file Geometry3D.cpp.

References v_GetEid().

Referenced by Nektar::SpatialDomains::PyrGeom::SetUpLocalEdges(), Nektar::SpatialDomains::TetGeom::SetUpLocalEdges(), Nektar::SpatialDomains::PrismGeom::SetUpLocalEdges(), and Nektar::SpatialDomains::HexGeom::SetUpLocalEdges().

     {
       return v_GetEid(i);
     }

Geometry2DSharedPtr Nektar::SpatialDomains::Geometry3D::GetFace ( int i )

Return face i in this element.

Definition at line 82 of file Geometry3D.cpp.

References v_GetFace().

Referenced by Nektar::SolverUtils::UpdateGeometry().

     {
       return v_GetFace(i);
     }

void Nektar::SpatialDomains::Geometry3D::NewtonIterationForLocCoord	(	const Array< OneD, const NekDouble > &	coords,
		const Array< OneD, const NekDouble > &	ptsx,
		const Array< OneD, const NekDouble > &	ptsy,
		const Array< OneD, const NekDouble > &	ptsz,
		Array< OneD, NekDouble > &	Lcoords,
		NekDouble &	resid
	)

protected

Definition at line 99 of file Geometry3D.cpp.

References Nektar::SpatialDomains::Geometry::m_geomFactors, Nektar::SpatialDomains::Geometry::m_xmap, Vmath::Vsum(), and WARNINGL1.

Referenced by Nektar::SpatialDomains::TetGeom::v_GetLocCoords(), Nektar::SpatialDomains::PrismGeom::v_GetLocCoords(), and Nektar::SpatialDomains::HexGeom::v_GetLocCoords().

     {
         // maximum iterations for convergence
         const int MaxIterations   = 51;
         // |x-xp|^2 < EPSILON  error tolerance
         const NekDouble Tol       = 1.e-8;
         // |r,s|    > LcoordDIV stop the search
         const NekDouble LcoordDiv = 15.0;
 
         Array<OneD, const NekDouble > Jac =
                             m_geomFactors->GetJac(m_xmap->GetPointsKeys());
 
         NekDouble ScaledTol = Vmath::Vsum(Jac.num_elements(), Jac, 1) /
                                         ((NekDouble)Jac.num_elements());
         ScaledTol *= Tol;
 
         NekDouble xmap,ymap,zmap, F1,F2, F3;
 
         NekDouble derx_1, derx_2, derx_3, dery_1, dery_2, dery_3,
                   derz_1, derz_2, derz_3, jac;
 
         // save intiial guess for later reference if required.
         NekDouble init0 = Lcoords[0], init1 = Lcoords[1], init2 = Lcoords[2];
 
         Array<OneD, NekDouble> DxD1(ptsx.num_elements());
         Array<OneD, NekDouble> DxD2(ptsx.num_elements());
         Array<OneD, NekDouble> DxD3(ptsx.num_elements());
         Array<OneD, NekDouble> DyD1(ptsx.num_elements());
         Array<OneD, NekDouble> DyD2(ptsx.num_elements());
         Array<OneD, NekDouble> DyD3(ptsx.num_elements());
         Array<OneD, NekDouble> DzD1(ptsx.num_elements());
         Array<OneD, NekDouble> DzD2(ptsx.num_elements());
         Array<OneD, NekDouble> DzD3(ptsx.num_elements());
 
         // Ideally this will be stored in m_geomfactors
         m_xmap->PhysDeriv(ptsx,DxD1,DxD2,DxD3);
         m_xmap->PhysDeriv(ptsy,DyD1,DyD2,DyD3);
         m_xmap->PhysDeriv(ptsz,DzD1,DzD2,DzD3);
 
         int cnt=0;
         Array<OneD, DNekMatSharedPtr > I(3);
         Array<OneD, NekDouble>       eta(3);
 
         F1 = F2 = F3 = 2000; // Starting value of Function
 
         while(cnt++ < MaxIterations)
         {
             //  evaluate lagrange interpolant at Lcoords
             m_xmap->LocCoordToLocCollapsed(Lcoords,eta);
             I[0] = m_xmap->GetBasis(0)->GetI(eta  );
             I[1] = m_xmap->GetBasis(1)->GetI(eta+1);
             I[2] = m_xmap->GetBasis(2)->GetI(eta+2);
 
             //calculate the global point `corresponding to Lcoords
             xmap = m_xmap->PhysEvaluate(I, ptsx);
             ymap = m_xmap->PhysEvaluate(I, ptsy);
             zmap = m_xmap->PhysEvaluate(I, ptsz);
 
             F1 = coords[0] - xmap;
             F2 = coords[1] - ymap;
             F3 = coords[2] - zmap;
 
             if(F1*F1 + F2*F2 + F3*F3 < ScaledTol)
             {
               resid = sqrt(F1*F1 + F2*F2 + F3*F3);
               break;
             }
 
             //Interpolate derivative metric at Lcoords
             derx_1 = m_xmap->PhysEvaluate(I, DxD1);
             derx_2 = m_xmap->PhysEvaluate(I, DxD2);
             derx_3 = m_xmap->PhysEvaluate(I, DxD3);
             dery_1 = m_xmap->PhysEvaluate(I, DyD1);
             dery_2 = m_xmap->PhysEvaluate(I, DyD2);
             dery_3 = m_xmap->PhysEvaluate(I, DyD3);
             derz_1 = m_xmap->PhysEvaluate(I, DzD1);
             derz_2 = m_xmap->PhysEvaluate(I, DzD2);
             derz_3 = m_xmap->PhysEvaluate(I, DzD3);
 
             jac = derx_1*(dery_2*derz_3 - dery_3*derz_2)
                 - derx_2*(dery_1*derz_3 - dery_3*derz_1)
                 + derx_3*(dery_1*derz_2 - dery_2*derz_1);
 
             // use analytical inverse of derivitives which are also similar to
             // those of metric factors.
             Lcoords[0] = Lcoords[0]
                              +((dery_2*derz_3 - dery_3*derz_2)*(coords[0]-xmap)
                              - (derx_2*derz_3 - derx_3*derz_2)*(coords[1]-ymap)
                              + (derx_2*dery_3 - derx_3*dery_2)*(coords[2]-zmap)
                                )/jac;
 
             Lcoords[1] = Lcoords[1]
                              -((dery_1*derz_3 - dery_3*derz_1)*(coords[0]-xmap)
                              - (derx_1*derz_3 - derx_3*derz_1)*(coords[1]-ymap)
                              + (derx_1*dery_3 - derx_3*dery_1)*(coords[2]-zmap)
                                )/jac;
 
             Lcoords[2] = Lcoords[2]
                              +((dery_1*derz_2 - dery_2*derz_1)*(coords[0]-xmap)
                              - (derx_1*derz_2 - derx_2*derz_1)*(coords[1]-ymap)
                              + (derx_1*dery_2 - derx_2*dery_1)*(coords[2]-zmap)
                                )/jac;
 
             if (fabs(Lcoords[0]) > LcoordDiv || fabs(Lcoords[1]) > LcoordDiv ||
                 fabs(Lcoords[0]) > LcoordDiv)
             {
                 break; // lcoords have diverged so stop iteration
             }
         }
 
         resid = sqrt(F1*F1 + F2*F2 + F3*F3);
 
         if(cnt >= MaxIterations)
         {
             Array<OneD, NekDouble> collCoords(3);
             m_xmap->LocCoordToLocCollapsed(Lcoords,collCoords);
 
             // if coordinate is inside element dump error!
             if((collCoords[0] >=  -1.0 && collCoords[0] <= 1.0)&&
                (collCoords[1] >=  -1.0 && collCoords[1] <= 1.0)&&
                (collCoords[2] >=  -1.0 && collCoords[2] <= 1.0))
             {
                 std::ostringstream ss;
 
                 ss << "Reached MaxIterations (" << MaxIterations << ") in Newton iteration ";
                 ss << "Init value ("<< setprecision(4) << init0 << "," << init1<< "," << init2 <<") ";
                 ss << "Fin  value ("<<Lcoords[0] << "," << Lcoords[1]<< "," << Lcoords[2] << ") ";
                 ss << "Resid = " << resid << " Tolerance = " << sqrt(ScaledTol) ;
 
                 WARNINGL1(cnt < MaxIterations,ss.str());
             }
         }
     }

void Nektar::SpatialDomains::Geometry3D::v_AddElmtConnected	(	int	gvo_id,
		int	locid
	)

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 502 of file Geometry3D.cpp.

References m_elmtmap.

     {
       CompToElmt ee(gvo_id,locid);
       m_elmtmap.push_back(ee);
     }

void Nektar::SpatialDomains::Geometry3D::v_FillGeom ( )

protectedvirtual

Put all quadrature information into face/edge structure and backward transform.

Note verts, edges, and faces are listed according to anticlockwise convention but points in _coeffs have to be in array format from left to right.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 247 of file Geometry3D.cpp.

References Nektar::SpatialDomains::ePtsFilled, Nektar::SpatialDomains::Geometry::GetXmap(), Nektar::SpatialDomains::Geometry::m_coeffs, Nektar::SpatialDomains::Geometry::m_coordim, m_faces, m_forient, Nektar::SpatialDomains::Geometry::m_state, and Nektar::SpatialDomains::Geometry::m_xmap.

     {
       if (m_state == ePtsFilled)
           return;
 
       int i,j,k;
 
       for(i = 0; i < m_forient.size(); i++)
       {
           m_faces[i]->FillGeom();
 
           int nFaceCoeffs = m_faces[i]->GetXmap()->GetNcoeffs();
 
           Array<OneD, unsigned int> mapArray (nFaceCoeffs);
           Array<OneD,          int> signArray(nFaceCoeffs);
 
           if (m_forient[i] < 9)
           {
               m_xmap->GetFaceToElementMap(
                   i, m_forient[i], mapArray, signArray,
                   m_faces[i]->GetXmap()->GetEdgeNcoeffs(0),
                   m_faces[i]->GetXmap()->GetEdgeNcoeffs(1));
           }
           else
           {
               m_xmap->GetFaceToElementMap(
                   i, m_forient[i], mapArray, signArray,
                   m_faces[i]->GetXmap()->GetEdgeNcoeffs(1),
                   m_faces[i]->GetXmap()->GetEdgeNcoeffs(0));
           }
 
           for (j = 0; j < m_coordim; j++)
           {
               const Array<OneD, const NekDouble> &coeffs =
                   m_faces[i]->GetCoeffs(j);
 
               for (k = 0; k < nFaceCoeffs; k++)
               {
                   NekDouble v = signArray[k] * coeffs[k];
                   m_coeffs[j][mapArray[k]] = v;
               }
           }
       }
 
       m_state = ePtsFilled;
     }

void Nektar::SpatialDomains::Geometry3D::v_GenGeomFactors ( )

protectedvirtual

Generate the geometry factors for this element.

Implements Nektar::SpatialDomains::Geometry.

Reimplemented in Nektar::SpatialDomains::HexGeom, Nektar::SpatialDomains::PrismGeom, and Nektar::SpatialDomains::PyrGeom.

Definition at line 297 of file Geometry3D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::SpatialDomains::eDeformed, Nektar::SpatialDomains::ePtsFilled, Nektar::SpatialDomains::eRegular, Nektar::SpatialDomains::Geometry::m_coeffs, Nektar::SpatialDomains::Geometry::m_coordim, Nektar::SpatialDomains::Geometry::m_geomFactors, Nektar::SpatialDomains::Geometry::m_geomFactorsState, Nektar::SpatialDomains::Geometry::m_xmap, and v_FillGeom().

     {
         if (m_geomFactorsState != ePtsFilled)
         {
             GeomType      Gtype  = eRegular;
 
             v_FillGeom();
 
             // check to see if expansions are linear
             for(int i = 0; i < m_coordim; ++i)
             {
                 if (m_xmap->GetBasisNumModes(0) != 2 ||
                     m_xmap->GetBasisNumModes(1) != 2 ||
                     m_xmap->GetBasisNumModes(2) != 2)
                 {
                     Gtype = eDeformed;
                 }
             }
 
             m_geomFactors = MemoryManager<GeomFactors>::AllocateSharedPtr(
                 Gtype, m_coordim, m_xmap, m_coeffs);
             m_geomFactorsState = ePtsFilled;
         }
     }

const LibUtilities::BasisSharedPtr Nektar::SpatialDomains::Geometry3D::v_GetBasis ( const int i )

protectedvirtual

Return the j-th basis of the i-th co-ordinate dimension.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 446 of file Geometry3D.cpp.

References Nektar::SpatialDomains::Geometry::m_xmap.

     {
       return m_xmap->GetBasis(i);
     }

NekDouble Nektar::SpatialDomains::Geometry3D::v_GetCoord	(	const int	i,
		const Array< OneD, const NekDouble > &	Lcoord
	)

protectedvirtual

Given local collapsed coordinate Lcoord return the value of physical coordinate in direction i.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 326 of file Geometry3D.cpp.

References ASSERTL1, Nektar::SpatialDomains::ePtsFilled, Nektar::SpatialDomains::Geometry::m_coeffs, Nektar::SpatialDomains::Geometry::m_state, and Nektar::SpatialDomains::Geometry::m_xmap.

     {
       ASSERTL1(m_state == ePtsFilled,
                "Geometry is not in physical space");
 
       Array<OneD, NekDouble> tmp(m_xmap->GetTotPoints());
       m_xmap->BwdTrans(m_coeffs[i], tmp);
 
       return m_xmap->PhysEvaluate(Lcoord, tmp);
     }

virtual int Nektar::SpatialDomains::Geometry3D::v_GetDir	(	const int	faceidx,
		const int	facedir
	)		const

protectedpure virtual

Implemented in Nektar::SpatialDomains::HexGeom, Nektar::SpatialDomains::PrismGeom, Nektar::SpatialDomains::TetGeom, and Nektar::SpatialDomains::PyrGeom.

Referenced by GetDir().

const SegGeomSharedPtr Nektar::SpatialDomains::Geometry3D::v_GetEdge ( int i ) const

protectedvirtual

Return edge i of this element.

Definition at line 373 of file Geometry3D.cpp.

References ASSERTL2, and m_edges.

     {
       ASSERTL2(i >= 0 && i <= m_edges.size() - 1,
                "Edge ID must be between 0 and "+
                boost::lexical_cast<string>(m_edges.size() - 1));
       return m_edges[i];
     }

int Nektar::SpatialDomains::Geometry3D::v_GetEid ( int i ) const

protectedvirtual

Return the ID of edge i in this element.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 396 of file Geometry3D.cpp.

References ASSERTL2, and m_edges.

     {
       ASSERTL2(i >= 0 && i <= m_edges.size() - 1,
                "Edge ID must be between 0 and "+
                boost::lexical_cast<string>(m_edges.size() - 1));
       return m_edges[i]->GetEid();
     }

int Nektar::SpatialDomains::Geometry3D::v_GetEid ( ) const

protectedvirtual

Return the ID of this element.

Definition at line 438 of file Geometry3D.cpp.

References m_eid.

Referenced by GetEid().

     {
       return m_eid;
     }

StdRegions::Orientation Nektar::SpatialDomains::Geometry3D::v_GetEorient ( const int i ) const

inlineprotectedvirtual

Return the orientation of edge i in this element.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 384 of file Geometry3D.cpp.

References ASSERTL2, m_edges, and m_eorient.

     {
       ASSERTL2(i >= 0 && i <= m_edges.size() - 1,
                "Edge ID must be between 0 and "+
                boost::lexical_cast<string>(m_edges.size() - 1));
       return m_eorient[i];
     }

const Geometry2DSharedPtr Nektar::SpatialDomains::Geometry3D::v_GetFace ( int i ) const

protectedvirtual

Return face i in this element.

Definition at line 407 of file Geometry3D.cpp.

References ASSERTL2, and m_faces.

Referenced by GetFace().

     {
       ASSERTL2((i >=0) && (i <= 5),"Edge id must be between 0 and 4");
       return m_faces[i];
     }

int Nektar::SpatialDomains::Geometry3D::v_GetFid ( int i ) const

protectedvirtual

Return the ID of face i in this element.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 427 of file Geometry3D.cpp.

References ASSERTL2, and m_faces.

     {
       ASSERTL2(i >= 0 && i <= m_faces.size() - 1,
                "Face ID must be between 0 and "+
                boost::lexical_cast<string>(m_faces.size() - 1));
       return m_faces[i]->GetFid();
     }

StdRegions::Orientation Nektar::SpatialDomains::Geometry3D::v_GetForient ( const int i ) const

protectedvirtual

Return the orientation of face i in this element.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 416 of file Geometry3D.cpp.

References ASSERTL2, m_faces, and m_forient.

     {
       ASSERTL2(i >= 0 && i <= m_faces.size() - 1,
                "Face ID must be between 0 and "+
                boost::lexical_cast<string>(m_faces.size() - 1));
       return m_forient[i];
     }

int Nektar::SpatialDomains::Geometry3D::v_GetShapeDim ( ) const

protectedvirtual

Return the dimension of this element.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 346 of file Geometry3D.cpp.

     {
       return 3;
     }

PointGeomSharedPtr Nektar::SpatialDomains::Geometry3D::v_GetVertex ( int i ) const

protectedvirtual

Return vertex i in this element.

Implements Nektar::SpatialDomains::Geometry.

Definition at line 365 of file Geometry3D.cpp.

References m_verts.

     {
         return m_verts[i];
     }

int Nektar::SpatialDomains::Geometry3D::v_GetVid ( int i ) const

protectedvirtual

Return the vertex ID of vertex i.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 354 of file Geometry3D.cpp.

References ASSERTL2, and m_verts.

     {
       ASSERTL2(i >= 0 && i <= m_verts.size() - 1,
                "Vertex ID must be between 0 and "+
                boost::lexical_cast<string>(m_verts.size() - 1));
       return m_verts[i]->GetVid();
     }

bool Nektar::SpatialDomains::Geometry3D::v_IsElmtConnected	(	int	gvo_id,
		int	locid
	)		const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 513 of file Geometry3D.cpp.

References Nektar::StdRegions::find(), and m_elmtmap.

     {
       std::list<CompToElmt>::const_iterator def;
       CompToElmt ee(gvo_id,locid);
 
       def = find(m_elmtmap.begin(),m_elmtmap.end(),ee);
 
       // Found the element connectivity object in the list
       return (def != m_elmtmap.end());
     }

int Nektar::SpatialDomains::Geometry3D::v_NumElmtConnected ( ) const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 508 of file Geometry3D.cpp.

References m_elmtmap.

     {
       return int(m_elmtmap.size());
     }

void Nektar::SpatialDomains::Geometry3D::v_SetOwnData ( )

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 524 of file Geometry3D.cpp.

References m_owndata.

     {
       m_owndata = true;
     }

int Nektar::SpatialDomains::Geometry3D::v_WhichEdge ( SegGeomSharedPtr edge )

protectedvirtual

Return the local ID of a given edge.

The local ID of an edge is a number between 0 and the number of edges in this element. If the edge is not found, this function returns -1.

Definition at line 458 of file Geometry3D.cpp.

References Nektar::iterator, and m_edges.

     {
       int returnval = -1;
 
       SegGeomVector::iterator edgeIter;
       int i;
 
       for (i=0,edgeIter = m_edges.begin(); edgeIter != m_edges.end(); ++edgeIter,++i)
       {
           if (*edgeIter == edge)
           {
               returnval = i;
               break;
           }
       }
 
       return returnval;
     }

int Nektar::SpatialDomains::Geometry3D::v_WhichFace ( Geometry2DSharedPtr face )

protectedvirtual

Return the local ID of a given face.

The local ID of a face is a number between 0 and the number of faces in this element. If the face is not found, this function returns -1.

Definition at line 483 of file Geometry3D.cpp.

References Nektar::iterator, and m_faces.

     {
       int i = 0;
 
       Geometry2DVector::iterator f;
       for (i = 0, f = m_faces.begin(); f != m_faces.end(); ++f,++i)
       {
           if (*f == face)
           {
               break;
           }
       }
       return i;
     }

Member Data Documentation

SegGeomVector Nektar::SpatialDomains::Geometry3D::m_edges

protected

Definition at line 90 of file Geometry3D.h.

int Nektar::SpatialDomains::Geometry3D::m_eid

protected

Definition at line 96 of file Geometry3D.h.

Referenced by v_GetEid().

std::list<CompToElmt> Nektar::SpatialDomains::Geometry3D::m_elmtmap

protected

Definition at line 94 of file Geometry3D.h.

Referenced by v_AddElmtConnected(), v_IsElmtConnected(), and v_NumElmtConnected().

std::vector<StdRegions::Orientation> Nektar::SpatialDomains::Geometry3D::m_eorient

protected

Definition at line 92 of file Geometry3D.h.

Referenced by Nektar::SpatialDomains::HexGeom::HexGeom(), Nektar::SpatialDomains::PrismGeom::PrismGeom(), Nektar::SpatialDomains::PyrGeom::PyrGeom(), Nektar::SpatialDomains::PyrGeom::SetUpEdgeOrientation(), Nektar::SpatialDomains::TetGeom::SetUpEdgeOrientation(), Nektar::SpatialDomains::PrismGeom::SetUpEdgeOrientation(), Nektar::SpatialDomains::HexGeom::SetUpEdgeOrientation(), Nektar::SpatialDomains::TetGeom::TetGeom(), and v_GetEorient().

Geometry2DVector Nektar::SpatialDomains::Geometry3D::m_faces

protected

Definition at line 91 of file Geometry3D.h.

std::vector<StdRegions::Orientation> Nektar::SpatialDomains::Geometry3D::m_forient

protected

Definition at line 93 of file Geometry3D.h.

bool Nektar::SpatialDomains::Geometry3D::m_owndata

protected

Definition at line 95 of file Geometry3D.h.

Referenced by v_SetOwnData().

bool Nektar::SpatialDomains::Geometry3D::m_ownverts

protected

Definition at line 97 of file Geometry3D.h.

PointGeomVector Nektar::SpatialDomains::Geometry3D::m_verts

protected

Definition at line 89 of file Geometry3D.h.

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation