#include <NodalTriExp.h>

Inheritance diagram for Nektar::LocalRegions::NodalTriExp:

Public Member Functions
	NodalTriExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::PointsType Ntype, const SpatialDomains::TriGeomSharedPtr &geom)
	Constructor using BasisKey class for quadrature points and order definition. More...

	NodalTriExp (const NodalTriExp &T)
	Copy Constructor. More...

	~NodalTriExp ()
	Destructor. More...

void	GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3=NullNekDouble1DArray)

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

NekDouble	Integral (const Array< OneD, const NekDouble > &inarray)
	Integrate the physical point list inarray over region. More...

void	IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Inner product of inarray over region with respect to the expansion basis (this)->_Base[0] and return in outarray. More...

void	IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
	Differentiation Methods. More...

void	FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Forward transform from physical quadrature space stored in inarray and evaluate the expansion coefficients and store in (this)->_coeffs. More...

NekDouble	PhysEvaluate (const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals)

void	MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

void	LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

void	LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

void	WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

void	HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

Public Member Functions inherited from Nektar::StdRegions::StdNodalTriExp
	StdNodalTriExp ()

	StdNodalTriExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::PointsType Ntype)

	StdNodalTriExp (const StdNodalTriExp &T)

	~StdNodalTriExp ()

void	NodalToModal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	NodalToModalTranspose (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	ModalToNodal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	GetNodalPoints (Array< OneD, const NekDouble > &x, Array< OneD, const NekDouble > &y)

DNekMatSharedPtr	GenNBasisTransMatrix ()

Public Member Functions inherited from Nektar::StdRegions::StdTriExp
	StdTriExp ()

	StdTriExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)

	StdTriExp (const StdTriExp &T)

	~StdTriExp ()

Public Member Functions inherited from Nektar::StdRegions::StdExpansion2D
	StdExpansion2D ()

	StdExpansion2D (int numcoeffs, const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)

	StdExpansion2D (const StdExpansion2D &T)

virtual	~StdExpansion2D ()

void	PhysTensorDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray_d0, Array< OneD, NekDouble > &outarray_d1)
	Calculate the 2D derivative in the local tensor/collapsed coordinate at the physical points. More...

NekDouble	Integral (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &w0, const Array< OneD, const NekDouble > &w1)

void	BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true)

void	IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true)

Public Member Functions inherited from Nektar::StdRegions::StdExpansion
	StdExpansion ()
	Default Constructor. More...

	StdExpansion (const int numcoeffs, const int numbases, const LibUtilities::BasisKey &Ba=LibUtilities::NullBasisKey, const LibUtilities::BasisKey &Bb=LibUtilities::NullBasisKey, const LibUtilities::BasisKey &Bc=LibUtilities::NullBasisKey)
	Constructor. More...

	StdExpansion (const StdExpansion &T)
	Copy Constructor. More...

virtual	~StdExpansion ()
	Destructor. More...

int	GetNumBases () const
	This function returns the number of 1D bases used in the expansion. More...

const Array< OneD, const LibUtilities::BasisSharedPtr > &	GetBase () const
	This function gets the shared point to basis. More...

const LibUtilities::BasisSharedPtr &	GetBasis (int dir) const
	This function gets the shared point to basis in the dir direction. More...

int	GetNcoeffs (void) const
	This function returns the total number of coefficients used in the expansion. More...

int	GetTotPoints () const
	This function returns the total number of quadrature points used in the element. More...

LibUtilities::BasisType	GetBasisType (const int dir) const
	This function returns the type of basis used in the dir direction. More...

int	GetBasisNumModes (const int dir) const
	This function returns the number of expansion modes in the dir direction. More...

int	EvalBasisNumModesMax (void) const
	This function returns the maximum number of expansion modes over all local directions. More...

LibUtilities::PointsType	GetPointsType (const int dir) const
	This function returns the type of quadrature points used in the dir direction. More...

int	GetNumPoints (const int dir) const
	This function returns the number of quadrature points in the dir direction. More...

const Array< OneD, const NekDouble > &	GetPoints (const int dir) const
	This function returns a pointer to the array containing the quadrature points in dir direction. More...

int	GetNverts () const
	This function returns the number of vertices of the expansion domain. More...

int	GetNedges () const
	This function returns the number of edges of the expansion domain. More...

int	GetEdgeNcoeffs (const int i) const
	This function returns the number of expansion coefficients belonging to the i-th edge. More...

int	GetTotalEdgeIntNcoeffs () const

int	GetEdgeNumPoints (const int i) const
	This function returns the number of quadrature points belonging to the i-th edge. More...

int	DetCartesianDirOfEdge (const int edge)

const LibUtilities::BasisKey	DetEdgeBasisKey (const int i) const

const LibUtilities::BasisKey	DetFaceBasisKey (const int i, const int k) const

int	GetFaceNumPoints (const int i) const
	This function returns the number of quadrature points belonging to the i-th face. More...

int	GetFaceNcoeffs (const int i) const
	This function returns the number of expansion coefficients belonging to the i-th face. More...

int	GetFaceIntNcoeffs (const int i) const

int	GetTotalFaceIntNcoeffs () const

int	GetTraceNcoeffs (const int i) const
	This function returns the number of expansion coefficients belonging to the i-th edge/face. More...

LibUtilities::PointsKey	GetFacePointsKey (const int i, const int j) const

int	NumBndryCoeffs (void) const

int	NumDGBndryCoeffs (void) const

LibUtilities::BasisType	GetEdgeBasisType (const int i) const
	This function returns the type of expansion basis on the i-th edge. More...

const LibUtilities::PointsKey	GetNodalPointsKey () const
	This function returns the type of expansion Nodal point type if defined. More...

int	GetNfaces () const
	This function returns the number of faces of the expansion domain. More...

int	GetNtrace () const
	Returns the number of trace elements connected to this element. More...

LibUtilities::ShapeType	DetShapeType () const
	This function returns the shape of the expansion domain. More...

std::shared_ptr< StdExpansion >	GetStdExp (void) const

std::shared_ptr< StdExpansion >	GetLinStdExp (void) const

int	GetShapeDimension () const

bool	IsBoundaryInteriorExpansion ()

bool	IsNodalNonTensorialExp ()

void	BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	This function performs the Backward transformation from coefficient space to physical space. More...

void	FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	This function performs the Forward transformation from physical space to coefficient space. More...

void	FwdTrans_BndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

NekDouble	Integral (const Array< OneD, const NekDouble > &inarray)
	This function integrates the specified function over the domain. More...

void	FillMode (const int mode, Array< OneD, NekDouble > &outarray)
	This function fills the array outarray with the mode-th mode of the expansion. More...

void	IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	this function calculates the inner product of a given function f with the different modes of the expansion More...

void	IProductWRTBase (const Array< OneD, const NekDouble > &base, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int coll_check)

void	IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	IProductWRTDirectionalDerivBase (const Array< OneD, const NekDouble > &direction, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

int	GetElmtId ()
	Get the element id of this expansion when used in a list by returning value of m_elmt_id. More...

void	SetElmtId (const int id)
	Set the element id of this expansion when used in a list by returning value of m_elmt_id. More...

void	GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2=NullNekDouble1DArray, Array< OneD, NekDouble > &coords_3=NullNekDouble1DArray)
	this function returns the physical coordinates of the quadrature points of the expansion More...

void	GetCoord (const Array< OneD, const NekDouble > &Lcoord, Array< OneD, NekDouble > &coord)
	given the coordinates of a point of the element in the local collapsed coordinate system, this function calculates the physical coordinates of the point More...

DNekMatSharedPtr	GetStdMatrix (const StdMatrixKey &mkey)

DNekBlkMatSharedPtr	GetStdStaticCondMatrix (const StdMatrixKey &mkey)

IndexMapValuesSharedPtr	GetIndexMap (const IndexMapKey &ikey)

const Array< OneD, const NekDouble > &	GetPhysNormals (void)

void	SetPhysNormals (Array< OneD, const NekDouble > &normal)

virtual void	SetUpPhysNormals (const int edge)

void	NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)

void	NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, NekDouble > &Fy, Array< OneD, NekDouble > &outarray)

void	NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, const Array< OneD, const NekDouble > &Fz, Array< OneD, NekDouble > &outarray)

void	NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray)

DNekScalBlkMatSharedPtr	GetLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)

void	DropLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)

StdRegions::Orientation	GetForient (int face)

StdRegions::Orientation	GetEorient (int edge)

void	SetCoeffsToOrientation (Array< OneD, NekDouble > &coeffs, StdRegions::Orientation dir)

void	SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

int	CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset)

NekDouble	StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals)

int	GetCoordim ()

void	GetBoundaryMap (Array< OneD, unsigned int > &outarray)

void	GetInteriorMap (Array< OneD, unsigned int > &outarray)

int	GetVertexMap (const int localVertexId, bool useCoeffPacking=false)

void	GetEdgeInteriorMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)

void	GetFaceNumModes (const int fid, const Orientation faceOrient, int &numModes0, int &numModes1)

void	GetFaceInteriorMap (const int fid, const Orientation faceOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)

void	GetEdgeToElementMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P=-1)

void	GetFaceToElementMap (const int fid, const Orientation faceOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int nummodesA=-1, int nummodesB=-1)

void	GetEdgePhysVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Extract the physical values along edge edge from inarray into outarray following the local edge orientation and point distribution defined by defined in EdgeExp. More...

void	GetEdgePhysVals (const int edge, const std::shared_ptr< StdExpansion > &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	GetTracePhysVals (const int edge, const std::shared_ptr< StdExpansion > &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	GetVertexPhysVals (const int vertex, const Array< OneD, const NekDouble > &inarray, NekDouble &outarray)

void	GetEdgeInterpVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	GetEdgeQFactors (const int edge, Array< OneD, NekDouble > &outarray)
	Extract the metric factors to compute the contravariant fluxes along edge edge and stores them into outarray following the local edge orientation (i.e. anticlockwise convention). More...

void	GetFacePhysVals (const int face, const std::shared_ptr< StdExpansion > &FaceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient=eNoOrientation)

void	GetEdgePhysMap (const int edge, Array< OneD, int > &outarray)

void	GetFacePhysMap (const int face, Array< OneD, int > &outarray)

void	MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

DNekMatSharedPtr	CreateGeneralMatrix (const StdMatrixKey &mkey)
	this function generates the mass matrix \(\mathbf{M}[i][j] = \int \phi_i(\mathbf{x}) \phi_j(\mathbf{x}) d\mathbf{x}\) More...

void	GeneralMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey)

void	ExponentialFilter (Array< OneD, NekDouble > &array, const NekDouble alpha, const NekDouble exponent, const NekDouble cutoff)

void	LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	WeakDirectionalDerivMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	MassLevelCurvatureMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	LinearAdvectionDiffusionReactionMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true)

void	HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

DNekMatSharedPtr	GenMatrix (const StdMatrixKey &mkey)

void	PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)

void	PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	PhysDeriv_s (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_ds)

void	PhysDeriv_n (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_dn)

void	PhysDirectionalDeriv (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &direction, Array< OneD, NekDouble > &outarray)

void	StdPhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)

void	StdPhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	AddRobinMassMatrix (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)

void	AddRobinEdgeContribution (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, Array< OneD, NekDouble > &coeffs)

NekDouble	PhysEvaluate (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals)
	This function evaluates the expansion at a single (arbitrary) point of the domain. More...

NekDouble	PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals)
	This function evaluates the expansion at a single (arbitrary) point of the domain. More...

void	LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
	Convert local cartesian coordinate xi into local collapsed coordinates eta. More...

virtual int	v_GetElmtId ()
	Get the element id of this expansion when used in a list by returning value of m_elmt_id. More...

virtual const Array< OneD, const NekDouble > &	v_GetPhysNormals (void)

virtual void	v_SetPhysNormals (Array< OneD, const NekDouble > &normal)

virtual void	v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)

virtual void	v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)

virtual void	v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, const Array< OneD, const NekDouble > &Fz, Array< OneD, NekDouble > &outarray)

virtual void	v_NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray)

virtual void	v_DropLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)

virtual StdRegions::Orientation	v_GetForient (int face)

virtual StdRegions::Orientation	v_GetEorient (int edge)

NekDouble	Linf (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray)
	Function to evaluate the discrete \( L_\infty\) error \( \|\epsilon\|_\infty = \max \|u - u_{exact}\|\) where \( u_{exact}\) is given by the array sol. More...

NekDouble	L2 (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray)
	Function to evaluate the discrete \( L_2\) error, \( \| \epsilon \|_{2} = \left [ \int^1_{-1} [u - u_{exact}]^2 dx \right]^{1/2} d\xi_1 \) where \( u_{exact}\) is given by the array sol. More...

NekDouble	H1 (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray)
	Function to evaluate the discrete \( H^1\) error, \( \| \epsilon \|^1_{2} = \left [ \int^1_{-1} [u - u_{exact}]^2 + \nabla(u - u_{exact})\cdot\nabla(u - u_{exact})\cdot dx \right]^{1/2} d\xi_1 \) where \( u_{exact}\) is given by the array sol. More...

const NormalVector &	GetEdgeNormal (const int edge) const

void	ComputeEdgeNormal (const int edge)

void	NegateEdgeNormal (const int edge)

bool	EdgeNormalNegated (const int edge)

void	ComputeFaceNormal (const int face)

void	NegateFaceNormal (const int face)

bool	FaceNormalNegated (const int face)

void	ComputeVertexNormal (const int vertex)

void	NegateVertexNormal (const int vertex)

bool	VertexNormalNegated (const int vertex)

const NormalVector &	GetFaceNormal (const int face) const

const NormalVector &	GetVertexNormal (const int vertex) const

const NormalVector &	GetSurfaceNormal (const int id) const

const LibUtilities::PointsKeyVector	GetPointsKeys () const

Array< OneD, unsigned int >	GetEdgeInverseBoundaryMap (int eid)

Array< OneD, unsigned int >	GetFaceInverseBoundaryMap (int fid, StdRegions::Orientation faceOrient=eNoOrientation, int P1=-1, int P2=-1)

void	GetInverseBoundaryMaps (Array< OneD, unsigned int > &vmap, Array< OneD, Array< OneD, unsigned int > > &emap, Array< OneD, Array< OneD, unsigned int > > &fmap)

DNekMatSharedPtr	BuildInverseTransformationMatrix (const DNekScalMatSharedPtr &m_transformationmatrix)

void	PhysInterpToSimplexEquiSpaced (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int npset=-1)
	This function performs an interpolation from the physical space points provided at input into an array of equispaced points which are not the collapsed coordinate. So for a tetrahedron you will only get a tetrahedral number of values. More...

void	GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true)
	This function provides the connectivity of local simplices (triangles or tets) to connect the equispaced data points provided by PhysInterpToSimplexEquiSpaced. More...

void	EquiSpacedToCoeffs (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	This function performs a projection/interpolation from the equispaced points sometimes used in post-processing onto the coefficient space. More...

template<class T >
std::shared_ptr< T >	as ()

void	IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)

Public Member Functions inherited from Nektar::LocalRegions::Expansion2D
	Expansion2D (SpatialDomains::Geometry2DSharedPtr pGeom)

virtual	~Expansion2D ()

void	SetTraceToGeomOrientation (Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &inout)

Expansion1DSharedPtr	GetEdgeExp (int edge, bool SetUpNormal=true)

void	SetEdgeExp (const int edge, Expansion1DSharedPtr &e)

void	AddNormTraceInt (const int dir, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &edgeCoeffs, Array< OneD, NekDouble > &outarray)

void	AddNormTraceInt (const int dir, Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs)

void	AddEdgeBoundaryInt (const int edge, ExpansionSharedPtr &EdgeExp, Array< OneD, NekDouble > &edgePhys, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap)

void	AddHDGHelmholtzEdgeTerms (const NekDouble tau, const int edge, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &edgePhys, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray)

void	AddHDGHelmholtzTraceTerms (const NekDouble tau, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray)

Expansion3DSharedPtr	GetLeftAdjacentElementExp () const

Expansion3DSharedPtr	GetRightAdjacentElementExp () const

int	GetLeftAdjacentElementFace () const

int	GetRightAdjacentElementFace () const

void	SetAdjacentElementExp (int face, Expansion3DSharedPtr &f)

SpatialDomains::Geometry2DSharedPtr	GetGeom2D () const

void	ReOrientEdgePhysMap (const int nvert, const StdRegions::Orientation orient, const int nq0, Array< OneD, int > &idmap)

Public Member Functions inherited from Nektar::LocalRegions::Expansion
	Expansion (SpatialDomains::GeometrySharedPtr pGeom)

	Expansion (const Expansion &pSrc)

virtual	~Expansion ()

DNekScalMatSharedPtr	GetLocMatrix (const LocalRegions::MatrixKey &mkey)

DNekScalMatSharedPtr	GetLocMatrix (const StdRegions::MatrixType mtype, const StdRegions::ConstFactorMap &factors=StdRegions::NullConstFactorMap, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap)

SpatialDomains::GeometrySharedPtr	GetGeom () const

void	Reset ()

const SpatialDomains::GeomFactorsSharedPtr &	GetMetricInfo () const

DNekMatSharedPtr	BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)

DNekMatSharedPtr	BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd)

void	ExtractDataToCoeffs (const NekDouble data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble coeffs, std::vector< LibUtilities::BasisType > &fromType)

void	AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)

void	AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)

void	AddFaceNormBoundaryInt (const int face, const std::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)

void	DGDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &coeffs, Array< OneD, NekDouble > &outarray)

NekDouble	VectorFlux (const Array< OneD, Array< OneD, NekDouble > > &vec)

Protected Member Functions
DNekMatSharedPtr	CreateStdMatrix (const StdRegions::StdMatrixKey &mkey)

DNekScalMatSharedPtr	CreateMatrix (const MatrixKey &mkey)

DNekScalBlkMatSharedPtr	CreateStaticCondMatrix (const MatrixKey &mkey)

void	IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)

void	IProductWRTBase_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	IProductWRTDerivBase_MatOp (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	GeneralMatrixOp_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual StdRegions::StdExpansionSharedPtr	v_GetStdExp (void) const

virtual StdRegions::StdExpansionSharedPtr	v_GetLinStdExp (void) const

virtual DNekMatSharedPtr	v_GenMatrix (const StdRegions::StdMatrixKey &mkey)

Protected Member Functions inherited from Nektar::StdRegions::StdNodalTriExp
virtual const LibUtilities::PointsKey	v_GetNodalPointsKey () const

virtual bool	v_IsNodalNonTensorialExp ()

virtual void	v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Backward tranform for triangular elements. More...

virtual void	v_FillMode (const int mode, Array< OneD, NekDouble > &outarray)

virtual int	v_NumBndryCoeffs () const

virtual void	v_GetEdgeToElementMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P)

virtual int	v_GetVertexMap (int localVertexId, bool useCoeffPacking=false)

virtual void	v_GetEdgeInteriorMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)

virtual void	v_GetInteriorMap (Array< OneD, unsigned int > &outarray)

virtual void	v_GetBoundaryMap (Array< OneD, unsigned int > &outarray)

Protected Member Functions inherited from Nektar::StdRegions::StdTriExp
virtual void	v_StdPhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1)

virtual void	v_FwdTrans_BndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_IProductWRTBase_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1)

virtual void	v_IProductWRTDerivBase_MatOp (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)

virtual int	v_GetNverts () const

virtual int	v_GetNedges () const

virtual LibUtilities::ShapeType	v_DetShapeType () const

virtual int	v_NumDGBndryCoeffs () const

virtual int	v_GetEdgeNcoeffs (const int i) const

virtual int	v_GetEdgeNumPoints (const int i) const

virtual int	v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset)

virtual LibUtilities::BasisType	v_GetEdgeBasisType (const int i) const

virtual bool	v_IsBoundaryInteriorExpansion ()

virtual int	v_DetCartesianDirOfEdge (const int edge)

virtual const LibUtilities::BasisKey	v_DetEdgeBasisKey (const int edge) const

virtual void	v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey)

virtual void	v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_GeneralMatrixOp_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

virtual void	v_MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true)

Protected Member Functions inherited from Nektar::StdRegions::StdExpansion2D
virtual NekDouble	v_PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals)

virtual void	v_LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual void	v_HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual int	v_GetTraceNcoeffs (const int i) const

Protected Member Functions inherited from Nektar::StdRegions::StdExpansion
DNekMatSharedPtr	CreateStdMatrix (const StdMatrixKey &mkey)

DNekBlkMatSharedPtr	CreateStdStaticCondMatrix (const StdMatrixKey &mkey)
	Create the static condensation of a matrix when using a boundary interior decomposition. More...

IndexMapValuesSharedPtr	CreateIndexMap (const IndexMapKey &ikey)
	Create an IndexMap which contains mapping information linking any specific element shape with either its boundaries, edges, faces, verteces, etc. More...

void	BwdTrans_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	IProductWRTDirectionalDerivBase_SumFac (const Array< OneD, const NekDouble > &direction, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

void	GeneralMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	MassMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp)

void	LaplacianMatrixOp_MatFree_GenericImpl (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	LaplacianMatrixOp_MatFree (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	WeakDerivMatrixOp_MatFree (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	WeakDirectionalDerivMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	MassLevelCurvatureMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	LinearAdvectionDiffusionReactionMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true)

void	HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

void	HelmholtzMatrixOp_MatFree_GenericImpl (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)

virtual void	v_SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_SetCoeffsToOrientation (Array< OneD, NekDouble > &coeffs, StdRegions::Orientation dir)

virtual NekDouble	v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals)

Protected Member Functions inherited from Nektar::LocalRegions::Expansion2D
virtual Array< OneD, NekDouble >	v_GetMF (const int dir, const int shapedim, const StdRegions::VarCoeffMap &varcoeffs)

virtual Array< OneD, NekDouble >	v_GetMFDiv (const int dir, const StdRegions::VarCoeffMap &varcoeffs)

virtual Array< OneD, NekDouble >	v_GetMFMag (const int dir, const StdRegions::VarCoeffMap &varcoeffs)

virtual void	v_DGDeriv (const int dir, const Array< OneD, const NekDouble > &incoeffs, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &edgeCoeffs, Array< OneD, NekDouble > &out_d)

virtual void	v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)

virtual void	v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)

virtual void	v_AddRobinMassMatrix (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)

virtual void	v_AddRobinEdgeContribution (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, Array< OneD, NekDouble > &coeffs)

virtual DNekMatSharedPtr	v_BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd)

void	GetPhysEdgeVarCoeffsFromElement (const int edge, ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &varcoeff, Array< OneD, NekDouble > &outarray)

Array< OneD, NekDouble >	v_GetnEdgecdotMF (const int dir, const int edge, ExpansionSharedPtr &EdgeExp_e, const Array< OneD, const Array< OneD, NekDouble > > &normals, const StdRegions::VarCoeffMap &varcoeffs)

void	ReOrientQuadEdgePhysMap (const StdRegions::Orientation orient, const int nq0, Array< OneD, int > &idmap)

Array< OneD, unsigned int >	v_GetEdgeInverseBoundaryMap (int eid)

virtual void	v_NegateEdgeNormal (const int edge)

virtual bool	v_EdgeNormalNegated (const int edge)

virtual void	v_SetUpPhysNormals (const int edge)

const StdRegions::NormalVector &	v_GetEdgeNormal (const int edge) const

const StdRegions::NormalVector &	v_GetSurfaceNormal (const int id) const

virtual NekDouble	v_VectorFlux (const Array< OneD, Array< OneD, NekDouble > > &vec)

Protected Member Functions inherited from Nektar::LocalRegions::Expansion
void	ComputeLaplacianMetric ()

void	ComputeQuadratureMetric ()

void	ComputeGmatcdotMF (const Array< TwoD, const NekDouble > &df, const Array< OneD, const NekDouble > &direction, Array< OneD, Array< OneD, NekDouble > > &dfdir)

virtual void	v_MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_ComputeLaplacianMetric ()

Array< OneD, NekDouble >	v_GetMF (const int dir, const int shapedim, const StdRegions::VarCoeffMap &varcoeffs)

Array< OneD, NekDouble >	v_GetMFDiv (const int dir, const StdRegions::VarCoeffMap &varcoeffs)

Array< OneD, NekDouble >	v_GetMFMag (const int dir, const StdRegions::VarCoeffMap &varcoeffs)

virtual DNekMatSharedPtr	v_BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)

virtual void	v_ExtractDataToCoeffs (const NekDouble data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble coeffs, std::vector< LibUtilities::BasisType > &fromType)

virtual void	v_AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)

virtual void	v_AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)

virtual void	v_AddFaceNormBoundaryInt (const int face, const std::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)

Private Member Functions
virtual DNekMatSharedPtr	v_GenNBasisTransMatrix ()

virtual void	v_GetCoords (Array< OneD, NekDouble > &coords_0, Array< OneD, NekDouble > &coords_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coords_2=NullNekDouble1DArray)

virtual void	v_GetCoord (const Array< OneD, const NekDouble > &lcoord, Array< OneD, NekDouble > &coord)

virtual void	v_GetNodalPoints (Array< OneD, const NekDouble > &x, Array< OneD, const NekDouble > &y)

virtual NekDouble	v_Integral (const Array< OneD, const NekDouble > &inarray)
	Virtual call to integrate the physical point list inarray over region (see SegExp::Integral) More...

virtual void	v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Virtual call to TriExp::IProduct_WRT_B. More...

virtual void	v_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_StdPhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)

virtual void	v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
	Calculate the derivative of the physical points. More...

virtual void	v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Calculate the derivative of the physical points in a given direction. More...

virtual void	v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Virtual call to SegExp::FwdTrans. More...

virtual NekDouble	v_PhysEvaluate (const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals)
	Virtual call to TriExp::Evaluate. More...

virtual DNekMatSharedPtr	v_CreateStdMatrix (const StdRegions::StdMatrixKey &mkey)

virtual DNekScalMatSharedPtr	v_GetLocMatrix (const MatrixKey &mkey)

virtual DNekScalBlkMatSharedPtr	v_GetLocStaticCondMatrix (const MatrixKey &mkey)

virtual void	v_BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)

virtual void	v_IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

virtual void	v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual void	v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual void	v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual void	v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

virtual void	v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)

void	v_ComputeEdgeNormal (const int edge)

Private Attributes
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess >	m_matrixManager

LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess >	m_staticCondMatrixManager

Additional Inherited Members
Protected Attributes inherited from Nektar::StdRegions::StdNodalTriExp
LibUtilities::PointsKey	m_nodalPointsKey

Protected Attributes inherited from Nektar::StdRegions::StdExpansion
Array< OneD, LibUtilities::BasisSharedPtr >	m_base

int	m_elmt_id

int	m_ncoeffs

LibUtilities::NekManager< StdMatrixKey, DNekMat, StdMatrixKey::opLess >	m_stdMatrixManager

LibUtilities::NekManager< StdMatrixKey, DNekBlkMat, StdMatrixKey::opLess >	m_stdStaticCondMatrixManager

LibUtilities::NekManager< IndexMapKey, IndexMapValues, IndexMapKey::opLess >	m_IndexMapManager

Protected Attributes inherited from Nektar::LocalRegions::Expansion2D
std::vector< Expansion1DWeakPtr >	m_edgeExp

std::vector< bool >	m_requireNeg

std::map< int, StdRegions::NormalVector >	m_edgeNormals

std::map< int, bool >	m_negatedNormals

Expansion3DWeakPtr	m_elementLeft

Expansion3DWeakPtr	m_elementRight

int	m_elementFaceLeft

int	m_elementFaceRight

Protected Attributes inherited from Nektar::LocalRegions::Expansion
SpatialDomains::GeometrySharedPtr	m_geom

SpatialDomains::GeomFactorsSharedPtr	m_metricinfo

MetricMap	m_metrics

Detailed Description

Definition at line 52 of file NodalTriExp.h.

Constructor & Destructor Documentation

◆ NodalTriExp() [1/2]

Nektar::LocalRegions::NodalTriExp::NodalTriExp	(	const LibUtilities::BasisKey &	Ba,
		const LibUtilities::BasisKey &	Bb,
		const LibUtilities::PointsType	Ntype,
		const SpatialDomains::TriGeomSharedPtr &	geom
	)

Constructor using BasisKey class for quadrature points and order definition.

Definition at line 44 of file NodalTriExp.cpp.

                                                                             :
             StdExpansion  (LibUtilities::StdTriData::getNumberOfCoefficients(Ba.GetNumModes(),(Bb.GetNumModes())),2,Ba,Bb),
             StdExpansion2D(LibUtilities::StdTriData::getNumberOfCoefficients(Ba.GetNumModes(),(Bb.GetNumModes())),Ba,Bb),
             StdNodalTriExp(Ba,Bb,Ntype),
             Expansion     (geom),
             Expansion2D   (geom),
             m_matrixManager(
                 std::bind(&NodalTriExp::CreateMatrix, this, std::placeholders::_1),
                 std::string("NodalTriExpMatrix")),
             m_staticCondMatrixManager(
                 std::bind(&NodalTriExp::CreateStaticCondMatrix, this, std::placeholders::_1),
                 std::string("NodalTriExpStaticCondMatrix"))
         {
         }

◆ NodalTriExp() [2/2]

Nektar::LocalRegions::NodalTriExp::NodalTriExp ( const NodalTriExp & T )

Copy Constructor.

Definition at line 62 of file NodalTriExp.cpp.

                                                     :
             StdExpansion(T),
             StdExpansion2D(T),
             StdRegions::StdNodalTriExp(T),
             Expansion   (T),
             Expansion2D (T),
             m_matrixManager(T.m_matrixManager),
             m_staticCondMatrixManager(T.m_staticCondMatrixManager)
         {
         }        

◆ ~NodalTriExp()

Nektar::LocalRegions::NodalTriExp::~NodalTriExp ( )

Destructor.

Definition at line 73 of file NodalTriExp.cpp.

74 {

75 }

Member Function Documentation

◆ CreateMatrix()

DNekScalMatSharedPtr Nektar::LocalRegions::NodalTriExp::CreateMatrix ( const MatrixKey & mkey )

protected

Definition at line 432 of file NodalTriExp.cpp.

Referenced by HelmholtzMatrixOp().

         {
             DNekScalMatSharedPtr returnval;
             LibUtilities::PointsKeyVector ptsKeys = GetPointsKeys();
 
             ASSERTL2(m_metricinfo->GetGtype() != SpatialDomains::eNoGeomType,"Geometric information is not set up");
 
             StdRegions::MatrixType mtype = mkey.GetMatrixType();
 
             switch(mtype)
             {
             case StdRegions::eMass:
                 {
                     if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
                     {
                         NekDouble one = 1.0;
                         DNekMatSharedPtr mat = GenMatrix(mkey);
                         returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,mat);
                     }
                     else
                     {
                         NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
                         DNekMatSharedPtr mat = GetStdMatrix(mkey);
                         returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(jac,mat);
                     }
                 }
                 break;
             case StdRegions::eInvMass:
                 {
                     if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
                     {
                         NekDouble one = 1.0;
                         StdRegions::StdMatrixKey masskey(StdRegions::eMass,DetShapeType(),
                                                          *this);
                         DNekMatSharedPtr mat = GenMatrix(masskey);
                         mat->Invert();
 
                         returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,mat);
                     }
                     else
                     {                       
                         NekDouble fac = 1.0/(m_metricinfo->GetJac(ptsKeys))[0];
                         DNekMatSharedPtr mat = GetStdMatrix(mkey);
                         returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(fac,mat);
                     }
                 }
                 break;
             case StdRegions::eLaplacian:
                 {
                     if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
                     {
                         NekDouble one = 1.0;
                         DNekMatSharedPtr mat = GenMatrix(mkey);
                         
                         returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,mat);
                     }
                     else
                     { 
                         ASSERTL1(m_geom->GetCoordim() == 2,"Standard Region Laplacian is only set up for Quads in two-dimensional");
                         MatrixKey lap00key(StdRegions::eLaplacian00,
                                            mkey.GetShapeType(), *this);  
                         MatrixKey lap01key(StdRegions::eLaplacian01,
                                            mkey.GetShapeType(), *this);  
                         MatrixKey lap11key(StdRegions::eLaplacian11,
                                            mkey.GetShapeType(), *this);  
                         
                         DNekMat &lap00 = *GetStdMatrix(lap00key);
                         DNekMat &lap01 = *GetStdMatrix(lap01key);
                         DNekMat &lap11 = *GetStdMatrix(lap11key);
 
                         NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
                         Array<TwoD, const NekDouble> gmat =
                                                 m_metricinfo->GetGmat(ptsKeys);
 
                         int rows = lap00.GetRows();
                         int cols = lap00.GetColumns();
                         
                         DNekMatSharedPtr lap = MemoryManager<DNekMat>::AllocateSharedPtr(rows,cols);
                         
                         (*lap) = gmat[0][0] * lap00 +
                                  gmat[1][0] * (lap01 + Transpose(lap01)) +
                                  gmat[3][0] * lap11;
                         
                         returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(jac,lap);
                     }
                 }
                 break;
             case StdRegions::eHelmholtz:
                 {
                     NekDouble factor = mkey.GetConstFactor(StdRegions::eFactorLambda);
                     MatrixKey masskey(StdRegions::eMass,
                                       mkey.GetShapeType(), *this);    
                     DNekScalMat &MassMat = *(this->m_matrixManager[masskey]);
                     MatrixKey lapkey(StdRegions::eLaplacian,
                                      mkey.GetShapeType(), *this);
                     DNekScalMat &LapMat = *(this->m_matrixManager[lapkey]);
 
                     int rows = LapMat.GetRows();
                     int cols = LapMat.GetColumns();
                     
                     DNekMatSharedPtr helm = MemoryManager<DNekMat>::AllocateSharedPtr(rows,cols);
                     
                     NekDouble one = 1.0;
                     (*helm) = LapMat + factor*MassMat;
                     
                     returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,helm);            
                 }
                 break;
             default:
                 NEKERROR(ErrorUtil::efatal, "Matrix creation not defined");
                 break;
             }
 
             return returnval;
         }

◆ CreateStaticCondMatrix()

DNekScalBlkMatSharedPtr Nektar::LocalRegions::NodalTriExp::CreateStaticCondMatrix ( const MatrixKey & mkey )

protected

Definition at line 548 of file NodalTriExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL2, Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::eHelmholtz, Nektar::StdRegions::eLaplacian, Nektar::SpatialDomains::eNoGeomType, Nektar::StdRegions::StdExpansion::GetBoundaryMap(), Nektar::StdRegions::StdExpansion::GetInteriorMap(), Nektar::LocalRegions::Expansion::GetLocMatrix(), Nektar::StdRegions::StdMatrixKey::GetMatrixType(), Nektar::StdRegions::StdExpansion::GetStdStaticCondMatrix(), Nektar::LocalRegions::Expansion::m_metricinfo, Nektar::StdRegions::StdExpansion::m_ncoeffs, and Nektar::StdRegions::StdExpansion::NumBndryCoeffs().

Referenced by HelmholtzMatrixOp().

         {
             DNekScalBlkMatSharedPtr returnval;
 
             ASSERTL2(m_metricinfo->GetGtype() != SpatialDomains::eNoGeomType,"Geometric information is not set up");
 
             // set up block matrix system
             unsigned int nbdry = NumBndryCoeffs();
             unsigned int nint = (unsigned int)(m_ncoeffs - nbdry);
             unsigned int exp_size[] = {nbdry,nint};
             unsigned int nblks = 2;
             returnval = MemoryManager<DNekScalBlkMat>::AllocateSharedPtr(nblks,nblks,exp_size,exp_size); //Really need a constructor which takes Arrays
             NekDouble factor = 1.0;
 
             switch(mkey.GetMatrixType())
             {
             case StdRegions::eLaplacian: 
             case StdRegions::eHelmholtz: // special case since Helmholtz not defined in StdRegions
 
                 // use Deformed case for both regular and deformed geometries 
                 factor = 1.0;
                 goto UseLocRegionsMatrix;
                 break;
             default:
                 if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
                 {
                     factor = 1.0;
                     goto UseLocRegionsMatrix;
                 }
                 else
                 {
                     DNekScalMatSharedPtr mat = GetLocMatrix(mkey);
                     factor = mat->Scale();
                     goto UseStdRegionsMatrix;
                 }
                 break;
             UseStdRegionsMatrix:
                 {
                     NekDouble            invfactor = 1.0/factor;
                     NekDouble            one = 1.0;
                     DNekBlkMatSharedPtr  mat = GetStdStaticCondMatrix(mkey);
                     DNekScalMatSharedPtr Atmp;
                     DNekMatSharedPtr     Asubmat;
 
                     returnval->SetBlock(0,0,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(factor,Asubmat = mat->GetBlock(0,0)));
                     returnval->SetBlock(0,1,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,Asubmat = mat->GetBlock(0,1)));
                     returnval->SetBlock(1,0,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(factor,Asubmat = mat->GetBlock(1,0)));
                     returnval->SetBlock(1,1,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(invfactor,Asubmat = mat->GetBlock(1,1)));
                 }
                 break;
             UseLocRegionsMatrix:
                 {
                     int i,j;
                     NekDouble            invfactor = 1.0/factor;
                     NekDouble            one = 1.0;
                     DNekScalMat &mat = *GetLocMatrix(mkey);
                     DNekMatSharedPtr A = MemoryManager<DNekMat>::AllocateSharedPtr(nbdry,nbdry);
                     DNekMatSharedPtr B = MemoryManager<DNekMat>::AllocateSharedPtr(nbdry,nint);
                     DNekMatSharedPtr C = MemoryManager<DNekMat>::AllocateSharedPtr(nint,nbdry);
                     DNekMatSharedPtr D = MemoryManager<DNekMat>::AllocateSharedPtr(nint,nint);
                     
                     Array<OneD,unsigned int> bmap(nbdry);
                     Array<OneD,unsigned int> imap(nint);
                     GetBoundaryMap(bmap);
                     GetInteriorMap(imap);
                     
                     for(i = 0; i < nbdry; ++i)
                     {
                         for(j = 0; j < nbdry; ++j)
                         {
                             (*A)(i,j) = mat(bmap[i],bmap[j]);
                         }
                         
                         for(j = 0; j < nint; ++j)
                         {
                             (*B)(i,j) = mat(bmap[i],imap[j]);
                         }
                     }                    
                     
                     for(i = 0; i < nint; ++i)
                     {
                         for(j = 0; j < nbdry; ++j)
                         {
                             (*C)(i,j) = mat(imap[i],bmap[j]);
                         }
                         
                         for(j = 0; j < nint; ++j)
                         {
                             (*D)(i,j) = mat(imap[i],imap[j]);
                         }
                     }
                     
                     // Calculate static condensed system 
                     if(nint)
                     {
                         D->Invert();
                         (*B) = (*B)*(*D);
                         (*A) = (*A) - (*B)*(*C);
                     }
                     
                     DNekScalMatSharedPtr     Atmp;
 
                     returnval->SetBlock(0,0,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(factor,A));
                     returnval->SetBlock(0,1,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(one,B));
                     returnval->SetBlock(1,0,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(factor,C));
                     returnval->SetBlock(1,1,Atmp = MemoryManager<DNekScalMat>::AllocateSharedPtr(invfactor,D));
 
                 }
             }
             
             return returnval;
         }

◆ CreateStdMatrix()

DNekMatSharedPtr Nektar::LocalRegions::NodalTriExp::CreateStdMatrix ( const StdRegions::StdMatrixKey & mkey )

protected

Definition at line 408 of file NodalTriExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdNodalTriExp::m_nodalPointsKey.

Referenced by HelmholtzMatrixOp(), and v_CreateStdMatrix().

         {
             LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
             LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
             LibUtilities::PointsType ntype = m_nodalPointsKey.GetPointsType();
             StdRegions::StdNodalTriExpSharedPtr tmp = MemoryManager<StdNodalTriExp>::
                 AllocateSharedPtr(bkey0,bkey1,ntype);
             
             return tmp->GetStdMatrix(mkey);  
         }

◆ FwdTrans()

void Nektar::LocalRegions::NodalTriExp::FwdTrans	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

Forward transform from physical quadrature space stored in inarray and evaluate the expansion coefficients and store in (this)->_coeffs.

Forward transform from physical quadrature space stored in inarray and evaluate the expansion coefficients and store in (this)->m_coeffs.

Inputs:

inarray: array of physical quadrature points to be transformed

Outputs:

(this)->_coeffs: updated array of expansion coefficients.

Definition at line 346 of file NodalTriExp.cpp.

References Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::eCopy, Nektar::StdRegions::eInvMass, Nektar::eWrapper, Nektar::LibUtilities::PointsKey::GetPointsType(), IProductWRTBase(), m_matrixManager, Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdNodalTriExp::m_nodalPointsKey, Nektar::StdRegions::NullConstFactorMap, and Nektar::StdRegions::NullVarCoeffMap.

Referenced by IProductWRTDerivBase(), and v_FwdTrans().

         {
             IProductWRTBase(inarray,outarray); 
 
             // get Mass matrix inverse
             MatrixKey  masskey(StdRegions::eInvMass, DetShapeType(),*this,StdRegions::NullConstFactorMap,StdRegions::NullVarCoeffMap,m_nodalPointsKey.GetPointsType());
             DNekScalMatSharedPtr  matsys = m_matrixManager[masskey];
             
             // copy inarray in case inarray == outarray
             NekVector<NekDouble> in(m_ncoeffs,outarray,eCopy);
             NekVector<NekDouble> out(m_ncoeffs,outarray,eWrapper);
             
             out = (*matsys)*in;
         }

◆ GeneralMatrixOp_MatOp()

void Nektar::LocalRegions::NodalTriExp::GeneralMatrixOp_MatOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

protected

Definition at line 362 of file NodalTriExp.cpp.

References Blas::Dgemv(), Nektar::LocalRegions::Expansion::GetLocMatrix(), Nektar::StdRegions::StdExpansion::m_ncoeffs, and Vmath::Vcopy().

Referenced by HelmholtzMatrixOp().

         {
             DNekScalMatSharedPtr   mat = GetLocMatrix(mkey);
 
             if(inarray.get() == outarray.get())
             {
                 Array<OneD,NekDouble> tmp(m_ncoeffs);
                 Vmath::Vcopy(m_ncoeffs,inarray.get(),1,tmp.get(),1);
                 
                 Blas::Dgemv('N',m_ncoeffs,m_ncoeffs,mat->Scale(),(mat->GetOwnedMatrix())->GetPtr().get(),
                             m_ncoeffs, tmp.get(), 1, 0.0, outarray.get(), 1);
             }
             else
             {                
                 Blas::Dgemv('N',m_ncoeffs,m_ncoeffs,mat->Scale(),(mat->GetOwnedMatrix())->GetPtr().get(),
                             m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
             }
         }

◆ GetCoord()

void Nektar::LocalRegions::NodalTriExp::GetCoord	(	const Array< OneD, const NekDouble > &	Lcoords,
		Array< OneD, NekDouble > &	coords
	)

Definition at line 391 of file NodalTriExp.cpp.

References ASSERTL1, and Nektar::LocalRegions::Expansion::m_geom.

Referenced by v_GetCoord().

         {
             int  i;
             
             ASSERTL1(Lcoords[0] >= -1.0 && Lcoords[1] <= 1.0 && 
                      Lcoords[1] >= -1.0 && Lcoords[1]  <=1.0,
                      "Local coordinates are not in region [-1,1]");
             
             m_geom->FillGeom();
             
             for(i = 0; i < m_geom->GetCoordim(); ++i)
             {
                 coords[i] = m_geom->GetCoord(i,Lcoords);
             }
         }

◆ GetCoords()

void Nektar::LocalRegions::NodalTriExp::GetCoords	(	Array< OneD, NekDouble > &	coords_1,
		Array< OneD, NekDouble > &	coords_2,
		Array< OneD, NekDouble > &	coords_3 = `NullNekDouble1DArray`
	)

Definition at line 383 of file NodalTriExp.cpp.

References Nektar::LocalRegions::Expansion::v_GetCoords().

Referenced by v_GetCoords().

         {
             Expansion::v_GetCoords(coords_0, coords_1, coords_2);
         }

◆ HelmholtzMatrixOp()

void Nektar::LocalRegions::NodalTriExp::HelmholtzMatrixOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inline

Definition at line 150 of file NodalTriExp.h.

References CreateMatrix(), CreateStaticCondMatrix(), CreateStdMatrix(), GeneralMatrixOp_MatOp(), IProductWRTBase_MatOp(), IProductWRTBase_SumFac(), IProductWRTDerivBase_MatOp(), IProductWRTDerivBase_SumFac(), v_GenMatrix(), v_GetLinStdExp(), and v_GetStdExp().

Referenced by v_HelmholtzMatrixOp().

             {
                 StdExpansion::HelmholtzMatrixOp_MatFree_GenericImpl(inarray,outarray,mkey);
             }  

◆ Integral()

NekDouble Nektar::LocalRegions::NodalTriExp::Integral ( const Array< OneD, const NekDouble > & inarray )

Integrate the physical point list inarray over region.

Integrate the physical point list inarray over region and return the value.

Inputs:

inarray: definition of function to be returned at quadrature point of expansion.

Outputs:

returns \(\int^1_{-1}\int^1_{-1} u(\xi_1, \xi_2) J[i,j] d \xi_1 d \xi_2 \) where \(inarray[i,j] = u(\xi_{1i},\xi_{2j}) \) and \( J[i,j] \) is the Jacobian evaluated at the quadrature point.

Definition at line 98 of file NodalTriExp.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Smul(), and Vmath::Vmul().

Referenced by v_Integral().

         {
             int    nquad0 = m_base[0]->GetNumPoints();
             int    nquad1 = m_base[1]->GetNumPoints();
             Array<OneD, const NekDouble> jac = m_metricinfo->GetJac(GetPointsKeys());
             NekDouble ival;
             Array<OneD,NekDouble> tmp(nquad0*nquad1);
             
             // multiply inarray with Jacobian
             if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
             {
                 Vmath::Vmul(nquad0*nquad1, jac, 1, inarray, 1,tmp, 1);
             }
             else
             {
                 Vmath::Smul(nquad0*nquad1, jac[0], inarray, 1, tmp, 1);
             }
             
             // call StdQuadExp version;
             ival = StdNodalTriExp::v_Integral(tmp);            
             return ival; 
         }

◆ IProductWRTBase()

void Nektar::LocalRegions::NodalTriExp::IProductWRTBase	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inline

Inner product of inarray over region with respect to the expansion basis (this)->_Base[0] and return in outarray.

Definition at line 84 of file NodalTriExp.h.

References IProductWRTBase_SumFac().

Referenced by FwdTrans(), and v_IProductWRTBase().

             {
                 NodalTriExp::IProductWRTBase_SumFac(inarray,outarray);
             }

◆ IProductWRTBase_MatOp()

void Nektar::LocalRegions::NodalTriExp::IProductWRTBase_MatOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

protected

Definition at line 148 of file NodalTriExp.cpp.

References Nektar::StdRegions::StdExpansion::DetShapeType(), Blas::Dgemv(), Nektar::StdRegions::eIProductWRTBase, Nektar::StdRegions::StdExpansion::GetTotPoints(), m_matrixManager, and Nektar::StdRegions::StdExpansion::m_ncoeffs.

Referenced by HelmholtzMatrixOp().

         {
             int nq = GetTotPoints();
             MatrixKey      iprodmatkey(StdRegions::eIProductWRTBase,DetShapeType(),*this);
             DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
             
             Blas::Dgemv('N',m_ncoeffs,nq,iprodmat->Scale(),(iprodmat->GetOwnedMatrix())->GetPtr().get(),
                         m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
         }

◆ IProductWRTBase_SumFac()

void Nektar::LocalRegions::NodalTriExp::IProductWRTBase_SumFac	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		bool	multiplybyweights = `true`
	)

protected

Definition at line 122 of file NodalTriExp.cpp.

References Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), and Nektar::StdRegions::StdNodalTriExp::NodalToModalTranspose().

Referenced by HelmholtzMatrixOp(), IProductWRTBase(), and v_IProductWRTBase_SumFac().

         { 
             int    nquad0 = m_base[0]->GetNumPoints();
             int    nquad1 = m_base[1]->GetNumPoints();
             int    order1 = m_base[1]->GetNumModes();
 
             if(multiplybyweights)
             {
                 Array<OneD,NekDouble> tmp(nquad0*nquad1+nquad0*order1);
                 Array<OneD,NekDouble> wsp(tmp+nquad0*nquad1);
                 
                 MultiplyByQuadratureMetric(inarray,tmp);
                 StdTriExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),m_base[1]->GetBdata(),tmp,outarray,wsp);
                 NodalToModalTranspose(outarray,outarray);  
             }
             else
             {
                 Array<OneD,NekDouble> wsp(nquad0*order1);
                 
                 StdTriExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),m_base[1]->GetBdata(),inarray,outarray,wsp);
                 NodalToModalTranspose(outarray,outarray);  
             }
         }

◆ IProductWRTDerivBase()

void Nektar::LocalRegions::NodalTriExp::IProductWRTDerivBase	(	const int	dir,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inline

Definition at line 90 of file NodalTriExp.h.

References FwdTrans(), IProductWRTDerivBase_SumFac(), LOCAL_REGIONS_EXPORT, Nektar::NullNekDouble1DArray, PhysDeriv(), and PhysEvaluate().

Referenced by v_IProductWRTDerivBase().

             {
                 NodalTriExp::IProductWRTDerivBase_SumFac(dir,inarray,outarray);
             }

◆ IProductWRTDerivBase_MatOp()

void Nektar::LocalRegions::NodalTriExp::IProductWRTDerivBase_MatOp	(	const int	dir,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

protected

Definition at line 229 of file NodalTriExp.cpp.

References ASSERTL1, Nektar::StdRegions::StdExpansion::DetShapeType(), Blas::Dgemv(), Nektar::StdRegions::eIProductWRTDerivBase0, Nektar::StdRegions::eIProductWRTDerivBase1, Nektar::StdRegions::eIProductWRTDerivBase2, Nektar::StdRegions::StdExpansion::GetTotPoints(), m_matrixManager, and Nektar::StdRegions::StdExpansion::m_ncoeffs.

Referenced by HelmholtzMatrixOp().

         { 
             int nq = GetTotPoints();            
             StdRegions::MatrixType mtype = StdRegions::eIProductWRTDerivBase0;
 
             switch(dir)
             {
             case 0:
                 {
                     mtype = StdRegions::eIProductWRTDerivBase0;
                 }
                 break;
             case 1:
                 {
                     mtype = StdRegions::eIProductWRTDerivBase1;
                 }
                 break;
             case 2:
                 {
                     mtype = StdRegions::eIProductWRTDerivBase2;
                 }
                 break;
             default:
                 {
                     ASSERTL1(false,"input dir is out of range");
                 }
                 break;
             }  
 
             MatrixKey      iprodmatkey(mtype,DetShapeType(),*this);
             DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
             
             Blas::Dgemv('N',m_ncoeffs,nq,iprodmat->Scale(),(iprodmat->GetOwnedMatrix())->GetPtr().get(),
                         m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
 
         }

◆ IProductWRTDerivBase_SumFac()

void Nektar::LocalRegions::NodalTriExp::IProductWRTDerivBase_SumFac	(	const int	dir,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

protected

Definition at line 159 of file NodalTriExp.cpp.

References ASSERTL1, Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion2D::IProductWRTBase_SumFacKernel(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_geom, Nektar::LocalRegions::Expansion::m_metricinfo, Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), Nektar::StdRegions::StdNodalTriExp::NodalToModalTranspose(), Vmath::Smul(), Vmath::Vadd(), and Vmath::Vmul().

Referenced by HelmholtzMatrixOp(), IProductWRTDerivBase(), and v_IProductWRTDerivBase_SumFac().

         {   
             ASSERTL1((dir==0)||(dir==1)||(dir==2),"Invalid direction.");
             ASSERTL1((dir==2)?(m_geom->GetCoordim()==3):true,"Invalid direction.");
 
             int    i;
             int    nquad0 = m_base[0]->GetNumPoints();
             int    nquad1 = m_base[1]->GetNumPoints();
             int    nqtot  = nquad0*nquad1; 
             int    wspsize = max(nqtot,m_ncoeffs);
 
             const Array<TwoD, const NekDouble>& df =
                                 m_metricinfo->GetDerivFactors(GetPointsKeys());
             
             Array<OneD, NekDouble> tmp0 (6*wspsize);
             Array<OneD, NekDouble> tmp1 (tmp0 +   wspsize);
             Array<OneD, NekDouble> tmp2 (tmp0 + 2*wspsize);
             Array<OneD, NekDouble> tmp3 (tmp0 + 3*wspsize);
             Array<OneD, NekDouble> gfac0(tmp0 + 4*wspsize);
             Array<OneD, NekDouble> gfac1(tmp0 + 5*wspsize);
 
             const Array<OneD, const NekDouble>& z0 = m_base[0]->GetZ();
             const Array<OneD, const NekDouble>& z1 = m_base[1]->GetZ();
 
             // set up geometric factor: 2/(1-z1)
             for(i = 0; i < nquad1; ++i)
             {
                 gfac0[i] = 2.0/(1-z1[i]);
             }
             for(i = 0; i < nquad0; ++i)
             {
                 gfac1[i] = 0.5*(1+z0[i]);
             }
 
             for(i = 0; i < nquad1; ++i)  
             {
                 Vmath::Smul(nquad0,gfac0[i],&inarray[0]+i*nquad0,1,&tmp0[0]+i*nquad0,1);
             }
             
             for(i = 0; i < nquad1; ++i) 
             {
                 Vmath::Vmul(nquad0,&gfac1[0],1,&tmp0[0]+i*nquad0,1,&tmp1[0]+i*nquad0,1);
             }
                                
             if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
             {
                 Vmath::Vmul(nqtot,&df[2*dir][0],  1,&tmp0[0],   1,&tmp0[0],1);
                 Vmath::Vmul(nqtot,&df[2*dir+1][0],1,&tmp1[0],   1,&tmp1[0],1);
                 Vmath::Vmul(nqtot,&df[2*dir+1][0],1,&inarray[0],1,&tmp2[0],1);
             }
             else
             {
                 Vmath::Smul(nqtot, df[2*dir][0],   tmp0,    1, tmp0, 1);
                 Vmath::Smul(nqtot, df[2*dir+1][0], tmp1,    1, tmp1, 1);
                 Vmath::Smul(nqtot, df[2*dir+1][0], inarray, 1, tmp2, 1);
             }
             Vmath::Vadd(nqtot, tmp0, 1, tmp1, 1, tmp1, 1); 
 
             MultiplyByQuadratureMetric(tmp1,tmp1);
             MultiplyByQuadratureMetric(tmp2,tmp2);
 
             IProductWRTBase_SumFacKernel(m_base[0]->GetDbdata(),m_base[1]->GetBdata() ,tmp1,tmp3    ,tmp0);
             IProductWRTBase_SumFacKernel(m_base[0]->GetBdata() ,m_base[1]->GetDbdata(),tmp2,outarray,tmp0);
             Vmath::Vadd(m_ncoeffs, tmp3, 1, outarray, 1, outarray, 1);
 
             NodalToModalTranspose(outarray,outarray);              
         }

◆ LaplacianMatrixOp() [1/2]

void Nektar::LocalRegions::NodalTriExp::LaplacianMatrixOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inline

Definition at line 127 of file NodalTriExp.h.

Referenced by v_LaplacianMatrixOp().

             {           
                 StdExpansion::LaplacianMatrixOp_MatFree_GenericImpl(inarray,outarray,mkey);
             }

◆ LaplacianMatrixOp() [2/2]

void Nektar::LocalRegions::NodalTriExp::LaplacianMatrixOp	(	const int	k1,
		const int	k2,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inline

Definition at line 134 of file NodalTriExp.h.

             {           
                 StdExpansion::LaplacianMatrixOp_MatFree(k1,k2,inarray,outarray,mkey);
             }

◆ MassMatrixOp()

void Nektar::LocalRegions::NodalTriExp::MassMatrixOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inline

Definition at line 120 of file NodalTriExp.h.

Referenced by v_MassMatrixOp().

             {              
                 StdExpansion::MassMatrixOp_MatFree(inarray,outarray,mkey);
             }

◆ PhysDeriv()

void Nektar::LocalRegions::NodalTriExp::PhysDeriv	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	out_d0,
		Array< OneD, NekDouble > &	out_d1,
		Array< OneD, NekDouble > &	out_d2 = `NullNekDouble1DArray`
	)

Differentiation Methods.

Calculate the deritive of the physical points

Definition at line 275 of file NodalTriExp.cpp.

References Blas::Daxpy(), Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Smul(), Vmath::Vmul(), and Vmath::Vvtvp().

Referenced by IProductWRTDerivBase(), and v_PhysDeriv().

         {
             int    nquad0 = m_base[0]->GetNumPoints();
             int    nquad1 = m_base[1]->GetNumPoints();
             int     nqtot = nquad0*nquad1;
             const Array<TwoD, const NekDouble>& df
                             = m_metricinfo->GetDerivFactors(GetPointsKeys());
 
             Array<OneD,NekDouble> diff0(2*nqtot);
             Array<OneD,NekDouble> diff1(diff0+nqtot);
 
             StdNodalTriExp::v_PhysDeriv(inarray, diff0, diff1);
 
             if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
             {
                 if(out_d0.num_elements())
                 {
                     Vmath::Vmul  (nqtot,df[0],1,diff0,1, out_d0, 1);
                     Vmath::Vvtvp (nqtot,df[1],1,diff1,1, out_d0, 1, out_d0,1);
                 }
 
                 if(out_d1.num_elements())
                 {
                     Vmath::Vmul  (nqtot,df[2],1,diff0,1, out_d1, 1);
                     Vmath::Vvtvp (nqtot,df[3],1,diff1,1, out_d1, 1, out_d1,1);
                 }
 
                 if(out_d2.num_elements())
                 {
                     Vmath::Vmul  (nqtot,df[4],1,diff0,1, out_d2, 1);
                     Vmath::Vvtvp (nqtot,df[5],1,diff1,1, out_d2, 1, out_d2,1);
                 }
             }
             else // regular geometry
             {
                 if(out_d0.num_elements())
                 {
                     Vmath::Smul (nqtot, df[0][0], diff0, 1, out_d0, 1);
                     Blas::Daxpy (nqtot, df[1][0], diff1, 1, out_d0, 1);
                 }
 
                 if(out_d1.num_elements())
                 {
                     Vmath::Smul (nqtot, df[2][0], diff0, 1, out_d1, 1);
                     Blas::Daxpy (nqtot, df[3][0], diff1, 1, out_d1, 1);
                 }
 
                 if(out_d2.num_elements())
                 {
                     Vmath::Smul (nqtot, df[4][0], diff0, 1, out_d2, 1);
                     Blas::Daxpy (nqtot, df[5][0], diff1, 1, out_d2, 1);
                 }
             }
         }    

◆ PhysEvaluate()

NekDouble Nektar::LocalRegions::NodalTriExp::PhysEvaluate	(	const Array< OneD, const NekDouble > &	coord,
		const Array< OneD, const NekDouble > &	physvals
	)

Definition at line 419 of file NodalTriExp.cpp.

References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.

Referenced by IProductWRTDerivBase(), and v_PhysEvaluate().

         {
             Array<OneD,NekDouble> Lcoord = Array<OneD,NekDouble>(2);
             
             ASSERTL0(m_geom,"m_geom not defined");
             m_geom->GetLocCoords(coord,Lcoord);
             
             return StdNodalTriExp::v_PhysEvaluate(Lcoord, physvals);
         }

◆ v_BwdTrans_SumFac()

virtual void Nektar::LocalRegions::NodalTriExp::v_BwdTrans_SumFac	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 322 of file NodalTriExp.h.

             {
                 StdNodalTriExp::v_BwdTrans_SumFac(inarray,outarray);
             }    

◆ v_ComputeEdgeNormal()

void Nektar::LocalRegions::NodalTriExp::v_ComputeEdgeNormal ( const int edge )

privatevirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 703 of file NodalTriExp.cpp.

References ASSERTL0, Nektar::SpatialDomains::eMovingRegular, Nektar::SpatialDomains::eRegular, Vmath::Fill(), Nektar::StdRegions::StdExpansion::GetCoordim(), Nektar::LocalRegions::Expansion::GetGeom(), Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::LibUtilities::Interp1D(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion2D::m_edgeNormals, Vmath::Reverse(), Vmath::Sdiv(), Vmath::Smul(), Vmath::Vmul(), Vmath::Vsqrt(), Vmath::Vvtvp(), and Vmath::Zero().

Referenced by v_HelmholtzMatrixOp().

         {
             int i;
             const SpatialDomains::GeomFactorsSharedPtr & geomFactors = GetGeom()->GetMetricInfo();
             const SpatialDomains::GeomType type = geomFactors->GetGtype();
 
             LibUtilities::PointsKeyVector ptsKeys = GetPointsKeys();
             const Array<TwoD, const NekDouble> & df = geomFactors->GetDerivFactors(ptsKeys);
             const Array<OneD, const NekDouble> & jac  = geomFactors->GetJac(ptsKeys);
             int nqe = m_base[0]->GetNumPoints();
             int dim = GetCoordim();
 
             m_edgeNormals[edge] = Array<OneD, Array<OneD, NekDouble> >(dim);
             Array<OneD, Array<OneD, NekDouble> > &normal = m_edgeNormals[edge];
             for (i = 0; i < dim; ++i)
             {
                 normal[i] = Array<OneD, NekDouble>(nqe);
             }
 
             // Regular geometry case
             if((type == SpatialDomains::eRegular)||(type == SpatialDomains::eMovingRegular))
             {
                 NekDouble fac;
                 // Set up normals
                 switch(edge)
                 {
                 case 0:
                     for(i = 0; i < GetCoordim(); ++i)
                     {
                         Vmath::Fill(nqe,-df[2*i+1][0],normal[i],1);
                     }
                     break;
                 case 1:
                     for(i = 0; i < GetCoordim(); ++i)
                     {
                         Vmath::Fill(nqe,df[2*i+1][0] + df[2*i][0],normal[i],1);
                     }
                         break;
                 case 2:
                     for(i = 0; i < GetCoordim(); ++i)
                     {
                         Vmath::Fill(nqe,-df[2*i][0],normal[i],1);
                     }
                     break;
                 default:
                     ASSERTL0(false,"Edge is out of range (edge < 3)");
                 }
 
                 // normalise
                 fac = 0.0;
                 for(i =0 ; i < GetCoordim(); ++i)
                 {
                     fac += normal[i][0]*normal[i][0];
                 }
                 fac = 1.0/sqrt(fac);
                 for (i = 0; i < GetCoordim(); ++i)
                 {
                     Vmath::Smul(nqe,fac,normal[i],1,normal[i],1);
                 }
             }
             else   // Set up deformed normals
             {
                 int j;
 
                 int nquad0 = ptsKeys[0].GetNumPoints();
                 int nquad1 = ptsKeys[1].GetNumPoints();
 
                 LibUtilities::PointsKey from_key;
 
                 Array<OneD,NekDouble> normals(GetCoordim()*max(nquad0,nquad1),0.0);
                 Array<OneD,NekDouble> edgejac(GetCoordim()*max(nquad0,nquad1),0.0);
 
                 // Extract Jacobian along edges and recover local
                 // derivates (dx/dr) for polynomial interpolation by
                 // multiplying m_gmat by jacobian
                 switch(edge)
                 {
                 case 0:
                     for(j = 0; j < nquad0; ++j)
                     {
                         edgejac[j] = jac[j];
                         for(i = 0; i < GetCoordim(); ++i)
                         {
                             normals[i*nquad0+j] = -df[2*i+1][j]*edgejac[j];
                         }
                     }
                     from_key = ptsKeys[0];
                     break;
                 case 1:
                     for(j = 0; j < nquad1; ++j)
                     {
                         edgejac[j] = jac[nquad0*j+nquad0-1];
                         for(i = 0; i < GetCoordim(); ++i)
                         {
                             normals[i*nquad1+j] = (df[2*i][nquad0*j + nquad0-1] +  df[2*i+1][nquad0*j + nquad0-1])*edgejac[j];
                         }
                     }
                     from_key = ptsKeys[1];
                     break;
                 case 2:
                     for(j = 0; j < nquad1; ++j)
                     {
                         edgejac[j] = jac[nquad0*j];
                         for(i = 0; i < GetCoordim(); ++i)
                         {
                             normals[i*nquad1+j] = -df[2*i][nquad0*j]*edgejac[j];
                         }
                     }
                     from_key = ptsKeys[1];
                     break;
                 default:
                     ASSERTL0(false,"edge is out of range (edge < 3)");
 
                 }
 
                 int nq  = from_key.GetNumPoints();
                 Array<OneD,NekDouble> work(nqe,0.0);
 
                 // interpolate Jacobian and invert
                 LibUtilities::Interp1D(from_key,jac,m_base[0]->GetPointsKey(),work);
                 Vmath::Sdiv(nq,1.0,&work[0],1,&work[0],1);
 
                 // interpolate
                 for(i = 0; i < GetCoordim(); ++i)
                 {
                     LibUtilities::Interp1D(from_key,&normals[i*nq],m_base[0]->GetPointsKey(),&normal[i][0]);
                     Vmath::Vmul(nqe,work,1,normal[i],1,normal[i],1);
                 }
 
                 //normalise normal vectors
                 Vmath::Zero(nqe,work,1);
                 for(i = 0; i < GetCoordim(); ++i)
                 {
                     Vmath::Vvtvp(nqe,normal[i],1, normal[i],1,work,1,work,1);
                 }
 
                 Vmath::Vsqrt(nqe,work,1,work,1);
                 Vmath::Sdiv(nqe,1.0,work,1,work,1);
 
                 for(i = 0; i < GetCoordim(); ++i)
                 {
                     Vmath::Vmul(nqe,normal[i],1,work,1,normal[i],1);
                 }
 
                 // Reverse direction so that points are in
                 // anticlockwise direction if edge >=2
                 if(edge >= 2)
                 {
                     for(i = 0; i < GetCoordim(); ++i)
                     {
                         Vmath::Reverse(nqe,normal[i],1, normal[i],1);
                     }
                 }
             }
         }

◆ v_CreateStdMatrix()

virtual DNekMatSharedPtr Nektar::LocalRegions::NodalTriExp::v_CreateStdMatrix ( const StdRegions::StdMatrixKey & mkey )

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 301 of file NodalTriExp.h.

References CreateStdMatrix().

             {
                 return CreateStdMatrix(mkey);
             }

◆ v_FwdTrans()

virtual void Nektar::LocalRegions::NodalTriExp::v_FwdTrans	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inlineprivatevirtual

Virtual call to SegExp::FwdTrans.

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 286 of file NodalTriExp.h.

References FwdTrans().

             {
                 FwdTrans(inarray,outarray);
             }

◆ v_GenMatrix()

DNekMatSharedPtr Nektar::LocalRegions::NodalTriExp::v_GenMatrix ( const StdRegions::StdMatrixKey & mkey )

protectedvirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 682 of file NodalTriExp.cpp.

References Nektar::StdRegions::eHybridDGHelmBndLam, Nektar::StdRegions::eHybridDGHelmholtz, Nektar::StdRegions::eHybridDGLamToQ0, Nektar::StdRegions::eHybridDGLamToQ1, Nektar::StdRegions::eHybridDGLamToQ2, Nektar::StdRegions::eHybridDGLamToU, Nektar::StdRegions::StdMatrixKey::GetMatrixType(), and Nektar::LocalRegions::Expansion2D::v_GenMatrix().

Referenced by HelmholtzMatrixOp().

         {
             DNekMatSharedPtr returnval;
 
             switch(mkey.GetMatrixType())
             {
             case StdRegions::eHybridDGHelmholtz:
             case StdRegions::eHybridDGLamToU:
             case StdRegions::eHybridDGLamToQ0:
             case StdRegions::eHybridDGLamToQ1:
             case StdRegions::eHybridDGLamToQ2:
             case StdRegions::eHybridDGHelmBndLam:
                 returnval = Expansion2D::v_GenMatrix(mkey);
                 break;
             default:
                 returnval = StdNodalTriExp::v_GenMatrix(mkey);
                 break;
             }
             return returnval;
         }

◆ v_GenNBasisTransMatrix()

virtual DNekMatSharedPtr Nektar::LocalRegions::NodalTriExp::v_GenNBasisTransMatrix ( )

inlineprivatevirtual

Definition at line 197 of file NodalTriExp.h.

             {
                 return StdNodalTriExp::GenNBasisTransMatrix();
             }

◆ v_GetCoord()

virtual void Nektar::LocalRegions::NodalTriExp::v_GetCoord	(	const Array< OneD, const NekDouble > &	lcoord,
		Array< OneD, NekDouble > &	coord
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 209 of file NodalTriExp.h.

References GetCoord().

             {
                 GetCoord(lcoord, coord);
             }

◆ v_GetCoords()

virtual void Nektar::LocalRegions::NodalTriExp::v_GetCoords	(	Array< OneD, NekDouble > &	coords_0,
		Array< OneD, NekDouble > &	coords_1 = `NullNekDouble1DArray`,
		Array< OneD, NekDouble > &	coords_2 = `NullNekDouble1DArray`
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 202 of file NodalTriExp.h.

References GetCoords().

             {
                 GetCoords(coords_0, coords_1, coords_2);
             }

◆ v_GetLinStdExp()

StdRegions::StdExpansionSharedPtr Nektar::LocalRegions::NodalTriExp::v_GetLinStdExp ( void ) const

protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 671 of file NodalTriExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::StdRegions::StdExpansion::GetBasisType(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdNodalTriExp::m_nodalPointsKey.

Referenced by HelmholtzMatrixOp().

         {
             LibUtilities::BasisKey bkey0(m_base[0]->GetBasisType(),
                            2, m_base[0]->GetPointsKey());
             LibUtilities::BasisKey bkey1(m_base[1]->GetBasisType(),
                            2, m_base[1]->GetPointsKey());
 
             return MemoryManager<StdRegions::StdNodalTriExp>
                 ::AllocateSharedPtr( bkey0, bkey1, m_nodalPointsKey.GetPointsType());
         }

◆ v_GetLocMatrix()

virtual DNekScalMatSharedPtr Nektar::LocalRegions::NodalTriExp::v_GetLocMatrix ( const MatrixKey & mkey )

inlineprivatevirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 306 of file NodalTriExp.h.

             {
                 return m_matrixManager[mkey];
             }

◆ v_GetLocStaticCondMatrix()

virtual DNekScalBlkMatSharedPtr Nektar::LocalRegions::NodalTriExp::v_GetLocStaticCondMatrix ( const MatrixKey & mkey )

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 317 of file NodalTriExp.h.

             {
                 return m_staticCondMatrixManager[mkey];
             }

◆ v_GetNodalPoints()

virtual void Nektar::LocalRegions::NodalTriExp::v_GetNodalPoints	(	Array< OneD, const NekDouble > &	x,
		Array< OneD, const NekDouble > &	y
	)

inlineprivatevirtual

Definition at line 215 of file NodalTriExp.h.

             {
                 return StdNodalTriExp::GetNodalPoints(x,y);
             }

◆ v_GetStdExp()

StdRegions::StdExpansionSharedPtr Nektar::LocalRegions::NodalTriExp::v_GetStdExp ( void ) const

protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 662 of file NodalTriExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdNodalTriExp::m_nodalPointsKey.

Referenced by HelmholtzMatrixOp().

         {
             
             return MemoryManager<StdRegions::StdNodalTriExp>
                     ::AllocateSharedPtr(m_base[0]->GetBasisKey(),
                                         m_base[1]->GetBasisKey(),
                                         m_nodalPointsKey.GetPointsType());
         }

◆ v_HelmholtzMatrixOp()

virtual void Nektar::LocalRegions::NodalTriExp::v_HelmholtzMatrixOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 373 of file NodalTriExp.h.

References HelmholtzMatrixOp(), and v_ComputeEdgeNormal().

             {
                 HelmholtzMatrixOp(inarray,outarray,mkey);
             }  

◆ v_Integral()

virtual NekDouble Nektar::LocalRegions::NodalTriExp::v_Integral ( const Array< OneD, const NekDouble > & inarray )

inlineprivatevirtual

Virtual call to integrate the physical point list inarray over region (see SegExp::Integral)

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 223 of file NodalTriExp.h.

References Integral().

             {
                 return Integral(inarray);
             }

◆ v_IProductWRTBase()

virtual void Nektar::LocalRegions::NodalTriExp::v_IProductWRTBase	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inlineprivatevirtual

Virtual call to TriExp::IProduct_WRT_B.

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 229 of file NodalTriExp.h.

References IProductWRTBase().

             {
                 IProductWRTBase(inarray,outarray);
             }

◆ v_IProductWRTBase_SumFac()

virtual void Nektar::LocalRegions::NodalTriExp::v_IProductWRTBase_SumFac	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		bool	multiplybyweights = `true`
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 328 of file NodalTriExp.h.

References IProductWRTBase_SumFac().

             {
                 boost::ignore_unused(multiplybyweights);
                 IProductWRTBase_SumFac(inarray,outarray);
             }            

◆ v_IProductWRTDerivBase()

virtual void Nektar::LocalRegions::NodalTriExp::v_IProductWRTDerivBase	(	const int	dir,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 235 of file NodalTriExp.h.

References IProductWRTDerivBase().

             {
                 IProductWRTDerivBase(dir,inarray,outarray);
             }

◆ v_IProductWRTDerivBase_SumFac()

virtual void Nektar::LocalRegions::NodalTriExp::v_IProductWRTDerivBase_SumFac	(	const int	dir,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 336 of file NodalTriExp.h.

References IProductWRTDerivBase_SumFac().

             {
                 IProductWRTDerivBase_SumFac(dir,inarray,outarray);
             }    

◆ v_LaplacianMatrixOp() [1/2]

virtual void Nektar::LocalRegions::NodalTriExp::v_LaplacianMatrixOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 350 of file NodalTriExp.h.

References LaplacianMatrixOp().

             {
                 LaplacianMatrixOp(inarray,outarray,mkey);
             }

◆ v_LaplacianMatrixOp() [2/2]

virtual void Nektar::LocalRegions::NodalTriExp::v_LaplacianMatrixOp	(	const int	k1,
		const int	k2,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 357 of file NodalTriExp.h.

References LaplacianMatrixOp().

             {
                 LaplacianMatrixOp(k1,k2,inarray,outarray,mkey);
             }

◆ v_MassMatrixOp()

virtual void Nektar::LocalRegions::NodalTriExp::v_MassMatrixOp	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 343 of file NodalTriExp.h.

References MassMatrixOp().

             {
                 MassMatrixOp(inarray,outarray,mkey);
             }  

◆ v_PhysDeriv() [1/2]

virtual void Nektar::LocalRegions::NodalTriExp::v_PhysDeriv	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	out_d0,
		Array< OneD, NekDouble > &	out_d1,
		Array< OneD, NekDouble > &	out_d2 = `NullNekDouble1DArray`
	)

inlineprivatevirtual

Calculate the derivative of the physical points.

\( \frac{\partial u}{\partial x_1} = \left . \frac{2.0}{1-\eta_2} \frac{\partial u}{\partial d\eta_1} \right |_{\eta_2}\)

\( \frac{\partial u}{\partial x_2} = \left . \frac{1+\eta_1}{1-\eta_2} \frac{\partial u}{\partial d\eta_1} \right |_{\eta_2} + \left . \frac{\partial u}{\partial d\eta_2} \right |_{\eta_1} \)

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 251 of file NodalTriExp.h.

References PhysDeriv().

             {
                 boost::ignore_unused(out_d2);
                 PhysDeriv(inarray, out_d0, out_d1);
             }

◆ v_PhysDeriv() [2/2]

virtual void Nektar::LocalRegions::NodalTriExp::v_PhysDeriv	(	const int	dir,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	out_d0
	)

inlineprivatevirtual

Calculate the derivative of the physical points in a given direction.

See also: StdRegions::StdExpansion::PhysDeriv

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 260 of file NodalTriExp.h.

References ASSERTL1, and PhysDeriv().

             {
                 Array<OneD,NekDouble> tmp;
                 switch(dir)
                 {
                 case 0:
                     {
                         PhysDeriv(inarray, outarray, tmp);   
                     }
                     break;
                 case 1:
                     {
                         PhysDeriv(inarray, tmp, outarray);   
                     }
                     break;
                 default:
                     {
                         ASSERTL1(dir >= 0 &&dir < 2,"input dir is out of range");
                     }
                     break;
                 }             
             }

◆ v_PhysEvaluate()

virtual NekDouble Nektar::LocalRegions::NodalTriExp::v_PhysEvaluate	(	const Array< OneD, const NekDouble > &	coord,
		const Array< OneD, const NekDouble > &	physvals
	)

inlineprivatevirtual

Virtual call to TriExp::Evaluate.

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 293 of file NodalTriExp.h.

References PhysEvaluate().

             {
                 return PhysEvaluate(coord, physvals);
             }

◆ v_StdPhysDeriv()

virtual void Nektar::LocalRegions::NodalTriExp::v_StdPhysDeriv	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	out_d0,
		Array< OneD, NekDouble > &	out_d1,
		Array< OneD, NekDouble > &	out_d2 = `NullNekDouble1DArray`
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 242 of file NodalTriExp.h.

             {
                 StdTriExp::v_PhysDeriv(inarray, out_d0, out_d1, out_d2);
             }

◆ v_WeakDerivMatrixOp()

virtual void Nektar::LocalRegions::NodalTriExp::v_WeakDerivMatrixOp	(	const int	i,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inlineprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 365 of file NodalTriExp.h.

References WeakDerivMatrixOp().

             {
                 WeakDerivMatrixOp(i,inarray,outarray,mkey);
             }

◆ WeakDerivMatrixOp()

void Nektar::LocalRegions::NodalTriExp::WeakDerivMatrixOp	(	const int	i,
		const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const StdRegions::StdMatrixKey &	mkey
	)

inline

Definition at line 142 of file NodalTriExp.h.

Referenced by v_WeakDerivMatrixOp().

             {
                 StdExpansion::WeakDerivMatrixOp_MatFree(i,inarray,outarray,mkey);
             }

Member Data Documentation

◆ m_matrixManager

LibUtilities::NekManager<MatrixKey, DNekScalMat, MatrixKey::opLess> Nektar::LocalRegions::NodalTriExp::m_matrixManager

private

Definition at line 188 of file NodalTriExp.h.

Referenced by CreateMatrix(), FwdTrans(), IProductWRTBase_MatOp(), and IProductWRTDerivBase_MatOp().

◆ m_staticCondMatrixManager

LibUtilities::NekManager<MatrixKey, DNekScalBlkMat, MatrixKey::opLess> Nektar::LocalRegions::NodalTriExp::m_staticCondMatrixManager

private

Definition at line 189 of file NodalTriExp.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ NodalTriExp() [1/2]

◆ NodalTriExp() [2/2]

◆ ~NodalTriExp()

Member Function Documentation

◆ CreateMatrix()

◆ CreateStaticCondMatrix()

◆ CreateStdMatrix()

◆ FwdTrans()

◆ GeneralMatrixOp_MatOp()

◆ GetCoord()

◆ GetCoords()

◆ HelmholtzMatrixOp()

◆ Integral()

◆ IProductWRTBase()

◆ IProductWRTBase_MatOp()

◆ IProductWRTBase_SumFac()

◆ IProductWRTDerivBase()

◆ IProductWRTDerivBase_MatOp()

◆ IProductWRTDerivBase_SumFac()

◆ LaplacianMatrixOp() [1/2]

◆ LaplacianMatrixOp() [2/2]

◆ MassMatrixOp()

◆ PhysDeriv()

◆ PhysEvaluate()

◆ v_BwdTrans_SumFac()

◆ v_ComputeEdgeNormal()

◆ v_CreateStdMatrix()

◆ v_FwdTrans()

◆ v_GenMatrix()

◆ v_GenNBasisTransMatrix()

◆ v_GetCoord()

◆ v_GetCoords()

◆ v_GetLinStdExp()

◆ v_GetLocMatrix()

◆ v_GetLocStaticCondMatrix()

◆ v_GetNodalPoints()

◆ v_GetStdExp()

◆ v_HelmholtzMatrixOp()

◆ v_Integral()

◆ v_IProductWRTBase()

◆ v_IProductWRTBase_SumFac()

◆ v_IProductWRTDerivBase()

◆ v_IProductWRTDerivBase_SumFac()

◆ v_LaplacianMatrixOp() [1/2]

◆ v_LaplacianMatrixOp() [2/2]

◆ v_MassMatrixOp()

◆ v_PhysDeriv() [1/2]

◆ v_PhysDeriv() [2/2]

◆ v_PhysEvaluate()

◆ v_StdPhysDeriv()

◆ v_WeakDerivMatrixOp()

◆ WeakDerivMatrixOp()

Member Data Documentation

◆ m_matrixManager

◆ m_staticCondMatrixManager