#include <SegGeom.h>

Inheritance diagram for Nektar::SpatialDomains::SegGeom:

Collaboration diagram for Nektar::SpatialDomains::SegGeom:

Public Member Functions
	SegGeom ()

	SegGeom (int id, const int coordim)

	SegGeom (int id, const int coordim, const PointGeomSharedPtr vertex[])

	SegGeom (int id, const int coordim, const PointGeomSharedPtr vertex[], const CurveSharedPtr &curve)

	SegGeom (const int id, const PointGeomSharedPtr &vert1, const PointGeomSharedPtr &vert2)

	SegGeom (const SegGeom &in)

SegGeomSharedPtr	GenerateOneSpaceDimGeom (void)
	Generate a one dimensional space segment geometry where the vert[0] has the same x value and vert[1] is set to vert[0] plus the length of the original segment. More...

	~SegGeom ()

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

StdRegions::StdExpansion1DSharedPtr	operator[] (const int i) const

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

	Geometry1D (const int coordim)

virtual	~Geometry1D ()

PointGeomSharedPtr	GetVertex (const int i) const

LibUtilities::ShapeType	DetShapeType () const

int	GetEid () const

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)

Static Public Member Functions
static StdRegions::Orientation	GetEdgeOrientation (const SegGeom &edge1, const SegGeom &edge2)
	Get the orientation of edge1. More...

Static Public Attributes
static const int	kNverts = 2

static const int	kNedges = 1

Protected Member Functions
virtual int	v_GetShapeDim () const

virtual int	v_GetVid (int i) const

virtual PointGeomSharedPtr	v_GetVertex (const int i) const

virtual int	v_GetEid () const

virtual const LibUtilities::BasisSharedPtr	v_GetBasis (const int i)

virtual StdRegions::StdExpansionSharedPtr	v_GetXmap () const

virtual void	v_SetOwnData ()

virtual void	v_AddElmtConnected (int gvo_id, int locid)

virtual int	v_NumElmtConnected () const

virtual bool	v_IsElmtConnected (int gvoId, int locId) const

virtual LibUtilities::ShapeType	v_DetShapeType () const

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

virtual void	v_GenGeomFactors ()
	Set up GeoFac for this geometry using Coord quadrature distribution. More...

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

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

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

virtual NekDouble	v_GetCoord (const int i, const Array< OneD, const NekDouble > &Lcoord)

virtual int	v_GetNumVerts () const

virtual int	v_GetNumEdges () const

virtual bool	v_ContainsPoint (const Array< OneD, const NekDouble > &gloCoord, NekDouble tol=0.0)
	Determines if a point specified in global coordinates is located within this tetrahedral geometry. More...

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)

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

virtual int	v_GetEid (int i) const

virtual int	v_GetFid (int i) const

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

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

virtual int	v_GetNumFaces () const

virtual int	v_GetCoordim () const

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

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

Protected Attributes
int	m_eid

std::list< CompToElmt >	m_elmtMap

SpatialDomains::PointGeomSharedPtr	m_verts [kNverts]

StdRegions::Orientation	m_porient [kNverts]

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

Private Member Functions
void	SetUpXmap ()

Private Attributes
bool	m_ownData
	Boolean indicating whether object owns the data. More...

CurveSharedPtr	m_curve

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

Definition at line 57 of file SegGeom.h.

Constructor & Destructor Documentation

Nektar::SpatialDomains::SegGeom::SegGeom ( )

Definition at line 48 of file SegGeom.cpp.

References Nektar::LibUtilities::eSegment, and Nektar::SpatialDomains::Geometry::m_shapeType.

         {
             m_shapeType = LibUtilities::eSegment;
         }

Nektar::SpatialDomains::SegGeom::SegGeom	(	int	id,
		const int	coordim
	)

Definition at line 53 of file SegGeom.cpp.

References Nektar::LibUtilities::eSegment, m_eid, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_shapeType, Nektar::SpatialDomains::Geometry::m_xmap, Nektar::SpatialDomains::Geometry::SetUpCoeffs(), and SetUpXmap().

                                                  :
             Geometry1D(coordim)
         {
 
             m_shapeType = LibUtilities::eSegment;
             m_eid = id;
             m_globalID = id;
             SetUpXmap();
             SetUpCoeffs(m_xmap->GetNcoeffs());
         }

Nektar::SpatialDomains::SegGeom::SegGeom	(	int	id,
		const int	coordim,
		const PointGeomSharedPtr	vertex[]
	)

Definition at line 64 of file SegGeom.cpp.

References Nektar::SpatialDomains::eNotFilled, Nektar::LibUtilities::eSegment, m_eid, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_shapeType, Nektar::SpatialDomains::Geometry::m_state, m_verts, Nektar::SpatialDomains::Geometry::m_xmap, Nektar::SpatialDomains::Geometry::SetUpCoeffs(), and SetUpXmap().

                                                   :
             Geometry1D(coordim)
         {
             m_shapeType = LibUtilities::eSegment;
             m_eid   = id;
             m_globalID = id;
             m_state = eNotFilled;
 
             if (coordim > 0)
             {
                 SetUpXmap();
                 SetUpCoeffs(m_xmap->GetNcoeffs());
             }
 
             m_verts[0] = vertex[0];
             m_verts[1] = vertex[1];
         }

Nektar::SpatialDomains::SegGeom::SegGeom	(	int	id,
		const int	coordim,
		const PointGeomSharedPtr	vertex[],
		const CurveSharedPtr &	curve
	)

Definition at line 85 of file SegGeom.cpp.

References Nektar::NekMeshUtils::curve, Nektar::SpatialDomains::eNotFilled, Nektar::LibUtilities::eSegment, m_curve, m_eid, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_shapeType, Nektar::SpatialDomains::Geometry::m_state, m_verts, Nektar::SpatialDomains::Geometry::m_xmap, Nektar::SpatialDomains::Geometry::SetUpCoeffs(), and SetUpXmap().

                                             :
             Geometry1D(coordim)
         {
             m_shapeType = LibUtilities::eSegment;
             m_eid = id;
             m_globalID = id; 
             m_state = eNotFilled;
             m_curve = curve;
 
             SetUpXmap();
             SetUpCoeffs(m_xmap->GetNcoeffs());
 
             m_verts[0] = vertex[0];
             m_verts[1] = vertex[1];
         }

Nektar::SpatialDomains::SegGeom::SegGeom	(	const int	id,
		const PointGeomSharedPtr &	vert1,
		const PointGeomSharedPtr &	vert2
	)

Definition at line 106 of file SegGeom.cpp.

References Nektar::SpatialDomains::eNotFilled, Nektar::LibUtilities::eSegment, m_eid, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_shapeType, Nektar::SpatialDomains::Geometry::m_state, m_verts, Nektar::SpatialDomains::Geometry::m_xmap, Nektar::SpatialDomains::Geometry::SetUpCoeffs(), and SetUpXmap().

                                                 :
             Geometry1D(vert1->GetCoordim())
         {
             m_shapeType = LibUtilities::eSegment;
 
             m_verts[0] = vert1;
             m_verts[1] = vert2;
 
             m_state = eNotFilled;
 
             m_eid = id;
             m_globalID = id;
 
             SetUpXmap();
             SetUpCoeffs(m_xmap->GetNcoeffs());
         }

Nektar::SpatialDomains::SegGeom::SegGeom ( const SegGeom & in )

Definition at line 126 of file SegGeom.cpp.

References Nektar::SpatialDomains::Geometry::m_coordim, m_eid, m_elmtMap, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_shapeType, Nektar::SpatialDomains::Geometry::m_state, m_verts, Nektar::SpatialDomains::Geometry::m_xmap, and Nektar::SpatialDomains::Geometry::SetUpCoeffs().

         {
             // From Geometry class
             m_shapeType = in.m_shapeType;
 
             // info from EdgeComponent class
             m_eid     = in.m_eid;
             m_globalID = in.m_globalID;
 
             std::list<CompToElmt>::const_iterator def;
             for(def = in.m_elmtMap.begin(); def != in.m_elmtMap.end(); def++)
             {
                 m_elmtMap.push_back(*def);
             }
             m_xmap = in.m_xmap;
             SetUpCoeffs(m_xmap->GetNcoeffs());
 
             // info from SegGeom class
             m_coordim  = in.m_coordim;
             m_verts[0] = in.m_verts[0];
             m_verts[1] = in.m_verts[1];
 
             m_state = in.m_state;
         }

Nektar::SpatialDomains::SegGeom::~SegGeom ( )

Definition at line 230 of file SegGeom.cpp.

231 {

232 }

Member Function Documentation

SegGeomSharedPtr Nektar::SpatialDomains::SegGeom::GenerateOneSpaceDimGeom ( void )

Generate a one dimensional space segment geometry where the vert[0] has the same x value and vert[1] is set to vert[0] plus the length of the original segment.

Definition at line 175 of file SegGeom.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::SpatialDomains::eNotFilled, Nektar::SpatialDomains::Geometry::GetVid(), m_eid, m_elmtMap, Nektar::SpatialDomains::Geometry::m_geomFactors, Nektar::SpatialDomains::Geometry::m_globalID, m_verts, Nektar::SpatialDomains::Geometry::m_xmap, and v_GenGeomFactors().

         {
             SegGeomSharedPtr returnval = MemoryManager<SegGeom>::AllocateSharedPtr();
             
             // info about numbering 
             returnval->m_eid       = m_eid;
             returnval->m_globalID  = m_globalID;
             returnval->m_elmtMap   = m_elmtMap; 
             
 
             // geometric information. 
             returnval->m_coordim = 1;
             NekDouble x0 = (*m_verts[0])[0];
             PointGeomSharedPtr vert0 = MemoryManager<PointGeom>::AllocateSharedPtr(1,m_verts[0]->GetVid(),x0,0.0,0.0);
             vert0->SetGlobalID(vert0->GetVid());
             returnval->m_verts[0] = vert0;
             
             // Get information to calculate length. 
             const Array<OneD, const LibUtilities::BasisSharedPtr> base =  m_xmap->GetBase();
             LibUtilities::PointsKeyVector v;
             v.push_back(base[0]->GetPointsKey());
             v_GenGeomFactors();
 
             const Array<OneD, const NekDouble> jac = m_geomFactors->GetJac(v);
             
             NekDouble len = 0.0;
             if(jac.num_elements() == 1)
             {
                 len = jac[0]*2.0;
             }
             else
             {
                 Array<OneD, const NekDouble> w0 = base[0]->GetW();
                 len = 0.0;
 
                 for(int i = 0; i < jac.num_elements(); ++i)
                 {
                     len += jac[i]*w0[i];
                 }
             }
             // Set up second vertex. 
             PointGeomSharedPtr vert1 = MemoryManager<PointGeom>::AllocateSharedPtr(1,m_verts[1]->GetVid(),x0+len,0.0,0.0);
             vert1->SetGlobalID(vert1->GetVid());
 
             returnval->m_verts[1] = vert1;
             
             // at present just use previous m_xmap[0]; 
             returnval->m_xmap    = m_xmap;
             returnval->SetUpCoeffs(m_xmap->GetNcoeffs());
             returnval->m_state   = eNotFilled;
 
             return returnval;
         }

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

StdRegions::Orientation Nektar::SpatialDomains::SegGeom::GetEdgeOrientation	(	const SegGeom &	edge1,
		const SegGeom &	edge2
	)

static

Get the orientation of edge1.

If edge1 is connected to edge2 in the same direction as the points comprising edge1 then it is forward, otherwise it is backward.

For example, assume edge1 is comprised of points 1 and 2, and edge2 is comprised of points 2 and 3, then edge1 is forward.

If edge1 is comprised of points 2 and 1 and edge2 is comprised of points 3 and 2, then edge1 is backward.

Since both edges are passed, it does not need any information from the EdgeComponent instance.

Backward direction. Vertex 0 is connected to edge 2.

Definition at line 293 of file SegGeom.cpp.

References ASSERTL0, Nektar::StdRegions::eBackwards, Nektar::StdRegions::eForwards, Nektar::SpatialDomains::Geometry1D::GetEid(), and Nektar::SpatialDomains::Geometry1D::GetVertex().

Referenced by Nektar::NekMeshUtils::Quadrilateral::GetGeom(), Nektar::NekMeshUtils::Triangle::GetGeom(), Nektar::SpatialDomains::MeshGraph2D::ReadElements(), and Nektar::SpatialDomains::MeshGraph3D::ReadFaces().

         {
             StdRegions::Orientation returnval = StdRegions::eForwards;
 
             /// Backward direction.  Vertex 0 is connected to edge 2.
             if ((*edge1.GetVertex(0) == *edge2.GetVertex(0)) ||
                 (*edge1.GetVertex(0) == *edge2.GetVertex(1)))
             {
                 returnval = StdRegions::eBackwards;
             }
             // Not forward either, then we have a problem.
             else if ((*edge1.GetVertex(1) != *edge2.GetVertex(0)) &&
                 (*edge1.GetVertex(1) != *edge2.GetVertex(1)))
             {
                 std::ostringstream errstrm;
                 errstrm << "Connected edges do not share a vertex. Edges ";
                 errstrm << edge1.GetEid() << ", " << edge2.GetEid();
                 ASSERTL0(false, errstrm.str());
             }
 
             return returnval;
         }

StdRegions::StdExpansion1DSharedPtr Nektar::SpatialDomains::SegGeom::operator[] ( const int i ) const

void Nektar::SpatialDomains::SegGeom::SetUpXmap ( )

private

Definition at line 151 of file SegGeom.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eModified_A, m_curve, Nektar::SpatialDomains::Geometry::m_xmap, and npts.

Referenced by SegGeom(), and v_Reset().

         {
             if (m_curve)
             {
                 int npts = m_curve->m_points.size();
                 LibUtilities::PointsKey pkey(npts+1,LibUtilities::eGaussLobattoLegendre);
                 const LibUtilities::BasisKey B(LibUtilities::eModified_A, npts, pkey);
                 m_xmap = MemoryManager<StdRegions::StdSegExp>::AllocateSharedPtr(B);
             }
             else
             {
                 const LibUtilities::BasisKey B(
                     LibUtilities::eModified_A, 2,
                     LibUtilities::PointsKey(3, LibUtilities::eGaussLobattoLegendre));
                 m_xmap = MemoryManager<StdRegions::StdSegExp>::AllocateSharedPtr(B);
             }
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 249 of file SegGeom.cpp.

References m_elmtMap.

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

bool Nektar::SpatialDomains::SegGeom::v_ContainsPoint	(	const Array< OneD, const NekDouble > &	gloCoord,
		NekDouble	tol = `0.0`
	)

protectedvirtual

Determines if a point specified in global coordinates is located within this tetrahedral geometry.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 460 of file SegGeom.cpp.

References Nektar::SpatialDomains::Geometry::GetCoordim().

Referenced by v_ContainsPoint().

         {
             Array<OneD,NekDouble> locCoord(GetCoordim(),0.0);
             return v_ContainsPoint(gloCoord,locCoord,tol);
 
         }

bool Nektar::SpatialDomains::SegGeom::v_ContainsPoint	(	const Array< OneD, const NekDouble > &	gloCoord,
		Array< OneD, NekDouble > &	locCoord,
		NekDouble	tol
	)

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 468 of file SegGeom.cpp.

References v_ContainsPoint().

         {
             NekDouble resid; 
             return v_ContainsPoint(gloCoord,stdCoord,tol,resid);
         }

bool Nektar::SpatialDomains::SegGeom::v_ContainsPoint	(	const Array< OneD, const NekDouble > &	gloCoord,
		Array< OneD, NekDouble > &	locCoord,
		NekDouble	tol,
		NekDouble &	resid
	)

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 477 of file SegGeom.cpp.

References Nektar::SpatialDomains::Geometry::GetLocCoords().

         {
             resid = GetLocCoords(gloCoord, stdCoord);
             if (stdCoord[0] >= -(1+tol) && stdCoord[0] <= 1+tol)
             {
                 return true;
             }
             return false;
         }

LibUtilities::ShapeType Nektar::SpatialDomains::SegGeom::v_DetShapeType ( ) const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry1D.

Definition at line 551 of file SegGeom.cpp.

References Nektar::LibUtilities::eSegment.

         {
             return  LibUtilities::eSegment;
         }

void Nektar::SpatialDomains::SegGeom::v_FillGeom ( )

protectedvirtual

put all quadrature information into edge structure and backward transform

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 341 of file SegGeom.cpp.

References Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::SpatialDomains::ePtsFilled, ErrorUtil::ewarning, Nektar::eWrapper, Nektar::NekVector< DataType >::GetPtr(), Nektar::NekConstants::kVertexTheSameDouble, Nektar::SpatialDomains::Geometry::m_coeffs, Nektar::SpatialDomains::Geometry::m_coordim, m_curve, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_state, m_verts, Nektar::SpatialDomains::Geometry::m_xmap, NEKERROR, npts, and Nektar::LibUtilities::PointsManager().

Referenced by v_GenGeomFactors(), and v_GetLocCoords().

         {
             if(m_state != ePtsFilled)
             {
                 int i;
 
                 if (m_coordim > 0 && m_curve)
                 {
                     int npts = m_curve->m_points.size();
                     LibUtilities::PointsKey pkey(npts+1,LibUtilities::eGaussLobattoLegendre);
                     Array<OneD,NekDouble> tmp(npts);
 
                     if(m_verts[0]->dist(*(m_curve->m_points[0])) > NekConstants::kVertexTheSameDouble)
                     { 
                         std::string err = "Vertex 0 is separated from first point by more than ";
                         std::stringstream strstrm;
                         strstrm << NekConstants::kVertexTheSameDouble
                                 << " in edge " << m_globalID;
                         err += strstrm.str();
                         NEKERROR(ErrorUtil::ewarning, err.c_str());
                     }
 
                     if(m_verts[1]->dist(*(m_curve->m_points[npts-1])) > NekConstants::kVertexTheSameDouble)
                     {
                         std::string err = "Vertex 1 is separated from last point by more than ";
                         std::stringstream strstrm;
                         strstrm << NekConstants::kVertexTheSameDouble
                                 << " in edge " << m_globalID;
                         err += strstrm.str();
                         NEKERROR(ErrorUtil::ewarning, err.c_str());
                     }
 
                     LibUtilities::PointsKey fkey(npts,m_curve->m_ptype);
                     DNekMatSharedPtr I0 =
                         LibUtilities::PointsManager()[fkey]->GetI(pkey);
                     NekVector<NekDouble> out(npts+1);
 
                     for(int i = 0; i < m_coordim; ++i)
                     {
                         // Load up coordinate values into tmp
                         for(int j = 0; j < npts; ++j)
                         {
                             tmp[j] = (m_curve->m_points[j]->GetPtr())[i];
                         }
 
                         // Interpolate to GLL points
                         NekVector<NekDouble> in (npts, tmp, eWrapper);
                         out = (*I0)*in;
 
                         m_xmap->FwdTrans(out.GetPtr(), m_coeffs[i]);
                     }
                 }
 
                 for (i = 0; i < m_coordim; ++i)
                 {
                     m_coeffs[i][0] = (*m_verts[0])[i];
                     m_coeffs[i][1] = (*m_verts[1])[i];
                 }
 
                 m_state = ePtsFilled;
             }
         }

void Nektar::SpatialDomains::SegGeom::v_GenGeomFactors ( )

protectedvirtual

Set up GeoFac for this geometry using Coord quadrature distribution.

Implements Nektar::SpatialDomains::Geometry.

Definition at line 319 of file SegGeom.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().

Referenced by GenerateOneSpaceDimGeom().

         {
             if (m_geomFactorsState != ePtsFilled)
             {
                 SpatialDomains::GeomType gType = eRegular;
 
                 SegGeom::v_FillGeom();
 
                 if(m_xmap->GetBasisNumModes(0)!=2)
                 {
                     gType = eDeformed;
                 }
 
                 m_geomFactors = MemoryManager<GeomFactors>::AllocateSharedPtr(
                     gType, m_coordim, m_xmap, m_coeffs);
                 m_geomFactorsState = ePtsFilled;
             }
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 520 of file SegGeom.cpp.

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

         {
             return m_xmap->GetBasis(i);
         }

NekDouble Nektar::SpatialDomains::SegGeom::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 236 of file SegGeom.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);
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry1D.

Definition at line 509 of file SegGeom.cpp.

References m_eid.

         {
             return m_eid;
         }

NekDouble Nektar::SpatialDomains::SegGeom::v_GetLocCoords	(	const Array< OneD, const NekDouble > &	coords,
		Array< OneD, NekDouble > &	Lcoords
	)

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 419 of file SegGeom.cpp.

References ErrorUtil::efatal, Nektar::SpatialDomains::eRegular, Nektar::SpatialDomains::Geometry::GetMetricInfo(), Nektar::SpatialDomains::Geometry::m_coeffs, Nektar::SpatialDomains::Geometry::m_coordim, Nektar::SpatialDomains::Geometry::m_xmap, NEKERROR, npts, and v_FillGeom().

         {
             int i;
             NekDouble resid = 0.0;
             SegGeom::v_FillGeom();
 
             // calculate local coordinate for coord
             if(GetMetricInfo()->GetGtype() == eRegular)
             {
                 NekDouble len = 0.0;
                 NekDouble xi  = 0.0;
 
                 const int npts = m_xmap->GetTotPoints();
                 Array<OneD, NekDouble> pts(npts);
 
                 for(i = 0; i < m_coordim; ++i)
                 {
                     m_xmap->BwdTrans(m_coeffs[i], pts);
                     len  += (pts[npts-1]-pts[0])*(pts[npts-1]-pts[0]);
                     xi   += (coords[i]-pts[0])*(coords[i]-pts[0]);
                 }
 
                 len = sqrt(len);
                 xi  = sqrt(xi);
 
                 Lcoords[0] = 2*xi/len-1.0;
             }
             else
             {
                 NEKERROR(ErrorUtil::efatal,
                          "inverse mapping must be set up to use this call");
             }
             return resid;
         }

int Nektar::SpatialDomains::SegGeom::v_GetNumEdges ( ) const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 561 of file SegGeom.cpp.

References kNedges.

         {
             return kNedges;
         }

int Nektar::SpatialDomains::SegGeom::v_GetNumVerts ( ) const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 556 of file SegGeom.cpp.

References kNverts.

         {
             return kNverts;
         }

StdRegions::Orientation Nektar::SpatialDomains::SegGeom::v_GetPorient ( const int i ) const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 535 of file SegGeom.cpp.

References ASSERTL2, Nektar::StdRegions::eBwd, and Nektar::StdRegions::eFwd.

         {
             //cout << "StdRegions::PointOrientation GetPorient"<<endl;
             ASSERTL2((i >=0) && (i <= 1),"Point id must be between 0 and 1");
 
             if (i % 2 == 0)
             {
                 return StdRegions::eBwd;
             }
             else 
             {
                 return StdRegions::eFwd;
             }
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry1D.

Definition at line 515 of file SegGeom.cpp.

         {
             return 1;
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry1D.

Definition at line 497 of file SegGeom.cpp.

References kNverts, and m_verts.

         {
             PointGeomSharedPtr returnval;
 
             if (i >= 0 && i < kNverts)
             {
                 returnval = m_verts[i];
             }
 
             return returnval;
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry1D.

Definition at line 491 of file SegGeom.cpp.

References ASSERTL2, and m_verts.

         {
             ASSERTL2((i >=0) && (i <= 1),"Verted id must be between 0 and 1");
             return m_verts[i]->GetVid();
         }

StdRegions::StdExpansionSharedPtr Nektar::SpatialDomains::SegGeom::v_GetXmap ( ) const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 525 of file SegGeom.cpp.

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

         {
             return m_xmap;
         }

bool Nektar::SpatialDomains::SegGeom::v_IsElmtConnected	(	int	gvoId,
		int	locId
	)		const

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 260 of file SegGeom.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
             if(def != m_elmtMap.end())
             {
                 return(true);
             }
 
             return(false);
         }

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

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 255 of file SegGeom.cpp.

References m_elmtMap.

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

void Nektar::SpatialDomains::SegGeom::v_Reset	(	CurveMap &	curvedEdges,
		CurveMap &	curvedFaces
	)

protectedvirtual

Reset this geometry object: unset the current state and remove allocated GeomFactors.

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 404 of file SegGeom.cpp.

References Nektar::iterator, m_curve, Nektar::SpatialDomains::Geometry::m_globalID, Nektar::SpatialDomains::Geometry::m_xmap, Nektar::SpatialDomains::Geometry::SetUpCoeffs(), SetUpXmap(), and Nektar::SpatialDomains::Geometry::v_Reset().

         {
             Geometry::v_Reset(curvedEdges, curvedFaces);
             CurveMap::iterator it = curvedEdges.find(m_globalID);
 
             if (it != curvedEdges.end())
             {
                 m_curve = it->second;
             }
 
             SetUpXmap();
             SetUpCoeffs(m_xmap->GetNcoeffs());
         }

void Nektar::SpatialDomains::SegGeom::v_SetOwnData ( )

protectedvirtual

Reimplemented from Nektar::SpatialDomains::Geometry.

Definition at line 530 of file SegGeom.cpp.

References m_ownData.

         {
             m_ownData = true;
         }

Member Data Documentation

const int Nektar::SpatialDomains::SegGeom::kNedges = 1

static

Definition at line 99 of file SegGeom.h.

Referenced by v_GetNumEdges().

const int Nektar::SpatialDomains::SegGeom::kNverts = 2

static

Definition at line 98 of file SegGeom.h.

Referenced by v_GetNumVerts(), and v_GetVertex().

CurveSharedPtr Nektar::SpatialDomains::SegGeom::m_curve

private

Definition at line 176 of file SegGeom.h.

Referenced by SegGeom(), SetUpXmap(), v_FillGeom(), and v_Reset().

int Nektar::SpatialDomains::SegGeom::m_eid

protected

Definition at line 102 of file SegGeom.h.

Referenced by GenerateOneSpaceDimGeom(), SegGeom(), and v_GetEid().

std::list<CompToElmt> Nektar::SpatialDomains::SegGeom::m_elmtMap

protected

Definition at line 103 of file SegGeom.h.

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

bool Nektar::SpatialDomains::SegGeom::m_ownData

private

Boolean indicating whether object owns the data.

Definition at line 175 of file SegGeom.h.

Referenced by v_SetOwnData().

StdRegions::Orientation Nektar::SpatialDomains::SegGeom::m_porient[kNverts]

protected

Definition at line 105 of file SegGeom.h.

SpatialDomains::PointGeomSharedPtr Nektar::SpatialDomains::SegGeom::m_verts[kNverts]

protected

Definition at line 104 of file SegGeom.h.

Referenced by GenerateOneSpaceDimGeom(), SegGeom(), v_FillGeom(), v_GetVertex(), and v_GetVid().

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation