Nektar++
Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
Nektar::LocalRegions::NodalTriExp Class Reference

#include <NodalTriExp.h>

Inheritance diagram for Nektar::LocalRegions::NodalTriExp:
[legend]

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::BasisSharedPtrGetBasis (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 GetTraceNcoeffs (const int i) const
 This function returns the number of expansion coefficients belonging to the i-th trace. More...
 
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. More...
 
const LibUtilities::BasisKey GetTraceBasisKey (const int i, int k=-1) const
 This function returns the basis key belonging to the i-th trace. More...
 
LibUtilities::PointsKey GetTracePointsKey (const int i, int k=-1) const
 This function returns the basis key belonging to the i-th trace. More...
 
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. More...
 
int GetNtraces () 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< StdExpansionGetStdExp (void) const
 
std::shared_ptr< StdExpansionGetLinStdExp (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 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. 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)
 
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}\) 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)
 
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...
 
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. More...
 
void LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
 Convert local cartesian coordinate xi into local collapsed coordinates eta. More...
 
void LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi)
 Convert local collapsed coordinates eta into local cartesian coordinate xi. More...
 
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)
 
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 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. 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)
 
void GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat)
 
- Public Member Functions inherited from Nektar::LocalRegions::Expansion2D
 Expansion2D (SpatialDomains::Geometry2DSharedPtr pGeom)
 
virtual ~Expansion2D ()
 
DNekScalMatSharedPtr CreateMatrix (const MatrixKey &mkey)
 
void SetTraceToGeomOrientation (Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &inout)
 
Array< OneD, unsigned int > GetTraceInverseBoundaryMap (int eid)
 
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)
 
SpatialDomains::Geometry2DSharedPtr GetGeom2D () const
 
void ReOrientEdgePhysMap (const int nvert, const StdRegions::Orientation orient, const int nq0, Array< OneD, int > &idmap)
 
virtual void v_GenTraceExp (const int traceid, ExpansionSharedPtr &exp)
 
- Public Member Functions inherited from Nektar::LocalRegions::Expansion
 Expansion (SpatialDomains::GeometrySharedPtr pGeom)
 
 Expansion (const Expansion &pSrc)
 
virtual ~Expansion ()
 
void SetTraceExp (const int traceid, ExpansionSharedPtr &f)
 
ExpansionSharedPtr GetTraceExp (const int traceid)
 
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 ()
 
IndexMapValuesSharedPtr CreateIndexMap (const IndexMapKey &ikey)
 
DNekScalBlkMatSharedPtr CreateStaticCondMatrix (const MatrixKey &mkey)
 
const SpatialDomains::GeomFactorsSharedPtrGetMetricInfo () 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)
 
void NormalTraceDerivFactors (Array< OneD, Array< OneD, NekDouble >> &factors, Array< OneD, Array< OneD, NekDouble >> &d0factors, Array< OneD, Array< OneD, NekDouble >> &d1factors)
 
IndexMapValuesSharedPtr GetIndexMap (const IndexMapKey &ikey)
 
void AlignVectorToCollapsedDir (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
 
ExpansionSharedPtr GetLeftAdjacentElementExp () const
 
ExpansionSharedPtr GetRightAdjacentElementExp () const
 
int GetLeftAdjacentElementTrace () const
 
int GetRightAdjacentElementTrace () const
 
void SetAdjacentElementExp (int traceid, ExpansionSharedPtr &e)
 
StdRegions::Orientation GetTraceOrient (int trace)
 
void SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void DivideByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 Divided by the metric jacobi and quadrature weights. More...
 
void GetTraceQFactors (const int trace, 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 GetTracePhysVals (const int trace, const StdRegions::StdExpansionSharedPtr &TraceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient=StdRegions::eNoOrientation)
 
void GetTracePhysMap (const int edge, Array< OneD, int > &outarray)
 
void ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1)
 
const NormalVectorGetTraceNormal (const int id)
 
void ComputeTraceNormal (const int id)
 
const Array< OneD, const NekDouble > & GetPhysNormals (void)
 
void SetPhysNormals (Array< OneD, const NekDouble > &normal)
 
void SetUpPhysNormals (const int trace)
 
void AddRobinMassMatrix (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)
 
void TraceNormLen (const int traceid, NekDouble &h, NekDouble &p)
 
virtual void AddRobinTraceContribution (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, const Array< OneD, NekDouble > &incoeffs, Array< OneD, NekDouble > &coeffs)
 
const Array< OneD, const NekDouble > & GetElmtBndNormDirElmtLen (const int nbnd) const
 
void StdDerivBaseOnTraceMat (Array< OneD, DNekMatSharedPtr > &DerivMat)
 

Protected Member Functions

DNekMatSharedPtr CreateStdMatrix (const StdRegions::StdMatrixKey &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 int v_GetVertexMap (int localVertexId, bool useCoeffPacking=false)
 
virtual void v_GetTraceToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient=eForwards, int P=-1, int Q=-1)
 
virtual void v_GetTraceInteriorToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation edgeOrient=eForwards)
 
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_FwdTransBndConstrained (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 void v_LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi)
 
NekDouble v_PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode) final
 
virtual int v_GetNverts () const
 
virtual int v_GetNtraces () const
 
virtual LibUtilities::ShapeType v_DetShapeType () const
 
virtual int v_NumDGBndryCoeffs () const
 
virtual int v_GetTraceNcoeffs (const int i) const
 
virtual int v_GetTraceNumPoints (const int i) const
 
virtual int v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset)
 
virtual bool v_IsBoundaryInteriorExpansion ()
 
virtual const LibUtilities::BasisKey v_GetTraceBasisKey (const int i, const int j) const
 
virtual void v_GetTraceCoeffMap (const unsigned int traceid, Array< OneD, unsigned int > &maparray)
 
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 void v_GetElmtTraceToTraceMap (const unsigned int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient, int P, int Q)
 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. More...
 
virtual void v_GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat)
 
- 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...
 
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 (Array< OneD, NekDouble > &coeffs, StdRegions::Orientation dir)
 
virtual NekDouble v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals)
 
template<int DIR>
NekDouble BaryEvaluate (const NekDouble &coord, const NekDouble *physvals)
 This function performs the barycentric interpolation of the polynomial stored in coord at a point physvals using barycentric interpolation weights in direction. More...
 
template<int DIR>
NekDouble BaryEvaluateBasis (const NekDouble &coord, const int &mode)
 
- Protected Member Functions inherited from Nektar::LocalRegions::Expansion2D
virtual Array< OneD, NekDoublev_GetMF (const int dir, const int shapedim, const StdRegions::VarCoeffMap &varcoeffs)
 
virtual Array< OneD, NekDoublev_GetMFDiv (const int dir, const StdRegions::VarCoeffMap &varcoeffs)
 
virtual Array< OneD, NekDoublev_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_AddRobinTraceContribution (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, const Array< OneD, NekDouble > &incoeffs, 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, NekDoublev_GetnEdgecdotMF (const int dir, const int edge, ExpansionSharedPtr &EdgeExp_e, const Array< OneD, const Array< OneD, NekDouble >> &normals, const StdRegions::VarCoeffMap &varcoeffs)
 
void v_ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1)
 
virtual void v_SetUpPhysNormals (const int edge)
 
virtual NekDouble v_VectorFlux (const Array< OneD, Array< OneD, NekDouble >> &vec)
 
virtual void v_TraceNormLen (const int traceid, NekDouble &h, NekDouble &p)
 
- 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_DivideByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_ComputeLaplacianMetric ()
 
Array< OneD, NekDoublev_GetMF (const int dir, const int shapedim, const StdRegions::VarCoeffMap &varcoeffs)
 
Array< OneD, NekDoublev_GetMFDiv (const int dir, const StdRegions::VarCoeffMap &varcoeffs)
 
Array< OneD, NekDoublev_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)
 
virtual void v_NormalTraceDerivFactors (Array< OneD, Array< OneD, NekDouble >> &factors, Array< OneD, Array< OneD, NekDouble >> &d0factors, Array< OneD, Array< OneD, NekDouble >> &d1factors)
 
virtual StdRegions::Orientation v_GetTraceOrient (int trace)
 
virtual void v_SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_GetTraceQFactors (const int trace, Array< OneD, NekDouble > &outarray)
 
virtual void v_GetTracePhysVals (const int trace, const StdRegions::StdExpansionSharedPtr &TraceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient)
 
virtual void v_GetTracePhysMap (const int edge, Array< OneD, int > &outarray)
 
virtual const Array< OneD, const NekDouble > & v_GetPhysNormals (void)
 
virtual void v_SetPhysNormals (Array< OneD, const NekDouble > &normal)
 

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_AlignVectorToCollapsedDir (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, 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_ComputeTraceNormal (const int edge)
 

Private Attributes

LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLessm_matrixManager
 
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLessm_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::BasisSharedPtrm_base
 
int m_elmt_id
 
int m_ncoeffs
 
LibUtilities::NekManager< StdMatrixKey, DNekMat, StdMatrixKey::opLessm_stdMatrixManager
 
LibUtilities::NekManager< StdMatrixKey, DNekBlkMat, StdMatrixKey::opLessm_stdStaticCondMatrixManager
 
- Protected Attributes inherited from Nektar::LocalRegions::Expansion2D
std::vector< bool > m_requireNeg
 
- Protected Attributes inherited from Nektar::LocalRegions::Expansion
LibUtilities::NekManager< IndexMapKey, IndexMapValues, IndexMapKey::opLessm_indexMapManager
 
std::map< int, ExpansionWeakPtrm_traceExp
 
SpatialDomains::GeometrySharedPtr m_geom
 
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
 
MetricMap m_metrics
 
std::map< int, NormalVectorm_traceNormals
 
ExpansionWeakPtr m_elementLeft
 
ExpansionWeakPtr m_elementRight
 
int m_elementTraceLeft = -1
 
int m_elementTraceRight = -1
 
std::map< int, Array< OneD, NekDouble > > m_elmtBndNormDirElmtLen
 the element length in each element boundary(Vertex, edge or face) normal direction calculated based on the local m_metricinfo times the standard element length (which is 2.0) More...
 

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.

49  Ba.GetNumModes(), (Bb.GetNumModes())),
50  2, Ba, Bb),
52  Ba.GetNumModes(), (Bb.GetNumModes())),
53  Ba, Bb),
54  StdNodalTriExp(Ba, Bb, Ntype), Expansion(geom), Expansion2D(geom),
56  std::bind(&Expansion2D::CreateMatrix, this, std::placeholders::_1),
57  std::string("NodalTriExpMatrix")),
59  this, std::placeholders::_1),
60  std::string("NodalTriExpStaticCondMatrix"))
61 {
62 }
DNekScalMatSharedPtr CreateMatrix(const MatrixKey &mkey)
Definition: Expansion2D.cpp:59
Expansion2D(SpatialDomains::Geometry2DSharedPtr pGeom)
Definition: Expansion2D.cpp:54
DNekScalBlkMatSharedPtr CreateStaticCondMatrix(const MatrixKey &mkey)
Definition: Expansion.cpp:272
Expansion(SpatialDomains::GeometrySharedPtr pGeom)
Definition: Expansion.cpp:47
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: NodalTriExp.h:191
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: NodalTriExp.h:193
StdExpansion()
Default Constructor.
int getNumberOfCoefficients(int Na, int Nb)
Definition: ShapeType.hpp:112

◆ NodalTriExp() [2/2]

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

Copy Constructor.

Definition at line 64 of file NodalTriExp.cpp.

65  : StdExpansion(T), StdExpansion2D(T), StdRegions::StdNodalTriExp(T),
66  Expansion(T), Expansion2D(T), m_matrixManager(T.m_matrixManager),
67  m_staticCondMatrixManager(T.m_staticCondMatrixManager)
68 {
69 }

◆ ~NodalTriExp()

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

Destructor.

Definition at line 71 of file NodalTriExp.cpp.

72 {
73 }

Member Function Documentation

◆ CreateStdMatrix()

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

Definition at line 440 of file NodalTriExp.cpp.

442 {
443  LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
444  LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
448 
449  return tmp->GetStdMatrix(mkey);
450 }
PointsType GetPointsType() const
Definition: Points.h:109
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Array< OneD, LibUtilities::BasisSharedPtr > m_base
LibUtilities::PointsKey m_nodalPointsKey
std::shared_ptr< StdNodalTriExp > StdNodalTriExpSharedPtr

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

Referenced by v_CreateStdMatrix().

◆ 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 373 of file NodalTriExp.cpp.

375 {
376  IProductWRTBase(inarray, outarray);
377 
378  // get Mass matrix inverse
379  MatrixKey masskey(StdRegions::eInvMass, DetShapeType(), *this,
383  DNekScalMatSharedPtr matsys = m_matrixManager[masskey];
384 
385  // copy inarray in case inarray == outarray
386  NekVector<NekDouble> in(m_ncoeffs, outarray, eCopy);
387  NekVector<NekDouble> out(m_ncoeffs, outarray, eWrapper);
388 
389  out = (*matsys) * in;
390 }
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 ...
Definition: NodalTriExp.h:86
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:375
static ConstFactorMap NullConstFactorMap
Definition: StdRegions.hpp:283
static VarCoeffMap NullVarCoeffMap
Definition: StdRegions.hpp:241
std::shared_ptr< DNekScalMat > DNekScalMatSharedPtr

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 v_FwdTrans().

◆ 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 392 of file NodalTriExp.cpp.

395 {
397 
398  if (inarray.get() == outarray.get())
399  {
400  Array<OneD, NekDouble> tmp(m_ncoeffs);
401  Vmath::Vcopy(m_ncoeffs, inarray.get(), 1, tmp.get(), 1);
402 
403  Blas::Dgemv('N', m_ncoeffs, m_ncoeffs, mat->Scale(),
404  (mat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
405  tmp.get(), 1, 0.0, outarray.get(), 1);
406  }
407  else
408  {
409  Blas::Dgemv('N', m_ncoeffs, m_ncoeffs, mat->Scale(),
410  (mat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
411  inarray.get(), 1, 0.0, outarray.get(), 1);
412  }
413 }
DNekScalMatSharedPtr GetLocMatrix(const LocalRegions::MatrixKey &mkey)
Definition: Expansion.cpp:88
static void Dgemv(const char &trans, const int &m, const int &n, const double &alpha, const double *a, const int &lda, const double *x, const int &incx, const double &beta, double *y, const int &incy)
BLAS level 2: Matrix vector multiply y = A x where A[m x n].
Definition: Blas.hpp:246
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1255

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

◆ GetCoord()

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

Definition at line 423 of file NodalTriExp.cpp.

425 {
426  int i;
427 
428  ASSERTL1(Lcoords[0] >= -1.0 && Lcoords[1] <= 1.0 && Lcoords[1] >= -1.0 &&
429  Lcoords[1] <= 1.0,
430  "Local coordinates are not in region [-1,1]");
431 
432  m_geom->FillGeom();
433 
434  for (i = 0; i < m_geom->GetCoordim(); ++i)
435  {
436  coords[i] = m_geom->GetCoord(i, Lcoords);
437  }
438 }
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode....
Definition: ErrorUtil.hpp:249
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:272

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

Referenced by v_GetCoord().

◆ 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 415 of file NodalTriExp.cpp.

418 {
419  Expansion::v_GetCoords(coords_0, coords_1, coords_2);
420 }
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3)
Definition: Expansion.cpp:524

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

Referenced by v_GetCoords().

◆ HelmholtzMatrixOp()

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

Definition at line 155 of file NodalTriExp.h.

158  {
159  StdExpansion::HelmholtzMatrixOp_MatFree_GenericImpl(inarray, outarray,
160  mkey);
161  }

Referenced by v_HelmholtzMatrixOp().

◆ 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 95 of file NodalTriExp.cpp.

96 {
97  int nquad0 = m_base[0]->GetNumPoints();
98  int nquad1 = m_base[1]->GetNumPoints();
99  Array<OneD, const NekDouble> jac = m_metricinfo->GetJac(GetPointsKeys());
100  NekDouble ival;
101  Array<OneD, NekDouble> tmp(nquad0 * nquad1);
102 
103  // multiply inarray with Jacobian
104  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
105  {
106  Vmath::Vmul(nquad0 * nquad1, jac, 1, inarray, 1, tmp, 1);
107  }
108  else
109  {
110  Vmath::Smul(nquad0 * nquad1, jac[0], inarray, 1, tmp, 1);
111  }
112 
113  // call StdQuadExp version;
114  ival = StdNodalTriExp::v_Integral(tmp);
115  return ival;
116 }
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:273
const LibUtilities::PointsKeyVector GetPointsKeys() const
@ eDeformed
Geometry is curved or has non-constant factors.
double NekDouble
void Vmul(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x*y.
Definition: Vmath.cpp:209
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*x.
Definition: Vmath.cpp:248

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().

◆ 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 86 of file NodalTriExp.h.

88  {
89  NodalTriExp::IProductWRTBase_SumFac(inarray, outarray);
90  }
void IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)

References IProductWRTBase_SumFac().

Referenced by FwdTrans(), and v_IProductWRTBase().

◆ IProductWRTBase_MatOp()

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

Definition at line 147 of file NodalTriExp.cpp.

150 {
151  int nq = GetTotPoints();
152  MatrixKey iprodmatkey(StdRegions::eIProductWRTBase, DetShapeType(), *this);
153  DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
154 
155  Blas::Dgemv('N', m_ncoeffs, nq, iprodmat->Scale(),
156  (iprodmat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
157  inarray.get(), 1, 0.0, outarray.get(), 1);
158 }
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
Definition: StdExpansion.h:140

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

◆ 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 118 of file NodalTriExp.cpp.

121 {
122  int nquad0 = m_base[0]->GetNumPoints();
123  int nquad1 = m_base[1]->GetNumPoints();
124  int order1 = m_base[1]->GetNumModes();
125 
126  if (multiplybyweights)
127  {
128  Array<OneD, NekDouble> tmp(nquad0 * nquad1 + nquad0 * order1);
129  Array<OneD, NekDouble> wsp(tmp + nquad0 * nquad1);
130 
131  MultiplyByQuadratureMetric(inarray, tmp);
132  StdTriExp::IProductWRTBase_SumFacKernel(
133  m_base[0]->GetBdata(), m_base[1]->GetBdata(), tmp, outarray, wsp);
134  NodalToModalTranspose(outarray, outarray);
135  }
136  else
137  {
138  Array<OneD, NekDouble> wsp(nquad0 * order1);
139 
140  StdTriExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),
141  m_base[1]->GetBdata(), inarray,
142  outarray, wsp);
143  NodalToModalTranspose(outarray, outarray);
144  }
145 }
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:731
void NodalToModalTranspose(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

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

Referenced by IProductWRTBase(), and v_IProductWRTBase_SumFac().

◆ IProductWRTDerivBase()

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

Definition at line 92 of file NodalTriExp.h.

95  {
96  NodalTriExp::IProductWRTDerivBase_SumFac(dir, inarray, outarray);
97  }
void IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)

References IProductWRTDerivBase_SumFac().

Referenced by v_IProductWRTDerivBase().

◆ 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 256 of file NodalTriExp.cpp.

259 {
260  int nq = GetTotPoints();
262 
263  switch (dir)
264  {
265  case 0:
266  {
268  }
269  break;
270  case 1:
271  {
273  }
274  break;
275  case 2:
276  {
278  }
279  break;
280  default:
281  {
282  ASSERTL1(false, "input dir is out of range");
283  }
284  break;
285  }
286 
287  MatrixKey iprodmatkey(mtype, DetShapeType(), *this);
288  DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
289 
290  Blas::Dgemv('N', m_ncoeffs, nq, iprodmat->Scale(),
291  (iprodmat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
292  inarray.get(), 1, 0.0, outarray.get(), 1);
293 }

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.

◆ 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 160 of file NodalTriExp.cpp.

163 {
164  int nquad0 = m_base[0]->GetNumPoints();
165  int nquad1 = m_base[1]->GetNumPoints();
166  int nqtot = nquad0 * nquad1;
167  int wspsize = max(nqtot, m_ncoeffs);
168 
169  Array<OneD, NekDouble> tmp0(4 * wspsize);
170  Array<OneD, NekDouble> tmp1(tmp0 + wspsize);
171  Array<OneD, NekDouble> tmp2(tmp0 + 2 * wspsize);
172  Array<OneD, NekDouble> tmp3(tmp0 + 3 * wspsize);
173 
174  Array<OneD, Array<OneD, NekDouble>> tmp2D{2};
175  tmp2D[0] = tmp1;
176  tmp2D[1] = tmp2;
177 
178  AlignVectorToCollapsedDir(dir, inarray, tmp2D);
179 
180  MultiplyByQuadratureMetric(tmp1, tmp1);
181  MultiplyByQuadratureMetric(tmp2, tmp2);
182 
183  IProductWRTBase_SumFacKernel(m_base[0]->GetDbdata(), m_base[1]->GetBdata(),
184  tmp1, tmp3, tmp0);
185  IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(), m_base[1]->GetDbdata(),
186  tmp2, outarray, tmp0);
187  Vmath::Vadd(m_ncoeffs, tmp3, 1, outarray, 1, outarray, 1);
188 
189  NodalToModalTranspose(outarray, outarray);
190 }
void AlignVectorToCollapsedDir(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
Definition: Expansion.h:150
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)
void Vadd(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Add vector z = x+y.
Definition: Vmath.cpp:359

References Nektar::LocalRegions::Expansion::AlignVectorToCollapsedDir(), Nektar::StdRegions::StdExpansion2D::IProductWRTBase_SumFacKernel(), Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), Nektar::StdRegions::StdNodalTriExp::NodalToModalTranspose(), and Vmath::Vadd().

Referenced by IProductWRTDerivBase(), and v_IProductWRTDerivBase_SumFac().

◆ 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 130 of file NodalTriExp.h.

133  {
134  StdExpansion::LaplacianMatrixOp_MatFree_GenericImpl(inarray, outarray,
135  mkey);
136  }

Referenced by v_LaplacianMatrixOp().

◆ 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 138 of file NodalTriExp.h.

142  {
143  StdExpansion::LaplacianMatrixOp_MatFree(k1, k2, inarray, outarray,
144  mkey);
145  }

◆ MassMatrixOp()

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

Definition at line 123 of file NodalTriExp.h.

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

Referenced by v_MassMatrixOp().

◆ 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 302 of file NodalTriExp.cpp.

306 {
307  int nquad0 = m_base[0]->GetNumPoints();
308  int nquad1 = m_base[1]->GetNumPoints();
309  int nqtot = nquad0 * nquad1;
310  const Array<TwoD, const NekDouble> &df =
311  m_metricinfo->GetDerivFactors(GetPointsKeys());
312 
313  Array<OneD, NekDouble> diff0(2 * nqtot);
314  Array<OneD, NekDouble> diff1(diff0 + nqtot);
315 
316  StdNodalTriExp::v_PhysDeriv(inarray, diff0, diff1);
317 
318  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
319  {
320  if (out_d0.size())
321  {
322  Vmath::Vmul(nqtot, df[0], 1, diff0, 1, out_d0, 1);
323  Vmath::Vvtvp(nqtot, df[1], 1, diff1, 1, out_d0, 1, out_d0, 1);
324  }
325 
326  if (out_d1.size())
327  {
328  Vmath::Vmul(nqtot, df[2], 1, diff0, 1, out_d1, 1);
329  Vmath::Vvtvp(nqtot, df[3], 1, diff1, 1, out_d1, 1, out_d1, 1);
330  }
331 
332  if (out_d2.size())
333  {
334  Vmath::Vmul(nqtot, df[4], 1, diff0, 1, out_d2, 1);
335  Vmath::Vvtvp(nqtot, df[5], 1, diff1, 1, out_d2, 1, out_d2, 1);
336  }
337  }
338  else // regular geometry
339  {
340  if (out_d0.size())
341  {
342  Vmath::Smul(nqtot, df[0][0], diff0, 1, out_d0, 1);
343  Blas::Daxpy(nqtot, df[1][0], diff1, 1, out_d0, 1);
344  }
345 
346  if (out_d1.size())
347  {
348  Vmath::Smul(nqtot, df[2][0], diff0, 1, out_d1, 1);
349  Blas::Daxpy(nqtot, df[3][0], diff1, 1, out_d1, 1);
350  }
351 
352  if (out_d2.size())
353  {
354  Vmath::Smul(nqtot, df[4][0], diff0, 1, out_d2, 1);
355  Blas::Daxpy(nqtot, df[5][0], diff1, 1, out_d2, 1);
356  }
357  }
358 }
static void Daxpy(const int &n, const double &alpha, const double *x, const int &incx, const double *y, const int &incy)
BLAS level 1: y = alpha x plus y.
Definition: Blas.hpp:154
void Vvtvp(int n, const T *w, const int incw, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
vvtvp (vector times vector plus vector): z = w*x + y
Definition: Vmath.cpp:574

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 v_PhysDeriv().

◆ PhysEvaluate()

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

Definition at line 452 of file NodalTriExp.cpp.

456 {
457  Array<OneD, NekDouble> Lcoord = Array<OneD, NekDouble>(2);
458 
459  ASSERTL0(m_geom, "m_geom not defined");
460  m_geom->GetLocCoords(coord, Lcoord);
461 
462  return StdNodalTriExp::v_PhysEvaluate(Lcoord, physvals);
463 }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:215

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

Referenced by v_PhysEvaluate().

◆ v_AlignVectorToCollapsedDir()

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

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 192 of file NodalTriExp.cpp.

195 {
196  ASSERTL1((dir == 0) || (dir == 1) || (dir == 2), "Invalid direction.");
197  ASSERTL1((dir == 2) ? (m_geom->GetCoordim() == 3) : true,
198  "Invalid direction.");
199 
200  int nquad0 = m_base[0]->GetNumPoints();
201  int nquad1 = m_base[1]->GetNumPoints();
202  int nqtot = nquad0 * nquad1;
203  int wspsize = max(nqtot, m_ncoeffs);
204 
205  const Array<TwoD, const NekDouble> &df =
206  m_metricinfo->GetDerivFactors(GetPointsKeys());
207 
208  Array<OneD, NekDouble> tmp0(4 * wspsize);
209  Array<OneD, NekDouble> tmp3(tmp0 + wspsize);
210  Array<OneD, NekDouble> gfac0(tmp0 + 2 * wspsize);
211  Array<OneD, NekDouble> gfac1(tmp0 + 3 * wspsize);
212 
213  Array<OneD, NekDouble> tmp1 = outarray[0];
214  Array<OneD, NekDouble> tmp2 = outarray[1];
215 
216  const Array<OneD, const NekDouble> &z0 = m_base[0]->GetZ();
217  const Array<OneD, const NekDouble> &z1 = m_base[1]->GetZ();
218 
219  // set up geometric factor: 2/(1-z1)
220  for (int i = 0; i < nquad1; ++i)
221  {
222  gfac0[i] = 2.0 / (1 - z1[i]);
223  }
224  for (int i = 0; i < nquad0; ++i)
225  {
226  gfac1[i] = 0.5 * (1 + z0[i]);
227  }
228 
229  for (int i = 0; i < nquad1; ++i)
230  {
231  Vmath::Smul(nquad0, gfac0[i], &inarray[0] + i * nquad0, 1,
232  &tmp0[0] + i * nquad0, 1);
233  }
234 
235  for (int i = 0; i < nquad1; ++i)
236  {
237  Vmath::Vmul(nquad0, &gfac1[0], 1, &tmp0[0] + i * nquad0, 1,
238  &tmp1[0] + i * nquad0, 1);
239  }
240 
241  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
242  {
243  Vmath::Vmul(nqtot, &df[2 * dir][0], 1, &tmp0[0], 1, &tmp0[0], 1);
244  Vmath::Vmul(nqtot, &df[2 * dir + 1][0], 1, &tmp1[0], 1, &tmp1[0], 1);
245  Vmath::Vmul(nqtot, &df[2 * dir + 1][0], 1, &inarray[0], 1, &tmp2[0], 1);
246  }
247  else
248  {
249  Vmath::Smul(nqtot, df[2 * dir][0], tmp0, 1, tmp0, 1);
250  Vmath::Smul(nqtot, df[2 * dir + 1][0], tmp1, 1, tmp1, 1);
251  Vmath::Smul(nqtot, df[2 * dir + 1][0], inarray, 1, tmp2, 1);
252  }
253  Vmath::Vadd(nqtot, tmp0, 1, tmp1, 1, tmp1, 1);
254 }

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

◆ 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 330 of file NodalTriExp.h.

332  {
333  StdNodalTriExp::v_BwdTrans_SumFac(inarray, outarray);
334  }

◆ v_ComputeTraceNormal()

void Nektar::LocalRegions::NodalTriExp::v_ComputeTraceNormal ( const int  edge)
privatevirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 505 of file NodalTriExp.cpp.

506 {
507  int i;
508  const SpatialDomains::GeomFactorsSharedPtr &geomFactors =
509  GetGeom()->GetMetricInfo();
510  const SpatialDomains::GeomType type = geomFactors->GetGtype();
511 
513  const Array<TwoD, const NekDouble> &df =
514  geomFactors->GetDerivFactors(ptsKeys);
515  const Array<OneD, const NekDouble> &jac = geomFactors->GetJac(ptsKeys);
516  int nqe = m_base[0]->GetNumPoints();
517  int dim = GetCoordim();
518 
519  m_traceNormals[edge] = Array<OneD, Array<OneD, NekDouble>>(dim);
520  Array<OneD, Array<OneD, NekDouble>> &normal = m_traceNormals[edge];
521  for (i = 0; i < dim; ++i)
522  {
523  normal[i] = Array<OneD, NekDouble>(nqe);
524  }
525 
526  size_t nqb = nqe;
527  size_t nbnd = edge;
528  m_elmtBndNormDirElmtLen[nbnd] = Array<OneD, NekDouble>{nqb, 0.0};
529  Array<OneD, NekDouble> &length = m_elmtBndNormDirElmtLen[nbnd];
530 
531  // Regular geometry case
532  if ((type == SpatialDomains::eRegular) ||
534  {
535  NekDouble fac;
536  // Set up normals
537  switch (edge)
538  {
539  case 0:
540  for (i = 0; i < GetCoordim(); ++i)
541  {
542  Vmath::Fill(nqe, -df[2 * i + 1][0], normal[i], 1);
543  }
544  break;
545  case 1:
546  for (i = 0; i < GetCoordim(); ++i)
547  {
548  Vmath::Fill(nqe, df[2 * i + 1][0] + df[2 * i][0], normal[i],
549  1);
550  }
551  break;
552  case 2:
553  for (i = 0; i < GetCoordim(); ++i)
554  {
555  Vmath::Fill(nqe, -df[2 * i][0], normal[i], 1);
556  }
557  break;
558  default:
559  ASSERTL0(false, "Edge is out of range (edge < 3)");
560  }
561 
562  // normalise
563  fac = 0.0;
564  for (i = 0; i < GetCoordim(); ++i)
565  {
566  fac += normal[i][0] * normal[i][0];
567  }
568  fac = 1.0 / sqrt(fac);
569 
570  Vmath::Fill(nqb, fac, length, 1);
571 
572  for (i = 0; i < GetCoordim(); ++i)
573  {
574  Vmath::Smul(nqe, fac, normal[i], 1, normal[i], 1);
575  }
576  }
577  else // Set up deformed normals
578  {
579  int j;
580 
581  int nquad0 = ptsKeys[0].GetNumPoints();
582  int nquad1 = ptsKeys[1].GetNumPoints();
583 
584  LibUtilities::PointsKey from_key;
585 
586  Array<OneD, NekDouble> normals(GetCoordim() * max(nquad0, nquad1), 0.0);
587  Array<OneD, NekDouble> edgejac(GetCoordim() * max(nquad0, nquad1), 0.0);
588 
589  // Extract Jacobian along edges and recover local
590  // derivates (dx/dr) for polynomial interpolation by
591  // multiplying m_gmat by jacobian
592  switch (edge)
593  {
594  case 0:
595  for (j = 0; j < nquad0; ++j)
596  {
597  edgejac[j] = jac[j];
598  for (i = 0; i < GetCoordim(); ++i)
599  {
600  normals[i * nquad0 + j] =
601  -df[2 * i + 1][j] * edgejac[j];
602  }
603  }
604  from_key = ptsKeys[0];
605  break;
606  case 1:
607  for (j = 0; j < nquad1; ++j)
608  {
609  edgejac[j] = jac[nquad0 * j + nquad0 - 1];
610  for (i = 0; i < GetCoordim(); ++i)
611  {
612  normals[i * nquad1 + j] =
613  (df[2 * i][nquad0 * j + nquad0 - 1] +
614  df[2 * i + 1][nquad0 * j + nquad0 - 1]) *
615  edgejac[j];
616  }
617  }
618  from_key = ptsKeys[1];
619  break;
620  case 2:
621  for (j = 0; j < nquad1; ++j)
622  {
623  edgejac[j] = jac[nquad0 * j];
624  for (i = 0; i < GetCoordim(); ++i)
625  {
626  normals[i * nquad1 + j] =
627  -df[2 * i][nquad0 * j] * edgejac[j];
628  }
629  }
630  from_key = ptsKeys[1];
631  break;
632  default:
633  ASSERTL0(false, "edge is out of range (edge < 3)");
634  }
635 
636  int nq = from_key.GetNumPoints();
637  Array<OneD, NekDouble> work(nqe, 0.0);
638 
639  // interpolate Jacobian and invert
640  LibUtilities::Interp1D(from_key, jac, m_base[0]->GetPointsKey(), work);
641  Vmath::Sdiv(nq, 1.0, &work[0], 1, &work[0], 1);
642 
643  // interpolate
644  for (i = 0; i < GetCoordim(); ++i)
645  {
646  LibUtilities::Interp1D(from_key, &normals[i * nq],
647  m_base[0]->GetPointsKey(), &normal[i][0]);
648  Vmath::Vmul(nqe, work, 1, normal[i], 1, normal[i], 1);
649  }
650 
651  // normalise normal vectors
652  Vmath::Zero(nqe, work, 1);
653  for (i = 0; i < GetCoordim(); ++i)
654  {
655  Vmath::Vvtvp(nqe, normal[i], 1, normal[i], 1, work, 1, work, 1);
656  }
657 
658  Vmath::Vsqrt(nqe, work, 1, work, 1);
659  Vmath::Sdiv(nqe, 1.0, work, 1, work, 1);
660 
661  Vmath::Vcopy(nqb, work, 1, length, 1);
662 
663  for (i = 0; i < GetCoordim(); ++i)
664  {
665  Vmath::Vmul(nqe, normal[i], 1, work, 1, normal[i], 1);
666  }
667 
668  // Reverse direction so that points are in
669  // anticlockwise direction if edge >=2
670  if (edge >= 2)
671  {
672  for (i = 0; i < GetCoordim(); ++i)
673  {
674  Vmath::Reverse(nqe, normal[i], 1, normal[i], 1);
675  }
676  }
677  }
678 }
std::map< int, NormalVector > m_traceNormals
Definition: Expansion.h:275
std::map< int, Array< OneD, NekDouble > > m_elmtBndNormDirElmtLen
the element length in each element boundary(Vertex, edge or face) normal direction calculated based o...
Definition: Expansion.h:285
SpatialDomains::GeometrySharedPtr GetGeom() const
Definition: Expansion.cpp:166
void Interp1D(const BasisKey &fbasis0, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, Array< OneD, NekDouble > &to)
this function interpolates a 1D function evaluated at the quadrature points of the basis fbasis0 to ...
Definition: Interp.cpp:52
std::vector< PointsKey > PointsKeyVector
Definition: Points.h:250
std::shared_ptr< GeomFactors > GeomFactorsSharedPtr
Pointer to a GeomFactors object.
Definition: GeomFactors.h:62
GeomType
Indicates the type of element geometry.
@ eRegular
Geometry is straight-sided with constant geometric factors.
@ eMovingRegular
Currently unused.
void Vsqrt(int n, const T *x, const int incx, T *y, const int incy)
sqrt y = sqrt(x)
Definition: Vmath.cpp:534
void Sdiv(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha/y.
Definition: Vmath.cpp:324
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:492
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:45
void Reverse(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1286
scalarT< T > sqrt(scalarT< T > in)
Definition: scalar.hpp:291

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::Expansion::m_elmtBndNormDirElmtLen, Nektar::LocalRegions::Expansion::m_traceNormals, Vmath::Reverse(), Vmath::Sdiv(), Vmath::Smul(), tinysimd::sqrt(), Vmath::Vcopy(), Vmath::Vmul(), Vmath::Vsqrt(), Vmath::Vvtvp(), and Vmath::Zero().

◆ v_CreateStdMatrix()

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

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 304 of file NodalTriExp.h.

306  {
307  return CreateStdMatrix(mkey);
308  }
DNekMatSharedPtr CreateStdMatrix(const StdRegions::StdMatrixKey &mkey)

References CreateStdMatrix().

◆ 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 289 of file NodalTriExp.h.

291  {
292  FwdTrans(inarray, outarray);
293  }
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 coeffic...

References FwdTrans().

◆ v_GenMatrix()

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

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 484 of file NodalTriExp.cpp.

485 {
486  DNekMatSharedPtr returnval;
487 
488  switch (mkey.GetMatrixType())
489  {
496  returnval = Expansion2D::v_GenMatrix(mkey);
497  break;
498  default:
499  returnval = StdNodalTriExp::v_GenMatrix(mkey);
500  break;
501  }
502  return returnval;
503 }
virtual DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey)
std::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:75

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().

◆ v_GenNBasisTransMatrix()

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

Definition at line 200 of file NodalTriExp.h.

201  {
202  return StdNodalTriExp::GenNBasisTransMatrix();
203  }

◆ 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 213 of file NodalTriExp.h.

215  {
216  GetCoord(lcoord, coord);
217  }
void GetCoord(const Array< OneD, const NekDouble > &Lcoords, Array< OneD, NekDouble > &coords)

References GetCoord().

◆ 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 205 of file NodalTriExp.h.

209  {
210  GetCoords(coords_0, coords_1, coords_2);
211  }
void GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3=NullNekDouble1DArray)

References GetCoords().

◆ v_GetLinStdExp()

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

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 473 of file NodalTriExp.cpp.

474 {
475  LibUtilities::BasisKey bkey0(m_base[0]->GetBasisType(), 2,
476  m_base[0]->GetPointsKey());
477  LibUtilities::BasisKey bkey1(m_base[1]->GetBasisType(), 2,
478  m_base[1]->GetPointsKey());
479 
481  bkey0, bkey1, m_nodalPointsKey.GetPointsType());
482 }
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
Definition: StdExpansion.h:163

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

◆ v_GetLocMatrix()

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

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 310 of file NodalTriExp.h.

311  {
312  return m_matrixManager[mkey];
313  }

References m_matrixManager.

◆ v_GetLocStaticCondMatrix()

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

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 324 of file NodalTriExp.h.

326  {
327  return m_staticCondMatrixManager[mkey];
328  }

References m_staticCondMatrixManager.

◆ v_GetNodalPoints()

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

Definition at line 219 of file NodalTriExp.h.

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

◆ v_GetStdExp()

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

◆ 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 383 of file NodalTriExp.h.

386  {
387  HelmholtzMatrixOp(inarray, outarray, mkey);
388  }
void HelmholtzMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:155

References HelmholtzMatrixOp().

◆ 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 227 of file NodalTriExp.h.

228  {
229  return Integral(inarray);
230  }
NekDouble Integral(const Array< OneD, const NekDouble > &inarray)
Integrate the physical point list inarray over region.
Definition: NodalTriExp.cpp:95

References Integral().

◆ 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 233 of file NodalTriExp.h.

235  {
236  IProductWRTBase(inarray, outarray);
237  }

References IProductWRTBase().

◆ 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 336 of file NodalTriExp.h.

339  {
340  boost::ignore_unused(multiplybyweights);
341  IProductWRTBase_SumFac(inarray, outarray);
342  }

References IProductWRTBase_SumFac().

◆ 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 239 of file NodalTriExp.h.

242  {
243  IProductWRTDerivBase(dir, inarray, outarray);
244  }
void IProductWRTDerivBase(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: NodalTriExp.h:92

References IProductWRTDerivBase().

◆ 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 344 of file NodalTriExp.h.

347  {
348  IProductWRTDerivBase_SumFac(dir, inarray, outarray);
349  }

References IProductWRTDerivBase_SumFac().

◆ 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 362 of file NodalTriExp.h.

365  {
366  LaplacianMatrixOp(inarray, outarray, mkey);
367  }
void LaplacianMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:130

References LaplacianMatrixOp().

◆ 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 369 of file NodalTriExp.h.

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

References LaplacianMatrixOp().

◆ 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 355 of file NodalTriExp.h.

358  {
359  MassMatrixOp(inarray, outarray, mkey);
360  }
void MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:123

References MassMatrixOp().

◆ 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 254 of file NodalTriExp.h.

258  {
259  boost::ignore_unused(out_d2);
260  PhysDeriv(inarray, out_d0, out_d1);
261  }
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.

References PhysDeriv().

◆ 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 263 of file NodalTriExp.h.

266  {
267  Array<OneD, NekDouble> tmp;
268  switch (dir)
269  {
270  case 0:
271  {
272  PhysDeriv(inarray, outarray, tmp);
273  }
274  break;
275  case 1:
276  {
277  PhysDeriv(inarray, tmp, outarray);
278  }
279  break;
280  default:
281  {
282  ASSERTL1(dir >= 0 && dir < 2, "input dir is out of range");
283  }
284  break;
285  }
286  }

References ASSERTL1, and PhysDeriv().

◆ 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 296 of file NodalTriExp.h.

300  {
301  return PhysEvaluate(coord, physvals);
302  }
NekDouble PhysEvaluate(const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals)

References PhysEvaluate().

◆ 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 246 of file NodalTriExp.h.

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

◆ 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 376 of file NodalTriExp.h.

379  {
380  WeakDerivMatrixOp(i, inarray, outarray, mkey);
381  }
void WeakDerivMatrixOp(const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:147

References WeakDerivMatrixOp().

◆ 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 147 of file NodalTriExp.h.

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

Referenced by v_WeakDerivMatrixOp().

Member Data Documentation

◆ m_matrixManager

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

◆ m_staticCondMatrixManager

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

Definition at line 193 of file NodalTriExp.h.

Referenced by v_GetLocStaticCondMatrix().