Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Protected Member Functions | Private Member Functions | Private Attributes | List of all members
Nektar::LocalRegions::QuadExp Class Reference

#include <QuadExp.h>

Inheritance diagram for Nektar::LocalRegions::QuadExp:
Inheritance graph
[legend]
Collaboration diagram for Nektar::LocalRegions::QuadExp:
Collaboration graph
[legend]

Public Member Functions

 QuadExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const SpatialDomains::QuadGeomSharedPtr &geom)
 Constructor using BasisKey class for quadrature points and order definition. More...
 
 QuadExp (const QuadExp &T)
 
virtual ~QuadExp ()
 
- Public Member Functions inherited from Nektar::StdRegions::StdQuadExp
 StdQuadExp ()
 
 StdQuadExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)
 Constructor using BasisKey class for quadrature points and order definition. More...
 
 StdQuadExp (const StdQuadExp &T)
 Copy Constructor. More...
 
 ~StdQuadExp ()
 Destructor. More...
 
- Public Member Functions inherited from Nektar::StdRegions::StdExpansion2D
 StdExpansion2D ()
 
 StdExpansion2D (int numcoeffs, const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)
 
 StdExpansion2D (const StdExpansion2D &T)
 
virtual ~StdExpansion2D ()
 
void PhysTensorDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray_d0, Array< OneD, NekDouble > &outarray_d1)
 Calculate the 2D derivative in the local tensor/collapsed coordinate at the physical points. More...
 
NekDouble Integral (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &w0, const Array< OneD, const NekDouble > &w1)
 
void BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true)
 
void IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true)
 
- Public Member Functions inherited from Nektar::StdRegions::StdExpansion
 StdExpansion ()
 Default Constructor. More...
 
 StdExpansion (const int numcoeffs, const int numbases, const LibUtilities::BasisKey &Ba=LibUtilities::NullBasisKey, const LibUtilities::BasisKey &Bb=LibUtilities::NullBasisKey, const LibUtilities::BasisKey &Bc=LibUtilities::NullBasisKey)
 Constructor. More...
 
 StdExpansion (const StdExpansion &T)
 Copy Constructor. More...
 
virtual ~StdExpansion ()
 Destructor. More...
 
int GetNumBases () const
 This function returns the number of 1D bases used in the expansion. More...
 
const Array< OneD, const
LibUtilities::BasisSharedPtr > & 		GetBase () const
 This function gets the shared point to basis. More...
 
const
LibUtilities::BasisSharedPtr & 		GetBasis (int dir) const
 This function gets the shared point to basis in the dir direction. More...
 
int GetNcoeffs (void) const
 This function returns the total number of coefficients used in the expansion. More...
 
int GetTotPoints () const
 This function returns the total number of quadrature points used in the element. More...
 
LibUtilities::BasisType GetBasisType (const int dir) const
 This function returns the type of basis used in the dir direction. More...
 
int GetBasisNumModes (const int dir) const
 This function returns the number of expansion modes in the dir direction. More...
 
int EvalBasisNumModesMax (void) const
 This function returns the maximum number of expansion modes over all local directions. More...
 
LibUtilities::PointsType GetPointsType (const int dir) const
 This function returns the type of quadrature points used in the dir direction. More...
 
int GetNumPoints (const int dir) const
 This function returns the number of quadrature points in the dir direction. More...
 
const Array< OneD, const
NekDouble > & 		GetPoints (const int dir) const
 This function returns a pointer to the array containing the quadrature points in dir direction. More...
 
int GetNverts () const
 This function returns the number of vertices of the expansion domain. More...
 
int GetNedges () const
 This function returns the number of edges of the expansion domain. More...
 
int GetEdgeNcoeffs (const int i) const
 This function returns the number of expansion coefficients belonging to the i-th edge. More...
 
int GetTotalEdgeIntNcoeffs () const
 
int GetEdgeNumPoints (const int i) const
 This function returns the number of quadrature points belonging to the i-th edge. More...
 
int DetCartesianDirOfEdge (const int edge)
 
const LibUtilities::BasisKey DetEdgeBasisKey (const int i) const
 
const LibUtilities::BasisKey DetFaceBasisKey (const int i, const int k) const
 
int GetFaceNumPoints (const int i) const
 This function returns the number of quadrature points belonging to the i-th face. More...
 
int GetFaceNcoeffs (const int i) const
 This function returns the number of expansion coefficients belonging to the i-th face. More...
 
int GetFaceIntNcoeffs (const int i) const
 
int GetTotalFaceIntNcoeffs () const
 
int GetTraceNcoeffs (const int i) const
 This function returns the number of expansion coefficients belonging to the i-th edge/face. More...
 
LibUtilities::PointsKey GetFacePointsKey (const int i, const int j) const
 
int NumBndryCoeffs (void) const
 
int NumDGBndryCoeffs (void) const
 
LibUtilities::BasisType GetEdgeBasisType (const int i) const
 This function returns the type of expansion basis on the i-th edge. More...
 
const LibUtilities::PointsKey GetNodalPointsKey () const
 This function returns the type of expansion Nodal point type if defined. More...
 
int GetNfaces () const
 This function returns the number of faces of the expansion domain. More...
 
int GetNtrace () const
 Returns the number of trace elements connected to this element. More...
 
LibUtilities::ShapeType DetShapeType () const
 This function returns the shape of the expansion domain. More...
 
boost::shared_ptr< StdExpansion > GetStdExp (void) const
 
int GetShapeDimension () const
 
bool IsBoundaryInteriorExpansion ()
 
bool IsNodalNonTensorialExp ()
 
void BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function performs the Backward transformation from coefficient space to physical space. More...
 
void FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function performs the Forward transformation from physical space to coefficient space. More...
 
void FwdTrans_BndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
NekDouble Integral (const Array< OneD, const NekDouble > &inarray)
 This function integrates the specified function over the domain. More...
 
void FillMode (const int mode, Array< OneD, NekDouble > &outarray)
 This function fills the array outarray with the mode-th mode of the expansion. More...
 
void IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 this function calculates the inner product of a given function f with the different modes of the expansion More...
 
void IProductWRTBase (const Array< OneD, const NekDouble > &base, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int coll_check)
 
void IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
int GetElmtId ()
 Get the element id of this expansion when used in a list by returning value of m_elmt_id. More...
 
void SetElmtId (const int id)
 Set the element id of this expansion when used in a list by returning value of m_elmt_id. More...
 
void GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2=NullNekDouble1DArray, Array< OneD, NekDouble > &coords_3=NullNekDouble1DArray)
 this function returns the physical coordinates of the quadrature points of the expansion More...
 
void GetCoord (const Array< OneD, const NekDouble > &Lcoord, Array< OneD, NekDouble > &coord)
 given the coordinates of a point of the element in the local collapsed coordinate system, this function calculates the physical coordinates of the point More...
 
DNekMatSharedPtr GetStdMatrix (const StdMatrixKey &mkey)
 
DNekBlkMatSharedPtr GetStdStaticCondMatrix (const StdMatrixKey &mkey)
 
IndexMapValuesSharedPtr GetIndexMap (const IndexMapKey &ikey)
 
const Array< OneD, const
NekDouble > & 		GetPhysNormals (void)
 
void SetPhysNormals (Array< OneD, const NekDouble > &normal)
 
virtual void SetUpPhysNormals (const int edge)
 
void NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)
 
void NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
 
void NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, const Array< OneD, const NekDouble > &Fz, Array< OneD, NekDouble > &outarray)
 
void NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray)
 
DNekScalBlkMatSharedPtr GetLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)
 
void DropLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey)
 
StdRegions::Orientation GetForient (int face)
 
StdRegions::Orientation GetEorient (int edge)
 
StdRegions::Orientation GetPorient (int point)
 
StdRegions::Orientation GetCartesianEorient (int edge)
 
void SetCoeffsToOrientation (Array< OneD, NekDouble > &coeffs, StdRegions::Orientation dir)
 
void SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
int CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset)
 
void ExtractDataToCoeffs (const NekDouble *data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble *coeffs)
 
NekDouble StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals)
 
int GetCoordim ()
 
void GetBoundaryMap (Array< OneD, unsigned int > &outarray)
 
void GetInteriorMap (Array< OneD, unsigned int > &outarray)
 
int GetVertexMap (const int localVertexId, bool useCoeffPacking=false)
 
void GetEdgeInteriorMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)
 
void GetFaceInteriorMap (const int fid, const Orientation faceOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)
 
void GetEdgeToElementMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P=-1)
 
void GetFaceToElementMap (const int fid, const Orientation faceOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int nummodesA=-1, int nummodesB=-1)
 
void GetEdgePhysVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 Extract the physical values along edge edge from inarray into outarray following the local edge orientation and point distribution defined by defined in EdgeExp. More...
 
void GetEdgePhysVals (const int edge, const boost::shared_ptr< StdExpansion > &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void GetTracePhysVals (const int edge, const boost::shared_ptr< StdExpansion > &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void GetVertexPhysVals (const int vertex, const Array< OneD, const NekDouble > &inarray, NekDouble &outarray)
 
void GetEdgeInterpVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void GetEdgeQFactors (const int edge, Array< OneD, NekDouble > &outarray)
 Extract the metric factors to compute the contravariant fluxes along edge edge and stores them into outarray following the local edge orientation (i.e. anticlockwise convention). More...
 
void GetFacePhysVals (const int face, const boost::shared_ptr< StdExpansion > &FaceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient=eNoOrientation)
 
void GetEdgePhysMap (const int edge, Array< OneD, int > &outarray)
 
void GetFacePhysMap (const int face, Array< OneD, int > &outarray)
 
void MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
DNekMatSharedPtr CreateGeneralMatrix (const StdMatrixKey &mkey)
 this function generates the mass matrix $\mathbf{M}[i][j] = \int \phi_i(\mathbf{x}) \phi_j(\mathbf{x}) d\mathbf{x}$ More...
 
void GeneralMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey)
 
void LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void WeakDirectionalDerivMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void MassLevelCurvatureMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void LinearAdvectionDiffusionReactionMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true)
 
void HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
DNekMatSharedPtr GenMatrix (const StdMatrixKey &mkey)
 
void PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
 
void PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void PhysDeriv_s (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_ds)
 
void PhysDeriv_n (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_dn)
 
void PhysDirectionalDeriv (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &direction, Array< OneD, NekDouble > &outarray)
 
void StdPhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
 
void StdPhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void AddRobinMassMatrix (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)
 
void AddRobinEdgeContribution (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, Array< OneD, NekDouble > &coeffs)
 
NekDouble PhysEvaluate (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals)
 This function evaluates the expansion at a single (arbitrary) point of the domain. More...
 
NekDouble PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals)
 This function evaluates the expansion at a single (arbitrary) point of the domain. More...
 
void LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
 Convert local cartesian coordinate xi into local collapsed coordinates eta. More...
 
const boost::shared_ptr
< SpatialDomains::GeomFactors > & 		GetMetricInfo (void) const
 
virtual int v_GetElmtId ()
 Get the element id of this expansion when used in a list by returning value of m_elmt_id. More...
 
virtual const Array< OneD,
const NekDouble > & 		v_GetPhysNormals (void)
 
virtual void v_SetPhysNormals (Array< OneD, const NekDouble > &normal)
 
virtual void v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)
 
virtual void v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
 
virtual StdRegions::Orientation v_GetForient (int face)
 
virtual StdRegions::Orientation v_GetPorient (int point)
 
NekDouble Linf (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray)
 Function to evaluate the discrete $ L_\infty$ error $ |\epsilon|_\infty = \max |u - u_{exact}|$ where $ u_{exact}$ is given by the array sol. More...
 
NekDouble L2 (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray)
 Function to evaluate the discrete $ L_2$ error, $ | \epsilon |_{2} = \left [ \int^1_{-1} [u - u_{exact}]^2 dx \right]^{1/2} d\xi_1 $ where $ u_{exact}$ is given by the array sol. More...
 
NekDouble H1 (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray)
 Function to evaluate the discrete $ H^1$ error, $ | \epsilon |^1_{2} = \left [ \int^1_{-1} [u - u_{exact}]^2 + \nabla(u - u_{exact})\cdot\nabla(u - u_{exact})\cdot dx \right]^{1/2} d\xi_1 $ where $ u_{exact}$ is given by the array sol. More...
 
const NormalVector & GetEdgeNormal (const int edge) const
 
void ComputeEdgeNormal (const int edge)
 
void NegateEdgeNormal (const int edge)
 
bool EdgeNormalNegated (const int edge)
 
void ComputeFaceNormal (const int face)
 
void NegateFaceNormal (const int face)
 
bool FaceNormalNegated (const int face)
 
void ComputeVertexNormal (const int vertex)
 
const NormalVector & GetFaceNormal (const int face) const
 
const NormalVector & GetVertexNormal (const int vertex) const
 
const NormalVector & GetSurfaceNormal (const int id) const
 
const LibUtilities::PointsKeyVector GetPointsKeys () const
 
Array< OneD, unsigned int > GetEdgeInverseBoundaryMap (int eid)
 
Array< OneD, unsigned int > GetFaceInverseBoundaryMap (int fid, StdRegions::Orientation faceOrient=eNoOrientation)
 
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 >
boost::shared_ptr< T > as ()
 
void IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
 
- Public Member Functions inherited from Nektar::LocalRegions::Expansion2D
 Expansion2D (SpatialDomains::Geometry2DSharedPtr pGeom)
 
virtual ~Expansion2D ()
 
void SetTraceToGeomOrientation (Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &inout)
 
Expansion1DSharedPtr GetEdgeExp (int edge, bool SetUpNormal=true)
 
void SetEdgeExp (const int edge, Expansion1DSharedPtr &e)
 
void AddNormTraceInt (const int dir, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &edgeCoeffs, Array< OneD, NekDouble > &outarray)
 
void AddNormTraceInt (const int dir, Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs)
 
void AddEdgeBoundaryInt (const int edge, ExpansionSharedPtr &EdgeExp, Array< OneD, NekDouble > &edgePhys, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap)
 
void AddHDGHelmholtzEdgeTerms (const NekDouble tau, const int edge, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &edgePhys, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray)
 
void AddHDGHelmholtzTraceTerms (const NekDouble tau, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray)
 
Expansion3DSharedPtr GetLeftAdjacentElementExp () const
 
Expansion3DSharedPtr GetRightAdjacentElementExp () const
 
int GetLeftAdjacentElementFace () const
 
int GetRightAdjacentElementFace () const
 
void SetAdjacentElementExp (int face, Expansion3DSharedPtr &f)
 
SpatialDomains::Geometry2DSharedPtr GetGeom2D () const
 
void ReOrientEdgePhysMap (const int nvert, const StdRegions::Orientation orient, const int nq0, Array< OneD, int > &idmap)
 
- Public Member Functions inherited from Nektar::LocalRegions::Expansion
 Expansion (SpatialDomains::GeometrySharedPtr pGeom)
 
 Expansion (const Expansion &pSrc)
 
virtual ~Expansion ()
 
DNekScalMatSharedPtr GetLocMatrix (const LocalRegions::MatrixKey &mkey)
 
DNekScalMatSharedPtr GetLocMatrix (const StdRegions::MatrixType mtype, const StdRegions::ConstFactorMap &factors=StdRegions::NullConstFactorMap, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap)
 
SpatialDomains::GeometrySharedPtr GetGeom () const
 
void Reset ()
 
DNekMatSharedPtr BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)
 
DNekMatSharedPtr BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd)
 
void AddEdgeNormBoundaryInt (const int edge, const boost::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 boost::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 
void AddFaceNormBoundaryInt (const int face, const boost::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)
 

Protected Member Functions

virtual NekDouble v_Integral (const Array< OneD, const NekDouble > &inarray)
 Integrates the specified function over the domain. More...
 
virtual NekDouble v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals)
 
virtual void v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
 Calculate the derivative of the physical points. More...
 
virtual void v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 Calculate the derivative of the physical points in a given direction. More...
 
virtual void v_PhysDirectionalDeriv (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &direction, Array< OneD, NekDouble > &out)
 Physical derivative along a direction vector. More...
 
virtual void v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 Transform a given function from physical quadrature space to coefficient space. More...
 
virtual void v_FwdTrans_BndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray. More...
 
virtual void v_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
 
virtual void v_IProductWRTBase_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_IProductWRTDerivBase_MatOp (const int dir, const Array< OneD, const NekDouble > &inarray, 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
StdRegions::StdExpansionSharedPtr 		v_GetStdExp (void) const
 
virtual void v_GetCoord (const Array< OneD, const NekDouble > &Lcoords, Array< OneD, NekDouble > &coords)
 
virtual void v_GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3)
 
virtual NekDouble v_PhysEvaluate (const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals)
 This function evaluates the expansion at a single (arbitrary) point of the domain. More...
 
virtual void v_GetEdgePhysVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 Extract the physical values along edge edge from inarray into outarray following the local edge orientation and point distribution defined by defined in EdgeExp. More...
 
virtual void v_GetEdgePhysVals (const int edge, const StdRegions::StdExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_GetTracePhysVals (const int edge, const StdRegions::StdExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient)
 
virtual void v_GetEdgeInterpVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_GetEdgeQFactors (const int edge, Array< OneD, NekDouble > &outarray)
 
virtual void v_ComputeEdgeNormal (const int edge)
 
virtual const
SpatialDomains::GeomFactorsSharedPtr & 		v_GetMetricInfo () const
 
virtual int v_GetCoordim ()
 
virtual void v_ExtractDataToCoeffs (const NekDouble *data, const std::vector< unsigned int > &nummodes, const int mode_offset, NekDouble *coeffs)
 Unpack data from input file assuming it comes from the same expansion type. More...
 
virtual StdRegions::Orientation v_GetEorient (int edge)
 
virtual StdRegions::Orientation v_GetCartesianEorient (int edge)
 
virtual const
LibUtilities::BasisSharedPtr & 		v_GetBasis (int dir) const
 
virtual int v_GetNumPoints (const int dir) const
 
virtual void v_GetEdgePhysMap (const int edge, Array< OneD, int > &outarray)
 
virtual DNekMatSharedPtr v_GenMatrix (const StdRegions::StdMatrixKey &mkey)
 
virtual DNekMatSharedPtr v_CreateStdMatrix (const StdRegions::StdMatrixKey &mkey)
 
DNekScalMatSharedPtr CreateMatrix (const MatrixKey &mkey)
 
DNekScalBlkMatSharedPtr CreateStaticCondMatrix (const MatrixKey &mkey)
 
virtual DNekScalMatSharedPtr v_GetLocMatrix (const MatrixKey &mkey)
 
virtual DNekScalBlkMatSharedPtr v_GetLocStaticCondMatrix (const MatrixKey &mkey)
 
void v_DropLocStaticCondMatrix (const MatrixKey &mkey)
 
virtual void v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_WeakDirectionalDerivMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_MassLevelCurvatureMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_GeneralMatrixOp_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp)
 
virtual void v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_ComputeLaplacianMetric ()
 
virtual void v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdRegions::StdMatrixKey &mkey)
 
- Protected Member Functions inherited from Nektar::StdRegions::StdQuadExp
virtual void v_StdPhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
 
virtual void v_StdPhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1)
 
virtual void v_IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1)
 
virtual void v_LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta)
 
virtual void v_FillMode (const int mode, Array< OneD, NekDouble > &array)
 Fill outarray with mode mode of expansion. More...
 
virtual int v_GetNverts () const
 
virtual int v_GetNedges () const
 
virtual int v_GetEdgeNcoeffs (const int i) const
 
virtual int v_GetEdgeNumPoints (const int i) const
 
virtual int v_NumBndryCoeffs () const
 
virtual int v_NumDGBndryCoeffs () const
 
virtual LibUtilities::BasisType v_GetEdgeBasisType (const int i) const
 
virtual int v_DetCartesianDirOfEdge (const int edge)
 
virtual int v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset)
 
virtual const
LibUtilities::BasisKey 		v_DetEdgeBasisKey (const int i) const
 
virtual LibUtilities::ShapeType v_DetShapeType () const
 
virtual bool v_IsBoundaryInteriorExpansion ()
 
virtual void v_GetBoundaryMap (Array< OneD, unsigned int > &outarray)
 
virtual void v_GetInteriorMap (Array< OneD, unsigned int > &outarray)
 
virtual int v_GetVertexMap (int localVertexId, bool useCoeffPacking=false)
 
virtual void v_GetEdgeInteriorMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray)
 
void v_GetEdgeToElementMap (const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P=-1)
 
void v_GeneralMatrixOp_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void v_MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true)
 
- Protected Member Functions inherited from Nektar::StdRegions::StdExpansion2D
virtual NekDouble v_PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals)
 
virtual void v_LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual void v_HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey)
 
virtual int v_GetTraceNcoeffs (const int i) const
 
- Protected Member Functions inherited from Nektar::StdRegions::StdExpansion
DNekMatSharedPtr CreateStdMatrix (const StdMatrixKey &mkey)
 
DNekBlkMatSharedPtr CreateStdStaticCondMatrix (const StdMatrixKey &mkey)
 Create the static condensation of a matrix when using a boundary interior decomposition. More...
 
IndexMapValuesSharedPtr CreateIndexMap (const IndexMapKey &ikey)
 Create an IndexMap which contains mapping information linking any specific element shape with either its boundaries, edges, faces, verteces, etc. More...
 
void BwdTrans_MatOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void GeneralMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void MassMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp)
 
void LaplacianMatrixOp_MatFree_GenericImpl (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void LaplacianMatrixOp_MatFree (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void WeakDerivMatrixOp_MatFree (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void WeakDirectionalDerivMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void MassLevelCurvatureMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void LinearAdvectionDiffusionReactionMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true)
 
void HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
void HelmholtzMatrixOp_MatFree_GenericImpl (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
 
virtual void v_SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_SetCoeffsToOrientation (Array< OneD, NekDouble > &coeffs, StdRegions::Orientation dir)
 
- 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)
 
virtual void v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
 
virtual void v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 
virtual void v_AddRobinMassMatrix (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat)
 
virtual void v_AddRobinEdgeContribution (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, Array< OneD, NekDouble > &coeffs)
 
virtual DNekMatSharedPtr v_BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd)
 
void GetPhysEdgeVarCoeffsFromElement (const int edge, ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &varcoeff, Array< OneD, NekDouble > &outarray)
 
void ReOrientQuadEdgePhysMap (const StdRegions::Orientation orient, const int nq0, Array< OneD, int > &idmap)
 
Array< OneD, unsigned int > v_GetEdgeInverseBoundaryMap (int eid)
 
virtual void v_NegateEdgeNormal (const int edge)
 
virtual bool v_EdgeNormalNegated (const int edge)
 
virtual void v_SetUpPhysNormals (const int edge)
 
const StdRegions::NormalVector & v_GetEdgeNormal (const int edge) const
 
const StdRegions::NormalVector & v_GetSurfaceNormal (const int id) const
 
- Protected Member Functions inherited from Nektar::LocalRegions::Expansion
void ComputeLaplacianMetric ()
 
void ComputeQuadratureMetric ()
 
virtual void v_MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual DNekMatSharedPtr v_BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType)
 
virtual void v_AddEdgeNormBoundaryInt (const int edge, const boost::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 boost::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 
virtual void v_AddFaceNormBoundaryInt (const int face, const boost::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 

Private Member Functions

 QuadExp ()
 

Private Attributes

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

Additional Inherited Members

- Protected Attributes inherited from Nektar::StdRegions::StdExpansion
Array< OneD,
LibUtilities::BasisSharedPtr > 		m_base
 
int m_elmt_id
 
int m_ncoeffs
 
LibUtilities::NekManager
< StdMatrixKey, DNekMat,
StdMatrixKey::opLess > 		m_stdMatrixManager
 
LibUtilities::NekManager
< StdMatrixKey, DNekBlkMat,
StdMatrixKey::opLess > 		m_stdStaticCondMatrixManager
 
LibUtilities::NekManager
< IndexMapKey, IndexMapValues,
IndexMapKey::opLess > 		m_IndexMapManager
 
- Protected Attributes inherited from Nektar::LocalRegions::Expansion2D
std::vector< Expansion1DWeakPtr > m_edgeExp
 
std::vector< bool > m_requireNeg
 
std::map< int,
StdRegions::NormalVector > 		m_edgeNormals
 
std::map< int, bool > m_negatedNormals
 
Expansion3DWeakPtr m_elementLeft
 
Expansion3DWeakPtr m_elementRight
 
int m_elementFaceLeft
 
int m_elementFaceRight
 
- Protected Attributes inherited from Nektar::LocalRegions::Expansion
SpatialDomains::GeometrySharedPtr m_geom
 
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
 
MetricMap m_metrics
 

Detailed Description

Definition at line 52 of file QuadExp.h.

Constructor & Destructor Documentation

Nektar::LocalRegions::QuadExp::QuadExp ( const LibUtilities::BasisKey &  Ba,
const LibUtilities::BasisKey &  Bb,
const SpatialDomains::QuadGeomSharedPtr &  geom 
)

Constructor using BasisKey class for quadrature points and order definition.

Definition at line 48 of file QuadExp.cpp.

50  :
51  StdExpansion (Ba.GetNumModes()*Bb.GetNumModes(),2,Ba,Bb),
52  StdExpansion2D(Ba.GetNumModes()*Bb.GetNumModes(),Ba,Bb),
53  StdQuadExp (Ba,Bb),
54  Expansion (geom),
55  Expansion2D (geom),
56  m_matrixManager(
57  boost::bind(&QuadExp::CreateMatrix, this, _1),
58  std::string("QuadExpMatrix")),
59  m_staticCondMatrixManager(
60  boost::bind(&QuadExp::CreateStaticCondMatrix, this, _1),
61  std::string("QuadExpStaticCondMatrix"))
62  {
63  }
Nektar::LocalRegions::QuadExp::m_staticCondMatrixManager
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: QuadExp.h:286
Nektar::LocalRegions::QuadExp::CreateStaticCondMatrix
DNekScalBlkMatSharedPtr CreateStaticCondMatrix(const MatrixKey &mkey)
Definition: QuadExp.cpp:1919
Nektar::LocalRegions::Expansion::Expansion
Expansion(SpatialDomains::GeometrySharedPtr pGeom)
Definition: Expansion.cpp:46
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
Nektar::LocalRegions::QuadExp::CreateMatrix
DNekScalMatSharedPtr CreateMatrix(const MatrixKey &mkey)
Definition: QuadExp.cpp:1605
Nektar::LocalRegions::Expansion2D::Expansion2D
Expansion2D(SpatialDomains::Geometry2DSharedPtr pGeom)
Definition: Expansion2D.cpp:47
Nektar::StdRegions::StdExpansion::StdExpansion
StdExpansion()
Default Constructor.
Definition: StdExpansion.cpp:44
Nektar::LocalRegions::QuadExp::QuadExp ( const QuadExp &  T)

Definition at line 66 of file QuadExp.cpp.

66  :
67  StdExpansion(T),
68  StdExpansion2D(T),
69  StdQuadExp(T),
70  Expansion (T),
71  Expansion2D (T),
72  m_matrixManager(T.m_matrixManager),
73  m_staticCondMatrixManager(T.m_staticCondMatrixManager)
74  {
75  }
Nektar::LocalRegions::QuadExp::m_staticCondMatrixManager
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: QuadExp.h:286
Nektar::LocalRegions::Expansion::Expansion
Expansion(SpatialDomains::GeometrySharedPtr pGeom)
Definition: Expansion.cpp:46
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
Nektar::LocalRegions::Expansion2D::Expansion2D
Expansion2D(SpatialDomains::Geometry2DSharedPtr pGeom)
Definition: Expansion2D.cpp:47
Nektar::StdRegions::StdExpansion::StdExpansion
StdExpansion()
Default Constructor.
Definition: StdExpansion.cpp:44
Nektar::LocalRegions::QuadExp::~QuadExp ( )
virtual

Definition at line 78 of file QuadExp.cpp.

79  {
80  }
Nektar::LocalRegions::QuadExp::QuadExp ( )
private

Member Function Documentation

DNekScalMatSharedPtr Nektar::LocalRegions::QuadExp::CreateMatrix ( const MatrixKey &  mkey)
protected

Definition at line 1605 of file QuadExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL2, Nektar::LocalRegions::Expansion::BuildVertexMatrix(), Nektar::StdRegions::StdMatrixKey::ConstFactorExists(), Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::eFactorGaussEdge, Nektar::StdRegions::eFactorLambda, Nektar::StdRegions::eFactorSVVCutoffRatio, Nektar::StdRegions::eHelmholtz, Nektar::StdRegions::eHybridDGHelmholtz, Nektar::StdRegions::eInterpGauss, Nektar::StdRegions::eInvHybridDGHelmholtz, Nektar::StdRegions::eInvLaplacianWithUnityMean, Nektar::StdRegions::eInvMass, Nektar::StdRegions::eIProductWRTBase, Nektar::StdRegions::eIProductWRTDerivBase0, Nektar::StdRegions::eIProductWRTDerivBase1, Nektar::StdRegions::eIProductWRTDerivBase2, Nektar::StdRegions::eLaplacian, Nektar::StdRegions::eLaplacian00, Nektar::StdRegions::eLaplacian01, Nektar::StdRegions::eLaplacian11, Nektar::StdRegions::eMass, Nektar::SpatialDomains::eNoGeomType, Nektar::StdRegions::ePreconLinearSpace, Nektar::StdRegions::eWeakDeriv0, Nektar::StdRegions::eWeakDeriv1, Nektar::StdRegions::eWeakDeriv2, Nektar::StdRegions::StdExpansion::GenMatrix(), Nektar::StdRegions::StdMatrixKey::GetConstFactor(), Nektar::StdRegions::StdMatrixKey::GetConstFactors(), Nektar::StdRegions::StdExpansion::GetLocStaticCondMatrix(), Nektar::StdRegions::StdMatrixKey::GetMatrixType(), Nektar::StdRegions::StdMatrixKey::GetNVarCoeff(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdMatrixKey::GetShapeType(), Nektar::StdRegions::StdExpansion::GetStdMatrix(), Nektar::StdRegions::StdMatrixKey::GetVarCoeffs(), Nektar::StdRegions::StdExpansion::m_base, m_matrixManager, Nektar::LocalRegions::Expansion::m_metricinfo, and Nektar::Transpose().

1606  {
1607  DNekScalMatSharedPtr returnval;
1608  LibUtilities::PointsKeyVector ptsKeys = GetPointsKeys();
1609 
1610  ASSERTL2(m_metricinfo->GetGtype() != SpatialDomains::eNoGeomType,
1611  "Geometric information is not set up");
1612 
1613  switch (mkey.GetMatrixType())
1614  {
1615  case StdRegions::eMass:
1616  {
1617  if ((m_metricinfo->GetGtype() ==
1618  SpatialDomains::eDeformed) || (mkey.GetNVarCoeff()))
1619  {
1620  NekDouble one = 1.0;
1621  DNekMatSharedPtr mat = GenMatrix(mkey);
1622  returnval = MemoryManager<DNekScalMat>::
1623  AllocateSharedPtr(one,mat);
1624  }
1625  else
1626  {
1627  NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
1628  DNekMatSharedPtr mat = GetStdMatrix(mkey);
1629  returnval = MemoryManager<DNekScalMat>::
1630  AllocateSharedPtr(jac,mat);
1631  }
1632  }
1633  break;
1634  case StdRegions::eInvMass:
1635  {
1636  if ((m_metricinfo->GetGtype() ==
1637  SpatialDomains::eDeformed) || (mkey.GetNVarCoeff()))
1638  {
1639  NekDouble one = 1.0;
1640  StdRegions::StdMatrixKey masskey(
1641  StdRegions::eMass, DetShapeType(), *this);
1642  DNekMatSharedPtr mat = GenMatrix(masskey);
1643  mat->Invert();
1644 
1645  returnval = MemoryManager<DNekScalMat>::
1646  AllocateSharedPtr(one,mat);
1647  }
1648  else
1649  {
1650  NekDouble fac = 1.0/(m_metricinfo->GetJac(ptsKeys))[0];
1651  DNekMatSharedPtr mat = GetStdMatrix(mkey);
1652  returnval = MemoryManager<DNekScalMat>::
1653  AllocateSharedPtr(fac,mat);
1654  }
1655  }
1656  break;
1657  case StdRegions::eWeakDeriv0:
1658  case StdRegions::eWeakDeriv1:
1659  case StdRegions::eWeakDeriv2:
1660  {
1661  if((m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
1662  || (mkey.GetNVarCoeff()))
1663  {
1664  NekDouble one = 1.0;
1665  DNekMatSharedPtr mat = GenMatrix(mkey);
1666 
1667  returnval = MemoryManager<DNekScalMat>::
1668  AllocateSharedPtr(one,mat);
1669  }
1670  else
1671  {
1672  NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
1673  Array<TwoD, const NekDouble> df =
1674  m_metricinfo->GetDerivFactors(ptsKeys);
1675  int dir = 0;
1676 
1677  switch(mkey.GetMatrixType())
1678  {
1679  case StdRegions::eWeakDeriv0:
1680  dir = 0;
1681  break;
1682  case StdRegions::eWeakDeriv1:
1683  dir = 1;
1684  break;
1685  case StdRegions::eWeakDeriv2:
1686  dir = 2;
1687  break;
1688  default:
1689  break;
1690  }
1691 
1692  MatrixKey deriv0key(StdRegions::eWeakDeriv0,
1693  mkey.GetShapeType(), *this);
1694  MatrixKey deriv1key(StdRegions::eWeakDeriv1,
1695  mkey.GetShapeType(), *this);
1696 
1697  DNekMat &deriv0 = *GetStdMatrix(deriv0key);
1698  DNekMat &deriv1 = *GetStdMatrix(deriv1key);
1699 
1700  int rows = deriv0.GetRows();
1701  int cols = deriv1.GetColumns();
1702 
1703  DNekMatSharedPtr WeakDeriv = MemoryManager<DNekMat>::
1704  AllocateSharedPtr(rows,cols);
1705  (*WeakDeriv) = df[2*dir][0]*deriv0 +
1706  df[2*dir+1][0]*deriv1;
1707  returnval = MemoryManager<DNekScalMat>::
1708  AllocateSharedPtr(jac,WeakDeriv);
1709  }
1710  }
1711  break;
1712  case StdRegions::eLaplacian:
1713  {
1714  if( (m_metricinfo->GetGtype() ==
1715  SpatialDomains::eDeformed) || (mkey.GetNVarCoeff() > 0)
1716  || (mkey.ConstFactorExists
1717  (StdRegions::eFactorSVVCutoffRatio)))
1718  {
1719  NekDouble one = 1.0;
1720  DNekMatSharedPtr mat = GenMatrix(mkey);
1721 
1722  returnval = MemoryManager<DNekScalMat>::
1723  AllocateSharedPtr(one,mat);
1724  }
1725  else
1726  {
1727  MatrixKey lap00key(StdRegions::eLaplacian00,
1728  mkey.GetShapeType(), *this);
1729  MatrixKey lap01key(StdRegions::eLaplacian01,
1730  mkey.GetShapeType(), *this);
1731  MatrixKey lap11key(StdRegions::eLaplacian11,
1732  mkey.GetShapeType(), *this);
1733 
1734  DNekMat &lap00 = *GetStdMatrix(lap00key);
1735  DNekMat &lap01 = *GetStdMatrix(lap01key);
1736  DNekMat &lap11 = *GetStdMatrix(lap11key);
1737 
1738  NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
1739  Array<TwoD, const NekDouble>
1740  gmat = m_metricinfo->GetGmat(ptsKeys);
1741 
1742  int rows = lap00.GetRows();
1743  int cols = lap00.GetColumns();
1744 
1745  DNekMatSharedPtr lap =
1746  MemoryManager<DNekMat>::AllocateSharedPtr(rows,cols);
1747 
1748  (*lap) = gmat[0][0] * lap00 +
1749  gmat[1][0] * (lap01 + Transpose(lap01)) +
1750  gmat[3][0] * lap11;
1751 
1752  returnval = MemoryManager<DNekScalMat>::
1753  AllocateSharedPtr(jac,lap);
1754  }
1755  }
1756  break;
1757  case StdRegions::eInvLaplacianWithUnityMean:
1758  {
1759  DNekMatSharedPtr mat = GenMatrix(mkey);
1760  returnval = MemoryManager<DNekScalMat>::AllocateSharedPtr(1.0,mat);
1761  }
1762  break;
1763  case StdRegions::eHelmholtz:
1764  {
1765  NekDouble lambda =
1766  mkey.GetConstFactor(StdRegions::eFactorLambda);
1767 
1768  MatrixKey masskey(mkey, StdRegions::eMass);
1769  DNekScalMat &MassMat = *(this->m_matrixManager[masskey]);
1770 
1771  MatrixKey lapkey(mkey, StdRegions::eLaplacian);
1772  DNekScalMat &LapMat = *(this->m_matrixManager[lapkey]);
1773 
1774  int rows = LapMat.GetRows();
1775  int cols = LapMat.GetColumns();
1776 
1777  DNekMatSharedPtr helm = MemoryManager<DNekMat>::
1778  AllocateSharedPtr(rows,cols);
1779 
1780  NekDouble one = 1.0;
1781  (*helm) = LapMat + lambda*MassMat;
1782 
1783  returnval =
1784  MemoryManager<DNekScalMat>::AllocateSharedPtr(one,helm);
1785  }
1786  break;
1787  case StdRegions::eIProductWRTBase:
1788  {
1789  if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
1790  {
1791  NekDouble one = 1.0;
1792  DNekMatSharedPtr mat = GenMatrix(mkey);
1793  returnval = MemoryManager<DNekScalMat>::
1794  AllocateSharedPtr(one,mat);
1795  }
1796  else
1797  {
1798  NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
1799  DNekMatSharedPtr mat = GetStdMatrix(mkey);
1800  returnval = MemoryManager<DNekScalMat>::
1801  AllocateSharedPtr(jac,mat);
1802  }
1803  }
1804  break;
1805  case StdRegions::eIProductWRTDerivBase0:
1806  case StdRegions::eIProductWRTDerivBase1:
1807  case StdRegions::eIProductWRTDerivBase2:
1808  {
1809  if(m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
1810  {
1811  NekDouble one = 1.0;
1812  DNekMatSharedPtr mat = GenMatrix(mkey);
1813  returnval = MemoryManager<DNekScalMat>::
1814  AllocateSharedPtr(one,mat);
1815  }
1816  else
1817  {
1818  NekDouble jac = (m_metricinfo->GetJac(ptsKeys))[0];
1819  const Array<TwoD, const NekDouble>& df =
1820  m_metricinfo->GetDerivFactors(ptsKeys);
1821  int dir = 0;
1822 
1823  switch(mkey.GetMatrixType())
1824  {
1825  case StdRegions::eIProductWRTDerivBase0:
1826  dir = 0;
1827  break;
1828  case StdRegions::eIProductWRTDerivBase1:
1829  dir = 1;
1830  break;
1831  case StdRegions::eIProductWRTDerivBase2:
1832  dir = 2;
1833  break;
1834  default:
1835  break;
1836  }
1837 
1838  MatrixKey iProdDeriv0Key(
1839  StdRegions::eIProductWRTDerivBase0,
1840  mkey.GetShapeType(), *this);
1841  MatrixKey iProdDeriv1Key(
1842  StdRegions::eIProductWRTDerivBase1,
1843  mkey.GetShapeType(), *this);
1844 
1845  DNekMat &stdiprod0 = *GetStdMatrix(iProdDeriv0Key);
1846  DNekMat &stdiprod1 = *GetStdMatrix(iProdDeriv0Key);
1847 
1848  int rows = stdiprod0.GetRows();
1849  int cols = stdiprod1.GetColumns();
1850 
1851  DNekMatSharedPtr mat = MemoryManager<DNekMat>::
1852  AllocateSharedPtr(rows,cols);
1853  (*mat) = df[2*dir][0]*stdiprod0 +
1854  df[2*dir+1][0]*stdiprod1;
1855 
1856  returnval = MemoryManager<DNekScalMat>::
1857  AllocateSharedPtr(jac,mat);
1858  }
1859  }
1860  break;
1861  case StdRegions::eInvHybridDGHelmholtz:
1862  {
1863  NekDouble one = 1.0;
1864 
1865  MatrixKey hkey(StdRegions::eHybridDGHelmholtz,
1866  DetShapeType(), *this,
1867  mkey.GetConstFactors(), mkey.GetVarCoeffs());
1868  DNekMatSharedPtr mat = GenMatrix(hkey);
1869 
1870  mat->Invert();
1871  returnval =
1872  MemoryManager<DNekScalMat>::AllocateSharedPtr(one,mat);
1873  }
1874  break;
1875  case StdRegions::eInterpGauss:
1876  {
1877  DNekMatSharedPtr m_Ix;
1878  Array<OneD, NekDouble> coords(1, 0.0);
1879  StdRegions::ConstFactorMap factors = mkey.GetConstFactors();
1880  int edge = (int)factors[StdRegions::eFactorGaussEdge];
1881 
1882  coords[0] = (edge == 0 || edge == 3) ? -1.0 : 1.0;
1883 
1884  m_Ix = m_base[(edge + 1) % 2]->GetI(coords);
1885  returnval =
1886  MemoryManager<DNekScalMat>::AllocateSharedPtr(1.0,m_Ix);
1887  }
1888  break;
1889  case StdRegions::ePreconLinearSpace:
1890  {
1891  NekDouble one = 1.0;
1892  MatrixKey helmkey(
1893  StdRegions::eHelmholtz, mkey.GetShapeType(), *this,
1894  mkey.GetConstFactors(), mkey.GetVarCoeffs());
1895  DNekScalBlkMatSharedPtr helmStatCond =
1896  GetLocStaticCondMatrix(helmkey);
1897  DNekScalMatSharedPtr A =helmStatCond->GetBlock(0,0);
1898  DNekMatSharedPtr R=BuildVertexMatrix(A);
1899 
1900  returnval =
1901  MemoryManager<DNekScalMat>::AllocateSharedPtr(one, R);
1902  }
1903  break;
1904  default:
1905  {
1906  NekDouble one = 1.0;
1907  DNekMatSharedPtr mat = GenMatrix(mkey);
1908 
1909  returnval =
1910  MemoryManager<DNekScalMat>::AllocateSharedPtr(one, mat);
1911  }
1912  break;
1913  }
1914 
1915  return returnval;
1916  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Nektar::StdRegions::StdExpansion::DetShapeType
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:470
Nektar::StdRegions::StdExpansion::GenMatrix
DNekMatSharedPtr GenMatrix(const StdMatrixKey &mkey)
Definition: StdExpansion.h:1049
Nektar::LibUtilities::PointsKeyVector
std::vector< PointsKey > PointsKeyVector
Definition: Points.h:220
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::LocalRegions::Expansion::BuildVertexMatrix
DNekMatSharedPtr BuildVertexMatrix(const DNekScalMatSharedPtr &r_bnd)
Definition: Expansion.cpp:96
Nektar::StdRegions::ConstFactorMap
std::map< ConstFactorType, NekDouble > ConstFactorMap
Definition: StdRegions.hpp:251
Nektar::StdRegions::StdExpansion::GetLocStaticCondMatrix
DNekScalBlkMatSharedPtr GetLocStaticCondMatrix(const LocalRegions::MatrixKey &mkey)
Definition: StdExpansion.h:747
Nektar::DNekMatSharedPtr
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
Nektar::StdRegions::StdExpansion::GetStdMatrix
DNekMatSharedPtr GetStdMatrix(const StdMatrixKey &mkey)
Definition: StdExpansion.h:700
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
Nektar::DNekScalBlkMatSharedPtr
boost::shared_ptr< DNekScalBlkMat > DNekScalBlkMatSharedPtr
Definition: NekTypeDefs.hpp:74
Nektar::Transpose
NekMatrix< InnerMatrixType, BlockMatrixTag > Transpose(NekMatrix< InnerMatrixType, BlockMatrixTag > &rhs)
Definition: BlockMatrix.hpp:221
Nektar::DNekMat
NekMatrix< NekDouble, StandardMatrixTag > DNekMat
Definition: NekTypeDefs.hpp:52
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
ASSERTL2
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Definition: ErrorUtil.hpp:213
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Nektar::DNekScalMat
NekMatrix< NekMatrix< NekDouble, StandardMatrixTag >, ScaledMatrixTag > DNekScalMat
Definition: NekMatrixFwd.hpp:72
DNekScalBlkMatSharedPtr Nektar::LocalRegions::QuadExp::CreateStaticCondMatrix ( const MatrixKey &  mkey)
protected

Definition at line 1919 of file QuadExp.cpp.

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

1921  {
1922  DNekScalBlkMatSharedPtr returnval;
1923 
1924  ASSERTL2(m_metricinfo->GetGtype()
1925  != SpatialDomains::eNoGeomType,
1926  "Geometric information is not set up");
1927 
1928  // set up block matrix system
1929  unsigned int nbdry = NumBndryCoeffs();
1930  unsigned int nint = (unsigned int)(m_ncoeffs - nbdry);
1931  unsigned int exp_size[] = {nbdry,nint};
1932  unsigned int nblks = 2;
1933  returnval = MemoryManager<DNekScalBlkMat>::
1934  AllocateSharedPtr(nblks,nblks,exp_size,exp_size);
1935  //Really need a constructor which takes Arrays
1936  NekDouble factor = 1.0;
1937 
1938  switch (mkey.GetMatrixType())
1939  {
1940  // this can only use stdregions statically condensed system
1941  // for mass matrix
1942  case StdRegions::eMass:
1943  if ((m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
1944  ||(mkey.GetNVarCoeff()))
1945  {
1946  factor = 1.0;
1947  goto UseLocRegionsMatrix;
1948  }
1949  else
1950  {
1951  factor = (m_metricinfo->GetJac(GetPointsKeys()))[0];
1952  goto UseStdRegionsMatrix;
1953  }
1954  break;
1955  default: // use Deformed case for both
1956  // regular and deformed geometries
1957  factor = 1.0;
1958  goto UseLocRegionsMatrix;
1959  break;
1960  UseStdRegionsMatrix:
1961  {
1962  NekDouble invfactor = 1.0/factor;
1963  NekDouble one = 1.0;
1964  DNekBlkMatSharedPtr mat = GetStdStaticCondMatrix(mkey);
1965  DNekScalMatSharedPtr Atmp;
1966  DNekMatSharedPtr Asubmat;
1967 
1968  returnval->SetBlock(0,0,Atmp = MemoryManager<DNekScalMat>::
1969  AllocateSharedPtr(factor,Asubmat = mat->GetBlock(0,0)));
1970  returnval->SetBlock(0,1,Atmp = MemoryManager<DNekScalMat>::
1971  AllocateSharedPtr(one,Asubmat = mat->GetBlock(0,1)));
1972  returnval->SetBlock(1,0,Atmp = MemoryManager<DNekScalMat>::
1973  AllocateSharedPtr(factor,Asubmat = mat->GetBlock(1,0)));
1974  returnval->SetBlock(1,1,Atmp = MemoryManager<DNekScalMat>::
1975  AllocateSharedPtr(invfactor,Asubmat = mat->GetBlock(1,1)));
1976  }
1977  break;
1978  UseLocRegionsMatrix:
1979  {
1980  int i,j;
1981  NekDouble invfactor = 1.0/factor;
1982  NekDouble one = 1.0;
1983  DNekScalMat &mat = *GetLocMatrix(mkey);
1984  DNekMatSharedPtr A = MemoryManager<DNekMat>::
1985  AllocateSharedPtr(nbdry,nbdry);
1986  DNekMatSharedPtr B = MemoryManager<DNekMat>::
1987  AllocateSharedPtr(nbdry,nint);
1988  DNekMatSharedPtr C = MemoryManager<DNekMat>::
1989  AllocateSharedPtr(nint,nbdry);
1990  DNekMatSharedPtr D = MemoryManager<DNekMat>::
1991  AllocateSharedPtr(nint,nint);
1992 
1993  Array<OneD,unsigned int> bmap(nbdry);
1994  Array<OneD,unsigned int> imap(nint);
1995  GetBoundaryMap(bmap);
1996  GetInteriorMap(imap);
1997 
1998  for (i = 0; i < nbdry; ++i)
1999  {
2000  for(j = 0; j < nbdry; ++j)
2001  {
2002  (*A)(i,j) = mat(bmap[i],bmap[j]);
2003  }
2004 
2005  for(j = 0; j < nint; ++j)
2006  {
2007  (*B)(i,j) = mat(bmap[i],imap[j]);
2008  }
2009  }
2010 
2011  for (i = 0; i < nint; ++i)
2012  {
2013  for(j = 0; j < nbdry; ++j)
2014  {
2015  (*C)(i,j) = mat(imap[i],bmap[j]);
2016  }
2017 
2018  for(j = 0; j < nint; ++j)
2019  {
2020  (*D)(i,j) = mat(imap[i],imap[j]);
2021  }
2022  }
2023 
2024  // Calculate static condensed system
2025  if(nint)
2026  {
2027  D->Invert();
2028  (*B) = (*B)*(*D);
2029  (*A) = (*A) - (*B)*(*C);
2030  }
2031 
2032  DNekScalMatSharedPtr Atmp;
2033 
2034  returnval->SetBlock(0,0,Atmp = MemoryManager<DNekScalMat>::
2035  AllocateSharedPtr(factor, A));
2036  returnval->SetBlock(0,1,Atmp = MemoryManager<DNekScalMat>::
2037  AllocateSharedPtr(one, B));
2038  returnval->SetBlock(1,0,Atmp = MemoryManager<DNekScalMat>::
2039  AllocateSharedPtr(factor, C));
2040  returnval->SetBlock(1,1,Atmp = MemoryManager<DNekScalMat>::
2041  AllocateSharedPtr(invfactor, D));
2042 
2043  }
2044  }
2045  return returnval;
2046  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::DNekMatSharedPtr
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::StdRegions::StdExpansion::GetStdStaticCondMatrix
DNekBlkMatSharedPtr GetStdStaticCondMatrix(const StdMatrixKey &mkey)
Definition: StdExpansion.h:705
Nektar::DNekScalBlkMatSharedPtr
boost::shared_ptr< DNekScalBlkMat > DNekScalBlkMatSharedPtr
Definition: NekTypeDefs.hpp:74
Nektar::StdRegions::StdExpansion::GetInteriorMap
void GetInteriorMap(Array< OneD, unsigned int > &outarray)
Definition: StdExpansion.h:821
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::DNekBlkMatSharedPtr
boost::shared_ptr< DNekBlkMat > DNekBlkMatSharedPtr
Definition: NekTypeDefs.hpp:72
Nektar::LocalRegions::Expansion::GetLocMatrix
DNekScalMatSharedPtr GetLocMatrix(const LocalRegions::MatrixKey &mkey)
Definition: Expansion.cpp:83
ASSERTL2
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Definition: ErrorUtil.hpp:213
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Nektar::DNekScalMat
NekMatrix< NekMatrix< NekDouble, StandardMatrixTag >, ScaledMatrixTag > DNekScalMat
Definition: NekMatrixFwd.hpp:72
Nektar::StdRegions::StdExpansion::GetBoundaryMap
void GetBoundaryMap(Array< OneD, unsigned int > &outarray)
Definition: StdExpansion.h:816
void Nektar::LocalRegions::QuadExp::v_ComputeEdgeNormal ( const int  edge)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 1188 of file QuadExp.cpp.

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

1189  {
1190  int i;
1191  const SpatialDomains::GeomFactorsSharedPtr & geomFactors =
1192  GetGeom()->GetMetricInfo();
1193  SpatialDomains::GeomType type = geomFactors->GetGtype();
1194  LibUtilities::PointsKeyVector ptsKeys = GetPointsKeys();
1195  const Array<TwoD, const NekDouble> & df = geomFactors->GetDerivFactors(ptsKeys);
1196  const Array<OneD, const NekDouble> & jac = geomFactors->GetJac(ptsKeys);
1197  int nqe = m_base[0]->GetNumPoints();
1198  int vCoordDim = GetCoordim();
1199 
1200  m_edgeNormals[edge] = Array<OneD, Array<OneD, NekDouble> >
1201  (vCoordDim);
1202  Array<OneD, Array<OneD, NekDouble> > &normal = m_edgeNormals[edge];
1203  for (i = 0; i < vCoordDim; ++i)
1204  {
1205  normal[i] = Array<OneD, NekDouble>(nqe);
1206  }
1207 
1208  // Regular geometry case
1209  if ((type == SpatialDomains::eRegular)||
1210  (type == SpatialDomains::eMovingRegular))
1211  {
1212  NekDouble fac;
1213  // Set up normals
1214  switch (edge)
1215  {
1216  case 0:
1217  for (i = 0; i < vCoordDim; ++i)
1218  {
1219  Vmath::Fill(nqe, -df[2*i+1][0], normal[i], 1);
1220  }
1221  break;
1222  case 1:
1223  for (i = 0; i < vCoordDim; ++i)
1224  {
1225  Vmath::Fill(nqe, df[2*i][0], normal[i], 1);
1226  }
1227  break;
1228  case 2:
1229  for (i = 0; i < vCoordDim; ++i)
1230  {
1231  Vmath::Fill(nqe, df[2*i+1][0], normal[i], 1);
1232  }
1233  break;
1234  case 3:
1235  for (i = 0; i < vCoordDim; ++i)
1236  {
1237  Vmath::Fill(nqe, -df[2*i][0], normal[i], 1);
1238  }
1239  break;
1240  default:
1241  ASSERTL0(false, "edge is out of range (edge < 4)");
1242  }
1243 
1244  // normalise
1245  fac = 0.0;
1246  for (i =0 ; i < vCoordDim; ++i)
1247  {
1248  fac += normal[i][0]*normal[i][0];
1249  }
1250  fac = 1.0/sqrt(fac);
1251  for (i = 0; i < vCoordDim; ++i)
1252  {
1253  Vmath::Smul(nqe, fac, normal[i], 1,normal[i], 1);
1254  }
1255  }
1256  else // Set up deformed normals
1257  {
1258  int j;
1259 
1260  int nquad0 = ptsKeys[0].GetNumPoints();
1261  int nquad1 = ptsKeys[1].GetNumPoints();
1262 
1263  LibUtilities::PointsKey from_key;
1264 
1265  Array<OneD,NekDouble> normals(vCoordDim*max(nquad0,nquad1),0.0);
1266  Array<OneD,NekDouble> edgejac(vCoordDim*max(nquad0,nquad1),0.0);
1267 
1268  // Extract Jacobian along edges and recover local
1269  // derivates (dx/dr) for polynomial interpolation by
1270  // multiplying m_gmat by jacobian
1271 
1272  // Implementation for all the basis except Gauss points
1273  if (m_base[0]->GetPointsType() !=
1274  LibUtilities::eGaussGaussLegendre
1275  && m_base[1]->GetPointsType() !=
1276  LibUtilities::eGaussGaussLegendre)
1277  {
1278  switch (edge)
1279  {
1280  case 0:
1281  for (j = 0; j < nquad0; ++j)
1282  {
1283  edgejac[j] = jac[j];
1284  for (i = 0; i < vCoordDim; ++i)
1285  {
1286  normals[i*nquad0+j] =
1287  -df[2*i+1][j]*edgejac[j];
1288  }
1289  }
1290  from_key = ptsKeys[0];
1291  break;
1292  case 1:
1293  for (j = 0; j < nquad1; ++j)
1294  {
1295  edgejac[j] = jac[nquad0*j+nquad0-1];
1296  for (i = 0; i < vCoordDim; ++i)
1297  {
1298  normals[i*nquad1+j] =
1299  df[2*i][nquad0*j + nquad0-1]
1300  *edgejac[j];
1301  }
1302  }
1303  from_key = ptsKeys[1];
1304  break;
1305  case 2:
1306  for (j = 0; j < nquad0; ++j)
1307  {
1308  edgejac[j] = jac[nquad0*(nquad1-1)+j];
1309  for (i = 0; i < vCoordDim; ++i)
1310  {
1311  normals[i*nquad0+j] =
1312  (df[2*i+1][nquad0*(nquad1-1)+j])
1313  *edgejac[j];
1314  }
1315  }
1316  from_key = ptsKeys[0];
1317  break;
1318  case 3:
1319  for (j = 0; j < nquad1; ++j)
1320  {
1321  edgejac[j] = jac[nquad0*j];
1322  for (i = 0; i < vCoordDim; ++i)
1323  {
1324  normals[i*nquad1+j] =
1325  -df[2*i][nquad0*j]*edgejac[j];
1326  }
1327  }
1328  from_key = ptsKeys[1];
1329  break;
1330  default:
1331  ASSERTL0(false,"edge is out of range (edge < 3)");
1332  }
1333  }
1334  else
1335  {
1336  int nqtot = nquad0 * nquad1;
1337  Array<OneD, NekDouble> tmp_gmat(nqtot, 0.0);
1338  Array<OneD, NekDouble> tmp_gmat_edge(nqe, 0.0);
1339 
1340  switch (edge)
1341  {
1342  case 0:
1343  for (j = 0; j < nquad0; ++j)
1344  {
1345  for (i = 0; i < vCoordDim; ++i)
1346  {
1347  Vmath::Vmul(nqtot,
1348  &(df[2*i+1][0]), 1,
1349  &jac[0], 1,
1350  &(tmp_gmat[0]), 1);
1351  QuadExp::v_GetEdgeInterpVals(
1352  edge, tmp_gmat, tmp_gmat_edge);
1353  normals[i*nquad0+j] = -tmp_gmat_edge[j];
1354  }
1355  }
1356  from_key = ptsKeys[0];
1357  break;
1358  case 1:
1359  for (j = 0; j < nquad1; ++j)
1360  {
1361  for (i = 0; i < vCoordDim; ++i)
1362  {
1363  Vmath::Vmul(nqtot,
1364  &(df[2*i][0]), 1,
1365  &jac[0], 1,
1366  &(tmp_gmat[0]), 1);
1367  QuadExp::v_GetEdgeInterpVals(
1368  edge, tmp_gmat, tmp_gmat_edge);
1369  normals[i*nquad1+j] = tmp_gmat_edge[j];
1370  }
1371  }
1372  from_key = ptsKeys[1];
1373  break;
1374  case 2:
1375  for (j = 0; j < nquad0; ++j)
1376  {
1377  for (i = 0; i < vCoordDim; ++i)
1378  {
1379  Vmath::Vmul(nqtot,
1380  &(df[2*i+1][0]), 1,
1381  &jac[0], 1,
1382  &(tmp_gmat[0]), 1);
1383  QuadExp::v_GetEdgeInterpVals(
1384  edge, tmp_gmat, tmp_gmat_edge);
1385  normals[i*nquad0+j] = tmp_gmat_edge[j];
1386  }
1387  }
1388  from_key = ptsKeys[0];
1389  break;
1390  case 3:
1391  for (j = 0; j < nquad1; ++j)
1392  {
1393  for (i = 0; i < vCoordDim; ++i)
1394  {
1395  Vmath::Vmul(nqtot,
1396  &(df[2*i][0]), 1,
1397  &jac[0], 1,
1398  &(tmp_gmat[0]) ,1);
1399  QuadExp::v_GetEdgeInterpVals(
1400  edge, tmp_gmat, tmp_gmat_edge);
1401  normals[i*nquad1+j] = -tmp_gmat_edge[j];
1402  }
1403  }
1404  from_key = ptsKeys[1];
1405  break;
1406  default:
1407  ASSERTL0(false,"edge is out of range (edge < 3)");
1408  }
1409  }
1410 
1411  int nq = from_key.GetNumPoints();
1412  Array<OneD,NekDouble> work(nqe,0.0);
1413 
1414  // interpolate Jacobian and invert
1415  LibUtilities::Interp1D(
1416  from_key,jac, m_base[0]->GetPointsKey(), work);
1417  Vmath::Sdiv(nq,1.0,&work[0],1,&work[0],1);
1418 
1419  // interpolate
1420  for (i = 0; i < GetCoordim(); ++i)
1421  {
1422  LibUtilities::Interp1D(
1423  from_key,&normals[i*nq],
1424  m_base[0]->GetPointsKey(),
1425  &normal[i][0]);
1426  Vmath::Vmul(nqe, work, 1, normal[i], 1, normal[i], 1);
1427  }
1428 
1429  //normalise normal vectors
1430  Vmath::Zero(nqe,work,1);
1431  for (i = 0; i < GetCoordim(); ++i)
1432  {
1433  Vmath::Vvtvp(nqe,
1434  normal[i], 1,
1435  normal[i],1 ,
1436  work, 1,
1437  work, 1);
1438  }
1439 
1440  Vmath::Vsqrt(nqe,work,1,work,1);
1441  Vmath::Sdiv(nqe,1.0,work,1,work,1);
1442 
1443  for (i = 0; i < GetCoordim(); ++i)
1444  {
1445  Vmath::Vmul(nqe, normal[i], 1, work, 1, normal[i], 1);
1446  }
1447 
1448  // Reverse direction so that points are in
1449  // anticlockwise direction if edge >=2
1450  if (edge >= 2)
1451  {
1452  for (i = 0; i < GetCoordim(); ++i)
1453  {
1454  Vmath::Reverse(nqe, normal[i], 1, normal[i], 1);
1455  }
1456  }
1457  }
1458  if (GetGeom()->GetEorient(edge) == StdRegions::eBackwards)
1459  {
1460  for (i = 0; i < vCoordDim; ++i)
1461  {
1462  if (geomFactors->GetGtype() == SpatialDomains::eDeformed)
1463  {
1464  Vmath::Reverse(nqe, normal[i], 1, normal[i],1);
1465  }
1466  }
1467  }
1468  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::LibUtilities::PointsKeyVector
std::vector< PointsKey > PointsKeyVector
Definition: Points.h:220
Vmath::Vsqrt
void Vsqrt(int n, const T *x, const int incx, T *y, const int incy)
sqrt y = sqrt(x)
Definition: Vmath.cpp:394
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:46
Vmath::Vvtvp
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:428
Nektar::LocalRegions::Expansion2D::m_edgeNormals
std::map< int, StdRegions::NormalVector > m_edgeNormals
Definition: Expansion2D.h:135
Nektar::StdRegions::StdExpansion::GetEorient
StdRegions::Orientation GetEorient(int edge)
Definition: StdExpansion.h:762
Nektar::LocalRegions::QuadExp::v_GetEdgeInterpVals
virtual void v_GetEdgeInterpVals(const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: QuadExp.cpp:814
Vmath::Sdiv
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:257
Nektar::LibUtilities::eGaussGaussLegendre
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:47
Vmath::Reverse
void Reverse(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1071
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:199
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::LocalRegions::Expansion::GetGeom
SpatialDomains::GeometrySharedPtr GetGeom() const
Definition: Expansion.cpp:148
Nektar::SpatialDomains::GeomFactorsSharedPtr
boost::shared_ptr< GeomFactors > GeomFactorsSharedPtr
Pointer to a GeomFactors object.
Definition: GeomFactors.h:62
Nektar::SpatialDomains::eRegular
Geometry is straight-sided with constant geometric factors.
Definition: SpatialDomains.hpp:55
Nektar::LibUtilities::Interp1D
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:54
Nektar::SpatialDomains::GeomType
GeomType
Indicates the type of element geometry.
Definition: SpatialDomains.hpp:52
Vmath::Zero
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:359
Nektar::StdRegions::StdExpansion::GetPointsType
LibUtilities::PointsType GetPointsType(const int dir) const
This function returns the type of quadrature points used in the dir direction.
Definition: StdExpansion.h:216
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vmul
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:169
void Nektar::LocalRegions::QuadExp::v_ComputeLaplacianMetric ( )
protectedvirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 2263 of file QuadExp.cpp.

References Nektar::LocalRegions::Expansion::ComputeQuadratureMetric(), Nektar::SpatialDomains::eDeformed, Nektar::LocalRegions::eMetricLaplacian00, Nektar::LocalRegions::eMetricLaplacian01, Nektar::LocalRegions::eMetricLaplacian02, Nektar::LocalRegions::eMetricLaplacian11, Nektar::LocalRegions::eMetricLaplacian12, Nektar::LocalRegions::eMetricLaplacian22, Nektar::LocalRegions::eMetricQuadrature, Vmath::Fill(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::LocalRegions::Expansion::m_metricinfo, Nektar::LocalRegions::Expansion::m_metrics, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), and Vmath::Vcopy().

2264  {
2265  if (m_metrics.count(eMetricQuadrature) == 0)
2266  {
2267  ComputeQuadratureMetric();
2268  }
2269 
2270  const SpatialDomains::GeomType type = m_metricinfo->GetGtype();
2271  const unsigned int nqtot = GetTotPoints();
2272  const unsigned int dim = 2;
2273  const MetricType m[3][3] = { {eMetricLaplacian00, eMetricLaplacian01, eMetricLaplacian02},
2274  {eMetricLaplacian01, eMetricLaplacian11, eMetricLaplacian12},
2275  {eMetricLaplacian02, eMetricLaplacian12, eMetricLaplacian22}
2276  };
2277 
2278  const Array<TwoD, const NekDouble> gmat =
2279  m_metricinfo->GetGmat(GetPointsKeys());
2280  for (unsigned int i = 0; i < dim; ++i)
2281  {
2282  for (unsigned int j = i; j < dim; ++j)
2283  {
2284  m_metrics[m[i][j]] = Array<OneD, NekDouble>(nqtot);
2285  if (type == SpatialDomains::eDeformed)
2286  {
2287  Vmath::Vcopy(nqtot, &gmat[i*dim+j][0], 1,
2288  &m_metrics[m[i][j]][0], 1);
2289  }
2290  else
2291  {
2292  Vmath::Fill(nqtot, gmat[i*dim+j][0],
2293  &m_metrics[m[i][j]][0], 1);
2294  }
2295  MultiplyByQuadratureMetric(m_metrics[m[i][j]],
2296  m_metrics[m[i][j]]);
2297 
2298  }
2299  }
2300  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:942
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:46
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::StdRegions::StdExpansion::GetTotPoints
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
Definition: StdExpansion.h:141
Nektar::SpatialDomains::GeomType
GeomType
Indicates the type of element geometry.
Definition: SpatialDomains.hpp:52
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
DNekMatSharedPtr Nektar::LocalRegions::QuadExp::v_CreateStdMatrix ( const StdRegions::StdMatrixKey &  mkey)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 1594 of file QuadExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.

1596  {
1597  LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
1598  LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
1599  StdRegions::StdQuadExpSharedPtr tmp =
1600  MemoryManager<StdQuadExp>::AllocateSharedPtr(bkey0,bkey1);
1601  return tmp->GetStdMatrix(mkey);
1602  }
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::StdRegions::StdQuadExpSharedPtr
boost::shared_ptr< StdQuadExp > StdQuadExpSharedPtr
Definition: StdQuadExp.h:262
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
void Nektar::LocalRegions::QuadExp::v_DropLocStaticCondMatrix ( const MatrixKey &  mkey)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 2061 of file QuadExp.cpp.

References m_staticCondMatrixManager.

2062  {
2063  m_staticCondMatrixManager.DeleteObject(mkey);
2064  }
Nektar::LocalRegions::QuadExp::m_staticCondMatrixManager
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: QuadExp.h:286
void Nektar::LocalRegions::QuadExp::v_ExtractDataToCoeffs ( const NekDouble *  data,
const std::vector< unsigned int > &  nummodes,
const int  nmode_offset,
NekDouble *  coeffs 
)
protectedvirtual

Unpack data from input file assuming it comes from the same expansion type.

See also
StdExpansion::ExtractDataToCoeffs

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 1482 of file QuadExp.cpp.

References ASSERTL0, ASSERTL1, Nektar::LibUtilities::eGauss_Lagrange, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGLL_Lagrange, Nektar::LibUtilities::eModified_A, Nektar::StdRegions::StdExpansion::GetBasisType(), Nektar::LibUtilities::Interp2D(), Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::m_ncoeffs, Vmath::Vcopy(), and Vmath::Zero().

1487  {
1488  int data_order0 = nummodes[mode_offset];
1489  int fillorder0 = std::min(m_base[0]->GetNumModes(),data_order0);
1490 
1491  int data_order1 = nummodes[mode_offset + 1];
1492  int order1 = m_base[1]->GetNumModes();
1493  int fillorder1 = min(order1,data_order1);
1494 
1495  switch (m_base[0]->GetBasisType())
1496  {
1497  case LibUtilities::eModified_A:
1498  {
1499  int i;
1500  int cnt = 0;
1501  int cnt1 = 0;
1502 
1503  ASSERTL1(m_base[1]->GetBasisType() ==
1504  LibUtilities::eModified_A,
1505  "Extraction routine not set up for this basis");
1506 
1507  Vmath::Zero(m_ncoeffs,coeffs,1);
1508  for (i = 0; i < fillorder0; ++i)
1509  {
1510  Vmath::Vcopy(fillorder1, data + cnt, 1, coeffs +cnt1, 1);
1511  cnt += data_order1;
1512  cnt1 += order1;
1513  }
1514  }
1515  break;
1516  case LibUtilities::eGLL_Lagrange:
1517  {
1518  // Assume that input is also Gll_Lagrange but no way to check;
1519  LibUtilities::PointsKey
1520  p0(nummodes[0], LibUtilities::eGaussLobattoLegendre);
1521  LibUtilities::PointsKey
1522  p1(nummodes[1], LibUtilities::eGaussLobattoLegendre);
1523  LibUtilities::Interp2D(p0, p1, data,
1524  m_base[0]->GetPointsKey(),
1525  m_base[1]->GetPointsKey(),
1526  coeffs);
1527  }
1528  break;
1529  case LibUtilities::eGauss_Lagrange:
1530  {
1531  // Assume that input is also Gll_Lagrange but no way to check;
1532  LibUtilities::PointsKey
1533  p0(nummodes[0],LibUtilities::eGaussGaussLegendre);
1534  LibUtilities::PointsKey
1535  p1(nummodes[1],LibUtilities::eGaussGaussLegendre);
1536  LibUtilities::Interp2D(p0, p1, data,
1537  m_base[0]->GetPointsKey(),
1538  m_base[1]->GetPointsKey(),
1539  coeffs);
1540  }
1541  break;
1542  default:
1543  ASSERTL0(false,
1544  "basis is either not set up or not hierarchicial");
1545  }
1546  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::LibUtilities::eModified_A
Principle Modified Functions .
Definition: BasisType.h:49
Nektar::LibUtilities::eGauss_Lagrange
Lagrange Polynomials using the Gauss points .
Definition: BasisType.h:54
Nektar::LibUtilities::eGaussGaussLegendre
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:47
Nektar::LibUtilities::Interp2D
void Interp2D(const BasisKey &fbasis0, const BasisKey &fbasis1, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, const BasisKey &tbasis1, Array< OneD, NekDouble > &to)
this function interpolates a 2D function evaluated at the quadrature points of the 2D basis...
Definition: Interp.cpp:116
Nektar::StdRegions::StdExpansion::GetBasisType
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
Definition: StdExpansion.h:165
Nektar::LibUtilities::eGLL_Lagrange
Lagrange for SEM basis .
Definition: BasisType.h:53
Vmath::Zero
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:359
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
Nektar::LibUtilities::eGaussLobattoLegendre
1D Gauss-Lobatto-Legendre quadrature points
Definition: PointsType.h:50
void Nektar::LocalRegions::QuadExp::v_FwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Transform a given function from physical quadrature space to coefficient space.

See also
StdExpansion::FwdTrans

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 259 of file QuadExp.cpp.

References Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::eCopy, Nektar::StdRegions::eInvMass, Nektar::eWrapper, Nektar::StdRegions::StdExpansion::IProductWRTBase(), Nektar::StdRegions::StdExpansion::m_base, m_matrixManager, Nektar::StdRegions::StdExpansion::m_ncoeffs, and Vmath::Vcopy().

262  {
263  if ((m_base[0]->Collocation())&&(m_base[1]->Collocation()))
264  {
265  Vmath::Vcopy(m_ncoeffs, inarray, 1, outarray, 1);
266  }
267  else
268  {
269  IProductWRTBase(inarray,outarray);
270 
271  // get Mass matrix inverse
272  MatrixKey masskey(StdRegions::eInvMass,
273  DetShapeType(),*this);
274  DNekScalMatSharedPtr matsys = m_matrixManager[masskey];
275 
276  // copy inarray in case inarray == outarray
277  NekVector<NekDouble> in(m_ncoeffs,outarray,eCopy);
278  NekVector<NekDouble> out(m_ncoeffs,outarray,eWrapper);
279 
280  out = (*matsys)*in;
281  }
282  }
Nektar::StdRegions::StdExpansion::DetShapeType
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:470
Nektar::StdRegions::StdExpansion::IProductWRTBase
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 expa...
Definition: StdExpansion.h:629
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Nektar::LocalRegions::QuadExp::v_FwdTrans_BndConstrained ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 285 of file QuadExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::StdRegions::eBackwards, Nektar::StdRegions::eMass, Nektar::StdRegions::StdExpansion::GetEdgeToElementMap(), Nektar::StdRegions::StdExpansion::GetEorient(), Nektar::LocalRegions::Expansion2D::GetGeom2D(), Nektar::StdRegions::StdExpansion::GetInteriorMap(), Nektar::StdRegions::StdExpansion::IProductWRTBase(), Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::m_ncoeffs, m_staticCondMatrixManager, Nektar::StdRegions::StdExpansion::MassMatrixOp(), Nektar::StdRegions::StdExpansion::NumBndryCoeffs(), sign, Vmath::Vcopy(), and Vmath::Vsub().

288  {
289  if ((m_base[0]->Collocation())&&(m_base[1]->Collocation()))
290  {
291  Vmath::Vcopy(m_ncoeffs, inarray, 1, outarray, 1);
292  }
293  else
294  {
295  int i,j;
296  int npoints[2] = {m_base[0]->GetNumPoints(),
297  m_base[1]->GetNumPoints()};
298  int nmodes[2] = {m_base[0]->GetNumModes(),
299  m_base[1]->GetNumModes()};
300 
301  fill(outarray.get(), outarray.get()+m_ncoeffs, 0.0 );
302 
303  Array<OneD, NekDouble> physEdge[4];
304  Array<OneD, NekDouble> coeffEdge[4];
305  StdRegions::Orientation orient[4];
306  for (i = 0; i < 4; i++)
307  {
308  physEdge[i] = Array<OneD, NekDouble>(npoints[i%2]);
309  coeffEdge[i] = Array<OneD, NekDouble>(nmodes[i%2]);
310  orient[i] = GetEorient(i);
311  }
312 
313  for (i = 0; i < npoints[0]; i++)
314  {
315  physEdge[0][i] = inarray[i];
316  physEdge[2][i] = inarray[npoints[0]*npoints[1]-1-i];
317  }
318 
319  for (i = 0; i < npoints[1]; i++)
320  {
321  physEdge[1][i] =
322  inarray[npoints[0]-1+i*npoints[0]];
323  physEdge[3][i] =
324  inarray[(npoints[1]-1)*npoints[0]-i*npoints[0]];
325  }
326 
327  for (i = 0; i < 4; i++)
328  {
329  if ( orient[i] == StdRegions::eBackwards )
330  {
331  reverse((physEdge[i]).get(),
332  (physEdge[i]).get() + npoints[i%2] );
333  }
334  }
335 
336  SegExpSharedPtr segexp[4];
337  for (i = 0; i < 4; i++)
338  {
339  segexp[i] = MemoryManager<LocalRegions::SegExp>::
340  AllocateSharedPtr(
341  m_base[i%2]->GetBasisKey(),GetGeom2D()->GetEdge(i));
342  }
343 
344  Array<OneD, unsigned int> mapArray;
345  Array<OneD, int> signArray;
346  NekDouble sign;
347 
348  for (i = 0; i < 4; i++)
349  {
350  segexp[i%2]->FwdTrans_BndConstrained(
351  physEdge[i],coeffEdge[i]);
352 
353  GetEdgeToElementMap(i,orient[i],mapArray,signArray);
354  for (j=0; j < nmodes[i%2]; j++)
355  {
356  sign = (NekDouble) signArray[j];
357  outarray[ mapArray[j] ] = sign * coeffEdge[i][j];
358  }
359  }
360 
361  int nBoundaryDofs = NumBndryCoeffs();
362  int nInteriorDofs = m_ncoeffs - nBoundaryDofs;
363 
364  if (nInteriorDofs > 0) {
365  Array<OneD, NekDouble> tmp0(m_ncoeffs);
366  Array<OneD, NekDouble> tmp1(m_ncoeffs);
367 
368  StdRegions::StdMatrixKey
369  stdmasskey(StdRegions::eMass,DetShapeType(),*this);
370  MassMatrixOp(outarray,tmp0,stdmasskey);
371  IProductWRTBase(inarray,tmp1);
372 
373  Vmath::Vsub(m_ncoeffs, tmp1, 1, tmp0, 1, tmp1, 1);
374 
375  // get Mass matrix inverse (only of interior DOF)
376  // use block (1,1) of the static condensed system
377  // note: this block alreay contains the inverse matrix
378  MatrixKey
379  masskey(StdRegions::eMass,DetShapeType(),*this);
380  DNekScalMatSharedPtr matsys =
381  (m_staticCondMatrixManager[masskey])->GetBlock(1,1);
382 
383  Array<OneD, NekDouble> rhs(nInteriorDofs);
384  Array<OneD, NekDouble> result(nInteriorDofs);
385 
386  GetInteriorMap(mapArray);
387 
388  for (i = 0; i < nInteriorDofs; i++)
389  {
390  rhs[i] = tmp1[ mapArray[i] ];
391  }
392 
393  Blas::Dgemv('N', nInteriorDofs, nInteriorDofs,
394  matsys->Scale(),
395  &((matsys->GetOwnedMatrix())->GetPtr())[0],
396  nInteriorDofs,rhs.get(),1,0.0,result.get(),1);
397 
398  for (i = 0; i < nInteriorDofs; i++)
399  {
400  outarray[ mapArray[i] ] = result[i];
401  }
402  }
403  }
404 
405  }
Nektar::StdRegions::StdExpansion::DetShapeType
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:470
Nektar::StdRegions::StdExpansion::MassMatrixOp
void MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:971
Nektar::LocalRegions::QuadExp::m_staticCondMatrixManager
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: QuadExp.h:286
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
sign
#define sign(a, b)
return the sign(b)*a
Definition: Polylib.cpp:22
Nektar::StdRegions::StdExpansion::IProductWRTBase
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 expa...
Definition: StdExpansion.h:629
Nektar::StdRegions::StdExpansion::GetEorient
StdRegions::Orientation GetEorient(int edge)
Definition: StdExpansion.h:762
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::LocalRegions::SegExpSharedPtr
boost::shared_ptr< SegExp > SegExpSharedPtr
Definition: SegExp.h:266
Nektar::LocalRegions::Expansion2D::GetGeom2D
SpatialDomains::Geometry2DSharedPtr GetGeom2D() const
Definition: Expansion2D.h:269
Nektar::StdRegions::StdExpansion::GetInteriorMap
void GetInteriorMap(Array< OneD, unsigned int > &outarray)
Definition: StdExpansion.h:821
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Vmath::Vsub
void Vsub(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Subtract vector z = x-y.
Definition: Vmath.cpp:329
Nektar::StdRegions::StdExpansion::GetEdgeToElementMap
void GetEdgeToElementMap(const int eid, const Orientation edgeOrient, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, int P=-1)
Definition: StdExpansion.h:846
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Nektar::LocalRegions::QuadExp::v_GeneralMatrixOp_MatOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Definition at line 2136 of file QuadExp.cpp.

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

2140  {
2141  MatrixKey newkey(mkey);
2142  DNekScalMatSharedPtr mat = GetLocMatrix(newkey);
2143 
2144  if (inarray.get() == outarray.get())
2145  {
2146  Array<OneD,NekDouble> tmp(m_ncoeffs);
2147  Vmath::Vcopy(m_ncoeffs,inarray.get(),1,tmp.get(),1);
2148 
2149  Blas::Dgemv('N',m_ncoeffs,m_ncoeffs, mat->Scale(),
2150  (mat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
2151  tmp.get(), 1, 0.0, outarray.get(), 1);
2152  }
2153  else
2154  {
2155  Blas::Dgemv('N',m_ncoeffs,m_ncoeffs,mat->Scale(),
2156  (mat->GetOwnedMatrix())->GetPtr().get(), m_ncoeffs,
2157  inarray.get(), 1, 0.0, outarray.get(), 1);
2158  }
2159  }
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::LocalRegions::Expansion::GetLocMatrix
DNekScalMatSharedPtr GetLocMatrix(const LocalRegions::MatrixKey &mkey)
Definition: Expansion.cpp:83
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
DNekMatSharedPtr Nektar::LocalRegions::QuadExp::v_GenMatrix ( const StdRegions::StdMatrixKey &  mkey)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 1573 of file QuadExp.cpp.

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

1575  {
1576  DNekMatSharedPtr returnval;
1577  switch (mkey.GetMatrixType())
1578  {
1579  case StdRegions::eHybridDGHelmholtz:
1580  case StdRegions::eHybridDGLamToU:
1581  case StdRegions::eHybridDGLamToQ0:
1582  case StdRegions::eHybridDGLamToQ1:
1583  case StdRegions::eHybridDGLamToQ2:
1584  case StdRegions::eHybridDGHelmBndLam:
1585  case StdRegions::eInvLaplacianWithUnityMean:
1586  returnval = Expansion2D::v_GenMatrix(mkey);
1587  break;
1588  default:
1589  returnval = StdQuadExp::v_GenMatrix(mkey);
1590  }
1591  return returnval;
1592  }
Nektar::DNekMatSharedPtr
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
Nektar::LocalRegions::Expansion2D::v_GenMatrix
virtual DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey)
Definition: Expansion2D.cpp:639
const LibUtilities::BasisSharedPtr & Nektar::LocalRegions::QuadExp::v_GetBasis ( int  dir) const
protectedvirtual

Definition at line 1561 of file QuadExp.cpp.

References ASSERTL1, and Nektar::StdRegions::StdExpansion::m_base.

1562  {
1563  ASSERTL1(dir >= 0 &&dir <= 1, "input dir is out of range");
1564  return m_base[dir];
1565  }
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
StdRegions::Orientation Nektar::LocalRegions::QuadExp::v_GetCartesianEorient ( int  edge)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 1555 of file QuadExp.cpp.

References Nektar::LocalRegions::Expansion2D::GetGeom2D().

1556  {
1557  return GetGeom2D()->GetCartesianEorient(edge);
1558  }
Nektar::LocalRegions::Expansion2D::GetGeom2D
SpatialDomains::Geometry2DSharedPtr GetGeom2D() const
Definition: Expansion2D.h:269
void Nektar::LocalRegions::QuadExp::v_GetCoord ( const Array< OneD, const NekDouble > &  Lcoords,
Array< OneD, NekDouble > &  coords 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 627 of file QuadExp.cpp.

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

629  {
630  int i;
631 
632  ASSERTL1(Lcoords[0] >= -1.0 && Lcoords[1] <= 1.0 &&
633  Lcoords[1] >= -1.0 && Lcoords[1] <=1.0,
634  "Local coordinates are not in region [-1,1]");
635 
636  m_geom->FillGeom();
637  for (i = 0; i < m_geom->GetCoordim(); ++i)
638  {
639  coords[i] = m_geom->GetCoord(i,Lcoords);
640  }
641  }
Nektar::LocalRegions::Expansion::m_geom
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
int Nektar::LocalRegions::QuadExp::v_GetCoordim ( void  )
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 1476 of file QuadExp.cpp.

References Nektar::LocalRegions::Expansion::m_geom.

1477  {
1478  return m_geom->GetCoordim();
1479  }
Nektar::LocalRegions::Expansion::m_geom
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
void Nektar::LocalRegions::QuadExp::v_GetCoords ( Array< OneD, NekDouble > &  coords_1,
Array< OneD, NekDouble > &  coords_2,
Array< OneD, NekDouble > &  coords_3 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 618 of file QuadExp.cpp.

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

622  {
623  Expansion::v_GetCoords(coords_0, coords_1, coords_2);
624  }
Nektar::LocalRegions::Expansion::v_GetCoords
virtual void v_GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3)
Definition: Expansion.cpp:211
void Nektar::LocalRegions::QuadExp::v_GetEdgeInterpVals ( const int  edge,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 814 of file QuadExp.cpp.

References ASSERTL0, Vmath::Ddot(), Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::StdRegions::eFactorGaussEdge, Nektar::StdRegions::eInterpGauss, Nektar::StdRegions::StdExpansion::m_base, and m_matrixManager.

Referenced by v_ComputeEdgeNormal(), v_GetEdgePhysVals(), and v_GetEdgeQFactors().

817  {
818  int i;
819  int nq0 = m_base[0]->GetNumPoints();
820  int nq1 = m_base[1]->GetNumPoints();
821 
822  StdRegions::ConstFactorMap factors;
823  factors[StdRegions::eFactorGaussEdge] = edge;
824 
825  StdRegions::StdMatrixKey key(
826  StdRegions::eInterpGauss,
827  DetShapeType(),*this,factors);
828 
829  DNekScalMatSharedPtr mat_gauss = m_matrixManager[key];
830 
831  switch (edge)
832  {
833  case 0:
834  {
835  for (i = 0; i < nq0; i++)
836  {
837  outarray[i] = Blas::Ddot(
838  nq1, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
839  1, &inarray[i], nq0);
840  }
841  break;
842  }
843  case 1:
844  {
845  for (i = 0; i < nq1; i++)
846  {
847  outarray[i] = Blas::Ddot(
848  nq0, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
849  1, &inarray[i * nq0], 1);
850  }
851  break;
852  }
853  case 2:
854  {
855  for (i = 0; i < nq0; i++)
856  {
857  outarray[i] = Blas::Ddot(
858  nq1, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
859  1, &inarray[i], nq0);
860  }
861  break;
862  }
863  case 3:
864  {
865  for (i = 0; i < nq1; i++)
866  {
867  outarray[i] = Blas::Ddot(
868  nq0, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
869  1, &inarray[i * nq0], 1);
870  }
871  break;
872  }
873  default:
874  ASSERTL0(false, "edge value (< 3) is out of range");
875  break;
876  }
877  }
Nektar::StdRegions::StdExpansion::DetShapeType
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:470
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::StdRegions::ConstFactorMap
std::map< ConstFactorType, NekDouble > ConstFactorMap
Definition: StdRegions.hpp:251
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
Vmath::Ddot
T Ddot(int n, const Array< OneD, const T > &w, const int incw, const Array< OneD, const T > &x, const int incx, const Array< OneD, const int > &y, const int incy)
Definition: VmathArray.hpp:434
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
void Nektar::LocalRegions::QuadExp::v_GetEdgePhysMap ( const int  edge,
Array< OneD, int > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 880 of file QuadExp.cpp.

References ASSERTL0, and Nektar::StdRegions::StdExpansion::m_base.

883  {
884  int nquad0 = m_base[0]->GetNumPoints();
885  int nquad1 = m_base[1]->GetNumPoints();
886 
887  // Get points in Cartesian orientation
888  switch (edge)
889  {
890  case 0:
891  outarray = Array<OneD, int>(nquad0);
892  for (int i = 0; i < nquad0; ++i)
893  {
894  outarray[i] = i;
895  }
896  break;
897  case 1:
898  outarray = Array<OneD, int>(nquad1);
899  for (int i = 0; i < nquad1; ++i)
900  {
901  outarray[i] = (nquad0-1) + i*nquad1;
902  }
903  break;
904  case 2:
905  outarray = Array<OneD, int>(nquad0);
906  for (int i = 0; i < nquad0; ++i)
907  {
908  outarray[i] = i + nquad0*(nquad1-1);
909  }
910  break;
911  case 3:
912  outarray = Array<OneD, int>(nquad1);
913  for (int i = 0; i < nquad1; ++i)
914  {
915  outarray[i] = i + i*(nquad0-1);
916  }
917  break;
918  default:
919  ASSERTL0(false, "edge value (< 3) is out of range");
920  break;
921  }
922 
923  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
void Nektar::LocalRegions::QuadExp::v_GetEdgePhysVals ( const int  edge,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Extract the physical values along edge edge from inarray into outarray following the local edge orientation and point distribution defined by defined in EdgeExp.

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 674 of file QuadExp.cpp.

References ASSERTL0, Nektar::StdRegions::eForwards, Nektar::StdRegions::StdExpansion::GetEorient(), Nektar::StdRegions::StdExpansion::m_base, and Vmath::Vcopy().

Referenced by v_GetTracePhysVals().

678  {
679  int nquad0 = m_base[0]->GetNumPoints();
680  int nquad1 = m_base[1]->GetNumPoints();
681 
682  StdRegions::Orientation edgedir = GetEorient(edge);
683  switch(edge)
684  {
685  case 0:
686  if (edgedir == StdRegions::eForwards)
687  {
688  Vmath::Vcopy(nquad0,&(inarray[0]),1,&(outarray[0]),1);
689  }
690  else
691  {
692  Vmath::Vcopy(nquad0,&(inarray[0])+(nquad0-1),-1,
693  &(outarray[0]),1);
694  }
695  break;
696  case 1:
697  if (edgedir == StdRegions::eForwards)
698  {
699  Vmath::Vcopy(nquad1,&(inarray[0])+(nquad0-1),nquad0,
700  &(outarray[0]),1);
701  }
702  else
703  {
704  Vmath::Vcopy(nquad1,&(inarray[0])+(nquad0*nquad1-1),
705  -nquad0, &(outarray[0]),1);
706  }
707  break;
708  case 2:
709  if (edgedir == StdRegions::eForwards)
710  {
711  Vmath::Vcopy(nquad0,&(inarray[0])+(nquad0*nquad1-1),-1,
712  &(outarray[0]),1);
713  }
714  else
715  {
716  Vmath::Vcopy(nquad0,&(inarray[0])+nquad0*(nquad1-1),1,
717  &(outarray[0]),1);
718  }
719  break;
720  case 3:
721  if (edgedir == StdRegions::eForwards)
722  {
723  Vmath::Vcopy(nquad1,&(inarray[0]) + nquad0*(nquad1-1),
724  -nquad0,&(outarray[0]),1);
725  }
726  else
727  {
728  Vmath::Vcopy(nquad1,&(inarray[0]),nquad0,
729  &(outarray[0]),1);
730  }
731  break;
732  default:
733  ASSERTL0(false,"edge value (< 3) is out of range");
734  break;
735  }
736  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::StdRegions::StdExpansion::GetEorient
StdRegions::Orientation GetEorient(int edge)
Definition: StdExpansion.h:762
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Nektar::LocalRegions::QuadExp::v_GetEdgePhysVals ( const int  edge,
const StdRegions::StdExpansionSharedPtr &  EdgeExp,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Definition at line 750 of file QuadExp.cpp.

References ASSERTL0, Nektar::StdRegions::eBackwards, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::StdRegions::StdExpansion::GetCartesianEorient(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::LibUtilities::Interp1D(), Nektar::StdRegions::StdExpansion::m_base, Vmath::Reverse(), v_GetEdgeInterpVals(), and Vmath::Vcopy().

755  {
756  int nquad0 = m_base[0]->GetNumPoints();
757  int nquad1 = m_base[1]->GetNumPoints();
758 
759  // Implementation for all the basis except Gauss points
760  if (m_base[0]->GetPointsType() !=
761  LibUtilities::eGaussGaussLegendre &&
762  m_base[1]->GetPointsType() !=
763  LibUtilities::eGaussGaussLegendre)
764  {
765  // get points in Cartesian orientation
766  switch (edge)
767  {
768  case 0:
769  Vmath::Vcopy(nquad0,&(inarray[0]),1,&(outarray[0]),1);
770  break;
771  case 1:
772  Vmath::Vcopy(nquad1,&(inarray[0])+(nquad0-1),
773  nquad0,&(outarray[0]),1);
774  break;
775  case 2:
776  Vmath::Vcopy(nquad0,&(inarray[0])+nquad0*(nquad1-1),1,
777  &(outarray[0]),1);
778  break;
779  case 3:
780  Vmath::Vcopy(nquad1,&(inarray[0]),nquad0,
781  &(outarray[0]),1);
782  break;
783  default:
784  ASSERTL0(false,"edge value (< 3) is out of range");
785  break;
786  }
787  }
788  else
789  {
790  QuadExp::v_GetEdgeInterpVals(edge, inarray, outarray);
791  }
792 
793  // Interpolate if required
794  if (m_base[edge%2]->GetPointsKey() !=
795  EdgeExp->GetBasis(0)->GetPointsKey())
796  {
797  Array<OneD,NekDouble> outtmp(max(nquad0,nquad1));
798 
799  outtmp = outarray;
800 
801  LibUtilities::Interp1D(
802  m_base[edge%2]->GetPointsKey(), outtmp,
803  EdgeExp->GetBasis(0)->GetPointsKey(), outarray);
804  }
805 
806  //Reverse data if necessary
807  if(GetCartesianEorient(edge) == StdRegions::eBackwards)
808  {
809  Vmath::Reverse(EdgeExp->GetNumPoints(0),&outarray[0], 1,
810  &outarray[0], 1);
811  }
812  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::StdRegions::StdExpansion::GetCartesianEorient
StdRegions::Orientation GetCartesianEorient(int edge)
Definition: StdExpansion.h:772
Nektar::LocalRegions::QuadExp::v_GetEdgeInterpVals
virtual void v_GetEdgeInterpVals(const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: QuadExp.cpp:814
Nektar::LibUtilities::eGaussGaussLegendre
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:47
Vmath::Reverse
void Reverse(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1071
Nektar::LibUtilities::Interp1D
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:54
Nektar::StdRegions::StdExpansion::GetPointsType
LibUtilities::PointsType GetPointsType(const int dir) const
This function returns the type of quadrature points used in the dir direction.
Definition: StdExpansion.h:216
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Nektar::LocalRegions::QuadExp::v_GetEdgeQFactors ( const int  edge,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 928 of file QuadExp.cpp.

References ASSERTL0, Nektar::SpatialDomains::eDeformed, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Reverse(), v_GetEdgeInterpVals(), Vmath::Vcopy(), and Vmath::Vmul().

931  {
932  int i;
933  int nquad0 = m_base[0]->GetNumPoints();
934  int nquad1 = m_base[1]->GetNumPoints();
935 
936  LibUtilities::PointsKeyVector ptsKeys = GetPointsKeys();
937  const Array<OneD, const NekDouble>& jac = m_metricinfo->GetJac(ptsKeys);
938  const Array<TwoD, const NekDouble>& df = m_metricinfo->GetDerivFactors(ptsKeys);
939 
940  Array<OneD, NekDouble> j (max(nquad0, nquad1), 0.0);
941  Array<OneD, NekDouble> g0(max(nquad0, nquad1), 0.0);
942  Array<OneD, NekDouble> g1(max(nquad0, nquad1), 0.0);
943  Array<OneD, NekDouble> g2(max(nquad0, nquad1), 0.0);
944  Array<OneD, NekDouble> g3(max(nquad0, nquad1), 0.0);
945 
946  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
947  {
948  // Implementation for all the basis except Gauss points
949  if (m_base[0]->GetPointsType()
950  != LibUtilities::eGaussGaussLegendre
951  && m_base[1]->GetPointsType() !=
952  LibUtilities::eGaussGaussLegendre)
953  {
954  switch (edge)
955  {
956  case 0:
957  Vmath::Vcopy(nquad0, &(df[1][0]),
958  1, &(g1[0]), 1);
959  Vmath::Vcopy(nquad0, &(df[3][0]),
960  1, &(g3[0]), 1);
961  Vmath::Vcopy(nquad0, &(jac[0]),1, &(j[0]), 1);
962 
963  for (i = 0; i < nquad0; ++i)
964  {
965  outarray[i] = j[i]*sqrt(g1[i]*g1[i]
966  + g3[i]*g3[i]);
967  }
968  break;
969  case 1:
970  Vmath::Vcopy(nquad1,
971  &(df[0][0])+(nquad0-1), nquad0,
972  &(g0[0]), 1);
973 
974  Vmath::Vcopy(nquad1,
975  &(df[2][0])+(nquad0-1), nquad0,
976  &(g2[0]), 1);
977 
978  Vmath::Vcopy(nquad1,
979  &(jac[0])+(nquad0-1), nquad0,
980  &(j[0]), 1);
981 
982  for (i = 0; i < nquad0; ++i)
983  {
984  outarray[i] = j[i]*sqrt(g0[i]*g0[i] +
985  g2[i]*g2[i]);
986  }
987  break;
988  case 2:
989 
990  Vmath::Vcopy(nquad0,
991  &(df[1][0])+(nquad0*nquad1-1), -1,
992  &(g1[0]), 1);
993 
994  Vmath::Vcopy(nquad0,
995  &(df[3][0])+(nquad0*nquad1-1), -1,
996  &(g3[0]), 1);
997 
998  Vmath::Vcopy(nquad0,
999  &(jac[0])+(nquad0*nquad1-1), -1,
1000  &(j[0]), 1);
1001 
1002  for (i = 0; i < nquad0; ++i)
1003  {
1004  outarray[i] = j[i]*sqrt(g1[i]*g1[i]
1005  + g3[i]*g3[i]);
1006  }
1007  break;
1008  case 3:
1009 
1010  Vmath::Vcopy(nquad1,
1011  &(df[0][0])+nquad0*(nquad1-1),
1012  -nquad0,&(g0[0]), 1);
1013 
1014  Vmath::Vcopy(nquad1,
1015  &(df[2][0])+nquad0*(nquad1-1),
1016  -nquad0,&(g2[0]), 1);
1017 
1018  Vmath::Vcopy(nquad1,
1019  &(jac[0])+nquad0*(nquad1-1), -nquad0,
1020  &(j[0]), 1);
1021 
1022  for (i = 0; i < nquad0; ++i)
1023  {
1024  outarray[i] = j[i]*sqrt(g0[i]*g0[i] +
1025  g2[i]*g2[i]);
1026  }
1027  break;
1028  default:
1029  ASSERTL0(false,"edge value (< 3) is out of range");
1030  break;
1031  }
1032  }
1033  else
1034  {
1035  int nqtot = nquad0 * nquad1;
1036  Array<OneD, NekDouble> tmp_gmat0(nqtot, 0.0);
1037  Array<OneD, NekDouble> tmp_gmat1(nqtot, 0.0);
1038  Array<OneD, NekDouble> tmp_gmat2(nqtot, 0.0);
1039  Array<OneD, NekDouble> tmp_gmat3(nqtot, 0.0);
1040  Array<OneD, NekDouble> g0_edge(max(nquad0, nquad1), 0.0);
1041  Array<OneD, NekDouble> g1_edge(max(nquad0, nquad1), 0.0);
1042  Array<OneD, NekDouble> g2_edge(max(nquad0, nquad1), 0.0);
1043  Array<OneD, NekDouble> g3_edge(max(nquad0, nquad1), 0.0);
1044  Array<OneD, NekDouble> jac_edge(max(nquad0, nquad1), 0.0);
1045 
1046  switch (edge)
1047  {
1048  case 0:
1049  Vmath::Vmul(nqtot,&(df[1][0]),1,&jac[0],1,
1050  &(tmp_gmat1[0]),1);
1051  Vmath::Vmul(nqtot,&(df[3][0]),1,&jac[0],1,
1052  &(tmp_gmat3[0]),1);
1053  QuadExp::v_GetEdgeInterpVals(
1054  edge, tmp_gmat1, g1_edge);
1055  QuadExp::v_GetEdgeInterpVals(
1056  edge, tmp_gmat3, g3_edge);
1057 
1058  for (i = 0; i < nquad0; ++i)
1059  {
1060  outarray[i] = sqrt(g1_edge[i]*g1_edge[i] +
1061  g3_edge[i]*g3_edge[i]);
1062  }
1063  break;
1064 
1065  case 1:
1066  Vmath::Vmul(nqtot,
1067  &(df[0][0]), 1,
1068  &jac[0], 1,
1069  &(tmp_gmat0[0]), 1);
1070  Vmath::Vmul(nqtot,
1071  &(df[2][0]), 1,
1072  &jac[0], 1,
1073  &(tmp_gmat2[0]),
1074  1);
1075  QuadExp::v_GetEdgeInterpVals(
1076  edge, tmp_gmat0, g0_edge);
1077  QuadExp::v_GetEdgeInterpVals(
1078  edge, tmp_gmat2, g2_edge);
1079 
1080  for (i = 0; i < nquad1; ++i)
1081  {
1082  outarray[i] = sqrt(g0_edge[i]*g0_edge[i]
1083  + g2_edge[i]*g2_edge[i]);
1084  }
1085 
1086  break;
1087  case 2:
1088 
1089  Vmath::Vmul(nqtot,
1090  &(df[1][0]), 1,
1091  &jac[0], 1,
1092  &(tmp_gmat1[0]), 1);
1093  Vmath::Vmul(nqtot,
1094  &(df[3][0]), 1,
1095  &jac[0], 1,
1096  &(tmp_gmat3[0]),1);
1097  QuadExp::v_GetEdgeInterpVals(
1098  edge, tmp_gmat1, g1_edge);
1099  QuadExp::v_GetEdgeInterpVals(
1100  edge, tmp_gmat3, g3_edge);
1101 
1102 
1103  for (i = 0; i < nquad0; ++i)
1104  {
1105  outarray[i] = sqrt(g1_edge[i]*g1_edge[i]
1106  + g3_edge[i]*g3_edge[i]);
1107  }
1108 
1109  Vmath::Reverse(nquad0,&outarray[0],1,&outarray[0],1);
1110 
1111  break;
1112  case 3:
1113  Vmath::Vmul(nqtot,
1114  &(df[0][0]), 1,
1115  &jac[0], 1,
1116  &(tmp_gmat0[0]), 1);
1117  Vmath::Vmul(nqtot,
1118  &(df[2][0]),1,
1119  &jac[0], 1,
1120  &(tmp_gmat2[0]),1);
1121  QuadExp::v_GetEdgeInterpVals(
1122  edge, tmp_gmat0, g0_edge);
1123  QuadExp::v_GetEdgeInterpVals(
1124  edge, tmp_gmat2, g2_edge);
1125 
1126 
1127  for (i = 0; i < nquad1; ++i)
1128  {
1129  outarray[i] = sqrt(g0_edge[i]*g0_edge[i] +
1130  g2_edge[i]*g2_edge[i]);
1131  }
1132 
1133  Vmath::Reverse(nquad1,
1134  &outarray[0], 1,
1135  &outarray[0], 1);
1136 
1137  break;
1138  default:
1139  ASSERTL0(false,"edge value (< 3) is out of range");
1140  break;
1141  }
1142  }
1143  }
1144  else
1145  {
1146 
1147  switch (edge)
1148  {
1149  case 0:
1150 
1151 
1152 
1153  for (i = 0; i < nquad0; ++i)
1154  {
1155  outarray[i] = jac[0]*sqrt(df[1][0]*df[1][0] +
1156  df[3][0]*df[3][0]);
1157  }
1158  break;
1159  case 1:
1160  for (i = 0; i < nquad1; ++i)
1161  {
1162  outarray[i] = jac[0]*sqrt(df[0][0]*df[0][0] +
1163  df[2][0]*df[2][0]);
1164  }
1165  break;
1166  case 2:
1167  for (i = 0; i < nquad0; ++i)
1168  {
1169  outarray[i] = jac[0]*sqrt(df[1][0]*df[1][0] +
1170  df[3][0]*df[3][0]);
1171  }
1172  break;
1173  case 3:
1174  for (i = 0; i < nquad1; ++i)
1175  {
1176  outarray[i] = jac[0]*sqrt(df[0][0]*df[0][0] +
1177  df[2][0]*df[2][0]);
1178  }
1179  break;
1180  default:
1181  ASSERTL0(false,"edge value (< 3) is out of range");
1182  break;
1183  }
1184  }
1185  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::LibUtilities::PointsKeyVector
std::vector< PointsKey > PointsKeyVector
Definition: Points.h:220
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::LocalRegions::QuadExp::v_GetEdgeInterpVals
virtual void v_GetEdgeInterpVals(const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: QuadExp.cpp:814
Nektar::LibUtilities::eGaussGaussLegendre
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:47
Vmath::Reverse
void Reverse(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1071
Nektar::StdRegions::StdExpansion::GetPointsType
LibUtilities::PointsType GetPointsType(const int dir) const
This function returns the type of quadrature points used in the dir direction.
Definition: StdExpansion.h:216
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vmul
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:169
StdRegions::Orientation Nektar::LocalRegions::QuadExp::v_GetEorient ( int  edge)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 1549 of file QuadExp.cpp.

References Nektar::LocalRegions::Expansion::m_geom.

1550  {
1551  return m_geom->GetEorient(edge);
1552  }
Nektar::LocalRegions::Expansion::m_geom
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
DNekScalMatSharedPtr Nektar::LocalRegions::QuadExp::v_GetLocMatrix ( const MatrixKey &  mkey)
protectedvirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 2049 of file QuadExp.cpp.

References m_matrixManager.

2050  {
2051  return m_matrixManager[mkey];
2052  }
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
DNekScalBlkMatSharedPtr Nektar::LocalRegions::QuadExp::v_GetLocStaticCondMatrix ( const MatrixKey &  mkey)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 2055 of file QuadExp.cpp.

References m_staticCondMatrixManager.

2057  {
2058  return m_staticCondMatrixManager[mkey];
2059  }
Nektar::LocalRegions::QuadExp::m_staticCondMatrixManager
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: QuadExp.h:286
const SpatialDomains::GeomFactorsSharedPtr & Nektar::LocalRegions::QuadExp::v_GetMetricInfo ( ) const
protectedvirtual

Reimplemented from Nektar::LocalRegions::Expansion.

Definition at line 1470 of file QuadExp.cpp.

References Nektar::LocalRegions::Expansion::m_metricinfo.

1471  {
1472  return m_metricinfo;
1473  }
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
int Nektar::LocalRegions::QuadExp::v_GetNumPoints ( const int  dir) const
protectedvirtual

Definition at line 1568 of file QuadExp.cpp.

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

1569  {
1570  return GetNumPoints(dir);
1571  }
Nektar::StdRegions::StdExpansion::GetNumPoints
int GetNumPoints(const int dir) const
This function returns the number of quadrature points in the dir direction.
Definition: StdExpansion.h:229
StdRegions::StdExpansionSharedPtr Nektar::LocalRegions::QuadExp::v_GetStdExp ( void  ) const
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 611 of file QuadExp.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.

612  {
613  return MemoryManager<StdRegions::StdQuadExp>
614  ::AllocateSharedPtr(m_base[0]->GetBasisKey(),
615  m_base[1]->GetBasisKey());
616  }
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
void Nektar::LocalRegions::QuadExp::v_GetTracePhysVals ( const int  edge,
const StdRegions::StdExpansionSharedPtr &  EdgeExp,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
StdRegions::Orientation  orient 
)
protectedvirtual

Definition at line 739 of file QuadExp.cpp.

References v_GetEdgePhysVals().

745  {
746  v_GetEdgePhysVals(edge,EdgeExp,inarray,outarray);
747  }
Nektar::LocalRegions::QuadExp::v_GetEdgePhysVals
virtual void v_GetEdgePhysVals(const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Extract the physical values along edge edge from inarray into outarray following the local edge orien...
Definition: QuadExp.cpp:674
void Nektar::LocalRegions::QuadExp::v_HelmholtzMatrixOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2127 of file QuadExp.cpp.

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

2131  {
2132  QuadExp::HelmholtzMatrixOp_MatFree(inarray, outarray, mkey);
2133  }
Nektar::StdRegions::StdExpansion::HelmholtzMatrixOp_MatFree
void HelmholtzMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:1531
NekDouble Nektar::LocalRegions::QuadExp::v_Integral ( const Array< OneD, const NekDouble > &  inarray)
protectedvirtual

Integrates the specified function over the domain.

See also
StdRegions::StdExpansion::Integral.

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 83 of file QuadExp.cpp.

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

86  {
87  int nquad0 = m_base[0]->GetNumPoints();
88  int nquad1 = m_base[1]->GetNumPoints();
89  Array<OneD, const NekDouble> jac = m_metricinfo->GetJac(GetPointsKeys());
90  NekDouble ival;
91  Array<OneD,NekDouble> tmp(nquad0*nquad1);
92 
93  // multiply inarray with Jacobian
94  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
95  {
96  Vmath::Vmul(nquad0*nquad1, jac, 1, inarray, 1, tmp, 1);
97  }
98  else
99  {
100  Vmath::Smul(nquad0*nquad1, jac[0], inarray, 1, tmp, 1);
101  }
102 
103  // call StdQuadExp version;
104  ival = StdQuadExp::v_Integral(tmp);
105  return ival;
106  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:199
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vmul
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:169
void Nektar::LocalRegions::QuadExp::v_IProductWRTBase ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

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

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

where

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

which can be implemented as

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

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 408 of file QuadExp.cpp.

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

411  {
412  if (m_base[0]->Collocation() && m_base[1]->Collocation())
413  {
414  MultiplyByQuadratureMetric(inarray,outarray);
415  }
416  else
417  {
418  IProductWRTBase_SumFac(inarray,outarray);
419  }
420  }
Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:942
Nektar::StdRegions::StdExpansion::IProductWRTBase_SumFac
void IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
Definition: StdExpansion.h:1414
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
void Nektar::LocalRegions::QuadExp::v_IProductWRTBase_MatOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 462 of file QuadExp.cpp.

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

465  {
466  int nq = GetTotPoints();
467  MatrixKey
468  iprodmatkey(StdRegions::eIProductWRTBase,DetShapeType(),*this);
469  DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
470 
471  Blas::Dgemv('N',m_ncoeffs,nq,iprodmat->Scale(),
472  (iprodmat->GetOwnedMatrix())->GetPtr().get(),
473  m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
474 
475  }
Nektar::StdRegions::StdExpansion::DetShapeType
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:470
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::StdRegions::StdExpansion::GetTotPoints
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
Definition: StdExpansion.h:141
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
void Nektar::LocalRegions::QuadExp::v_IProductWRTBase_SumFac ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
bool  multiplybyweights = true 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 432 of file QuadExp.cpp.

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

436  {
437  int nquad0 = m_base[0]->GetNumPoints();
438  int nquad1 = m_base[1]->GetNumPoints();
439  int order0 = m_base[0]->GetNumModes();
440 
441  if(multiplybyweights)
442  {
443  Array<OneD,NekDouble> tmp(nquad0*nquad1+nquad1*order0);
444  Array<OneD,NekDouble> wsp(tmp+nquad0*nquad1);
445 
446  MultiplyByQuadratureMetric(inarray,tmp);
447  StdQuadExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),
448  m_base[1]->GetBdata(),
449  tmp,outarray,wsp,true,true);
450  }
451  else
452  {
453  Array<OneD,NekDouble> wsp(nquad1*order0);
454 
455  StdQuadExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),
456  m_base[1]->GetBdata(),
457  inarray,outarray,wsp,true,true);
458  }
459  }
Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:942
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
void Nektar::LocalRegions::QuadExp::v_IProductWRTDerivBase ( const int  dir,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 423 of file QuadExp.cpp.

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

427  {
428  IProductWRTDerivBase_SumFac(dir,inarray,outarray);
429  }
Nektar::StdRegions::StdExpansion::IProductWRTDerivBase_SumFac
void IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:1467
void Nektar::LocalRegions::QuadExp::v_IProductWRTDerivBase_MatOp ( const int  dir,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 534 of file QuadExp.cpp.

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

538  {
539  int nq = GetTotPoints();
540  StdRegions::MatrixType mtype = StdRegions::eIProductWRTDerivBase0;
541 
542  switch (dir)
543  {
544  case 0:
545  {
546  mtype = StdRegions::eIProductWRTDerivBase0;
547  }
548  break;
549  case 1:
550  {
551  mtype = StdRegions::eIProductWRTDerivBase1;
552  }
553  break;
554  case 2:
555  {
556  mtype = StdRegions::eIProductWRTDerivBase2;
557  }
558  break;
559  default:
560  {
561  ASSERTL1(false,"input dir is out of range");
562  }
563  break;
564  }
565 
566  MatrixKey iprodmatkey(mtype,DetShapeType(),*this);
567  DNekScalMatSharedPtr iprodmat = m_matrixManager[iprodmatkey];
568 
569  Blas::Dgemv('N', m_ncoeffs, nq, iprodmat->Scale(),
570  (iprodmat->GetOwnedMatrix())->GetPtr().get(),
571  m_ncoeffs, inarray.get(), 1, 0.0, outarray.get(), 1);
572  }
Nektar::StdRegions::StdExpansion::DetShapeType
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:470
Nektar::DNekScalMatSharedPtr
boost::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:73
Nektar::StdRegions::StdExpansion::GetTotPoints
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
Definition: StdExpansion.h:141
Nektar::LocalRegions::QuadExp::m_matrixManager
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:285
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
void Nektar::LocalRegions::QuadExp::v_IProductWRTDerivBase_SumFac ( const int  dir,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 478 of file QuadExp.cpp.

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

482  {
483  ASSERTL1((dir==0) || (dir==1) || (dir==2),
484  "Invalid direction.");
485  ASSERTL1((dir==2) ? (m_geom->GetCoordim() ==3):true,
486  "Invalid direction.");
487 
488  int nquad0 = m_base[0]->GetNumPoints();
489  int nquad1 = m_base[1]->GetNumPoints();
490  int nqtot = nquad0*nquad1;
491  int nmodes0 = m_base[0]->GetNumModes();
492 
493  const Array<TwoD, const NekDouble>& df = m_metricinfo->GetDerivFactors(GetPointsKeys());
494 
495  Array<OneD, NekDouble> tmp1(2*nqtot+m_ncoeffs+nmodes0*nquad1);
496  Array<OneD, NekDouble> tmp2(tmp1 + nqtot);
497  Array<OneD, NekDouble> tmp3(tmp1 + 2*nqtot);
498  Array<OneD, NekDouble> tmp4(tmp1 + 2*nqtot+m_ncoeffs);
499 
500  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
501  {
502  Vmath::Vmul(nqtot,
503  &df[2*dir][0], 1,
504  inarray.get(), 1,
505  tmp1.get(), 1);
506  Vmath::Vmul(nqtot,
507  &df[2*dir+1][0], 1,
508  inarray.get(), 1,
509  tmp2.get(),1);
510  }
511  else
512  {
513  Vmath::Smul(nqtot,
514  df[2*dir][0], inarray.get(), 1,
515  tmp1.get(), 1);
516  Vmath::Smul(nqtot,
517  df[2*dir+1][0], inarray.get(), 1,
518  tmp2.get(), 1);
519  }
520 
521  MultiplyByQuadratureMetric(tmp1,tmp1);
522  MultiplyByQuadratureMetric(tmp2,tmp2);
523 
524  IProductWRTBase_SumFacKernel(
525  m_base[0]->GetDbdata(), m_base[1]->GetBdata(),
526  tmp1, tmp3, tmp4, false, true);
527  IProductWRTBase_SumFacKernel(
528  m_base[0]->GetBdata() , m_base[1]->GetDbdata(),
529  tmp2, outarray, tmp4, true, false);
530  Vmath::Vadd(m_ncoeffs, tmp3, 1, outarray, 1, outarray, 1);
531  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:942
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::LocalRegions::Expansion::m_geom
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:199
Nektar::StdRegions::StdExpansion2D::IProductWRTBase_SumFacKernel
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)
Definition: StdExpansion2D.cpp:200
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vadd
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:285
Vmath::Vmul
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:169
void Nektar::LocalRegions::QuadExp::v_LaplacianMatrixOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2076 of file QuadExp.cpp.

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

2080  {
2081  QuadExp::LaplacianMatrixOp_MatFree(inarray, outarray, mkey);
2082  }
Nektar::StdRegions::StdExpansion::LaplacianMatrixOp_MatFree
void LaplacianMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:1489
void Nektar::LocalRegions::QuadExp::v_LaplacianMatrixOp ( const int  k1,
const int  k2,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2085 of file QuadExp.cpp.

2091  {
2092  StdExpansion::LaplacianMatrixOp_MatFree(
2093  k1, k2, inarray, outarray, mkey);
2094  }
void Nektar::LocalRegions::QuadExp::v_LaplacianMatrixOp_MatFree_Kernel ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
Array< OneD, NekDouble > &  wsp 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 2211 of file QuadExp.cpp.

References ASSERTL1, Nektar::LocalRegions::Expansion::ComputeLaplacianMetric(), Nektar::LocalRegions::eMetricLaplacian00, Nektar::LocalRegions::eMetricLaplacian01, Nektar::LocalRegions::eMetricLaplacian11, Nektar::StdRegions::StdExpansion2D::IProductWRTBase_SumFacKernel(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metrics, Nektar::StdRegions::StdExpansion::m_ncoeffs, Vmath::Vadd(), and Vmath::Vvtvvtp().

2215  {
2216  if (m_metrics.count(eMetricLaplacian00) == 0)
2217  {
2218  ComputeLaplacianMetric();
2219  }
2220 
2221  int nquad0 = m_base[0]->GetNumPoints();
2222  int nquad1 = m_base[1]->GetNumPoints();
2223  int nqtot = nquad0*nquad1;
2224  int nmodes0 = m_base[0]->GetNumModes();
2225  int nmodes1 = m_base[1]->GetNumModes();
2226  int wspsize = max(max(max(nqtot,m_ncoeffs),nquad1*nmodes0),nquad0*nmodes1);
2227 
2228  ASSERTL1(wsp.num_elements() >= 3*wspsize,
2229  "Workspace is of insufficient size.");
2230 
2231  const Array<OneD, const NekDouble>& base0 = m_base[0]->GetBdata();
2232  const Array<OneD, const NekDouble>& base1 = m_base[1]->GetBdata();
2233  const Array<OneD, const NekDouble>& dbase0 = m_base[0]->GetDbdata();
2234  const Array<OneD, const NekDouble>& dbase1 = m_base[1]->GetDbdata();
2235  const Array<OneD, const NekDouble>& metric00 = m_metrics[eMetricLaplacian00];
2236  const Array<OneD, const NekDouble>& metric01 = m_metrics[eMetricLaplacian01];
2237  const Array<OneD, const NekDouble>& metric11 = m_metrics[eMetricLaplacian11];
2238 
2239  // Allocate temporary storage
2240  Array<OneD,NekDouble> wsp0(wsp);
2241  Array<OneD,NekDouble> wsp1(wsp+wspsize);
2242  Array<OneD,NekDouble> wsp2(wsp+2*wspsize);
2243 
2244  StdExpansion2D::PhysTensorDeriv(inarray,wsp1,wsp2);
2245 
2246  // wsp0 = k = g0 * wsp1 + g1 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
2247  // wsp2 = l = g1 * wsp1 + g2 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
2248  // where g0, g1 and g2 are the metric terms set up in the GeomFactors class
2249  // especially for this purpose
2250  Vmath::Vvtvvtp(nqtot,&metric00[0],1,&wsp1[0],1,&metric01[0],1,&wsp2[0],1,&wsp0[0],1);
2251  Vmath::Vvtvvtp(nqtot,&metric01[0],1,&wsp1[0],1,&metric11[0],1,&wsp2[0],1,&wsp2[0],1);
2252 
2253  // outarray = m = (D_xi1 * B)^T * k
2254  // wsp1 = n = (D_xi2 * B)^T * l
2255  IProductWRTBase_SumFacKernel(dbase0,base1,wsp0,outarray,wsp1,false,true);
2256  IProductWRTBase_SumFacKernel(base0,dbase1,wsp2,wsp1, wsp0,true,false);
2257 
2258  // outarray = outarray + wsp1
2259  // = L * u_hat
2260  Vmath::Vadd(m_ncoeffs,wsp1.get(),1,outarray.get(),1,outarray.get(),1);
2261  }
Vmath::Vvtvvtp
void Vvtvvtp(int n, const T *v, int incv, const T *w, int incw, const T *x, int incx, const T *y, int incy, T *z, int incz)
vvtvvtp (vector times vector plus vector times vector):
Definition: Vmath.cpp:523
Nektar::StdRegions::StdExpansion2D::IProductWRTBase_SumFacKernel
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)
Definition: StdExpansion2D.cpp:200
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Vmath::Vadd
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:285
void Nektar::LocalRegions::QuadExp::v_MassLevelCurvatureMatrixOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 2117 of file QuadExp.cpp.

2121  {
2122  StdExpansion::MassLevelCurvatureMatrixOp_MatFree(
2123  inarray, outarray, mkey);
2124  }
void Nektar::LocalRegions::QuadExp::v_MassMatrixOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2067 of file QuadExp.cpp.

2071  {
2072  StdExpansion::MassMatrixOp_MatFree(inarray, outarray, mkey);
2073  }
void Nektar::LocalRegions::QuadExp::v_NormVectorIProductWRTBase ( const Array< OneD, const NekDouble > &  Fx,
const Array< OneD, const NekDouble > &  Fy,
const Array< OneD, const NekDouble > &  Fz,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 575 of file QuadExp.cpp.

References Nektar::SpatialDomains::eDeformed, Nektar::LocalRegions::Expansion2D::GetLeftAdjacentElementExp(), Nektar::LocalRegions::Expansion2D::GetLeftAdjacentElementFace(), Nektar::StdRegions::StdExpansion::IProductWRTBase(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Svtsvtp(), Vmath::Svtvp(), Vmath::Vvtvp(), and Vmath::Vvtvvtp().

580  {
581  int nq = m_base[0]->GetNumPoints()*m_base[1]->GetNumPoints();
582  Array<OneD, NekDouble> Fn(nq);
583 
584  const Array<OneD, const Array<OneD, NekDouble> > &normals =
585  GetLeftAdjacentElementExp()->GetFaceNormal(
586  GetLeftAdjacentElementFace());
587 
588  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
589  {
590  Vmath::Vvtvvtp(nq,&normals[0][0],1,&Fx[0],1,
591  &normals[1][0],1,&Fy[0],1,&Fn[0],1);
592  Vmath::Vvtvp (nq,&normals[2][0],1,&Fz[0],1,&Fn[0],1,&Fn[0],1);
593  }
594  else
595  {
596  Vmath::Svtsvtp(nq,normals[0][0],&Fx[0],1,
597  normals[1][0],&Fy[0],1,&Fn[0],1);
598  Vmath::Svtvp (nq,normals[2][0],&Fz[0],1,&Fn[0],1,&Fn[0],1);
599  }
600 
601  IProductWRTBase(Fn,outarray);
602  }
Nektar::StdRegions::StdExpansion::IProductWRTBase
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 expa...
Definition: StdExpansion.h:629
Vmath::Svtvp
void Svtvp(int n, const T alpha, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
svtvp (scalar times vector plus vector): z = alpha*x + y
Definition: Vmath.cpp:471
Vmath::Vvtvp
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:428
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::LocalRegions::Expansion2D::GetLeftAdjacentElementExp
Expansion3DSharedPtr GetLeftAdjacentElementExp() const
Definition: Expansion2D.h:223
Vmath::Vvtvvtp
void Vvtvvtp(int n, const T *v, int incv, const T *w, int incw, const T *x, int incx, const T *y, int incy, T *z, int incz)
vvtvvtp (vector times vector plus vector times vector):
Definition: Vmath.cpp:523
Vmath::Svtsvtp
void Svtsvtp(int n, const T alpha, const T *x, int incx, const T beta, const T *y, int incy, T *z, int incz)
vvtvvtp (scalar times vector plus scalar times vector):
Definition: Vmath.cpp:577
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
void Nektar::LocalRegions::QuadExp::v_NormVectorIProductWRTBase ( const Array< OneD, const Array< OneD, NekDouble > > &  Fvec,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 604 of file QuadExp.cpp.

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

607  {
608  NormVectorIProductWRTBase(Fvec[0], Fvec[1], Fvec[2], outarray);
609  }
Nektar::StdRegions::StdExpansion::NormVectorIProductWRTBase
void NormVectorIProductWRTBase(const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:727
void Nektar::LocalRegions::QuadExp::v_PhysDeriv ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  out_d0,
Array< OneD, NekDouble > &  out_d1,
Array< OneD, NekDouble > &  out_d2 = NullNekDouble1DArray 
)
protectedvirtual

Calculate the derivative of the physical points.

For quadrilateral region can use the Tensor_Deriv function defined under StdExpansion.

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 109 of file QuadExp.cpp.

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

114  {
115  int nquad0 = m_base[0]->GetNumPoints();
116  int nquad1 = m_base[1]->GetNumPoints();
117  int nqtot = nquad0*nquad1;
118  const Array<TwoD, const NekDouble>& df = m_metricinfo->GetDerivFactors(GetPointsKeys());
119  Array<OneD,NekDouble> diff0(2*nqtot);
120  Array<OneD,NekDouble> diff1(diff0+nqtot);
121 
122  StdQuadExp::v_PhysDeriv(inarray, diff0, diff1);
123 
124  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
125  {
126  if (out_d0.num_elements())
127  {
128  Vmath::Vmul (nqtot, df[0], 1, diff0, 1, out_d0, 1);
129  Vmath::Vvtvp (nqtot, df[1], 1, diff1, 1, out_d0, 1,
130  out_d0,1);
131  }
132 
133  if(out_d1.num_elements())
134  {
135  Vmath::Vmul (nqtot,df[2],1,diff0,1, out_d1, 1);
136  Vmath::Vvtvp (nqtot,df[3],1,diff1,1, out_d1, 1, out_d1,1);
137  }
138 
139  if (out_d2.num_elements())
140  {
141  Vmath::Vmul (nqtot,df[4],1,diff0,1, out_d2, 1);
142  Vmath::Vvtvp (nqtot,df[5],1,diff1,1, out_d2, 1, out_d2,1);
143  }
144  }
145  else // regular geometry
146  {
147  if (out_d0.num_elements())
148  {
149  Vmath::Smul (nqtot, df[0][0], diff0, 1, out_d0, 1);
150  Blas::Daxpy (nqtot, df[1][0], diff1, 1, out_d0, 1);
151  }
152 
153  if (out_d1.num_elements())
154  {
155  Vmath::Smul (nqtot, df[2][0], diff0, 1, out_d1, 1);
156  Blas::Daxpy (nqtot, df[3][0], diff1, 1, out_d1, 1);
157  }
158 
159  if (out_d2.num_elements())
160  {
161  Vmath::Smul (nqtot, df[4][0], diff0, 1, out_d2, 1);
162  Blas::Daxpy (nqtot, df[5][0], diff1, 1, out_d2, 1);
163  }
164  }
165  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Vmath::Vvtvp
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:428
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:199
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vmul
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:169
void Nektar::LocalRegions::QuadExp::v_PhysDeriv ( const int  dir,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  out_d0 
)
protectedvirtual

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

See also
StdRegions::StdExpansion::PhysDeriv

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 168 of file QuadExp.cpp.

References ASSERTL1, Nektar::NullNekDouble1DArray, and v_PhysDeriv().

172  {
173  switch (dir)
174  {
175  case 0:
176  {
177  v_PhysDeriv(inarray, outarray, NullNekDouble1DArray,
178  NullNekDouble1DArray);
179  }
180  break;
181  case 1:
182  {
183  v_PhysDeriv(inarray, NullNekDouble1DArray, outarray,
184  NullNekDouble1DArray);
185  }
186  break;
187  case 2:
188  {
189  v_PhysDeriv(inarray, NullNekDouble1DArray,
190  NullNekDouble1DArray, outarray);
191  }
192  break;
193  default:
194  {
195  ASSERTL1(false,"input dir is out of range");
196  }
197  break;
198  }
199  }
Nektar::NullNekDouble1DArray
static Array< OneD, NekDouble > NullNekDouble1DArray
Definition: SharedArray.hpp:784
Nektar::LocalRegions::QuadExp::v_PhysDeriv
virtual void v_PhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray)
Calculate the derivative of the physical points.
Definition: QuadExp.cpp:109
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
void Nektar::LocalRegions::QuadExp::v_PhysDirectionalDeriv ( const Array< OneD, const NekDouble > &  inarray,
const Array< OneD, const NekDouble > &  direction,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Physical derivative along a direction vector.

See also
StdRegions::StdExpansion::PhysDirectionalDeriv

D_v = d/dx_v^s + d/dx_v^r

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 202 of file QuadExp.cpp.

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

206  {
207  int nquad0 = m_base[0]->GetNumPoints();
208  int nquad1 = m_base[1]->GetNumPoints();
209  int nqtot = nquad0*nquad1;
210 
211  const Array<TwoD, const NekDouble>& df = m_metricinfo->GetDerivFactors(GetPointsKeys());
212 
213  Array<OneD,NekDouble> diff0(2*nqtot);
214  Array<OneD,NekDouble> diff1(diff0+nqtot);
215 
216  StdQuadExp::v_PhysDeriv(inarray, diff0, diff1);
217 
218  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
219  {
220  Array<OneD, Array<OneD, NekDouble> > tangmat(2);
221 
222  // d/dx_v^s = v_x*ds/dx + v_y*ds/dy + v_z*dx/dz
223  for (int i=0; i< 2; ++i)
224  {
225  tangmat[i] = Array<OneD, NekDouble>(nqtot,0.0);
226  for (int k=0; k<(m_geom->GetCoordim()); ++k)
227  {
228  Vmath::Vvtvp(nqtot,
229  &df[2*k+i][0], 1,
230  &direction[k*nqtot], 1,
231  &tangmat[i][0], 1,
232  &tangmat[i][0], 1);
233  }
234  }
235 
236  /// D_v = d/dx_v^s + d/dx_v^r
237  if (out.num_elements())
238  {
239  Vmath::Vmul (nqtot,
240  &tangmat[0][0], 1,
241  &diff0[0], 1,
242  &out[0], 1);
243  Vmath::Vvtvp (nqtot,
244  &tangmat[1][0], 1,
245  &diff1[0], 1,
246  &out[0], 1,
247  &out[0], 1);
248  }
249 
250  }
251  else
252  {
253  ASSERTL1(m_metricinfo->GetGtype() ==
254  SpatialDomains::eDeformed,"Wrong route");
255  }
256  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Vmath::Vvtvp
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:428
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Nektar::LocalRegions::Expansion::m_geom
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
ASSERTL1
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vmul
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:169
NekDouble Nektar::LocalRegions::QuadExp::v_PhysEvaluate ( const Array< OneD, const NekDouble > &  coords,
const Array< OneD, const NekDouble > &  physvals 
)
protectedvirtual

This function evaluates the expansion at a single (arbitrary) point of the domain.

This function is a wrapper around the virtual function v_PhysEvaluate()

Based on the value of the expansion at the quadrature points, this function calculates the value of the expansion at an arbitrary single points (with coordinates $ \mathbf{x_c}$ given by the pointer coords). This operation, equivalent to

\[ u(\mathbf{x_c}) = \sum_p \phi_p(\mathbf{x_c}) \hat{u}_p \]

is evaluated using Lagrangian interpolants through the quadrature points:

\[ u(\mathbf{x_c}) = \sum_p h_p(\mathbf{x_c}) u_p\]

This function requires that the physical value array $\mathbf{u}$ (implemented as the attribute #m_phys) is set.

Parameters
coordsthe coordinates of the single point
Returns
returns the value of the expansion at the single point

Reimplemented from Nektar::StdRegions::StdExpansion2D.

Definition at line 658 of file QuadExp.cpp.

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

661  {
662  Array<OneD,NekDouble> Lcoord = Array<OneD, NekDouble>(2);
663 
664  ASSERTL0(m_geom,"m_geom not defined");
665  m_geom->GetLocCoords(coord,Lcoord);
666 
667  return StdQuadExp::v_PhysEvaluate(Lcoord, physvals);
668  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
Nektar::LocalRegions::Expansion::m_geom
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:125
void Nektar::LocalRegions::QuadExp::v_ReduceOrderCoeffs ( int  numMin,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2161 of file QuadExp.cpp.

References Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eOrtho_A, Nektar::StdRegions::StdExpansion::GetBasisType(), Nektar::LibUtilities::InterpCoeff2D(), Nektar::StdRegions::StdExpansion::m_base, Vmath::Vcopy(), and Vmath::Zero().

2165  {
2166  int n_coeffs = inarray.num_elements();
2167 
2168  Array<OneD, NekDouble> coeff (n_coeffs);
2169  Array<OneD, NekDouble> coeff_tmp(n_coeffs, 0.0);
2170  Array<OneD, NekDouble> tmp, tmp2;
2171 
2172  int nmodes0 = m_base[0]->GetNumModes();
2173  int nmodes1 = m_base[1]->GetNumModes();
2174  int numMax = nmodes0;
2175 
2176  Vmath::Vcopy(n_coeffs,inarray,1,coeff_tmp,1);
2177 
2178  const LibUtilities::PointsKey Pkey0(
2179  nmodes0, LibUtilities::eGaussLobattoLegendre);
2180  const LibUtilities::PointsKey Pkey1(
2181  nmodes1, LibUtilities::eGaussLobattoLegendre);
2182  LibUtilities::BasisKey b0(
2183  m_base[0]->GetBasisType(), nmodes0, Pkey0);
2184  LibUtilities::BasisKey b1(
2185  m_base[1]->GetBasisType(), nmodes1, Pkey1);
2186  LibUtilities::BasisKey bortho0(
2187  LibUtilities::eOrtho_A, nmodes0, Pkey0);
2188  LibUtilities::BasisKey bortho1(
2189  LibUtilities::eOrtho_A, nmodes1, Pkey1);
2190 
2191  LibUtilities::InterpCoeff2D(
2192  b0, b1, coeff_tmp, bortho0, bortho1, coeff);
2193 
2194  Vmath::Zero(n_coeffs, coeff_tmp, 1);
2195 
2196  int cnt = 0;
2197  for (int i = 0; i < numMin+1; ++i)
2198  {
2199  Vmath::Vcopy(numMin,
2200  tmp = coeff+cnt,1,
2201  tmp2 = coeff_tmp+cnt,1);
2202 
2203  cnt = i*numMax;
2204  }
2205 
2206  LibUtilities::InterpCoeff2D(
2207  bortho0, bortho1, coeff_tmp,
2208  b0, b1, outarray);
2209  }
Nektar::LibUtilities::eOrtho_A
Principle Orthogonal Functions .
Definition: BasisType.h:46
Nektar::StdRegions::StdExpansion::GetBasisType
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
Definition: StdExpansion.h:165
Vmath::Zero
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:359
Nektar::StdRegions::StdExpansion::m_base
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1422
Nektar::LibUtilities::InterpCoeff2D
void InterpCoeff2D(const BasisKey &fbasis0, const BasisKey &fbasis1, const Array< OneD, const NekDouble > &from, const BasisKey &tbasis0, const BasisKey &tbasis1, Array< OneD, NekDouble > &to)
Definition: InterpCoeff.cpp:73
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
Nektar::LibUtilities::eGaussLobattoLegendre
1D Gauss-Lobatto-Legendre quadrature points
Definition: PointsType.h:50
NekDouble Nektar::LocalRegions::QuadExp::v_StdPhysEvaluate ( const Array< OneD, const NekDouble > &  Lcoord,
const Array< OneD, const NekDouble > &  physvals 
)
protectedvirtual

Given the local cartesian coordinate Lcoord evaluate the value of physvals at this point by calling through to the StdExpansion method

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 650 of file QuadExp.cpp.

653  {
654  // Evaluate point in local (eta) coordinates.
655  return StdQuadExp::v_PhysEvaluate(Lcoord,physvals);
656  }
void Nektar::LocalRegions::QuadExp::v_SVVLaplacianFilter ( Array< OneD, NekDouble > &  array,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2302 of file QuadExp.cpp.

References Nektar::SpatialDomains::eDeformed, Vmath::Fill(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Vdiv(), Vmath::Vmul(), and Vmath::Vsqrt().

2305  {
2306  int nq = GetTotPoints();
2307 
2308  // Calculate sqrt of the Jacobian
2309  Array<OneD, const NekDouble> jac =
2310  m_metricinfo->GetJac(GetPointsKeys());
2311  Array<OneD, NekDouble> sqrt_jac(nq);
2312  if (m_metricinfo->GetGtype() == SpatialDomains::eDeformed)
2313  {
2314  Vmath::Vsqrt(nq,jac,1,sqrt_jac,1);
2315  }
2316  else
2317  {
2318  Vmath::Fill(nq,sqrt(jac[0]),sqrt_jac,1);
2319  }
2320 
2321  // Multiply array by sqrt(Jac)
2322  Vmath::Vmul(nq,sqrt_jac,1,array,1,array,1);
2323 
2324  // Apply std region filter
2325  StdQuadExp::v_SVVLaplacianFilter( array, mkey);
2326 
2327  // Divide by sqrt(Jac)
2328  Vmath::Vdiv(nq,array,1,sqrt_jac,1,array,1);
2329  }
Nektar::StdRegions::StdExpansion::GetPointsKeys
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1328
Vmath::Vsqrt
void Vsqrt(int n, const T *x, const int incx, T *y, const int incy)
sqrt y = sqrt(x)
Definition: Vmath.cpp:394
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:46
Nektar::LocalRegions::Expansion::m_metricinfo
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:126
Vmath::Vdiv
void Vdiv(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:227
Nektar::StdRegions::StdExpansion::GetTotPoints
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
Definition: StdExpansion.h:141
Nektar::SpatialDomains::eDeformed
Geometry is curved or has non-constant factors.
Definition: SpatialDomains.hpp:57
Vmath::Vmul
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:169
void Nektar::LocalRegions::QuadExp::v_WeakDerivMatrixOp ( const int  i,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdQuadExp.

Definition at line 2097 of file QuadExp.cpp.

2102  {
2103  StdExpansion::WeakDerivMatrixOp_MatFree(i, inarray, outarray, mkey);
2104  }
void Nektar::LocalRegions::QuadExp::v_WeakDirectionalDerivMatrixOp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const StdRegions::StdMatrixKey &  mkey 
)
protectedvirtual

Reimplemented from Nektar::StdRegions::StdExpansion.

Definition at line 2107 of file QuadExp.cpp.

2111  {
2112  StdExpansion::WeakDirectionalDerivMatrixOp_MatFree(
2113  inarray, outarray, mkey);
2114  }

Member Data Documentation

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

Definition at line 285 of file QuadExp.h.

Referenced by CreateMatrix(), v_FwdTrans(), v_GetEdgeInterpVals(), v_GetLocMatrix(), v_IProductWRTBase_MatOp(), and v_IProductWRTDerivBase_MatOp().

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

Definition at line 286 of file QuadExp.h.

Referenced by v_DropLocStaticCondMatrix(), v_FwdTrans_BndConstrained(), and v_GetLocStaticCondMatrix().

Generated on Wed Mar 23 2016 20:14:32 for Nektar++ by   doxygen 1.8.8