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

#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...
 
virtual ~NodalTriExp () override=default
 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)
 
NekDouble v_PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) final
 
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 ()=default
 
 StdNodalTriExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::PointsType Ntype)
 
 StdNodalTriExp (const StdNodalTriExp &T)
 
 ~StdNodalTriExp () override=default
 
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)
 
virtual ~StdTriExp () override
 
- 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 () override
 
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)
 
NekDouble BaryTensorDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs)
 
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 () const
 
std::shared_ptr< StdExpansionGetLinStdExp (void) const
 
int GetShapeDimension () const
 
bool IsBoundaryInteriorExpansion () const
 
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, 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. More...
 
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. 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 () override=default
 
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) override
 
- 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)
 
void DropLocMatrix (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)
 
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 override
 
virtual StdRegions::StdExpansionSharedPtr v_GetLinStdExp (void) const override
 
virtual DNekMatSharedPtr v_GenMatrix (const StdRegions::StdMatrixKey &mkey) override
 
- Protected Member Functions inherited from Nektar::StdRegions::StdNodalTriExp
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_FillMode (const int mode, Array< OneD, NekDouble > &outarray) override
 
int v_NumBndryCoeffs () 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
 
- 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) override
 
void v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 Backward tranform for triangular elements. More...
 
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) override
 
virtual void v_FwdTransBndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 
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) override
 
virtual void v_LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta) override
 
virtual void v_LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi) override
 
NekDouble v_PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode) final override
 
virtual int v_GetNverts () const final override
 
virtual int v_GetNtraces () const final override
 
virtual LibUtilities::ShapeType v_DetShapeType () const final override
 
virtual int v_NumDGBndryCoeffs () const override
 
virtual int v_GetTraceNcoeffs (const int i) const override
 
virtual int v_GetTraceIntNcoeffs (const int i) const override
 
virtual int v_GetTraceNumPoints (const int i) const override
 
virtual int v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset) override
 
virtual bool v_IsBoundaryInteriorExpansion () const override
 
virtual const LibUtilities::BasisKey v_GetTraceBasisKey (const int i, const int j) const override
 
virtual 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
 
virtual void v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey) override
 
virtual void v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 
virtual void v_MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 
virtual void v_GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true) override
 
- Protected Member Functions inherited from Nektar::StdRegions::StdExpansion2D
virtual NekDouble v_PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals) override
 
virtual void v_LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
virtual void v_HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
virtual 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. More...
 
virtual void v_GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat) override
 
- 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_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, 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. More...
 
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. More...
 
template<int DIR, bool DERIV = false, bool DERIV2 = false>
NekDouble BaryEvaluate (const NekDouble &coord, const NekDouble *physvals, NekDouble &deriv)
 
- Protected Member Functions inherited from Nektar::LocalRegions::Expansion2D
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) override
 
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) override
 
virtual void v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) override
 
virtual void v_AddRobinMassMatrix (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat) override
 
virtual void v_AddRobinTraceContribution (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, const Array< OneD, NekDouble > &incoeffs, Array< OneD, NekDouble > &coeffs) override
 
virtual DNekMatSharedPtr v_BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd) override
 
virtual void v_ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1) override
 
virtual void v_SetUpPhysNormals (const int edge) override
 
virtual NekDouble v_VectorFlux (const Array< OneD, Array< OneD, NekDouble >> &vec) override
 
virtual void v_TraceNormLen (const int traceid, NekDouble &h, NekDouble &p) override
 
- 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)
 
Array< OneD, NekDoubleGetMF (const int dir, const int shapedim, const StdRegions::VarCoeffMap &varcoeffs)
 
Array< OneD, NekDoubleGetMFDiv (const int dir, const StdRegions::VarCoeffMap &varcoeffs)
 
Array< OneD, NekDoubleGetMFMag (const int dir, const StdRegions::VarCoeffMap &varcoeffs)
 
virtual void v_MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 
virtual void v_DivideByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_ComputeLaplacianMetric ()
 
virtual int v_GetCoordim () const override
 
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) override
 
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 ()
 
virtual void v_SetPhysNormals (Array< OneD, const NekDouble > &normal)
 

Private Member Functions

virtual void v_GetCoords (Array< OneD, NekDouble > &coords_0, Array< OneD, NekDouble > &coords_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coords_2=NullNekDouble1DArray) override
 
virtual void v_GetCoord (const Array< OneD, const NekDouble > &lcoord, Array< OneD, NekDouble > &coord) override
 
virtual NekDouble v_Integral (const Array< OneD, const NekDouble > &inarray) override
 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) override
 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) override
 
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
 
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) override
 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) override
 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) override
 Virtual call to SegExp::FwdTrans. More...
 
virtual NekDouble v_PhysEvaluate (const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals) override
 Virtual call to TriExp::Evaluate. More...
 
virtual DNekMatSharedPtr v_CreateStdMatrix (const StdRegions::StdMatrixKey &mkey) override
 
virtual DNekScalMatSharedPtr v_GetLocMatrix (const MatrixKey &mkey) override
 
virtual void v_DropLocMatrix (const MatrixKey &mkey) override
 
virtual DNekScalBlkMatSharedPtr v_GetLocStaticCondMatrix (const MatrixKey &mkey) override
 
virtual void v_BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 
virtual void v_IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true) override
 
void v_IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
 
void v_AlignVectorToCollapsedDir (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray) override
 
virtual void v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
virtual void v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
virtual void v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
virtual void v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
virtual void v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override
 
void v_ComputeTraceNormal (const int edge) override
 

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:277
Expansion(SpatialDomains::GeometrySharedPtr pGeom)
Definition: Expansion.cpp:47
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: NodalTriExp.h:198
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: NodalTriExp.h:200
StdExpansion()
Default Constructor.
int getNumberOfCoefficients(int Na, int Nb)
Definition: ShapeType.hpp:114

◆ NodalTriExp() [2/2]

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

Copy Constructor.

Definition at line 64 of file NodalTriExp.cpp.

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

◆ ~NodalTriExp()

virtual Nektar::LocalRegions::NodalTriExp::~NodalTriExp ( )
overridevirtualdefault

Destructor.

Member Function Documentation

◆ CreateStdMatrix()

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

Definition at line 438 of file NodalTriExp.cpp.

440 {
441  LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
442  LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
446 
447  return tmp->GetStdMatrix(mkey);
448 }
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 371 of file NodalTriExp.cpp.

373 {
374  IProductWRTBase(inarray, outarray);
375 
376  // get Mass matrix inverse
377  MatrixKey masskey(StdRegions::eInvMass, DetShapeType(), *this,
381  DNekScalMatSharedPtr matsys = m_matrixManager[masskey];
382 
383  // copy inarray in case inarray == outarray
384  NekVector<NekDouble> in(m_ncoeffs, outarray, eCopy);
385  NekVector<NekDouble> out(m_ncoeffs, outarray, eWrapper);
386 
387  out = (*matsys) * in;
388 }
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:373
static ConstFactorMap NullConstFactorMap
Definition: StdRegions.hpp:400
static VarCoeffMap NullVarCoeffMap
Definition: StdRegions.hpp:344
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 390 of file NodalTriExp.cpp.

393 {
395 
396  if (inarray.get() == outarray.get())
397  {
398  Array<OneD, NekDouble> tmp(m_ncoeffs);
399  Vmath::Vcopy(m_ncoeffs, inarray.get(), 1, tmp.get(), 1);
400 
401  Blas::Dgemv('N', m_ncoeffs, m_ncoeffs, mat->Scale(),
402  (mat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
403  tmp.get(), 1, 0.0, outarray.get(), 1);
404  }
405  else
406  {
407  Blas::Dgemv('N', m_ncoeffs, m_ncoeffs, mat->Scale(),
408  (mat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
409  inarray.get(), 1, 0.0, outarray.get(), 1);
410  }
411 }
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 421 of file NodalTriExp.cpp.

423 {
424  int i;
425 
426  ASSERTL1(Lcoords[0] >= -1.0 && Lcoords[1] <= 1.0 && Lcoords[1] >= -1.0 &&
427  Lcoords[1] <= 1.0,
428  "Local coordinates are not in region [-1,1]");
429 
430  m_geom->FillGeom();
431 
432  for (i = 0; i < m_geom->GetCoordim(); ++i)
433  {
434  coords[i] = m_geom->GetCoord(i, Lcoords);
435  }
436 }
#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:275

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

416 {
417  Expansion::v_GetCoords(coords_0, coords_1, coords_2);
418 }
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3) override
Definition: Expansion.cpp:535

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

163  {
164  StdExpansion::HelmholtzMatrixOp_MatFree_GenericImpl(inarray, outarray,
165  mkey);
166  }

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

94 {
95  int nquad0 = m_base[0]->GetNumPoints();
96  int nquad1 = m_base[1]->GetNumPoints();
97  Array<OneD, const NekDouble> jac = m_metricinfo->GetJac(GetPointsKeys());
98  NekDouble ival;
99  Array<OneD, NekDouble> tmp(nquad0 * nquad1);
100 
101  // multiply inarray with Jacobian
102  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
103  {
104  Vmath::Vmul(nquad0 * nquad1, jac, 1, inarray, 1, tmp, 1);
105  }
106  else
107  {
108  Vmath::Smul(nquad0 * nquad1, jac[0], inarray, 1, tmp, 1);
109  }
110 
111  // call StdQuadExp version;
112  ival = StdNodalTriExp::v_Integral(tmp);
113  return ival;
114 }
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:276
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 145 of file NodalTriExp.cpp.

148 {
149  int nq = GetTotPoints();
150  MatrixKey iprodmatkey(StdRegions::eIProductWRTBase, DetShapeType(), *this);
151  DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
152 
153  Blas::Dgemv('N', m_ncoeffs, nq, iprodmat->Scale(),
154  (iprodmat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
155  inarray.get(), 1, 0.0, outarray.get(), 1);
156 }
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 116 of file NodalTriExp.cpp.

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

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

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

161 {
162  int nquad0 = m_base[0]->GetNumPoints();
163  int nquad1 = m_base[1]->GetNumPoints();
164  int nqtot = nquad0 * nquad1;
165  int wspsize = max(nqtot, m_ncoeffs);
166 
167  Array<OneD, NekDouble> tmp0(4 * wspsize);
168  Array<OneD, NekDouble> tmp1(tmp0 + wspsize);
169  Array<OneD, NekDouble> tmp2(tmp0 + 2 * wspsize);
170  Array<OneD, NekDouble> tmp3(tmp0 + 3 * wspsize);
171 
172  Array<OneD, Array<OneD, NekDouble>> tmp2D{2};
173  tmp2D[0] = tmp1;
174  tmp2D[1] = tmp2;
175 
176  AlignVectorToCollapsedDir(dir, inarray, tmp2D);
177 
178  MultiplyByQuadratureMetric(tmp1, tmp1);
179  MultiplyByQuadratureMetric(tmp2, tmp2);
180 
181  IProductWRTBase_SumFacKernel(m_base[0]->GetDbdata(), m_base[1]->GetBdata(),
182  tmp1, tmp3, tmp0);
183  IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(), m_base[1]->GetDbdata(),
184  tmp2, outarray, tmp0);
185  Vmath::Vadd(m_ncoeffs, tmp3, 1, outarray, 1, outarray, 1);
186 
187  NodalToModalTranspose(outarray, outarray);
188 }
void AlignVectorToCollapsedDir(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray)
Definition: Expansion.h:153
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 135 of file NodalTriExp.h.

138  {
139  StdExpansion::LaplacianMatrixOp_MatFree_GenericImpl(inarray, outarray,
140  mkey);
141  }

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

147  {
148  StdExpansion::LaplacianMatrixOp_MatFree(k1, k2, inarray, outarray,
149  mkey);
150  }

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

304 {
305  int nquad0 = m_base[0]->GetNumPoints();
306  int nquad1 = m_base[1]->GetNumPoints();
307  int nqtot = nquad0 * nquad1;
308  const Array<TwoD, const NekDouble> &df =
309  m_metricinfo->GetDerivFactors(GetPointsKeys());
310 
311  Array<OneD, NekDouble> diff0(2 * nqtot);
312  Array<OneD, NekDouble> diff1(diff0 + nqtot);
313 
314  StdNodalTriExp::v_PhysDeriv(inarray, diff0, diff1);
315 
316  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
317  {
318  if (out_d0.size())
319  {
320  Vmath::Vmul(nqtot, df[0], 1, diff0, 1, out_d0, 1);
321  Vmath::Vvtvp(nqtot, df[1], 1, diff1, 1, out_d0, 1, out_d0, 1);
322  }
323 
324  if (out_d1.size())
325  {
326  Vmath::Vmul(nqtot, df[2], 1, diff0, 1, out_d1, 1);
327  Vmath::Vvtvp(nqtot, df[3], 1, diff1, 1, out_d1, 1, out_d1, 1);
328  }
329 
330  if (out_d2.size())
331  {
332  Vmath::Vmul(nqtot, df[4], 1, diff0, 1, out_d2, 1);
333  Vmath::Vvtvp(nqtot, df[5], 1, diff1, 1, out_d2, 1, out_d2, 1);
334  }
335  }
336  else // regular geometry
337  {
338  if (out_d0.size())
339  {
340  Vmath::Smul(nqtot, df[0][0], diff0, 1, out_d0, 1);
341  Blas::Daxpy(nqtot, df[1][0], diff1, 1, out_d0, 1);
342  }
343 
344  if (out_d1.size())
345  {
346  Vmath::Smul(nqtot, df[2][0], diff0, 1, out_d1, 1);
347  Blas::Daxpy(nqtot, df[3][0], diff1, 1, out_d1, 1);
348  }
349 
350  if (out_d2.size())
351  {
352  Vmath::Smul(nqtot, df[4][0], diff0, 1, out_d2, 1);
353  Blas::Daxpy(nqtot, df[5][0], diff1, 1, out_d2, 1);
354  }
355  }
356 }
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 450 of file NodalTriExp.cpp.

454 {
455  Array<OneD, NekDouble> Lcoord = Array<OneD, NekDouble>(2);
456 
457  ASSERTL0(m_geom, "m_geom not defined");
458  m_geom->GetLocCoords(coord, Lcoord);
459 
460  return StdExpansion2D::v_PhysEvaluate(Lcoord, physvals);
461 }
#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 
)
overrideprivatevirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 190 of file NodalTriExp.cpp.

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

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 317 of file NodalTriExp.h.

319  {
320  StdNodalTriExp::v_BwdTrans_SumFac(inarray, outarray);
321  }

◆ v_ComputeTraceNormal()

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

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 514 of file NodalTriExp.cpp.

515 {
516  int i;
517  const SpatialDomains::GeomFactorsSharedPtr &geomFactors =
518  GetGeom()->GetMetricInfo();
519  const SpatialDomains::GeomType type = geomFactors->GetGtype();
520 
522  const Array<TwoD, const NekDouble> &df =
523  geomFactors->GetDerivFactors(ptsKeys);
524  const Array<OneD, const NekDouble> &jac = geomFactors->GetJac(ptsKeys);
525  int nqe = m_base[0]->GetNumPoints();
526  int dim = GetCoordim();
527 
528  m_traceNormals[edge] = Array<OneD, Array<OneD, NekDouble>>(dim);
529  Array<OneD, Array<OneD, NekDouble>> &normal = m_traceNormals[edge];
530  for (i = 0; i < dim; ++i)
531  {
532  normal[i] = Array<OneD, NekDouble>(nqe);
533  }
534 
535  size_t nqb = nqe;
536  size_t nbnd = edge;
537  m_elmtBndNormDirElmtLen[nbnd] = Array<OneD, NekDouble>{nqb, 0.0};
538  Array<OneD, NekDouble> &length = m_elmtBndNormDirElmtLen[nbnd];
539 
540  // Regular geometry case
541  if ((type == SpatialDomains::eRegular) ||
543  {
544  NekDouble fac;
545  // Set up normals
546  switch (edge)
547  {
548  case 0:
549  for (i = 0; i < GetCoordim(); ++i)
550  {
551  Vmath::Fill(nqe, -df[2 * i + 1][0], normal[i], 1);
552  }
553  break;
554  case 1:
555  for (i = 0; i < GetCoordim(); ++i)
556  {
557  Vmath::Fill(nqe, df[2 * i + 1][0] + df[2 * i][0], normal[i],
558  1);
559  }
560  break;
561  case 2:
562  for (i = 0; i < GetCoordim(); ++i)
563  {
564  Vmath::Fill(nqe, -df[2 * i][0], normal[i], 1);
565  }
566  break;
567  default:
568  ASSERTL0(false, "Edge is out of range (edge < 3)");
569  }
570 
571  // normalise
572  fac = 0.0;
573  for (i = 0; i < GetCoordim(); ++i)
574  {
575  fac += normal[i][0] * normal[i][0];
576  }
577  fac = 1.0 / sqrt(fac);
578 
579  Vmath::Fill(nqb, fac, length, 1);
580 
581  for (i = 0; i < GetCoordim(); ++i)
582  {
583  Vmath::Smul(nqe, fac, normal[i], 1, normal[i], 1);
584  }
585  }
586  else // Set up deformed normals
587  {
588  int j;
589 
590  int nquad0 = ptsKeys[0].GetNumPoints();
591  int nquad1 = ptsKeys[1].GetNumPoints();
592 
593  LibUtilities::PointsKey from_key;
594 
595  Array<OneD, NekDouble> normals(GetCoordim() * max(nquad0, nquad1), 0.0);
596  Array<OneD, NekDouble> edgejac(GetCoordim() * max(nquad0, nquad1), 0.0);
597 
598  // Extract Jacobian along edges and recover local
599  // derivates (dx/dr) for polynomial interpolation by
600  // multiplying m_gmat by jacobian
601  switch (edge)
602  {
603  case 0:
604  for (j = 0; j < nquad0; ++j)
605  {
606  edgejac[j] = jac[j];
607  for (i = 0; i < GetCoordim(); ++i)
608  {
609  normals[i * nquad0 + j] =
610  -df[2 * i + 1][j] * edgejac[j];
611  }
612  }
613  from_key = ptsKeys[0];
614  break;
615  case 1:
616  for (j = 0; j < nquad1; ++j)
617  {
618  edgejac[j] = jac[nquad0 * j + nquad0 - 1];
619  for (i = 0; i < GetCoordim(); ++i)
620  {
621  normals[i * nquad1 + j] =
622  (df[2 * i][nquad0 * j + nquad0 - 1] +
623  df[2 * i + 1][nquad0 * j + nquad0 - 1]) *
624  edgejac[j];
625  }
626  }
627  from_key = ptsKeys[1];
628  break;
629  case 2:
630  for (j = 0; j < nquad1; ++j)
631  {
632  edgejac[j] = jac[nquad0 * j];
633  for (i = 0; i < GetCoordim(); ++i)
634  {
635  normals[i * nquad1 + j] =
636  -df[2 * i][nquad0 * j] * edgejac[j];
637  }
638  }
639  from_key = ptsKeys[1];
640  break;
641  default:
642  ASSERTL0(false, "edge is out of range (edge < 3)");
643  }
644 
645  int nq = from_key.GetNumPoints();
646  Array<OneD, NekDouble> work(nqe, 0.0);
647 
648  // interpolate Jacobian and invert
649  LibUtilities::Interp1D(from_key, jac, m_base[0]->GetPointsKey(), work);
650  Vmath::Sdiv(nq, 1.0, &work[0], 1, &work[0], 1);
651 
652  // interpolate
653  for (i = 0; i < GetCoordim(); ++i)
654  {
655  LibUtilities::Interp1D(from_key, &normals[i * nq],
656  m_base[0]->GetPointsKey(), &normal[i][0]);
657  Vmath::Vmul(nqe, work, 1, normal[i], 1, normal[i], 1);
658  }
659 
660  // normalise normal vectors
661  Vmath::Zero(nqe, work, 1);
662  for (i = 0; i < GetCoordim(); ++i)
663  {
664  Vmath::Vvtvp(nqe, normal[i], 1, normal[i], 1, work, 1, work, 1);
665  }
666 
667  Vmath::Vsqrt(nqe, work, 1, work, 1);
668  Vmath::Sdiv(nqe, 1.0, work, 1, work, 1);
669 
670  Vmath::Vcopy(nqb, work, 1, length, 1);
671 
672  for (i = 0; i < GetCoordim(); ++i)
673  {
674  Vmath::Vmul(nqe, normal[i], 1, work, 1, normal[i], 1);
675  }
676 
677  // Reverse direction so that points are in
678  // anticlockwise direction if edge >=2
679  if (edge >= 2)
680  {
681  for (i = 0; i < GetCoordim(); ++i)
682  {
683  Vmath::Reverse(nqe, normal[i], 1, normal[i], 1);
684  }
685  }
686  }
687 }
std::map< int, NormalVector > m_traceNormals
Definition: Expansion.h:278
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:288
SpatialDomains::GeometrySharedPtr GetGeom() const
Definition: Expansion.cpp:171
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:294

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)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 295 of file NodalTriExp.h.

297  {
298  return CreateStdMatrix(mkey);
299  }
DNekMatSharedPtr CreateStdMatrix(const StdRegions::StdMatrixKey &mkey)

References CreateStdMatrix().

◆ v_DropLocMatrix()

virtual void Nektar::LocalRegions::NodalTriExp::v_DropLocMatrix ( const MatrixKey mkey)
inlineoverrideprivatevirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 306 of file NodalTriExp.h.

307  {
308  m_matrixManager.DeleteObject(mkey);
309  }

References m_matrixManager.

◆ v_FwdTrans()

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

Virtual call to SegExp::FwdTrans.

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 281 of file NodalTriExp.h.

283  {
284  FwdTrans(inarray, outarray);
285  }
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)
overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 493 of file NodalTriExp.cpp.

494 {
495  DNekMatSharedPtr returnval;
496 
497  switch (mkey.GetMatrixType())
498  {
505  returnval = Expansion2D::v_GenMatrix(mkey);
506  break;
507  default:
508  returnval = StdNodalTriExp::v_GenMatrix(mkey);
509  break;
510  }
511  return returnval;
512 }
virtual DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey) override
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_GetCoord()

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

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 210 of file NodalTriExp.h.

212  {
213  GetCoord(lcoord, coord);
214  }
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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 202 of file NodalTriExp.h.

206  {
207  GetCoords(coords_0, coords_1, coords_2);
208  }
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
overrideprotectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 482 of file NodalTriExp.cpp.

483 {
484  LibUtilities::BasisKey bkey0(m_base[0]->GetBasisType(), 2,
485  m_base[0]->GetPointsKey());
486  LibUtilities::BasisKey bkey1(m_base[1]->GetBasisType(), 2,
487  m_base[1]->GetPointsKey());
488 
490  bkey0, bkey1, m_nodalPointsKey.GetPointsType());
491 }
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
Definition: StdExpansion.h:162

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)
inlineoverrideprivatevirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 301 of file NodalTriExp.h.

302  {
303  return m_matrixManager[mkey];
304  }

References m_matrixManager.

◆ v_GetLocStaticCondMatrix()

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

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 311 of file NodalTriExp.h.

313  {
314  return m_staticCondMatrixManager[mkey];
315  }

References m_staticCondMatrixManager.

◆ v_GetStdExp()

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

◆ v_HelmholtzMatrixOp()

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

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 374 of file NodalTriExp.h.

378  {
379  HelmholtzMatrixOp(inarray, outarray, mkey);
380  }
void HelmholtzMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:160

References HelmholtzMatrixOp().

◆ v_Integral()

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

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

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 218 of file NodalTriExp.h.

220  {
221  return Integral(inarray);
222  }
NekDouble Integral(const Array< OneD, const NekDouble > &inarray)
Integrate the physical point list inarray over region.
Definition: NodalTriExp.cpp:93

References Integral().

◆ v_IProductWRTBase()

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

Virtual call to TriExp::IProduct_WRT_B.

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 225 of file NodalTriExp.h.

227  {
228  IProductWRTBase(inarray, outarray);
229  }

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 323 of file NodalTriExp.h.

327  {
328  boost::ignore_unused(multiplybyweights);
329  IProductWRTBase_SumFac(inarray, outarray);
330  }

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 231 of file NodalTriExp.h.

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

References IProductWRTDerivBase().

◆ v_IProductWRTDerivBase_SumFac()

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

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 332 of file NodalTriExp.h.

335  {
336  IProductWRTDerivBase_SumFac(dir, inarray, outarray);
337  }

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 350 of file NodalTriExp.h.

354  {
355  LaplacianMatrixOp(inarray, outarray, mkey);
356  }
void LaplacianMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:135

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 358 of file NodalTriExp.h.

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

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 343 of file NodalTriExp.h.

346  {
347  MassMatrixOp(inarray, outarray, mkey);
348  }
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]

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 
)
inlineoverrideprivatevirtual

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

250  {
251  boost::ignore_unused(out_d2);
252  PhysDeriv(inarray, out_d0, out_d1);
253  }
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 
)
inlineoverrideprivatevirtual

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

See also
StdRegions::StdExpansion::PhysDeriv

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 255 of file NodalTriExp.h.

258  {
259  Array<OneD, NekDouble> tmp;
260  switch (dir)
261  {
262  case 0:
263  {
264  PhysDeriv(inarray, outarray, tmp);
265  }
266  break;
267  case 1:
268  {
269  PhysDeriv(inarray, tmp, outarray);
270  }
271  break;
272  default:
273  {
274  ASSERTL1(dir >= 0 && dir < 2, "input dir is out of range");
275  }
276  break;
277  }
278  }

References ASSERTL1, and PhysDeriv().

◆ v_PhysEvaluate() [1/2]

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

Virtual call to TriExp::Evaluate.

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 288 of file NodalTriExp.h.

291  {
292  return PhysEvaluate(coord, physvals);
293  }
NekDouble PhysEvaluate(const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals)

References PhysEvaluate().

◆ v_PhysEvaluate() [2/2]

NekDouble Nektar::LocalRegions::NodalTriExp::v_PhysEvaluate ( const Array< OneD, NekDouble > &  coord,
const Array< OneD, const NekDouble > &  inarray,
std::array< NekDouble, 3 > &  firstOrderDerivs 
)
finalvirtual

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 463 of file NodalTriExp.cpp.

467 {
468  Array<OneD, NekDouble> Lcoord(2);
469  ASSERTL0(m_geom, "m_geom not defined");
470  m_geom->GetLocCoords(coord, Lcoord);
471  return StdExpansion2D::v_PhysEvaluate(Lcoord, inarray, firstOrderDerivs);
472 }

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

◆ v_StdPhysDeriv()

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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdTriExp.

Definition at line 238 of file NodalTriExp.h.

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

◆ 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 
)
inlineoverrideprivatevirtual

Reimplemented from Nektar::StdRegions::StdNodalTriExp.

Definition at line 366 of file NodalTriExp.h.

370  {
371  WeakDerivMatrixOp(i, inarray, outarray, mkey);
372  }
void WeakDerivMatrixOp(const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
Definition: NodalTriExp.h:152

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

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

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

Referenced by v_GetLocStaticCondMatrix().