Nektar::StdRegions::StdNodalTriExp Class Reference

#include <StdNodalTriExp.h>

Inheritance diagram for Nektar::StdRegions::StdNodalTriExp:

Public Member Functions
	StdNodalTriExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::PointsType Ntype)
	StdNodalTriExp ()=default
	StdNodalTriExp (const StdNodalTriExp &T)=default
	~StdNodalTriExp () override=default
void	v_NodalToModal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
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 (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)
	StdTriExp ()=default
	StdTriExp (const StdTriExp &T)=default
	~StdTriExp () override=default
Public Member Functions inherited from Nektar::StdRegions::StdExpansion2D
	StdExpansion2D (int numcoeffs, const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)
	StdExpansion2D ()=default
	StdExpansion2D (const StdExpansion2D &T)=default
	~StdExpansion2D () override=default
NekDouble	BaryTensorDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs)
void	IProductWRTBaseKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const Array< OneD, NekDouble > &jac, const bool Deformed, bool CollDir0=false, bool CollDir1=false)
Public Member Functions inherited from Nektar::StdRegions::StdExpansion
	StdExpansion ()
	Default Constructor.
	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.
	StdExpansion (const StdExpansion &T)
	Copy Constructor.
virtual	~StdExpansion ()
	Destructor.
int	GetNumBases () const
	This function returns the number of 1D bases used in the expansion.
const Array< OneD, const LibUtilities::BasisSharedPtr > &	GetBase () const
	This function gets the shared point to basis.
const LibUtilities::BasisSharedPtr &	GetBasis (int dir) const
	This function gets the shared point to basis in the dir direction.
int	GetNcoeffs (void) const
	This function returns the total number of coefficients used in the expansion.
int	GetTotPoints () const
	This function returns the total number of quadrature points used in the element.
LibUtilities::BasisType	GetBasisType (const int dir) const
	This function returns the type of basis used in the dir direction.
int	GetBasisNumModes (const int dir) const
	This function returns the number of expansion modes in the dir direction.
int	EvalBasisNumModesMax (void) const
	This function returns the maximum number of expansion modes over all local directions.
LibUtilities::PointsType	GetPointsType (const int dir) const
	This function returns the type of quadrature points used in the dir direction.
int	GetNumPoints (const int dir) const
	This function returns the number of quadrature points in the dir direction.
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.
int	GetNverts () const
	This function returns the number of vertices of the expansion domain.
int	GetTraceNcoeffs (const int i) const
	This function returns the number of expansion coefficients belonging to the i-th trace.
int	GetTraceIntNcoeffs (const int i) const
int	GetTraceNumPoints (const int i) const
	This function returns the number of quadrature points belonging to the i-th trace.
const LibUtilities::BasisKey	GetTraceBasisKey (const int i, int k=-1, bool UseGLL=false) const
	This function returns the basis key belonging to the i-th trace.
LibUtilities::PointsKey	GetTracePointsKey (const int i, int k=-1) const
	This function returns the basis key belonging to the i-th trace.
int	NumBndryCoeffs (void) const
int	NumDGBndryCoeffs (void) const
const LibUtilities::PointsKey	GetNodalPointsKey () const
	This function returns the type of expansion Nodal point type if defined.
int	GetNtraces () const
	Returns the number of trace elements connected to this element.
LibUtilities::ShapeType	DetShapeType () const
	This function returns the shape of the expansion domain.
int	GetShapeDimension () const
bool	IsBoundaryInteriorExpansion () const
bool	IsNodalNonTensorialExp ()
void	NodalToModal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
void	BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	This function performs the Backward transformation from coefficient space to physical space.
void	FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
void	FwdTransBndConstrained (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.
void	FillMode (const int mode, Array< OneD, NekDouble > &outarray)
	This function fills the array outarray with the mode-th mode of the expansion.
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
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.
void	SetElmtId (const int id)
	Set the element id of this expansion when used in a list by returning value of m_elmt_id.
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
Array< OneD, Array< OneD, NekDouble > >	GetCoords ()
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
DNekMatSharedPtr	GetStdMatrix (const StdMatrixKey &mkey)
DNekBlkMatSharedPtr	GetStdStaticCondMatrix (const StdMatrixKey &mkey)
Array< OneD, const NekDouble >	GetStdFac (const StdFacKey &mkey)
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)
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	GetTraceToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards, int P=-1, int Q=-1)
void	GetTraceCoeffMap (const unsigned int traceid, Array< OneD, unsigned int > &maparray)
void	GetElmtTraceToTraceMap (const unsigned int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards, int P=-1, int Q=-1)
void	GetTraceInteriorToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation traceOrient=eForwards)
void	GetTraceNumModes (const int tid, int &numModes0, int &numModes1, const Orientation traceOrient=eDir1FwdDir1_Dir2FwdDir2)
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}\)
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	LinearAdvectionMatrixOp (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	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)
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.
NekDouble	PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs)
	This function evaluates the first derivative of the expansion at a single (arbitrary) point of the domain.
NekDouble	PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs, std::array< NekDouble, 6 > &secondOrderDerivs)
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.
NekDouble	PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode)
	This function evaluates the basis function mode mode at a point coords of the domain.
void	ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1, bool Forwards=true)
void	LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
	Convert local cartesian coordinate xi into local collapsed coordinates eta.
void	LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi)
	Convert local collapsed coordinates eta into local cartesian coordinate xi.
void	PhysInterp (std::shared_ptr< StdExpansion > fromExp, const Array< OneD, const NekDouble > &fromData, Array< OneD, NekDouble > &toData, bool Transpose=false)
	interpolate from one set of quadrature points available from FromExp to the set of quadrature points in the current expansion. If the points are the same this routine will just copy the data
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 DNekScalBlkMatSharedPtr	v_GetLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)
virtual void	v_DropLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)
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.
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.
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.
const LibUtilities::PointsKeyVector	GetPointsKeys () const
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.
void	PhysInterpToGLL (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int npset=-1)
void	PhysInterpToPoints (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int npset, MatrixType distrib)
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.
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.
void	EquiSpacedToPhys (const int nequi, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
template<class T >
std::shared_ptr< T >	as ()
void	GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat)

Protected Member Functions
const LibUtilities::PointsKey	v_GetNodalPointsKey () const override
bool	v_IsNodalNonTensorialExp () override
void	v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
void	v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
	Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray.
void	v_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
void	v_FillMode (const int mode, Array< OneD, NekDouble > &outarray) override
LibUtilities::ShapeType	v_DetShapeType () const final
int	v_NumBndryCoeffs () const override
const LibUtilities::BasisKey	v_GetTraceBasisKey (const int i, const int j, bool UseGLL=false) const override
int	v_GetVertexMap (int localVertexId, bool useCoeffPacking=false) override
void	v_GetTraceToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient=eForwards, int P=-1, int Q=-1) override
void	v_GetTraceInteriorToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation edgeOrient=eForwards) override
void	v_GetInteriorMap (Array< OneD, unsigned int > &outarray) override
void	v_GetBoundaryMap (Array< OneD, unsigned int > &outarray) override
void	v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
DNekMatSharedPtr	v_GenMatrix (const StdMatrixKey &mkey) override
DNekMatSharedPtr	v_CreateStdMatrix (const StdMatrixKey &mkey) override
Protected Member Functions inherited from Nektar::StdRegions::StdTriExp
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) override
	Calculate the derivative of the physical points.
void	v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
	Backward tranform for triangular elements.
void	v_FwdTransBndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
void	v_IProductWRTBaseKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const Array< OneD, NekDouble > &jac, const bool Deformed, const bool CollDir0=false, const bool CollDir1=false) override
	Inner product of inarray over region with respect to the expansion basis (this)->m_base[0] and return in outarray.
void	v_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
void	v_LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta) override
void	v_LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi) override
void	v_FillMode (const int mode, Array< OneD, NekDouble > &outarray) override
NekDouble	v_PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode) final
NekDouble	v_PhysEvalFirstDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) override
int	v_GetNverts () const final
int	v_GetNtraces () const final
LibUtilities::ShapeType	v_DetShapeType () const override
int	v_NumBndryCoeffs () const override
int	v_NumDGBndryCoeffs () const override
int	v_GetTraceNcoeffs (const int i) const override
int	v_GetTraceIntNcoeffs (const int i) const override
int	v_GetTraceNumPoints (const int i) const override
int	v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset) override
void	v_GetCoords (Array< OneD, NekDouble > &coords_x, Array< OneD, NekDouble > &coords_y, Array< OneD, NekDouble > &coords_z) override
bool	v_IsBoundaryInteriorExpansion () const override
const LibUtilities::BasisKey	v_GetTraceBasisKey (const int i, const int j, bool UseGLL=false) const override
int	v_GetVertexMap (int localVertexId, bool useCoeffPacking=false) override
void	v_GetInteriorMap (Array< OneD, unsigned int > &outarray) override
void	v_GetBoundaryMap (Array< OneD, unsigned int > &outarray) override
void	v_GetTraceCoeffMap (const unsigned int traceid, Array< OneD, unsigned int > &maparray) override
void	v_GetTraceInteriorToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation edgeOrient=eForwards) override
DNekMatSharedPtr	v_GenMatrix (const StdMatrixKey &mkey) override
DNekMatSharedPtr	v_CreateStdMatrix (const StdMatrixKey &mkey) override
void	v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey) override
void	v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
void	v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override
void	v_GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true) override
Protected Member Functions inherited from Nektar::StdRegions::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.
void	v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
	Calculate the derivative of the physical points in a given direction.
NekDouble	v_StdPhysEvaluate (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals) override
	This function evaluates the expansion at a single (arbitrary) point of the domain.
NekDouble	v_PhysEvaluateInterp (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals) override
void	v_MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
void	v_LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
void	v_HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
void	v_GetElmtTraceToTraceMap (const unsigned int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient, int P, int Q) override
	Determine the mapping to re-orientate the coefficients along the element trace (assumed to align with the standard element) into the orientation of the local trace given by edgeOrient.
void	v_GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat) override
void	v_PhysInterp (std::shared_ptr< StdExpansion > fromExp, const Array< OneD, const NekDouble > &fromData, Array< OneD, NekDouble > &toData, bool Transpose) override
void	v_ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1, bool Forwards) override
int	v_GetShapeDimension () const final
bool	v_IsCollocatedBasis () const final
virtual void	v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2, Array< OneD, NekDouble > &out_d3)
	Calculate the derivative of the physical points.
virtual void	v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0)
	Calculate the derivative of the physical points in a given direction.
Protected Member Functions inherited from Nektar::StdRegions::StdExpansion
DNekMatSharedPtr	CreateStdMatrix (const StdMatrixKey &mkey)
std::shared_ptr< Array< OneD, const NekDouble > >	CreateStdFac (const StdFacKey &mkey)
DNekBlkMatSharedPtr	CreateStdStaticCondMatrix (const StdMatrixKey &mkey)
	Create the static condensation of a matrix when using a boundary interior decomposition.
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	LinearAdvectionMatrixOp_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)
template<int DIR, bool DERIV = false, bool DERIV2 = false>
NekDouble	BaryEvaluate (const NekDouble &coord, const NekDouble *physvals, NekDouble &deriv, NekDouble &deriv2)
	This function performs the barycentric interpolation of the polynomial stored in coord at a point physvals using barycentric interpolation weights in direction.
template<int DIR>
NekDouble	BaryEvaluateBasis (const NekDouble &coord, const int &mode)
template<int DIR, bool DERIV = false, bool DERIV2 = false>
NekDouble	BaryEvaluate (const NekDouble &coord, const NekDouble *physvals)
	Helper function to pass an unused value by reference into BaryEvaluate.
template<int DIR, bool DERIV = false, bool DERIV2 = false>
NekDouble	BaryEvaluate (const NekDouble &coord, const NekDouble *physvals, NekDouble &deriv)
virtual LibUtilities::PointsKey	v_GetTracePointsKey (const int i, const int j) const
virtual void	v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
	Transform a given function from physical quadrature space to coefficient space.
virtual void	v_IProductWRTDirectionalDerivBase (const Array< OneD, const NekDouble > &direction, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void	v_PhysDirectionalDeriv (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &direction, Array< OneD, NekDouble > &outarray)
	Physical derivative along a direction vector.
virtual NekDouble	v_PhysEvaluate (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals)
virtual NekDouble	v_PhysEvalFirstSecondDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs, std::array< NekDouble, 6 > &secondOrderDerivs)
virtual void	v_GetCoord (const Array< OneD, const NekDouble > &Lcoord, Array< OneD, NekDouble > &coord)
virtual int	v_GetCoordim () const
virtual void	v_GetTraceNumModes (const int fid, int &numModes0, int &numModes1, Orientation traceOrient=eDir1FwdDir1_Dir2FwdDir2)
virtual void	v_GetVertexPhysVals (const int vertex, const Array< OneD, const NekDouble > &inarray, NekDouble &outarray)
virtual void	v_MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
virtual void	v_ExponentialFilter (Array< OneD, NekDouble > &array, const NekDouble alpha, const NekDouble exponent, const NekDouble cutoff)
virtual void	v_WeakDirectionalDerivMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
virtual void	v_MassLevelCurvatureMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
virtual void	v_LinearAdvectionMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
virtual void	v_LinearAdvectionDiffusionReactionMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true)
virtual void	v_LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp)
virtual DNekMatSharedPtr	v_BuildInverseTransformationMatrix (const DNekScalMatSharedPtr &m_transformationmatrix)

Protected Attributes
LibUtilities::PointsKey	m_nodalPointsKey
Protected Attributes inherited from Nektar::StdRegions::StdExpansion
Array< OneD, LibUtilities::BasisSharedPtr >	m_base
int	m_elmt_id
int	m_ncoeffs
std::vector< Array< OneD, const NekDouble > >	m_weights
LibUtilities::NekManager< StdMatrixKey, DNekMat, StdMatrixKey::opLess >	m_stdMatrixManager
LibUtilities::NekManager< StdMatrixKey, DNekBlkMat, StdMatrixKey::opLess >	m_stdStaticCondMatrixManager
LibUtilities::NekManager< StdFacKey, Array< OneD, const NekDouble > >	m_stdFacManager

Detailed Description

Definition at line 43 of file StdNodalTriExp.h.

Constructor & Destructor Documentation

StdNodalTriExp() [1/3]

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

Definition at line 42 of file StdNodalTriExp.cpp.

    : StdExpansion(LibUtilities::StdTriData::getNumberOfCoefficients(
                       Ba.GetNumModes(), Bb.GetNumModes()),
                   2, Ba, Bb),
      StdExpansion2D(LibUtilities::StdTriData::getNumberOfCoefficients(
                         Ba.GetNumModes(), Bb.GetNumModes()),
                     Ba, Bb),
      StdTriExp(Ba, Bb), m_nodalPointsKey(Ba.GetNumModes(), Ntype)
{
    ASSERTL0(m_base[0]->GetNumModes() == m_base[1]->GetNumModes(),
             "Nodal basis initiated with different orders in the a "
             "and b directions");
 
    // cache integration weights for future use
    m_weights.push_back(m_base[0]->GetW());
 
    StdFacKey w1key(eWeights1, Bb);
    // get weights[1] from manager where points are rescaled
    m_weights.push_back(GetStdFac(w1key));
}

References ASSERTL0, Nektar::StdRegions::eWeights1, Nektar::StdRegions::StdExpansion::GetStdFac(), Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdExpansion::m_weights.

StdNodalTriExp() [2/3]

Nektar::StdRegions::StdNodalTriExp::StdNodalTriExp ( )

default

StdNodalTriExp() [3/3]

Nektar::StdRegions::StdNodalTriExp::StdNodalTriExp ( const StdNodalTriExp &  T )

default

~StdNodalTriExp()

Nektar::StdRegions::StdNodalTriExp::~StdNodalTriExp ( )

overridedefault

Member Function Documentation

GenNBasisTransMatrix()

DNekMatSharedPtr Nektar::StdRegions::StdNodalTriExp::GenNBasisTransMatrix

(

)

Definition at line 122 of file StdNodalTriExp.cpp.

{
    int i, j;
    Array<OneD, const NekDouble> r, s;
    Array<OneD, NekDouble> c(2);
    DNekMatSharedPtr Mat;
 
    Mat = MemoryManager<DNekMat>::AllocateSharedPtr(m_ncoeffs, m_ncoeffs);
    GetNodalPoints(r, s);
 
    // Store the values of m_phys in a temporary array
    int nqtot = GetTotPoints();
    Array<OneD, NekDouble> phys(nqtot);
 
    for (i = 0; i < m_ncoeffs; ++i)
    {
        // fill physical space with mode i
        StdTriExp::v_FillMode(i, phys);
 
        // interpolate mode i to the Nodal points 'j' and
        // store in outarray
        for (j = 0; j < m_ncoeffs; ++j)
        {
            c[0]         = r[j];
            c[1]         = s[j];
            (*Mat)(j, i) = StdExpansion2D::v_PhysEvaluate(c, phys);
        }
    }
    return Mat;
}

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), GetNodalPoints(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdTriExp::v_FillMode(), and Nektar::StdRegions::StdExpansion::v_PhysEvaluate().

Referenced by v_GenMatrix().

GetNodalPoints()

void Nektar::StdRegions::StdNodalTriExp::GetNodalPoints	(	Array< OneD, const NekDouble > &	x,
		Array< OneD, const NekDouble > &	y
	)

Definition at line 116 of file StdNodalTriExp.cpp.

{
    LibUtilities::PointsManager()[m_nodalPointsKey]->GetPoints(x, y);
}

References m_nodalPointsKey, and Nektar::LibUtilities::PointsManager().

Referenced by GenNBasisTransMatrix().

ModalToNodal()

void Nektar::StdRegions::StdNodalTriExp::ModalToNodal	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

Definition at line 102 of file StdNodalTriExp.cpp.

{
    StdMatrixKey Nkey(eNBasisTrans, DetShapeType(), *this, NullConstFactorMap,
                      NullVarCoeffMap, NullVarFactorsMap,
                      m_nodalPointsKey.GetPointsType());
    DNekMatSharedPtr vdm = GetStdMatrix(Nkey);
 
    // Multiply out matrix
    NekVector<NekDouble> modal(m_ncoeffs, inarray, eWrapper);
    NekVector<NekDouble> nodal(m_ncoeffs, outarray, eWrapper);
    nodal = (*vdm) * modal;
}

References Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::StdRegions::eNBasisTrans, Nektar::eWrapper, Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::StdRegions::StdExpansion::GetStdMatrix(), Nektar::StdRegions::StdExpansion::m_ncoeffs, m_nodalPointsKey, Nektar::StdRegions::NullConstFactorMap, Nektar::StdRegions::NullVarCoeffMap, and Nektar::StdRegions::NullVarFactorsMap.

Referenced by Nektar::GlobalMapping::UpdateGeometry(), and Nektar::LocalRegions::NodalTriExp::v_ExtractDataToCoeffs().

NodalToModalTranspose()

void Nektar::StdRegions::StdNodalTriExp::NodalToModalTranspose	(	const Array< OneD, const NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray
	)

Definition at line 88 of file StdNodalTriExp.cpp.

{
    StdMatrixKey Nkey(eInvNBasisTrans, DetShapeType(), *this,
                      NullConstFactorMap, NullVarCoeffMap, NullVarFactorsMap,
                      m_nodalPointsKey.GetPointsType());
    DNekMatSharedPtr inv_vdm = GetStdMatrix(Nkey);
 
    NekVector<NekDouble> nodal(m_ncoeffs, inarray, eCopy);
    NekVector<NekDouble> modal(m_ncoeffs, outarray, eWrapper);
    modal = Transpose(*inv_vdm) * nodal;
}

References Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::eCopy, Nektar::StdRegions::eInvNBasisTrans, Nektar::eWrapper, Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::StdRegions::StdExpansion::GetStdMatrix(), Nektar::StdRegions::StdExpansion::m_ncoeffs, m_nodalPointsKey, Nektar::StdRegions::NullConstFactorMap, Nektar::StdRegions::NullVarCoeffMap, Nektar::StdRegions::NullVarFactorsMap, and Nektar::Transpose().

Referenced by Nektar::LocalRegions::NodalTriExp::v_IProductWRTBase(), v_IProductWRTBase(), Nektar::LocalRegions::NodalTriExp::v_IProductWRTDerivBase(), and v_IProductWRTDerivBase().

v_BwdTrans()

void Nektar::StdRegions::StdNodalTriExp::v_BwdTrans ( const Array< OneD, const NekDouble > &  inarray, Array< OneD, NekDouble > &  outarray  )

overrideprotectedvirtual

Implements Nektar::StdRegions::StdExpansion.

Definition at line 157 of file StdNodalTriExp.cpp.

{
    Array<OneD, NekDouble> tmp(m_ncoeffs);
    v_NodalToModal(inarray, tmp);
    StdTriExp::v_BwdTrans(tmp, outarray);
}

References Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdTriExp::v_BwdTrans(), and v_NodalToModal().

Referenced by v_FillMode().

v_CreateStdMatrix()

DNekMatSharedPtr Nektar::StdRegions::StdNodalTriExp::v_CreateStdMatrix ( const StdMatrixKey &  mkey )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 440 of file StdNodalTriExp.cpp.

{
    return StdNodalTriExp::v_GenMatrix(mkey);
}

References v_GenMatrix().

v_DetShapeType()

LibUtilities::ShapeType Nektar::StdRegions::StdNodalTriExp::v_DetShapeType ( ) const

finalprotectedvirtual

Implements Nektar::StdRegions::StdExpansion.

Definition at line 204 of file StdNodalTriExp.cpp.

{
    return LibUtilities::eNodalTri;
}

References Nektar::LibUtilities::eNodalTri.

v_FillMode()

void Nektar::StdRegions::StdNodalTriExp::v_FillMode ( const int  mode, Array< OneD, NekDouble > &  outarray  )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 189 of file StdNodalTriExp.cpp.

{
    ASSERTL2(mode >= m_ncoeffs,
             "calling argument mode is larger than total expansion order");
 
    Vmath::Zero(m_ncoeffs, outarray, 1);
    outarray[mode] = 1.0;
    v_BwdTrans(outarray, outarray);
}

References ASSERTL2, Nektar::StdRegions::StdExpansion::m_ncoeffs, v_BwdTrans(), and Vmath::Zero().

v_GenMatrix()

DNekMatSharedPtr Nektar::StdRegions::StdNodalTriExp::v_GenMatrix ( const StdMatrixKey &  mkey )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 423 of file StdNodalTriExp.cpp.

{
    DNekMatSharedPtr Mat;
 
    switch (mkey.GetMatrixType())
    {
        case eNBasisTrans:
            Mat = GenNBasisTransMatrix();
            break;
        default:
            Mat = StdExpansion::CreateGeneralMatrix(mkey);
            break;
    }
 
    return Mat;
}

References Nektar::StdRegions::StdExpansion::CreateGeneralMatrix(), Nektar::StdRegions::eNBasisTrans, GenNBasisTransMatrix(), and Nektar::StdRegions::StdMatrixKey::GetMatrixType().

Referenced by v_CreateStdMatrix().

v_GetBoundaryMap()

void Nektar::StdRegions::StdNodalTriExp::v_GetBoundaryMap ( Array< OneD, unsigned int > &  outarray )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 405 of file StdNodalTriExp.cpp.

{
    unsigned int i;
    if (outarray.size() != NumBndryCoeffs())
    {
        outarray = Array<OneD, unsigned int>(NumBndryCoeffs());
    }
 
    for (i = 0; i < NumBndryCoeffs(); i++)
    {
        outarray[i] = i;
    }
}

References Nektar::StdRegions::StdExpansion::NumBndryCoeffs().

v_GetInteriorMap()

void Nektar::StdRegions::StdNodalTriExp::v_GetInteriorMap ( Array< OneD, unsigned int > &  outarray )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 391 of file StdNodalTriExp.cpp.

{
    unsigned int i;
    if (outarray.size() != GetNcoeffs() - NumBndryCoeffs())
    {
        outarray = Array<OneD, unsigned int>(GetNcoeffs() - NumBndryCoeffs());
    }
 
    for (i = NumBndryCoeffs(); i < GetNcoeffs(); i++)
    {
        outarray[i - NumBndryCoeffs()] = i;
    }
}

References Nektar::StdRegions::StdExpansion::GetNcoeffs(), and Nektar::StdRegions::StdExpansion::NumBndryCoeffs().

v_GetNodalPointsKey()

const LibUtilities::PointsKey Nektar::StdRegions::StdNodalTriExp::v_GetNodalPointsKey ( ) const

inlineoverrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 72 of file StdNodalTriExp.h.

    {
        return m_nodalPointsKey;
    };

References m_nodalPointsKey.

v_GetTraceBasisKey()

const LibUtilities::BasisKey Nektar::StdRegions::StdNodalTriExp::v_GetTraceBasisKey ( const int  i, const int  j, bool  UseGLL = false  ) const

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 214 of file StdNodalTriExp.cpp.

{
    ASSERTL2(i >= 0 && i <= 2, "edge id is out of range");
 
    // Get basiskey (0 or 1) according to edge id i
    int dir = (i != 0);
 
    switch (m_base[dir]->GetBasisType())
    {
        case LibUtilities::eOrtho_A:
        case LibUtilities::eModified_A:
        {
            switch (m_base[dir]->GetPointsType())
            {
                case LibUtilities::eGaussLobattoLegendre:
                {
                    LibUtilities::PointsKey pkey(
                        m_base[dir]
                            ->GetBasisKey()
                            .GetPointsKey()
                            .GetNumPoints(),
                        LibUtilities::eGaussLobattoLegendre);
                    return LibUtilities::BasisKey(LibUtilities::eGLL_Lagrange,
                                                  m_base[dir]->GetNumModes(),
                                                  pkey);
                }
                break;
                default:
                {
                    NEKERROR(ErrorUtil::efatal,
                             "Unexpected points distribution " +
                                 LibUtilities::kPointsTypeStr
                                     [m_base[dir]->GetPointsType()] +
                                 " in StdNodalTriExp::v_GetTraceBasisKey");
                }
            }
        }
        break;
        case LibUtilities::eModified_B:
        case LibUtilities::eOrtho_B:
        {
            switch (m_base[dir]->GetPointsType())
            {
                case LibUtilities::eGaussRadauMAlpha1Beta0:
                {
                    LibUtilities::PointsKey pkey(
                        m_base[dir]
                                ->GetBasisKey()
                                .GetPointsKey()
                                .GetNumPoints() +
                            1,
                        LibUtilities::eGaussLobattoLegendre);
                    return LibUtilities::BasisKey(LibUtilities::eGLL_Lagrange,
                                                  m_base[dir]->GetNumModes(),
                                                  pkey);
                }
                break;
                default:
                {
                    NEKERROR(ErrorUtil::efatal,
                             "Unexpected points distribution " +
                                 LibUtilities::kPointsTypeStr
                                     [m_base[dir]->GetPointsType()] +
                                 " in StdNodalTriExp::v_GetTraceBasisKey");
                }
            }
        }
        break;
        default:
        {
            NEKERROR(ErrorUtil::efatal,
                     "Information not available to set edge key");
        }
    }
    return LibUtilities::NullBasisKey;
}

References ASSERTL2, Nektar::ErrorUtil::efatal, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGLL_Lagrange, Nektar::LibUtilities::eModified_A, Nektar::LibUtilities::eModified_B, Nektar::LibUtilities::eOrtho_A, Nektar::LibUtilities::eOrtho_B, Nektar::StdRegions::StdExpansion::GetBasisType(), Nektar::StdRegions::StdExpansion::GetNumPoints(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::LibUtilities::kPointsTypeStr, Nektar::StdRegions::StdExpansion::m_base, NEKERROR, and Nektar::LibUtilities::NullBasisKey().

v_GetTraceInteriorToElementMap()

void Nektar::StdRegions::StdNodalTriExp::v_GetTraceInteriorToElementMap ( const int  eid, Array< OneD, unsigned int > &  maparray, Array< OneD, int > &  signarray, const Orientation  edgeOrient = eForwards  )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 352 of file StdNodalTriExp.cpp.

{
    ASSERTL0(eid >= 0 && eid <= 2, "Local Edge ID must be between 0 and 2");
 
    const int nEdgeIntCoeffs = GetTraceNcoeffs(eid) - 2;
 
    if (maparray.size() != nEdgeIntCoeffs)
    {
        maparray = Array<OneD, unsigned int>(nEdgeIntCoeffs);
    }
 
    if (signarray.size() != nEdgeIntCoeffs)
    {
        signarray = Array<OneD, int>(nEdgeIntCoeffs, 1);
    }
    else
    {
        fill(signarray.data(), signarray.data() + nEdgeIntCoeffs, 1);
    }
 
    Orientation orient = edgeOrient;
    if (eid == 2)
    {
        orient = orient == eForwards ? eBackwards : eForwards;
    }
 
    for (int i = 0; i < nEdgeIntCoeffs; i++)
    {
        maparray[i] = eid * nEdgeIntCoeffs + 3 + i;
    }
 
    if (orient == eBackwards)
    {
        reverse(maparray.data(), maparray.data() + nEdgeIntCoeffs);
    }
}

References ASSERTL0, Nektar::StdRegions::eBackwards, Nektar::StdRegions::eForwards, and Nektar::StdRegions::StdExpansion::GetTraceNcoeffs().

v_GetTraceToElementMap()

void Nektar::StdRegions::StdNodalTriExp::v_GetTraceToElementMap ( const int  eid, Array< OneD, unsigned int > &  maparray, Array< OneD, int > &  signarray, Orientation  edgeOrient = eForwards, int  P = -1, int  Q = -1  )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 305 of file StdNodalTriExp.cpp.

{
    ASSERTL0(eid >= 0 && eid <= 2, "Local Edge ID must be between 0 and 2");
 
    const int nEdgeCoeffs = GetTraceNcoeffs(eid);
 
    ASSERTL0(P == -1 || P == nEdgeCoeffs,
             "Nodal triangle not set up to deal with variable "
             "polynomial order.");
 
    if (maparray.size() != nEdgeCoeffs)
    {
        maparray = Array<OneD, unsigned int>(nEdgeCoeffs);
    }
 
    if (signarray.size() != nEdgeCoeffs)
    {
        signarray = Array<OneD, int>(nEdgeCoeffs, 1);
    }
    else
    {
        fill(signarray.data(), signarray.data() + nEdgeCoeffs, 1);
    }
 
    Orientation orient = edgeOrient;
    if (eid == 2)
    {
        orient = orient == eForwards ? eBackwards : eForwards;
    }
 
    maparray[0]               = eid;
    maparray[nEdgeCoeffs - 1] = eid == 2 ? 0 : eid + 1;
    for (int i = 2; i < nEdgeCoeffs; i++)
    {
        maparray[i - 1] = eid * (nEdgeCoeffs - 2) + 1 + i;
    }
 
    if (orient == eBackwards)
    {
        reverse(maparray.data(), maparray.data() + nEdgeCoeffs);
    }
}

References ASSERTL0, Nektar::StdRegions::eBackwards, Nektar::StdRegions::eForwards, Nektar::StdRegions::StdExpansion::GetTraceNcoeffs(), and Nektar::LibUtilities::P.

v_GetVertexMap()

int Nektar::StdRegions::StdNodalTriExp::v_GetVertexMap ( int  localVertexId, bool  useCoeffPacking = false  )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 297 of file StdNodalTriExp.cpp.

{
    ASSERTL0(localVertexId >= 0 && localVertexId <= 2,
             "Local Vertex ID must be between 0 and 2");
    return localVertexId;
}

References ASSERTL0.

v_HelmholtzMatrixOp()

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

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 479 of file StdNodalTriExp.cpp.

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

References Nektar::StdRegions::StdExpansion::HelmholtzMatrixOp_MatFree_GenericImpl().

v_IProductWRTBase()

void Nektar::StdRegions::StdNodalTriExp::v_IProductWRTBase ( const Array< OneD, const NekDouble > &  inarray, Array< OneD, NekDouble > &  outarray  )

overrideprotectedvirtual

Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray.

\( \begin{array}{rcl} I_{pq} = (\phi_p \phi_q, u) & = & \sum_{i=0}^{nq_0} \sum_{j=0}^{nq_1} \phi_p(\xi_{0,i}) \phi_q(\xi_{1,j}) w^0_i w^1_j u(\xi_{0,i} \xi_{1,j}) \\ & = & \sum_{i=0}^{nq_0} \phi_p(\xi_{0,i}) \sum_{j=0}^{nq_1} \phi_q(\xi_{1,j}) \tilde{u}_{i,j} \end{array} \)

where

\( \tilde{u}_{i,j} = w^0_i w^1_j u(\xi_{0,i},\xi_{1,j}) \)

which can be implemented as

\( f_{qi} = \sum_{j=0}^{nq_1} \phi_q(\xi_{1,j}) \tilde{u}_{i,j} = {\bf B_1 U} \) \( I_{pq} = \sum_{i=0}^{nq_0} \phi_p(\xi_{0,i}) f_{qi} = {\bf B_0 F} \)

This is a wrapper function around IProductWRTBaseKernel()

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 169 of file StdNodalTriExp.cpp.

{
    StdTriExp::v_IProductWRTBase(inarray, outarray);
    NodalToModalTranspose(outarray, outarray);
}

References NodalToModalTranspose(), and Nektar::StdRegions::StdExpansion::v_IProductWRTBase().

v_IProductWRTDerivBase()

void Nektar::StdRegions::StdNodalTriExp::v_IProductWRTDerivBase ( const int  dir, const Array< OneD, const NekDouble > &  inarray, Array< OneD, NekDouble > &  outarray  )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 177 of file StdNodalTriExp.cpp.

{
    StdTriExp::v_IProductWRTDerivBase(dir, inarray, outarray);
    NodalToModalTranspose(outarray, outarray);
}

References NodalToModalTranspose(), and Nektar::StdRegions::StdTriExp::v_IProductWRTDerivBase().

v_IsNodalNonTensorialExp()

bool Nektar::StdRegions::StdNodalTriExp::v_IsNodalNonTensorialExp ( )

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 65 of file StdNodalTriExp.cpp.

{
    return true;
}

v_LaplacianMatrixOp() [1/2]

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

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 456 of file StdNodalTriExp.cpp.

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

References Nektar::StdRegions::StdExpansion::LaplacianMatrixOp_MatFree_GenericImpl().

v_LaplacianMatrixOp() [2/2]

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

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 464 of file StdNodalTriExp.cpp.

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

References Nektar::StdRegions::StdExpansion::LaplacianMatrixOp_MatFree().

v_MassMatrixOp()

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

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 449 of file StdNodalTriExp.cpp.

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

References Nektar::StdRegions::StdExpansion::MassMatrixOp_MatFree().

v_NodalToModal()

void Nektar::StdRegions::StdNodalTriExp::v_NodalToModal ( const Array< OneD, const NekDouble > &  inarray, Array< OneD, NekDouble > &  outarray  )

overridevirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 74 of file StdNodalTriExp.cpp.

{
    StdMatrixKey Nkey(eInvNBasisTrans, DetShapeType(), *this,
                      NullConstFactorMap, NullVarCoeffMap, NullVarFactorsMap,
                      m_nodalPointsKey.GetPointsType());
    DNekMatSharedPtr inv_vdm = GetStdMatrix(Nkey);
 
    NekVector<NekDouble> nodal(m_ncoeffs, inarray, eWrapper);
    NekVector<NekDouble> modal(m_ncoeffs, outarray, eWrapper);
    modal = (*inv_vdm) * nodal;
}

References Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::StdRegions::eInvNBasisTrans, Nektar::eWrapper, Nektar::LibUtilities::PointsKey::GetPointsType(), Nektar::StdRegions::StdExpansion::GetStdMatrix(), Nektar::StdRegions::StdExpansion::m_ncoeffs, m_nodalPointsKey, Nektar::StdRegions::NullConstFactorMap, Nektar::StdRegions::NullVarCoeffMap, and Nektar::StdRegions::NullVarFactorsMap.

Referenced by Nektar::LocalRegions::NodalTriExp::v_BwdTrans(), and v_BwdTrans().

v_NumBndryCoeffs()

int Nektar::StdRegions::StdNodalTriExp::v_NumBndryCoeffs ( ) const

overrideprotectedvirtual

Implements Nektar::StdRegions::StdExpansion.

Definition at line 209 of file StdNodalTriExp.cpp.

{
    return 3 + (GetBasisNumModes(0) - 2) + 2 * (GetBasisNumModes(1) - 2);
}

References Nektar::StdRegions::StdExpansion::GetBasisNumModes().

v_WeakDerivMatrixOp()

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

overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 472 of file StdNodalTriExp.cpp.

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

References Nektar::StdRegions::StdExpansion::WeakDerivMatrixOp_MatFree().

Member Data Documentation

m_nodalPointsKey

LibUtilities::PointsKey Nektar::StdRegions::StdNodalTriExp::m_nodalPointsKey

protected

Definition at line 70 of file StdNodalTriExp.h.

Referenced by GetNodalPoints(), ModalToNodal(), NodalToModalTranspose(), Nektar::LocalRegions::NodalTriExp::v_CreateStdMatrix(), Nektar::LocalRegions::NodalTriExp::v_GetLinStdExp(), v_GetNodalPointsKey(), Nektar::LocalRegions::NodalTriExp::v_GetStdExp(), and v_NodalToModal().