#include <QuadExp.h>
Inheritance diagram for Nektar::LocalRegions::QuadExp:
Public Member Functions | |
| QuadExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const SpatialDomains::Geometry2DSharedPtr &geom) | |
| Constructor using BasisKey class for quadrature points and order definition. More... | |
| QuadExp (const QuadExp &T) | |
| ~QuadExp () override=default | |
 Public Member Functions inherited from Nektar::StdRegions::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)=default | |
| ~StdQuadExp () override=default | |
| Public Member Functions inherited from Nektar::StdRegions::StdExpansion2D | |
| StdExpansion2D (int numcoeffs, const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb) | |
| StdExpansion2D ()=default | |
| StdExpansion2D (const StdExpansion2D &T)=default | |
| ~StdExpansion2D () override=default | |
| void | PhysTensorDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray_d0, Array< OneD, NekDouble > &outarray_d1) |
| Calculate the 2D derivative in the local tensor/collapsed coordinate at the physical points. More... | |
| NekDouble | Integral (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &w0, const Array< OneD, const NekDouble > &w1) |
| NekDouble | BaryTensorDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) |
| void | BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true) |
| void | IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0=true, bool doCheckCollDir1=true) |
| Public Member Functions inherited from Nektar::StdRegions::StdExpansion | |
| StdExpansion () | |
| Default Constructor. More... | |
| StdExpansion (const int numcoeffs, const int numbases, const LibUtilities::BasisKey &Ba=LibUtilities::NullBasisKey, const LibUtilities::BasisKey &Bb=LibUtilities::NullBasisKey, const LibUtilities::BasisKey &Bc=LibUtilities::NullBasisKey) | |
| Constructor. More... | |
| StdExpansion (const StdExpansion &T) | |
| Copy Constructor. More... | |
| virtual | ~StdExpansion () |
| Destructor. More... | |
| int | GetNumBases () const |
| This function returns the number of 1D bases used in the expansion. More... | |
| const Array< OneD, const LibUtilities::BasisSharedPtr > & | GetBase () const |
| This function gets the shared point to basis. More... | |
| const LibUtilities::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 | GetTraceNcoeffs (const int i) const |
| This function returns the number of expansion coefficients belonging to the i-th trace. More... | |
| int | GetTraceIntNcoeffs (const int i) const |
| int | GetTraceNumPoints (const int i) const |
| This function returns the number of quadrature points belonging to the i-th trace. More... | |
| const LibUtilities::BasisKey | GetTraceBasisKey (const int i, int k=-1, bool UseGLL=false) const |
| This function returns the basis key belonging to the i-th trace. More... | |
| LibUtilities::PointsKey | GetTracePointsKey (const int i, int k=-1) const |
| This function returns the basis key belonging to the i-th trace. More... | |
| int | NumBndryCoeffs (void) const |
| int | NumDGBndryCoeffs (void) const |
| const LibUtilities::PointsKey | GetNodalPointsKey () const |
| This function returns the type of expansion Nodal point type if defined. More... | |
| int | GetNtraces () const |
| Returns the number of trace elements connected to this element. More... | |
| LibUtilities::ShapeType | DetShapeType () const |
| This function returns the shape of the expansion domain. More... | |
| std::shared_ptr< StdExpansion > | GetStdExp () const |
| std::shared_ptr< StdExpansion > | GetLinStdExp (void) const |
| int | GetShapeDimension () const |
| bool | IsBoundaryInteriorExpansion () const |
| bool | IsNodalNonTensorialExp () |
| void | NodalToModal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| This function performs the Backward transformation from coefficient space to physical space. More... | |
| void | FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| This function performs the Forward transformation from physical space to coefficient space. More... | |
| void | FwdTransBndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| NekDouble | Integral (const Array< OneD, const NekDouble > &inarray) |
| This function integrates the specified function over the domain. More... | |
| void | FillMode (const int mode, Array< OneD, NekDouble > &outarray) |
| This function fills the array outarray with the mode-th mode of the expansion. More... | |
| void | IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| this function calculates the inner product of a given function f with the different modes of the expansion More... | |
| void | IProductWRTBase (const Array< OneD, const NekDouble > &base, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int coll_check) |
| void | IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | IProductWRTDirectionalDerivBase (const Array< OneD, const NekDouble > &direction, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| int | GetElmtId () |
| Get the element id of this expansion when used in a list by returning value of m_elmt_id. More... | |
| void | SetElmtId (const int id) |
| Set the element id of this expansion when used in a list by returning value of m_elmt_id. More... | |
| void | GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2=NullNekDouble1DArray, Array< OneD, NekDouble > &coords_3=NullNekDouble1DArray) |
| this function returns the physical coordinates of the quadrature points of the expansion More... | |
| void | GetCoord (const Array< OneD, const NekDouble > &Lcoord, Array< OneD, NekDouble > &coord) |
| given the coordinates of a point of the element in the local collapsed coordinate system, this function calculates the physical coordinates of the point More... | |
| DNekMatSharedPtr | GetStdMatrix (const StdMatrixKey &mkey) |
| DNekBlkMatSharedPtr | GetStdStaticCondMatrix (const StdMatrixKey &mkey) |
| void | NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray) |
| void | NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, NekDouble > &Fy, Array< OneD, NekDouble > &outarray) |
| void | NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, const Array< OneD, const NekDouble > &Fz, Array< OneD, NekDouble > &outarray) |
| void | NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray) |
| DNekScalBlkMatSharedPtr | GetLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey) |
| void | DropLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey) |
| int | CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset) |
| NekDouble | StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals) |
| int | GetCoordim () |
| void | GetBoundaryMap (Array< OneD, unsigned int > &outarray) |
| void | GetInteriorMap (Array< OneD, unsigned int > &outarray) |
| int | GetVertexMap (const int localVertexId, bool useCoeffPacking=false) |
| void | GetTraceToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards, int P=-1, int Q=-1) |
| void | GetTraceCoeffMap (const unsigned int traceid, Array< OneD, unsigned int > &maparray) |
| void | GetElmtTraceToTraceMap (const unsigned int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards, int P=-1, int Q=-1) |
| void | GetTraceInteriorToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation traceOrient=eForwards) |
| void | GetTraceNumModes (const int tid, int &numModes0, int &numModes1, const Orientation traceOrient=eDir1FwdDir1_Dir2FwdDir2) |
| void | MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| DNekMatSharedPtr | CreateGeneralMatrix (const StdMatrixKey &mkey) |
| this function generates the mass matrix \(\mathbf{M}[i][j] =
\int \phi_i(\mathbf{x}) \phi_j(\mathbf{x}) d\mathbf{x}\) More... | |
| void | GeneralMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey) |
| void | ExponentialFilter (Array< OneD, NekDouble > &array, const NekDouble alpha, const NekDouble exponent, const NekDouble cutoff) |
| void | LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | WeakDirectionalDerivMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | MassLevelCurvatureMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LinearAdvectionMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LinearAdvectionDiffusionReactionMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true) |
| void | HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| DNekMatSharedPtr | GenMatrix (const StdMatrixKey &mkey) |
| void | PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) |
| void | PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | PhysDeriv_s (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_ds) |
| void | PhysDeriv_n (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_dn) |
| void | PhysDirectionalDeriv (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &direction, Array< OneD, NekDouble > &outarray) |
| void | StdPhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1=NullNekDouble1DArray, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) |
| void | StdPhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| NekDouble | PhysEvaluate (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals) |
| This function evaluates the expansion at a single (arbitrary) point of the domain. More... | |
| NekDouble | PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) |
| This function evaluates the first derivative of the expansion at a single (arbitrary) point of the domain. More... | |
| NekDouble | PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs, std::array< NekDouble, 6 > &secondOrderDerivs) |
| NekDouble | PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals) |
| This function evaluates the expansion at a single (arbitrary) point of the domain. More... | |
| NekDouble | PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode) |
This function evaluates the basis function mode mode at a point coords of the domain. More... | |
| void | LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta) |
| Convert local cartesian coordinate xi into local collapsed coordinates eta. More... | |
| void | LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi) |
| Convert local collapsed coordinates eta into local cartesian coordinate xi. More... | |
| void | PhysInterp (std::shared_ptr< StdExpansion > fromExp, const Array< OneD, const NekDouble > &fromData, Array< OneD, NekDouble > &toData) |
| interpolate from one set of quadrature points available from FromExp to the set of quadrature points in the current expansion. If the points are the same this routine will just copy the data More... | |
| virtual int | v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset) |
| virtual void | v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray) |
| virtual void | v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray) |
| virtual void | v_NormVectorIProductWRTBase (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, const Array< OneD, const NekDouble > &Fz, Array< OneD, NekDouble > &outarray) |
| virtual void | v_NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray) |
| virtual DNekScalBlkMatSharedPtr | v_GetLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey) |
| virtual void | v_DropLocStaticCondMatrix (const LocalRegions::MatrixKey &mkey) |
| NekDouble | Linf (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray) |
| Function to evaluate the discrete \( L_\infty\) error \( |\epsilon|_\infty = \max |u - u_{exact}|\) where \(
u_{exact}\) is given by the array sol. More... | |
| NekDouble | L2 (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray) |
| Function to evaluate the discrete \( L_2\) error, \( | \epsilon |_{2} = \left [ \int^1_{-1} [u - u_{exact}]^2
dx \right]^{1/2} d\xi_1 \) where \( u_{exact}\) is given by the array sol. More... | |
| NekDouble | H1 (const Array< OneD, const NekDouble > &phys, const Array< OneD, const NekDouble > &sol=NullNekDouble1DArray) |
| Function to evaluate the discrete \( H^1\) error, \( | \epsilon |^1_{2} = \left [ \int^1_{-1} [u -
u_{exact}]^2 + \nabla(u - u_{exact})\cdot\nabla(u -
u_{exact})\cdot dx \right]^{1/2} d\xi_1 \) where \(
u_{exact}\) is given by the array sol. More... | |
| const LibUtilities::PointsKeyVector | GetPointsKeys () const |
| DNekMatSharedPtr | BuildInverseTransformationMatrix (const DNekScalMatSharedPtr &m_transformationmatrix) |
| void | PhysInterpToSimplexEquiSpaced (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, int npset=-1) |
| This function performs an interpolation from the physical space points provided at input into an array of equispaced points which are not the collapsed coordinate. So for a tetrahedron you will only get a tetrahedral number of values. More... | |
| void | GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true) |
| This function provides the connectivity of local simplices (triangles or tets) to connect the equispaced data points provided by PhysInterpToSimplexEquiSpaced. More... | |
| void | EquiSpacedToCoeffs (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| This function performs a projection/interpolation from the equispaced points sometimes used in post-processing onto the coefficient space. More... | |
| template<class T > | |
| std::shared_ptr< T > | as () |
| void | IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true) |
| void | GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat) |
| Public Member Functions inherited from Nektar::LocalRegions::Expansion2D | |
| Expansion2D (SpatialDomains::Geometry2DSharedPtr pGeom) | |
| ~Expansion2D () override=default | |
| DNekScalMatSharedPtr | CreateMatrix (const MatrixKey &mkey) |
| void | SetTraceToGeomOrientation (Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &inout) |
| Array< OneD, unsigned int > | GetTraceInverseBoundaryMap (int eid) |
| void | AddNormTraceInt (const int dir, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &edgeCoeffs, Array< OneD, NekDouble > &outarray) |
| void | AddNormTraceInt (const int dir, Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs) |
| void | AddEdgeBoundaryInt (const int edge, ExpansionSharedPtr &EdgeExp, Array< OneD, NekDouble > &edgePhys, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap) |
| void | AddHDGHelmholtzEdgeTerms (const NekDouble tau, const int edge, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, NekDouble > &edgePhys, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray) |
| void | AddHDGHelmholtzTraceTerms (const NekDouble tau, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray) |
| SpatialDomains::Geometry2DSharedPtr | GetGeom2D () const |
| void | ReOrientEdgePhysMap (const int nvert, const StdRegions::Orientation orient, const int nq0, Array< OneD, int > &idmap) |
| DNekMatSharedPtr | v_GenMatrix (const StdRegions::StdMatrixKey &mkey) override |
| void | v_GenTraceExp (const int traceid, ExpansionSharedPtr &exp) override |
| Public Member Functions inherited from Nektar::LocalRegions::Expansion | |
| Expansion (SpatialDomains::GeometrySharedPtr pGeom) | |
| Expansion (const Expansion &pSrc) | |
| ~Expansion () override | |
| void | SetTraceExp (const int traceid, ExpansionSharedPtr &f) |
| ExpansionSharedPtr | GetTraceExp (const int traceid) |
| DNekScalMatSharedPtr | GetLocMatrix (const LocalRegions::MatrixKey &mkey) |
| void | DropLocMatrix (const LocalRegions::MatrixKey &mkey) |
| DNekScalMatSharedPtr | GetLocMatrix (const StdRegions::MatrixType mtype, const StdRegions::ConstFactorMap &factors=StdRegions::NullConstFactorMap, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap) |
| SpatialDomains::GeometrySharedPtr | GetGeom () const |
| void | Reset () |
| IndexMapValuesSharedPtr | CreateIndexMap (const IndexMapKey &ikey) |
| DNekScalBlkMatSharedPtr | CreateStaticCondMatrix (const MatrixKey &mkey) |
| const SpatialDomains::GeomFactorsSharedPtr & | GetMetricInfo () const |
| DNekMatSharedPtr | BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType) |
| DNekMatSharedPtr | BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd) |
| void | ExtractDataToCoeffs (const NekDouble *data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble *coeffs, std::vector< LibUtilities::BasisType > &fromType) |
| void | AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray) |
| void | AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) |
| void | AddFaceNormBoundaryInt (const int face, const std::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) |
| void | DGDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &coeffs, Array< OneD, NekDouble > &outarray) |
| NekDouble | VectorFlux (const Array< OneD, Array< OneD, NekDouble > > &vec) |
| void | NormalTraceDerivFactors (Array< OneD, Array< OneD, NekDouble > > &factors, Array< OneD, Array< OneD, NekDouble > > &d0factors, Array< OneD, Array< OneD, NekDouble > > &d1factors) |
| IndexMapValuesSharedPtr | GetIndexMap (const IndexMapKey &ikey) |
| void | AlignVectorToCollapsedDir (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) |
| ExpansionSharedPtr | GetLeftAdjacentElementExp () const |
| ExpansionSharedPtr | GetRightAdjacentElementExp () const |
| int | GetLeftAdjacentElementTrace () const |
| int | GetRightAdjacentElementTrace () const |
| void | SetAdjacentElementExp (int traceid, ExpansionSharedPtr &e) |
| StdRegions::Orientation | GetTraceOrient (int trace) |
| void | SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | DivideByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| Divided by the metric jacobi and quadrature weights. More... | |
| void | GetTraceQFactors (const int trace, Array< OneD, NekDouble > &outarray) |
| Extract the metric factors to compute the contravariant fluxes along edge edge and stores them into outarray following the local edge orientation (i.e. anticlockwise convention). More... | |
| void | GetTracePhysVals (const int trace, const StdRegions::StdExpansionSharedPtr &TraceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient=StdRegions::eNoOrientation) |
| void | GetTracePhysMap (const int edge, Array< OneD, int > &outarray) |
| void | ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1) |
| const NormalVector & | GetTraceNormal (const int id) |
| void | ComputeTraceNormal (const int id) |
| const Array< OneD, const NekDouble > & | GetPhysNormals (void) |
| void | SetPhysNormals (Array< OneD, const NekDouble > &normal) |
| void | SetUpPhysNormals (const int trace) |
| void | AddRobinMassMatrix (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat) |
| void | TraceNormLen (const int traceid, NekDouble &h, NekDouble &p) |
| void | AddRobinTraceContribution (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, const Array< OneD, NekDouble > &incoeffs, Array< OneD, NekDouble > &coeffs) |
| const Array< OneD, const NekDouble > & | GetElmtBndNormDirElmtLen (const int nbnd) const |
| void | StdDerivBaseOnTraceMat (Array< OneD, DNekMatSharedPtr > &DerivMat) |
Protected Member Functions | |
| NekDouble | v_Integral (const Array< OneD, const NekDouble > &inarray) override |
| Integrates the specified function over the domain. More... | |
| NekDouble | v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals) override |
| void | v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) override |
| Calculate the derivative of the physical points. More... | |
| void | v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Calculate the derivative of the physical points in a given direction. More... | |
| void | v_PhysDirectionalDeriv (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &direction, Array< OneD, NekDouble > &out) override |
| Physical derivative along a direction vector. More... | |
| void | v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Transform a given function from physical quadrature space to coefficient space. More... | |
| void | v_FwdTransBndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray. More... | |
| void | v_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true) override |
| void | v_IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_AlignVectorToCollapsedDir (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override |
| 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) override |
| void | v_NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray) override |
| StdRegions::StdExpansionSharedPtr | v_GetStdExp (void) const override |
| StdRegions::StdExpansionSharedPtr | v_GetLinStdExp (void) const override |
| void | v_GetCoord (const Array< OneD, const NekDouble > &Lcoords, Array< OneD, NekDouble > &coords) override |
| void | v_GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3) override |
| NekDouble | v_PhysEvaluate (const Array< OneD, const NekDouble > &coord, const Array< OneD, const NekDouble > &physvals) override |
| This function evaluates the expansion at a single (arbitrary) point of the domain. More... | |
| NekDouble | v_PhysEvalFirstDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) override |
| void | v_GetTracePhysVals (const int edge, const StdRegions::StdExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient) override |
| void | v_GetTraceQFactors (const int edge, Array< OneD, NekDouble > &outarray) override |
| void | v_ComputeTraceNormal (const int edge) override |
| void | v_ExtractDataToCoeffs (const NekDouble *data, const std::vector< unsigned int > &nummodes, const int mode_offset, NekDouble *coeffs, std::vector< LibUtilities::BasisType > &fromType) override |
| StdRegions::Orientation | v_GetTraceOrient (int edge) override |
| void | v_GetTracePhysMap (const int edge, Array< OneD, int > &outarray) override |
| DNekMatSharedPtr | v_GenMatrix (const StdRegions::StdMatrixKey &mkey) override |
| DNekMatSharedPtr | v_CreateStdMatrix (const StdRegions::StdMatrixKey &mkey) override |
| DNekScalMatSharedPtr | v_GetLocMatrix (const MatrixKey &mkey) override |
| void | v_DropLocMatrix (const MatrixKey &mkey) override |
| DNekScalBlkMatSharedPtr | v_GetLocStaticCondMatrix (const MatrixKey &mkey) override |
| void | v_DropLocStaticCondMatrix (const MatrixKey &mkey) override |
| void | v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_WeakDirectionalDerivMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_MassLevelCurvatureMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp) override |
| void | v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_ComputeLaplacianMetric () override |
| void | v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdRegions::StdMatrixKey &mkey) override |
| void | v_NormalTraceDerivFactors (Array< OneD, Array< OneD, NekDouble > > &factors, Array< OneD, Array< OneD, NekDouble > > &d0factors, Array< OneD, Array< OneD, NekDouble > > &d1factors) override |
| : This method gets all of the factors which are required as part of the Gradient Jump Penalty (GJP) stabilisation and involves the product of the normal and geometric factors along the element trace. More... | |
| Protected Member Functions inherited from Nektar::StdRegions::StdQuadExp | |
| NekDouble | v_Integral (const Array< OneD, const NekDouble > &inarray) override |
| Integrates the specified function over the domain. More... | |
| void | v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) override |
| Calculate the derivative of the physical points. More... | |
| void | v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Calculate the derivative of the physical points in a given direction. More... | |
| void | v_StdPhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) override |
| void | v_StdPhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1) override |
| void | v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Transform a given function from physical quadrature space to coefficient space. More... | |
| void | v_FwdTransBndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray. More... | |
| void | v_IProductWRTBase_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true) override |
| void | v_IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1) override |
| void | v_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_LocCoordToLocCollapsed (const Array< OneD, const NekDouble > &xi, Array< OneD, NekDouble > &eta) override |
| void | v_LocCollapsedToLocCoord (const Array< OneD, const NekDouble > &eta, Array< OneD, NekDouble > &xi) override |
| void | v_FillMode (const int mode, Array< OneD, NekDouble > &array) override |
| Fill outarray with mode mode of expansion. More... | |
| int | v_GetNverts () const final |
| int | v_GetNtraces () const final |
| int | v_GetTraceNcoeffs (const int i) const final |
| int | v_GetTraceIntNcoeffs (const int i) const final |
| int | v_GetTraceNumPoints (const int i) const final |
| int | v_NumBndryCoeffs () const final |
| int | v_NumDGBndryCoeffs () const final |
| int | v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset) override |
| const LibUtilities::BasisKey | v_GetTraceBasisKey (const int i, const int j, bool UseGLL=false) const final |
| LibUtilities::ShapeType | v_DetShapeType () const final |
| bool | v_IsBoundaryInteriorExpansion () const override |
| void | v_GetCoords (Array< OneD, NekDouble > &coords_0, Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2) override |
| NekDouble | v_PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode) override |
This function evaluates the basis function mode mode at a point coords of the domain. More... | |
| NekDouble | v_PhysEvalFirstDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) override |
| void | v_GetBoundaryMap (Array< OneD, unsigned int > &outarray) override |
| void | v_GetInteriorMap (Array< OneD, unsigned int > &outarray) override |
| int | v_GetVertexMap (int localVertexId, bool useCoeffPacking=false) override |
| void | v_GetTraceCoeffMap (const unsigned int traceid, Array< OneD, unsigned int > &maparray) override |
| Get the map of the coefficient location to teh local trace coefficients. More... | |
| void | v_GetTraceInteriorToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation edgeOrient=eForwards) override |
| DNekMatSharedPtr | v_GenMatrix (const StdMatrixKey &mkey) override |
| DNekMatSharedPtr | v_CreateStdMatrix (const StdMatrixKey &mkey) override |
| void | v_MassMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override |
| void | v_LaplacianMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override |
| void | v_LaplacianMatrixOp (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override |
| void | v_WeakDerivMatrixOp (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override |
| void | v_HelmholtzMatrixOp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) override |
| void | v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey) override |
| void | v_ExponentialFilter (Array< OneD, NekDouble > &array, const NekDouble alpha, const NekDouble exponent, const NekDouble cutoff) override |
| void | v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| void | v_GetSimplexEquiSpacedConnectivity (Array< OneD, int > &conn, bool standard=true) override |
| Protected Member Functions inherited from Nektar::StdRegions::StdExpansion2D | |
| NekDouble | v_PhysEvaluate (const Array< OneD, const NekDouble > &coords, const Array< OneD, const NekDouble > &physvals) override |
| This function evaluates the expansion at a single (arbitrary) point of the domain. More... | |
| NekDouble | v_PhysEvaluateInterp (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals) override |
| 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)=0 |
| 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)=0 |
| void | v_LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdRegions::StdMatrixKey &mkey) override |
| void | v_GetTraceCoeffMap (const unsigned int traceid, Array< OneD, unsigned int > &maparray) override |
| void | v_GetElmtTraceToTraceMap (const unsigned int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient, int P, int Q) override |
| Determine the mapping to re-orientate the coefficients along the element trace (assumed to align with the standard element) into the orientation of the local trace given by edgeOrient. More... | |
| void | v_GetTraceToElementMap (const int eid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation edgeOrient=eForwards, int P=-1, int Q=-1) override |
| void | v_GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat) override |
| void | v_PhysInterp (std::shared_ptr< StdExpansion > fromExp, const Array< OneD, const NekDouble > &fromData, Array< OneD, NekDouble > &toData) override |
| Protected Member Functions inherited from Nektar::StdRegions::StdExpansion | |
| DNekMatSharedPtr | CreateStdMatrix (const StdMatrixKey &mkey) |
| DNekBlkMatSharedPtr | CreateStdStaticCondMatrix (const StdMatrixKey &mkey) |
| Create the static condensation of a matrix when using a boundary interior decomposition. More... | |
| void | BwdTrans_SumFac (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | IProductWRTDerivBase_SumFac (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | IProductWRTDirectionalDerivBase_SumFac (const Array< OneD, const NekDouble > &direction, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| void | GeneralMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | MassMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LaplacianMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp) |
| void | LaplacianMatrixOp_MatFree_GenericImpl (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LaplacianMatrixOp_MatFree (const int k1, const int k2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | WeakDerivMatrixOp_MatFree (const int i, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | WeakDirectionalDerivMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | MassLevelCurvatureMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LinearAdvectionMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | LinearAdvectionDiffusionReactionMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey, bool addDiffusionTerm=true) |
| void | HelmholtzMatrixOp_MatFree (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| void | HelmholtzMatrixOp_MatFree_GenericImpl (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey) |
| virtual void | v_SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| virtual NekDouble | v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals) |
| virtual void | v_GenStdMatBwdDeriv (const int dir, DNekMatSharedPtr &mat) |
| virtual void | v_MultiplyByStdQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| template<int DIR, bool DERIV = false, bool DERIV2 = false> | |
| NekDouble | BaryEvaluate (const NekDouble &coord, const NekDouble *physvals, NekDouble &deriv, NekDouble &deriv2) |
This function performs the barycentric interpolation of the polynomial stored in coord at a point physvals using barycentric interpolation weights in direction. More... | |
| template<int DIR> | |
| NekDouble | BaryEvaluateBasis (const NekDouble &coord, const int &mode) |
| template<int DIR, bool DERIV = false, bool DERIV2 = false> | |
| NekDouble | BaryEvaluate (const NekDouble &coord, const NekDouble *physvals) |
| Helper function to pass an unused value by reference into BaryEvaluate. More... | |
| template<int DIR, bool DERIV = false, bool DERIV2 = false> | |
| NekDouble | BaryEvaluate (const NekDouble &coord, const NekDouble *physvals, NekDouble &deriv) |
| Protected Member Functions inherited from Nektar::LocalRegions::Expansion2D | |
| void | v_DGDeriv (const int dir, const Array< OneD, const NekDouble > &incoeffs, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &edgeCoeffs, Array< OneD, NekDouble > &out_d) override |
| void | v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray) override |
| void | v_AddEdgeNormBoundaryInt (const int edge, const ExpansionSharedPtr &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) override |
| void | v_AddRobinMassMatrix (const int edgeid, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat) override |
| void | v_AddRobinTraceContribution (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, const Array< OneD, NekDouble > &incoeffs, Array< OneD, NekDouble > &coeffs) override |
| DNekMatSharedPtr | v_BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd) override |
| void | v_ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1) override |
| void | v_SetUpPhysNormals (const int edge) override |
| NekDouble | v_VectorFlux (const Array< OneD, Array< OneD, NekDouble > > &vec) override |
| void | v_TraceNormLen (const int traceid, NekDouble &h, NekDouble &p) override |
| Protected Member Functions inherited from Nektar::LocalRegions::Expansion | |
| void | ComputeLaplacianMetric () |
| void | ComputeQuadratureMetric () |
| void | ComputeGmatcdotMF (const Array< TwoD, const NekDouble > &df, const Array< OneD, const NekDouble > &direction, Array< OneD, Array< OneD, NekDouble > > &dfdir) |
| Array< OneD, NekDouble > | GetMF (const int dir, const int shapedim, const StdRegions::VarCoeffMap &varcoeffs) |
| Array< OneD, NekDouble > | GetMFDiv (const int dir, const StdRegions::VarCoeffMap &varcoeffs) |
| Array< OneD, NekDouble > | GetMFMag (const int dir, const StdRegions::VarCoeffMap &varcoeffs) |
| void | v_MultiplyByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| virtual void | v_DivideByQuadratureMetric (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| virtual void | v_ComputeLaplacianMetric () |
| int | v_GetCoordim () const override |
| void | v_GetCoords (Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3) override |
| virtual DNekScalMatSharedPtr | v_GetLocMatrix (const LocalRegions::MatrixKey &mkey) |
| virtual void | v_DropLocMatrix (const LocalRegions::MatrixKey &mkey) |
| virtual DNekMatSharedPtr | v_BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType) |
| virtual DNekMatSharedPtr | v_BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd) |
| virtual void | v_ExtractDataToCoeffs (const NekDouble *data, const std::vector< unsigned int > &nummodes, const int nmodes_offset, NekDouble *coeffs, std::vector< LibUtilities::BasisType > &fromType) |
| virtual void | v_AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray) |
| virtual void | v_AddEdgeNormBoundaryInt (const int edge, const std::shared_ptr< Expansion > &EdgeExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) |
| virtual void | v_AddFaceNormBoundaryInt (const int face, const std::shared_ptr< Expansion > &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) |
| virtual void | v_DGDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &EdgeExp, Array< OneD, Array< OneD, NekDouble > > &coeffs, Array< OneD, NekDouble > &outarray) |
| virtual NekDouble | v_VectorFlux (const Array< OneD, Array< OneD, NekDouble > > &vec) |
| virtual void | v_NormalTraceDerivFactors (Array< OneD, Array< OneD, NekDouble > > &factors, Array< OneD, Array< OneD, NekDouble > > &d0factors, Array< OneD, Array< OneD, NekDouble > > &d1factors) |
| virtual void | v_AlignVectorToCollapsedDir (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) |
| virtual StdRegions::Orientation | v_GetTraceOrient (int trace) |
| void | v_SetCoeffsToOrientation (StdRegions::Orientation dir, Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| virtual void | v_GetTraceQFactors (const int trace, Array< OneD, NekDouble > &outarray) |
| virtual void | v_GetTracePhysVals (const int trace, const StdRegions::StdExpansionSharedPtr &TraceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient) |
| virtual void | v_GetTracePhysMap (const int edge, Array< OneD, int > &outarray) |
| virtual void | v_ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1=-1) |
| virtual void | v_ComputeTraceNormal (const int id) |
| virtual const Array< OneD, const NekDouble > & | v_GetPhysNormals () |
| virtual void | v_SetPhysNormals (Array< OneD, const NekDouble > &normal) |
| virtual void | v_SetUpPhysNormals (const int id) |
| virtual void | v_AddRobinMassMatrix (const int face, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat) |
| virtual void | v_AddRobinTraceContribution (const int traceid, const Array< OneD, const NekDouble > &primCoeffs, const Array< OneD, NekDouble > &incoeffs, Array< OneD, NekDouble > &coeffs) |
| virtual void | v_TraceNormLen (const int traceid, NekDouble &h, NekDouble &p) |
| virtual void | v_GenTraceExp (const int traceid, ExpansionSharedPtr &exp) |
Private Member Functions | |
| void | GetEdgeInterpVals (const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) |
| QuadExp () | |
Detailed Description
Constructor & Destructor Documentation
◆ QuadExp() [1/3]
| Nektar::LocalRegions::QuadExp::QuadExp | ( | const LibUtilities::BasisKey & | Ba, |
| const LibUtilities::BasisKey & | Bb, | ||
| const SpatialDomains::Geometry2DSharedPtr & | geom | ||
| ) |
Constructor using BasisKey class for quadrature points and order definition.
Definition at line 47 of file QuadExp.cpp.
50 : StdExpansion(Ba.GetNumModes() * Bb.GetNumModes(), 2, Ba, Bb),
51 StdExpansion2D(Ba.GetNumModes() * Bb.GetNumModes(), Ba, Bb),
53 m_matrixManager(
55 std::string("QuadExpMatrix")),
57 this, std::placeholders::_1),
58 std::string("QuadExpStaticCondMatrix"))
59{
60}
DNekScalMatSharedPtr CreateMatrix(const MatrixKey &mkey)
Definition: Expansion2D.cpp:55
Expansion2D(SpatialDomains::Geometry2DSharedPtr pGeom)
Definition: Expansion2D.cpp:50
DNekScalBlkMatSharedPtr CreateStaticCondMatrix(const MatrixKey &mkey)
Definition: Expansion.cpp:272
Expansion(SpatialDomains::GeometrySharedPtr pGeom)
Definition: Expansion.cpp:43
LibUtilities::NekManager< MatrixKey, DNekScalMat, MatrixKey::opLess > m_matrixManager
Definition: QuadExp.h:244
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: QuadExp.h:246
StdExpansion2D()=default
StdQuadExp(const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb)
Constructor using BasisKey class for quadrature points and order definition.
Definition: StdQuadExp.cpp:47
◆ QuadExp() [2/3]
| Nektar::LocalRegions::QuadExp::QuadExp | ( | const QuadExp & | T | ) |
Definition at line 62 of file QuadExp.cpp.
◆ ~QuadExp()
|
overridedefault |
◆ QuadExp() [3/3]
|
private |
Member Function Documentation
◆ GetEdgeInterpVals()
|
private |
Definition at line 660 of file QuadExp.cpp.
663{
664 int i;
667
670
673
675
676 switch (edge)
677 {
678 case 0:
679 {
680 for (i = 0; i < nq0; i++)
681 {
682 outarray[i] = Blas::Ddot(
683 nq1, mat_gauss->GetOwnedMatrix()->GetPtr().data(), 1,
684 &inarray[i], nq0);
685 }
686 break;
687 }
688 case 1:
689 {
690 for (i = 0; i < nq1; i++)
691 {
692 outarray[i] = Blas::Ddot(
693 nq0, mat_gauss->GetOwnedMatrix()->GetPtr().data(), 1,
694 &inarray[i * nq0], 1);
695 }
696 break;
697 }
698 case 2:
699 {
700 for (i = 0; i < nq0; i++)
701 {
702 outarray[i] = Blas::Ddot(
703 nq1, mat_gauss->GetOwnedMatrix()->GetPtr().data(), 1,
704 &inarray[i], nq0);
705 }
706 break;
707 }
708 case 3:
709 {
710 for (i = 0; i < nq1; i++)
711 {
712 outarray[i] = Blas::Ddot(
713 nq0, mat_gauss->GetOwnedMatrix()->GetPtr().data(), 1,
714 &inarray[i * nq0], 1);
715 }
716 break;
717 }
718 default:
720 break;
721 }
722}
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:370
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1193
static double Ddot(const int &n, const double *x, const int &incx, const double *y, const int &incy)
BLAS level 1: output = .
Definition: Blas.hpp:163
std::map< ConstFactorType, NekDouble > ConstFactorMap
Definition: StdRegions.hpp:430
StdRegions::ConstFactorMap factors
Definition: TestVarcoeffHashing.cpp:52
std::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:66
References ASSERTL0, Blas::Ddot(), Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::StdRegions::eFactorGaussEdge, Nektar::StdRegions::eInterpGauss, Nektar::VarcoeffHashingTest::factors, Nektar::StdRegions::StdExpansion::m_base, and m_matrixManager.
Referenced by v_ComputeTraceNormal(), v_GetTracePhysVals(), and v_GetTraceQFactors().
◆ v_AlignVectorToCollapsedDir()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 452 of file QuadExp.cpp.
455{
458 "Invalid direction.");
459
462 int nqtot = nquad0 * nquad1;
464
465 const Array<TwoD, const NekDouble> &df =
467
468 Array<OneD, NekDouble> tmp1 = outarray[0];
469 Array<OneD, NekDouble> tmp2 = outarray[1];
470 Array<OneD, NekDouble> tmp3(m_ncoeffs);
471 Array<OneD, NekDouble> tmp4(nmodes0 * nquad1);
472
474 {
475 Vmath::Vmul(nqtot, &df[2 * dir][0], 1, inarray.data(), 1, tmp1.data(),
476 1);
477 Vmath::Vmul(nqtot, &df[2 * dir + 1][0], 1, inarray.data(), 1,
478 tmp2.data(), 1);
479 }
480 else
481 {
482 Vmath::Smul(nqtot, df[2 * dir][0], inarray.data(), 1, tmp1.data(), 1);
483 Vmath::Smul(nqtot, df[2 * dir + 1][0], inarray.data(), 1, tmp2.data(),
484 1);
485 }
486}
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode....
Definition: ErrorUtil.hpp:242
SpatialDomains::GeometrySharedPtr m_geom
Definition: Expansion.h:273
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:274
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1115
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.hpp:72
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*x.
Definition: Vmath.hpp:100
References ASSERTL1, Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_geom, Nektar::LocalRegions::Expansion::m_metricinfo, Nektar::StdRegions::StdExpansion::m_ncoeffs, Vmath::Smul(), and Vmath::Vmul().
Referenced by v_IProductWRTDerivBase_SumFac().
◆ v_ComputeLaplacianMetric()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1553 of file QuadExp.cpp.
1554{
1556 {
1557 ComputeQuadratureMetric();
1558 }
1559
1562 const unsigned int dim = 2;
1567
1568 const Array<TwoD, const NekDouble> gmat =
1570 for (unsigned int i = 0; i < dim; ++i)
1571 {
1572 for (unsigned int j = i; j < dim; ++j)
1573 {
1574 m_metrics[m[i][j]] = Array<OneD, NekDouble>(nqtot);
1576 {
1577 Vmath::Vcopy(nqtot, &gmat[i * dim + j][0], 1,
1578 &m_metrics[m[i][j]][0], 1);
1579 }
1580 else
1581 {
1583 1);
1584 }
1586 }
1587 }
1588}
void ComputeQuadratureMetric()
Definition: Expansion.cpp:458
int GetTotPoints() const
This function returns the total number of quadrature points used in the element.
Definition: StdExpansion.h:134
void MultiplyByQuadratureMetric(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:732
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.hpp:54
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.hpp:825
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().
◆ v_ComputeTraceNormal()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 980 of file QuadExp.cpp.
981{
982 int i;
984 GetGeom()->GetMetricInfo();
985 SpatialDomains::GeomType type = geomFactors->GetGtype();
986
988 for (i = 0; i < ptsKeys.size(); ++i)
989 {
990 // Need at least 2 points for computing normals
992 {
993 LibUtilities::PointsKey pKey(2, ptsKeys[i].GetPointsType());
994 ptsKeys[i] = pKey;
995 }
996 }
997
998 const Array<TwoD, const NekDouble> &df =
999 geomFactors->GetDerivFactors(ptsKeys);
1000 const Array<OneD, const NekDouble> &jac = geomFactors->GetJac(ptsKeys);
1001
1002 // The points of normals should follow trace basis, not local basis.
1003 LibUtilities::BasisKey tobasis = GetTraceBasisKey(edge);
1004
1005 int nqe = tobasis.GetNumPoints();
1007
1008 m_traceNormals[edge] = Array<OneD, Array<OneD, NekDouble>>(vCoordDim);
1009 Array<OneD, Array<OneD, NekDouble>> &normal = m_traceNormals[edge];
1010 for (i = 0; i < vCoordDim; ++i)
1011 {
1012 normal[i] = Array<OneD, NekDouble>(nqe);
1013 }
1014
1015 size_t nqb = nqe;
1016 size_t nbnd = edge;
1017 m_elmtBndNormDirElmtLen[nbnd] = Array<OneD, NekDouble>{nqb, 0.0};
1018 Array<OneD, NekDouble> &length = m_elmtBndNormDirElmtLen[nbnd];
1019
1020 // Regular geometry case
1022 (type == SpatialDomains::eMovingRegular))
1023 {
1024 NekDouble fac;
1025 // Set up normals
1026 switch (edge)
1027 {
1028 case 0:
1029 for (i = 0; i < vCoordDim; ++i)
1030 {
1031 Vmath::Fill(nqe, -df[2 * i + 1][0], normal[i], 1);
1032 }
1033 break;
1034 case 1:
1035 for (i = 0; i < vCoordDim; ++i)
1036 {
1037 Vmath::Fill(nqe, df[2 * i][0], normal[i], 1);
1038 }
1039 break;
1040 case 2:
1041 for (i = 0; i < vCoordDim; ++i)
1042 {
1043 Vmath::Fill(nqe, df[2 * i + 1][0], normal[i], 1);
1044 }
1045 break;
1046 case 3:
1047 for (i = 0; i < vCoordDim; ++i)
1048 {
1049 Vmath::Fill(nqe, -df[2 * i][0], normal[i], 1);
1050 }
1051 break;
1052 default:
1054 }
1055
1056 // normalise
1057 fac = 0.0;
1058 for (i = 0; i < vCoordDim; ++i)
1059 {
1060 fac += normal[i][0] * normal[i][0];
1061 }
1062 fac = 1.0 / sqrt(fac);
1063
1064 Vmath::Fill(nqb, fac, length, 1);
1065
1066 for (i = 0; i < vCoordDim; ++i)
1067 {
1068 Vmath::Smul(nqe, fac, normal[i], 1, normal[i], 1);
1069 }
1070 }
1071 else // Set up deformed normals
1072 {
1073 int j;
1074
1075 int nquad0 = ptsKeys[0].GetNumPoints();
1076 int nquad1 = ptsKeys[1].GetNumPoints();
1077
1078 LibUtilities::PointsKey from_key;
1079
1080 Array<OneD, NekDouble> normals(vCoordDim * max(nquad0, nquad1), 0.0);
1081 Array<OneD, NekDouble> edgejac(vCoordDim * max(nquad0, nquad1), 0.0);
1082
1083 // Extract Jacobian along edges and recover local
1084 // derivates (dx/dr) for polynomial interpolation by
1085 // multiplying m_gmat by jacobian
1086
1087 // Implementation for all the basis except Gauss points
1090 {
1091 switch (edge)
1092 {
1093 case 0:
1094 for (j = 0; j < nquad0; ++j)
1095 {
1096 edgejac[j] = jac[j];
1097 for (i = 0; i < vCoordDim; ++i)
1098 {
1099 normals[i * nquad0 + j] =
1100 -df[2 * i + 1][j] * edgejac[j];
1101 }
1102 }
1103 from_key = ptsKeys[0];
1104 break;
1105 case 1:
1106 for (j = 0; j < nquad1; ++j)
1107 {
1108 edgejac[j] = jac[nquad0 * j + nquad0 - 1];
1109 for (i = 0; i < vCoordDim; ++i)
1110 {
1111 normals[i * nquad1 + j] =
1112 df[2 * i][nquad0 * j + nquad0 - 1] * edgejac[j];
1113 }
1114 }
1115 from_key = ptsKeys[1];
1116 break;
1117 case 2:
1118 for (j = 0; j < nquad0; ++j)
1119 {
1120 edgejac[j] = jac[nquad0 * (nquad1 - 1) + j];
1121 for (i = 0; i < vCoordDim; ++i)
1122 {
1123 normals[i * nquad0 + j] =
1124 (df[2 * i + 1][nquad0 * (nquad1 - 1) + j]) *
1125 edgejac[j];
1126 }
1127 }
1128 from_key = ptsKeys[0];
1129 break;
1130 case 3:
1131 for (j = 0; j < nquad1; ++j)
1132 {
1133 edgejac[j] = jac[nquad0 * j];
1134 for (i = 0; i < vCoordDim; ++i)
1135 {
1136 normals[i * nquad1 + j] =
1137 -df[2 * i][nquad0 * j] * edgejac[j];
1138 }
1139 }
1140 from_key = ptsKeys[1];
1141 break;
1142 default:
1144 }
1145 }
1146 else
1147 {
1148 int nqtot = nquad0 * nquad1;
1149 Array<OneD, NekDouble> tmp_gmat(nqtot, 0.0);
1150 Array<OneD, NekDouble> tmp_gmat_edge(nqe, 0.0);
1151
1152 switch (edge)
1153 {
1154 case 0:
1155 for (j = 0; j < nquad0; ++j)
1156 {
1157 for (i = 0; i < vCoordDim; ++i)
1158 {
1159 Vmath::Vmul(nqtot, &(df[2 * i + 1][0]), 1, &jac[0],
1160 1, &(tmp_gmat[0]), 1);
1161 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1162 tmp_gmat_edge);
1163 normals[i * nquad0 + j] = -tmp_gmat_edge[j];
1164 }
1165 }
1166 from_key = ptsKeys[0];
1167 break;
1168 case 1:
1169 for (j = 0; j < nquad1; ++j)
1170 {
1171 for (i = 0; i < vCoordDim; ++i)
1172 {
1173 Vmath::Vmul(nqtot, &(df[2 * i][0]), 1, &jac[0], 1,
1174 &(tmp_gmat[0]), 1);
1175 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1176 tmp_gmat_edge);
1177 normals[i * nquad1 + j] = tmp_gmat_edge[j];
1178 }
1179 }
1180 from_key = ptsKeys[1];
1181 break;
1182 case 2:
1183 for (j = 0; j < nquad0; ++j)
1184 {
1185 for (i = 0; i < vCoordDim; ++i)
1186 {
1187 Vmath::Vmul(nqtot, &(df[2 * i + 1][0]), 1, &jac[0],
1188 1, &(tmp_gmat[0]), 1);
1189 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1190 tmp_gmat_edge);
1191 normals[i * nquad0 + j] = tmp_gmat_edge[j];
1192 }
1193 }
1194 from_key = ptsKeys[0];
1195 break;
1196 case 3:
1197 for (j = 0; j < nquad1; ++j)
1198 {
1199 for (i = 0; i < vCoordDim; ++i)
1200 {
1201 Vmath::Vmul(nqtot, &(df[2 * i][0]), 1, &jac[0], 1,
1202 &(tmp_gmat[0]), 1);
1203 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1204 tmp_gmat_edge);
1205 normals[i * nquad1 + j] = -tmp_gmat_edge[j];
1206 }
1207 }
1208 from_key = ptsKeys[1];
1209 break;
1210 default:
1212 }
1213 }
1214
1215 int nq = from_key.GetNumPoints();
1216 Array<OneD, NekDouble> work(nqe, 0.0);
1217
1218 // interpolate Jacobian and invert
1219 LibUtilities::Interp1D(from_key, jac, tobasis.GetPointsKey(), work);
1220 Vmath::Sdiv(nqe, 1.0, &work[0], 1, &work[0], 1);
1221
1222 // interpolate
1224 {
1225 LibUtilities::Interp1D(from_key, &normals[i * nq],
1226 tobasis.GetPointsKey(), &normal[i][0]);
1227 Vmath::Vmul(nqe, work, 1, normal[i], 1, normal[i], 1);
1228 }
1229
1230 // normalise normal vectors
1231 Vmath::Zero(nqe, work, 1);
1233 {
1234 Vmath::Vvtvp(nqe, normal[i], 1, normal[i], 1, work, 1, work, 1);
1235 }
1236
1237 Vmath::Vsqrt(nqe, work, 1, work, 1);
1238 Vmath::Sdiv(nqe, 1.0, work, 1, work, 1);
1239
1240 Vmath::Vcopy(nqb, work, 1, length, 1);
1241
1243 {
1244 Vmath::Vmul(nqe, normal[i], 1, work, 1, normal[i], 1);
1245 }
1246 }
1248 {
1249 for (i = 0; i < vCoordDim; ++i)
1250 {
1252 {
1253 Vmath::Reverse(nqe, normal[i], 1, normal[i], 1);
1254 }
1255 }
1256 }
1257}
std::map< int, NormalVector > m_traceNormals
Definition: Expansion.h:276
std::map< int, Array< OneD, NekDouble > > m_elmtBndNormDirElmtLen
the element length in each element boundary(Vertex, edge or face) normal direction calculated based o...
Definition: Expansion.h:286
SpatialDomains::GeometrySharedPtr GetGeom() const
Definition: Expansion.cpp:167
void GetEdgeInterpVals(const int edge, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: QuadExp.cpp:660
LibUtilities::PointsType GetPointsType(const int dir) const
This function returns the type of quadrature points used in the dir direction.
Definition: StdExpansion.h:205
int GetNumPoints(const int dir) const
This function returns the number of quadrature points in the dir direction.
Definition: StdExpansion.h:218
const LibUtilities::BasisKey GetTraceBasisKey(const int i, int k=-1, bool UseGLL=false) const
This function returns the basis key belonging to the i-th trace.
Definition: StdExpansion.h:301
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:47
std::vector< PointsKey > PointsKeyVector
Definition: Points.h:231
std::shared_ptr< GeomFactors > GeomFactorsSharedPtr
Pointer to a GeomFactors object.
Definition: GeomFactors.h:60
@ eRegular
Geometry is straight-sided with constant geometric factors.
Definition: SpatialDomains.hpp:52
void Vsqrt(int n, const T *x, const int incx, T *y, const int incy)
sqrt y = sqrt(x)
Definition: Vmath.hpp:340
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.hpp:366
void Sdiv(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha/x.
Definition: Vmath.hpp:154
void Reverse(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.hpp:844
References ASSERTL0, Nektar::StdRegions::eBackwards, Nektar::SpatialDomains::eDeformed, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::SpatialDomains::eMovingRegular, Nektar::SpatialDomains::eRegular, Vmath::Fill(), Nektar::StdRegions::StdExpansion::GetCoordim(), GetEdgeInterpVals(), Nektar::LocalRegions::Expansion::GetGeom(), Nektar::LibUtilities::BasisKey::GetNumPoints(), Nektar::LibUtilities::PointsKey::GetNumPoints(), Nektar::StdRegions::StdExpansion::GetNumPoints(), Nektar::LibUtilities::BasisKey::GetPointsKey(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::StdRegions::StdExpansion::GetTraceBasisKey(), Nektar::LibUtilities::Interp1D(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_elmtBndNormDirElmtLen, Nektar::LocalRegions::Expansion::m_traceNormals, Vmath::Reverse(), Vmath::Sdiv(), Vmath::Smul(), tinysimd::sqrt(), Vmath::Vcopy(), Vmath::Vmul(), Vmath::Vsqrt(), Vmath::Vvtvp(), and Vmath::Zero().
◆ v_CreateStdMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1374 of file QuadExp.cpp.
1376{
1377 LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
1378 LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
1379 StdRegions::StdQuadExpSharedPtr tmp =
1380 MemoryManager<StdQuadExp>::AllocateSharedPtr(bkey0, bkey1);
1381 return tmp->GetStdMatrix(mkey);
1382}
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:166
std::shared_ptr< StdQuadExp > StdQuadExpSharedPtr
Definition: StdQuadExp.h:221
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_DropLocMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1389 of file QuadExp.cpp.
References m_matrixManager.
◆ v_DropLocStaticCondMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1399 of file QuadExp.cpp.
References m_staticCondMatrixManager.
◆ v_ExtractDataToCoeffs()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1259 of file QuadExp.cpp.
1263{
1264 int data_order0 = nummodes[mode_offset];
1266
1267 int data_order1 = nummodes[mode_offset + 1];
1269 int fillorder1 = min(order1, data_order1);
1270
1271 // Check if same basis
1274 {
1275 // Construct a quad with the appropriate basis type at our
1276 // quadrature points, and one more to do a forwards
1277 // transform. We can then copy the output to coeffs.
1278 StdRegions::StdQuadExp tmpQuad(
1279 LibUtilities::BasisKey(fromType[0], data_order0,
1280 m_base[0]->GetPointsKey()),
1281 LibUtilities::BasisKey(fromType[1], data_order1,
1282 m_base[1]->GetPointsKey()));
1283 StdRegions::StdQuadExp tmpQuad2(m_base[0]->GetBasisKey(),
1284 m_base[1]->GetBasisKey());
1285
1286 Array<OneD, const NekDouble> tmpData(tmpQuad.GetNcoeffs(), data);
1287 Array<OneD, NekDouble> tmpBwd(tmpQuad2.GetTotPoints());
1288 Array<OneD, NekDouble> tmpOut(tmpQuad2.GetNcoeffs());
1289
1290 tmpQuad.BwdTrans(tmpData, tmpBwd);
1291 tmpQuad2.FwdTrans(tmpBwd, tmpOut);
1292 Vmath::Vcopy(tmpOut.size(), &tmpOut[0], 1, coeffs, 1);
1293
1294 return;
1295 }
1296
1298 {
1300 {
1301 int i;
1302 int cnt = 0;
1303 int cnt1 = 0;
1304
1306 "Extraction routine not set up for this basis");
1307
1309 for (i = 0; i < fillorder0; ++i)
1310 {
1311 Vmath::Vcopy(fillorder1, data + cnt, 1, coeffs + cnt1, 1);
1312 cnt += data_order1;
1313 cnt1 += order1;
1314 }
1315 }
1316 break;
1318 {
1319 LibUtilities::PointsKey p0(nummodes[0],
1321 LibUtilities::PointsKey p1(nummodes[1],
1323 LibUtilities::PointsKey t0(m_base[0]->GetNumModes(),
1325 LibUtilities::PointsKey t1(m_base[1]->GetNumModes(),
1327 LibUtilities::Interp2D(p0, p1, data, t0, t1, coeffs);
1328 }
1329 break;
1331 {
1332 // Assume that input is also Gll_Lagrange but no way to check;
1333 LibUtilities::PointsKey p0(nummodes[0],
1335 LibUtilities::PointsKey p1(nummodes[1],
1337 LibUtilities::PointsKey t0(m_base[0]->GetNumModes(),
1339 LibUtilities::PointsKey t1(m_base[1]->GetNumModes(),
1341 LibUtilities::Interp2D(p0, p1, data, t0, t1, coeffs);
1342 }
1343 break;
1344 default:
1346 }
1347}
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
Definition: StdExpansion.h:156
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:101
References ASSERTL0, ASSERTL1, Nektar::StdRegions::StdExpansion::BwdTrans(), Nektar::LibUtilities::eGauss_Lagrange, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::LibUtilities::eGaussLobattoLegendre, Nektar::LibUtilities::eGLL_Lagrange, Nektar::LibUtilities::eModified_A, Nektar::StdRegions::StdExpansion::FwdTrans(), Nektar::StdRegions::StdExpansion::GetBasisType(), Nektar::StdRegions::StdExpansion::GetNcoeffs(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::LibUtilities::Interp2D(), Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::m_ncoeffs, Vmath::Vcopy(), and Vmath::Zero().
◆ v_FwdTrans()
|
overrideprotectedvirtual |
Transform a given function from physical quadrature space to coefficient space.
- See also
- StdExpansion::FwdTrans
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 227 of file QuadExp.cpp.
229{
231 {
233 }
234 else
235 {
236 IProductWRTBase(inarray, outarray);
237
238 // get Mass matrix inverse
241
242 // copy inarray in case inarray == outarray
245
246 out = (*matsys) * in;
247 }
248}
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:537
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().
◆ v_FwdTransBndConstrained()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 250 of file QuadExp.cpp.
253{
255 {
257 }
258 else
259 {
260 int i, j;
263
264 fill(outarray.data(), outarray.data() + m_ncoeffs, 0.0);
265
266 if (nmodes[0] == 1 && nmodes[1] == 1)
267 {
268 outarray[0] = inarray[0];
269 return;
270 }
271
272 Array<OneD, NekDouble> physEdge[4];
273 Array<OneD, NekDouble> coeffEdge[4];
274 StdRegions::Orientation orient[4];
275 for (i = 0; i < 4; i++)
276 {
277 physEdge[i] = Array<OneD, NekDouble>(npoints[i % 2]);
278 coeffEdge[i] = Array<OneD, NekDouble>(nmodes[i % 2]);
279 orient[i] = GetTraceOrient(i);
280 }
281
282 for (i = 0; i < npoints[0]; i++)
283 {
284 physEdge[0][i] = inarray[i];
285 physEdge[2][i] = inarray[npoints[0] * (npoints[1] - 1) + i];
286 }
287
288 for (i = 0; i < npoints[1]; i++)
289 {
290 physEdge[1][i] = inarray[npoints[0] - 1 + i * npoints[0]];
291 physEdge[3][i] = inarray[i * npoints[0]];
292 }
293
294 for (i = 0; i < 4; i++)
295 {
297 {
298 reverse((physEdge[i]).data(),
299 (physEdge[i]).data() + npoints[i % 2]);
300 }
301 }
302
303 SegExpSharedPtr segexp[4];
304 for (i = 0; i < 4; i++)
305 {
306 segexp[i] = MemoryManager<LocalRegions::SegExp>::AllocateSharedPtr(
308 }
309
310 Array<OneD, unsigned int> mapArray;
311 Array<OneD, int> signArray;
313
314 for (i = 0; i < 4; i++)
315 {
316 segexp[i % 2]->FwdTransBndConstrained(physEdge[i], coeffEdge[i]);
317
318 GetTraceToElementMap(i, mapArray, signArray, orient[i]);
319 for (j = 0; j < nmodes[i % 2]; j++)
320 {
322 outarray[mapArray[j]] = sign * coeffEdge[i][j];
323 }
324 }
325
328
329 if (nInteriorDofs > 0)
330 {
331 Array<OneD, NekDouble> tmp0(m_ncoeffs);
332 Array<OneD, NekDouble> tmp1(m_ncoeffs);
333
334 StdRegions::StdMatrixKey stdmasskey(StdRegions::eMass,
336 MassMatrixOp(outarray, tmp0, stdmasskey);
337 IProductWRTBase(inarray, tmp1);
338
340
341 // get Mass matrix inverse (only of interior DOF)
342 // use block (1,1) of the static condensed system
343 // note: this block alreay contains the inverse matrix
345 DNekScalMatSharedPtr matsys =
346 (m_staticCondMatrixManager[masskey])->GetBlock(1, 1);
347
348 Array<OneD, NekDouble> rhs(nInteriorDofs);
349 Array<OneD, NekDouble> result(nInteriorDofs);
350
351 GetInteriorMap(mapArray);
352
353 for (i = 0; i < nInteriorDofs; i++)
354 {
355 rhs[i] = tmp1[mapArray[i]];
356 }
357
359 &((matsys->GetOwnedMatrix())->GetPtr())[0],
360 nInteriorDofs, rhs.data(), 1, 0.0, result.data(), 1);
361
362 for (i = 0; i < nInteriorDofs; i++)
363 {
364 outarray[mapArray[i]] = result[i];
365 }
366 }
367 }
368}
SpatialDomains::Geometry2DSharedPtr GetGeom2D() const
Definition: Expansion2D.h:164
StdRegions::Orientation GetTraceOrient(int trace)
Definition: Expansion.h:168
int NumBndryCoeffs(void) const
Definition: StdExpansion.h:325
void MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:761
void GetTraceToElementMap(const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards, int P=-1, int Q=-1)
Definition: StdExpansion.h:693
void GetInteriorMap(Array< OneD, unsigned int > &outarray)
Definition: StdExpansion.h:683
static void Dgemv(const char &trans, const int &m, const int &n, const double &alpha, const double *a, const int &lda, const double *x, const int &incx, const double &beta, double *y, const int &incy)
BLAS level 2: Matrix vector multiply y = alpha A x plus beta y where A[m x n].
Definition: Blas.hpp:211
std::shared_ptr< SegExp > SegExpSharedPtr
Definition: SegExp.h:253
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.hpp:220
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::StdRegions::StdExpansion::DetShapeType(), Blas::Dgemv(), Nektar::StdRegions::eBackwards, Nektar::StdRegions::eMass, Nektar::LocalRegions::Expansion2D::GetGeom2D(), Nektar::StdRegions::StdExpansion::GetInteriorMap(), Nektar::LocalRegions::Expansion::GetTraceOrient(), Nektar::StdRegions::StdExpansion::GetTraceToElementMap(), 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().
◆ v_GenMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1354 of file QuadExp.cpp.
1355{
1356 DNekMatSharedPtr returnval;
1357 switch (mkey.GetMatrixType())
1358 {
1366 returnval = Expansion2D::v_GenMatrix(mkey);
1367 break;
1368 default:
1369 returnval = StdQuadExp::v_GenMatrix(mkey);
1370 }
1371 return returnval;
1372}
DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey) override
Definition: Expansion2D.cpp:1256
@ eInvLaplacianWithUnityMean
Definition: StdRegions.hpp:96
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().
◆ v_GetCoord()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 547 of file QuadExp.cpp.
549{
550 int i;
551
552 ASSERTL1(Lcoords[0] >= -1.0 && Lcoords[1] <= 1.0 && Lcoords[1] >= -1.0 &&
553 Lcoords[1] <= 1.0,
554 "Local coordinates are not in region [-1,1]");
555
556 m_geom->FillGeom();
558 {
559 coords[i] = m_geom->GetCoord(i, Lcoords);
560 }
561}
References ASSERTL1, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_GetCoords()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 540 of file QuadExp.cpp.
543{
544 Expansion::v_GetCoords(coords_0, coords_1, coords_2);
545}
void v_GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3) override
Definition: Expansion.cpp:534
References Nektar::LocalRegions::Expansion::v_GetCoords().
◆ v_GetLinStdExp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 529 of file QuadExp.cpp.
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::StdRegions::StdExpansion::GetBasisType(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_GetLocMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1384 of file QuadExp.cpp.
References m_matrixManager.
◆ v_GetLocStaticCondMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1394 of file QuadExp.cpp.
References m_staticCondMatrixManager.
◆ v_GetStdExp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 523 of file QuadExp.cpp.
524{
527}
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_GetTraceOrient()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1349 of file QuadExp.cpp.
References Nektar::LocalRegions::Expansion::m_geom.
◆ v_GetTracePhysMap()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 724 of file QuadExp.cpp.
725{
728
729 // Get points in Cartesian orientation
730 switch (edge)
731 {
732 case 0:
733 outarray = Array<OneD, int>(nquad0);
734 for (int i = 0; i < nquad0; ++i)
735 {
736 outarray[i] = i;
737 }
738 break;
739 case 1:
740 outarray = Array<OneD, int>(nquad1);
741 for (int i = 0; i < nquad1; ++i)
742 {
743 outarray[i] = (nquad0 - 1) + i * nquad0;
744 }
745 break;
746 case 2:
747 outarray = Array<OneD, int>(nquad0);
748 for (int i = 0; i < nquad0; ++i)
749 {
750 outarray[i] = i + nquad0 * (nquad1 - 1);
751 }
752 break;
753 case 3:
754 outarray = Array<OneD, int>(nquad1);
755 for (int i = 0; i < nquad1; ++i)
756 {
757 outarray[i] = i * nquad0;
758 }
759 break;
760 default:
762 break;
763 }
764}
References ASSERTL0, and Nektar::StdRegions::StdExpansion::m_base.
◆ v_GetTracePhysVals()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 598 of file QuadExp.cpp.
602{
605
606 // Implementation for all the basis except Gauss points
609 {
610 switch (edge)
611 {
612 case 0:
613 Vmath::Vcopy(nquad0, &(inarray[0]), 1, &(outarray[0]), 1);
614 break;
615 case 1:
616 Vmath::Vcopy(nquad1, &(inarray[0]) + (nquad0 - 1), nquad0,
617 &(outarray[0]), 1);
618 break;
619 case 2:
620 Vmath::Vcopy(nquad0, &(inarray[0]) + nquad0 * (nquad1 - 1), 1,
621 &(outarray[0]), 1);
622 break;
623 case 3:
624 Vmath::Vcopy(nquad1, &(inarray[0]), nquad0, &(outarray[0]), 1);
625 break;
626 default:
628 break;
629 }
630 }
631 else
632 {
633 QuadExp::GetEdgeInterpVals(edge, inarray, outarray);
634 }
635
636 // Interpolate if required
638 EdgeExp->GetBasis(0)->GetPointsKey())
639 {
640 Array<OneD, NekDouble> outtmp(max(nquad0, nquad1));
641
642 outtmp = outarray;
643
645 EdgeExp->GetBasis(0)->GetPointsKey(), outarray);
646 }
647
649 {
650 orient = GetTraceOrient(edge);
651 }
652 // Reverse data if necessary
654 {
655 Vmath::Reverse(EdgeExp->GetNumPoints(0), &outarray[0], 1, &outarray[0],
656 1);
657 }
658}
References ASSERTL0, Nektar::StdRegions::eBackwards, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::StdRegions::eNoOrientation, GetEdgeInterpVals(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::LocalRegions::Expansion::GetTraceOrient(), Nektar::LibUtilities::Interp1D(), Nektar::StdRegions::StdExpansion::m_base, Vmath::Reverse(), and Vmath::Vcopy().
◆ v_GetTraceQFactors()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 766 of file QuadExp.cpp.
768{
769 int i;
772
775 const Array<TwoD, const NekDouble> &df =
776 m_metricinfo->GetDerivFactors(ptsKeys);
777
778 Array<OneD, NekDouble> j(max(nquad0, nquad1), 0.0);
779 Array<OneD, NekDouble> g0(max(nquad0, nquad1), 0.0);
780 Array<OneD, NekDouble> g1(max(nquad0, nquad1), 0.0);
781 Array<OneD, NekDouble> g2(max(nquad0, nquad1), 0.0);
782 Array<OneD, NekDouble> g3(max(nquad0, nquad1), 0.0);
783
785 {
786 // Implementation for all the basis except Gauss points
789 {
790 switch (edge)
791 {
792 case 0:
793 Vmath::Vcopy(nquad0, &(df[1][0]), 1, &(g1[0]), 1);
794 Vmath::Vcopy(nquad0, &(df[3][0]), 1, &(g3[0]), 1);
795 Vmath::Vcopy(nquad0, &(jac[0]), 1, &(j[0]), 1);
796
797 for (i = 0; i < nquad0; ++i)
798 {
799 outarray[i] =
800 j[i] * sqrt(g1[i] * g1[i] + g3[i] * g3[i]);
801 }
802 break;
803 case 1:
804 Vmath::Vcopy(nquad1, &(df[0][0]) + (nquad0 - 1), nquad0,
805 &(g0[0]), 1);
806
807 Vmath::Vcopy(nquad1, &(df[2][0]) + (nquad0 - 1), nquad0,
808 &(g2[0]), 1);
809
810 Vmath::Vcopy(nquad1, &(jac[0]) + (nquad0 - 1), nquad0,
811 &(j[0]), 1);
812
813 for (i = 0; i < nquad1; ++i)
814 {
815 outarray[i] =
816 j[i] * sqrt(g0[i] * g0[i] + g2[i] * g2[i]);
817 }
818 break;
819 case 2:
820
821 Vmath::Vcopy(nquad0, &(df[1][0]) + (nquad0 * (nquad1 - 1)),
822 1, &(g1[0]), 1);
823
824 Vmath::Vcopy(nquad0, &(df[3][0]) + (nquad0 * (nquad1 - 1)),
825 1, &(g3[0]), 1);
826
827 Vmath::Vcopy(nquad0, &(jac[0]) + (nquad0 * (nquad1 - 1)), 1,
828 &(j[0]), 1);
829
830 for (i = 0; i < nquad0; ++i)
831 {
832 outarray[i] =
833 j[i] * sqrt(g1[i] * g1[i] + g3[i] * g3[i]);
834 }
835 break;
836 case 3:
837
838 Vmath::Vcopy(nquad1, &(df[0][0]), nquad0, &(g0[0]), 1);
839 Vmath::Vcopy(nquad1, &(df[2][0]), nquad0, &(g2[0]), 1);
840 Vmath::Vcopy(nquad1, &(jac[0]), nquad0, &(j[0]), 1);
841
842 for (i = 0; i < nquad1; ++i)
843 {
844 outarray[i] =
845 j[i] * sqrt(g0[i] * g0[i] + g2[i] * g2[i]);
846 }
847 break;
848 default:
850 break;
851 }
852 }
853 else
854 {
855 int nqtot = nquad0 * nquad1;
856 Array<OneD, NekDouble> tmp_gmat0(nqtot, 0.0);
857 Array<OneD, NekDouble> tmp_gmat1(nqtot, 0.0);
858 Array<OneD, NekDouble> tmp_gmat2(nqtot, 0.0);
859 Array<OneD, NekDouble> tmp_gmat3(nqtot, 0.0);
860 Array<OneD, NekDouble> g0_edge(max(nquad0, nquad1), 0.0);
861 Array<OneD, NekDouble> g1_edge(max(nquad0, nquad1), 0.0);
862 Array<OneD, NekDouble> g2_edge(max(nquad0, nquad1), 0.0);
863 Array<OneD, NekDouble> g3_edge(max(nquad0, nquad1), 0.0);
864 Array<OneD, NekDouble> jac_edge(max(nquad0, nquad1), 0.0);
865
866 switch (edge)
867 {
868 case 0:
869 Vmath::Vmul(nqtot, &(df[1][0]), 1, &jac[0], 1,
870 &(tmp_gmat1[0]), 1);
871 Vmath::Vmul(nqtot, &(df[3][0]), 1, &jac[0], 1,
872 &(tmp_gmat3[0]), 1);
873 QuadExp::GetEdgeInterpVals(edge, tmp_gmat1, g1_edge);
874 QuadExp::GetEdgeInterpVals(edge, tmp_gmat3, g3_edge);
875
876 for (i = 0; i < nquad0; ++i)
877 {
878 outarray[i] = sqrt(g1_edge[i] * g1_edge[i] +
879 g3_edge[i] * g3_edge[i]);
880 }
881 break;
882
883 case 1:
884 Vmath::Vmul(nqtot, &(df[0][0]), 1, &jac[0], 1,
885 &(tmp_gmat0[0]), 1);
886 Vmath::Vmul(nqtot, &(df[2][0]), 1, &jac[0], 1,
887 &(tmp_gmat2[0]), 1);
888 QuadExp::GetEdgeInterpVals(edge, tmp_gmat0, g0_edge);
889 QuadExp::GetEdgeInterpVals(edge, tmp_gmat2, g2_edge);
890
891 for (i = 0; i < nquad1; ++i)
892 {
893 outarray[i] = sqrt(g0_edge[i] * g0_edge[i] +
894 g2_edge[i] * g2_edge[i]);
895 }
896
897 break;
898 case 2:
899
900 Vmath::Vmul(nqtot, &(df[1][0]), 1, &jac[0], 1,
901 &(tmp_gmat1[0]), 1);
902 Vmath::Vmul(nqtot, &(df[3][0]), 1, &jac[0], 1,
903 &(tmp_gmat3[0]), 1);
904 QuadExp::GetEdgeInterpVals(edge, tmp_gmat1, g1_edge);
905 QuadExp::GetEdgeInterpVals(edge, tmp_gmat3, g3_edge);
906
907 for (i = 0; i < nquad0; ++i)
908 {
909 outarray[i] = sqrt(g1_edge[i] * g1_edge[i] +
910 g3_edge[i] * g3_edge[i]);
911 }
912
913 Vmath::Reverse(nquad0, &outarray[0], 1, &outarray[0], 1);
914
915 break;
916 case 3:
917 Vmath::Vmul(nqtot, &(df[0][0]), 1, &jac[0], 1,
918 &(tmp_gmat0[0]), 1);
919 Vmath::Vmul(nqtot, &(df[2][0]), 1, &jac[0], 1,
920 &(tmp_gmat2[0]), 1);
921 QuadExp::GetEdgeInterpVals(edge, tmp_gmat0, g0_edge);
922 QuadExp::GetEdgeInterpVals(edge, tmp_gmat2, g2_edge);
923
924 for (i = 0; i < nquad1; ++i)
925 {
926 outarray[i] = sqrt(g0_edge[i] * g0_edge[i] +
927 g2_edge[i] * g2_edge[i]);
928 }
929
930 Vmath::Reverse(nquad1, &outarray[0], 1, &outarray[0], 1);
931
932 break;
933 default:
935 break;
936 }
937 }
938 }
939 else
940 {
941
942 switch (edge)
943 {
944 case 0:
945
946 for (i = 0; i < nquad0; ++i)
947 {
948 outarray[i] = jac[0] * sqrt(df[1][0] * df[1][0] +
949 df[3][0] * df[3][0]);
950 }
951 break;
952 case 1:
953 for (i = 0; i < nquad1; ++i)
954 {
955 outarray[i] = jac[0] * sqrt(df[0][0] * df[0][0] +
956 df[2][0] * df[2][0]);
957 }
958 break;
959 case 2:
960 for (i = 0; i < nquad0; ++i)
961 {
962 outarray[i] = jac[0] * sqrt(df[1][0] * df[1][0] +
963 df[3][0] * df[3][0]);
964 }
965 break;
966 case 3:
967 for (i = 0; i < nquad1; ++i)
968 {
969 outarray[i] = jac[0] * sqrt(df[0][0] * df[0][0] +
970 df[2][0] * df[2][0]);
971 }
972 break;
973 default:
975 break;
976 }
977 }
978}
References ASSERTL0, Nektar::SpatialDomains::eDeformed, Nektar::LibUtilities::eGaussGaussLegendre, GetEdgeInterpVals(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Reverse(), tinysimd::sqrt(), Vmath::Vcopy(), and Vmath::Vmul().
◆ v_HelmholtzMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1448 of file QuadExp.cpp.
1451{
1452 QuadExp::HelmholtzMatrixOp_MatFree(inarray, outarray, mkey);
1453}
void HelmholtzMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:1304
References Nektar::StdRegions::StdExpansion::HelmholtzMatrixOp_MatFree().
◆ v_Integral()
|
overrideprotectedvirtual |
Integrates the specified function over the domain.
- See also
- StdRegions::StdExpansion::Integral.
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 69 of file QuadExp.cpp.
70{
74 NekDouble ival;
75 Array<OneD, NekDouble> tmp(nquad0 * nquad1);
76
77 // multiply inarray with Jacobian
79 {
80 Vmath::Vmul(nquad0 * nquad1, jac, 1, inarray, 1, tmp, 1);
81 }
82 else
83 {
84 Vmath::Smul(nquad0 * nquad1, jac[0], inarray, 1, tmp, 1);
85 }
86
87 // call StdQuadExp version;
88 ival = StdQuadExp::v_Integral(tmp);
89 return ival;
90}
References Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Smul(), and Vmath::Vmul().
◆ v_IProductWRTBase()
|
overrideprotectedvirtual |
Calculate the inner product of inarray with respect to the basis B=base0*base1 and put into outarray.
\( \begin{array}{rcl} I_{pq} = (\phi_p \phi_q, u) & = & \sum_{i=0}^{nq_0} \sum_{j=0}^{nq_1} \phi_p(\xi_{0,i}) \phi_q(\xi_{1,j}) w^0_i w^1_j u(\xi_{0,i} \xi_{1,j}) \\ & = & \sum_{i=0}^{nq_0} \phi_p(\xi_{0,i}) \sum_{j=0}^{nq_1} \phi_q(\xi_{1,j}) \tilde{u}_{i,j} \end{array} \)
where
\( \tilde{u}_{i,j} = w^0_i w^1_j u(\xi_{0,i},\xi_{1,j}) \)
which can be implemented as
\( f_{qi} = \sum_{j=0}^{nq_1} \phi_q(\xi_{1,j}) \tilde{u}_{i,j} = {\bf B_1 U} \) \( I_{pq} = \sum_{i=0}^{nq_0} \phi_p(\xi_{0,i}) f_{qi} = {\bf B_0 F} \)
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 370 of file QuadExp.cpp.
372{
374 {
375 MultiplyByQuadratureMetric(inarray, outarray);
376 }
377 else
378 {
379 IProductWRTBase_SumFac(inarray, outarray);
380 }
381}
void IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true)
Definition: StdExpansion.h:1178
References Nektar::StdRegions::StdExpansion::IProductWRTBase_SumFac(), Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric().
◆ v_IProductWRTBase_SumFac()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 390 of file QuadExp.cpp.
393{
397
398 if (multiplybyweights)
399 {
400 Array<OneD, NekDouble> tmp(nquad0 * nquad1 + nquad1 * order0);
401 Array<OneD, NekDouble> wsp(tmp + nquad0 * nquad1);
402
403 MultiplyByQuadratureMetric(inarray, tmp);
404 StdQuadExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),
405 m_base[1]->GetBdata(), tmp,
406 outarray, wsp, true, true);
407 }
408 else
409 {
410 Array<OneD, NekDouble> wsp(nquad1 * order0);
411
412 StdQuadExp::IProductWRTBase_SumFacKernel(m_base[0]->GetBdata(),
413 m_base[1]->GetBdata(), inarray,
414 outarray, wsp, true, true);
415 }
416}
References Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric().
◆ v_IProductWRTDerivBase()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 383 of file QuadExp.cpp.
386{
387 IProductWRTDerivBase_SumFac(dir, inarray, outarray);
388}
void IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:1229
References Nektar::StdRegions::StdExpansion::IProductWRTDerivBase_SumFac().
◆ v_IProductWRTDerivBase_SumFac()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 418 of file QuadExp.cpp.
421{
424 "Invalid direction.");
425
428 int nqtot = nquad0 * nquad1;
430
431 Array<OneD, NekDouble> tmp1(2 * nqtot + m_ncoeffs + nmodes0 * nquad1);
432 Array<OneD, NekDouble> tmp2(tmp1 + nqtot);
433 Array<OneD, NekDouble> tmp3(tmp1 + 2 * nqtot);
434 Array<OneD, NekDouble> tmp4(tmp1 + 2 * nqtot + m_ncoeffs);
435
436 Array<OneD, Array<OneD, NekDouble>> tmp2D{2};
437 tmp2D[0] = tmp1;
438 tmp2D[1] = tmp2;
439
440 QuadExp::v_AlignVectorToCollapsedDir(dir, inarray, tmp2D);
441
442 MultiplyByQuadratureMetric(tmp1, tmp1);
443 MultiplyByQuadratureMetric(tmp2, tmp2);
444
446 tmp1, tmp3, tmp4, false, true);
448 tmp2, outarray, tmp4, true, false);
450}
void v_AlignVectorToCollapsedDir(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
Definition: QuadExp.cpp:452
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:192
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.hpp:180
References ASSERTL1, Nektar::StdRegions::StdExpansion2D::IProductWRTBase_SumFacKernel(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_geom, Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), v_AlignVectorToCollapsedDir(), and Vmath::Vadd().
◆ v_LaplacianMatrixOp() [1/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1411 of file QuadExp.cpp.
1414{
1415 QuadExp::LaplacianMatrixOp_MatFree(inarray, outarray, mkey);
1416}
void LaplacianMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:1261
References Nektar::StdRegions::StdExpansion::LaplacianMatrixOp_MatFree().
◆ v_LaplacianMatrixOp() [2/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1418 of file QuadExp.cpp.
1422{
1423 StdExpansion::LaplacianMatrixOp_MatFree(k1, k2, inarray, outarray, mkey);
1424}
◆ v_LaplacianMatrixOp_MatFree_Kernel()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1495 of file QuadExp.cpp.
1498{
1500 {
1501 ComputeLaplacianMetric();
1502 }
1503
1506 int nqtot = nquad0 * nquad1;
1509 int wspsize =
1510 max(max(max(nqtot, m_ncoeffs), nquad1 * nmodes0), nquad0 * nmodes1);
1511
1513
1518 const Array<OneD, const NekDouble> &metric00 =
1520 const Array<OneD, const NekDouble> &metric01 =
1522 const Array<OneD, const NekDouble> &metric11 =
1524
1525 // Allocate temporary storage
1526 Array<OneD, NekDouble> wsp0(wsp);
1527 Array<OneD, NekDouble> wsp1(wsp + wspsize);
1528 Array<OneD, NekDouble> wsp2(wsp + 2 * wspsize);
1529
1530 StdExpansion2D::PhysTensorDeriv(inarray, wsp1, wsp2);
1531
1532 // wsp0 = k = g0 * wsp1 + g1 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
1533 // wsp2 = l = g1 * wsp1 + g2 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
1534 // where g0, g1 and g2 are the metric terms set up in the GeomFactors class
1535 // especially for this purpose
1536 Vmath::Vvtvvtp(nqtot, &metric00[0], 1, &wsp1[0], 1, &metric01[0], 1,
1537 &wsp2[0], 1, &wsp0[0], 1);
1538 Vmath::Vvtvvtp(nqtot, &metric01[0], 1, &wsp1[0], 1, &metric11[0], 1,
1539 &wsp2[0], 1, &wsp2[0], 1);
1540
1541 // outarray = m = (D_xi1 * B)^T * k
1542 // wsp1 = n = (D_xi2 * B)^T * l
1544 true);
1546
1547 // outarray = outarray + wsp1
1548 // = L * u_hat
1550 1);
1551}
void ComputeLaplacianMetric()
Definition: Expansion.cpp:453
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.hpp:439
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().
◆ v_MassLevelCurvatureMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1441 of file QuadExp.cpp.
1444{
1445 StdExpansion::MassLevelCurvatureMatrixOp_MatFree(inarray, outarray, mkey);
1446}
◆ v_MassMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1404 of file QuadExp.cpp.
1407{
1408 StdExpansion::MassMatrixOp_MatFree(inarray, outarray, mkey);
1409}
◆ v_NormalTraceDerivFactors()
|
overrideprotectedvirtual |
: This method gets all of the factors which are required as part of the Gradient Jump Penalty (GJP) stabilisation and involves the product of the normal and geometric factors along the element trace.
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1622 of file QuadExp.cpp.
1626{
1629
1630 const Array<TwoD, const NekDouble> &df =
1632
1633 if (d0factors.size() != 4)
1634 {
1635 d0factors = Array<OneD, Array<OneD, NekDouble>>(4);
1636 d1factors = Array<OneD, Array<OneD, NekDouble>>(4);
1637 }
1638
1639 if (d0factors[0].size() != nquad0)
1640 {
1641 d0factors[0] = Array<OneD, NekDouble>(nquad0);
1642 d0factors[2] = Array<OneD, NekDouble>(nquad0);
1643 d1factors[0] = Array<OneD, NekDouble>(nquad0);
1644 d1factors[2] = Array<OneD, NekDouble>(nquad0);
1645 }
1646
1647 if (d0factors[1].size() != nquad1)
1648 {
1649 d0factors[1] = Array<OneD, NekDouble>(nquad1);
1650 d0factors[3] = Array<OneD, NekDouble>(nquad1);
1651 d1factors[1] = Array<OneD, NekDouble>(nquad1);
1652 d1factors[3] = Array<OneD, NekDouble>(nquad1);
1653 }
1654
1655 // Outwards normals
1656 const Array<OneD, const Array<OneD, NekDouble>> &normal_0 =
1657 GetTraceNormal(0);
1658 const Array<OneD, const Array<OneD, NekDouble>> &normal_1 =
1659 GetTraceNormal(1);
1660 const Array<OneD, const Array<OneD, NekDouble>> &normal_2 =
1661 GetTraceNormal(2);
1662 const Array<OneD, const Array<OneD, NekDouble>> &normal_3 =
1663 GetTraceNormal(3);
1664
1665 int ncoords = normal_0.size();
1666
1668 {
1669 // needs checking for 3D coords
1670
1671 // factors 0 and 2
1672 for (int i = 0; i < nquad0; ++i)
1673 {
1674 d0factors[0][i] = df[0][i] * normal_0[0][i];
1675 d0factors[2][i] = df[0][nquad0 * (nquad1 - 1) + i] * normal_2[0][i];
1676
1677 d1factors[0][i] = df[1][i] * normal_0[0][i];
1678 d1factors[2][i] = df[1][nquad0 * (nquad1 - 1) + i] * normal_2[0][i];
1679 }
1680
1681 for (int n = 1; n < ncoords; ++n)
1682 {
1683 // d xi_1/dy n_y
1684 // needs checking for 3D coords
1685 for (int i = 0; i < nquad0; ++i)
1686 {
1687 d0factors[0][i] += df[2 * n][i] * normal_0[n][i];
1688 d0factors[2][i] +=
1689 df[2 * n][nquad0 * (nquad1 - 1) + i] * normal_2[n][i];
1690
1691 d1factors[0][i] += df[2 * n + 1][i] * normal_0[n][i];
1692 d1factors[2][i] +=
1693 df[2 * n + 1][nquad0 * (nquad1 - 1) + i] * normal_2[n][i];
1694 }
1695 }
1696
1697 // faces 1 and 3
1698 for (int i = 0; i < nquad1; ++i)
1699 {
1700 d0factors[1][i] = df[0][(i + 1) * nquad0 - 1] * normal_1[0][i];
1701 d0factors[3][i] = df[0][i * nquad0] * normal_3[0][i];
1702
1703 d1factors[1][i] = df[1][(i + 1) * nquad0 - 1] * normal_1[0][i];
1704 d1factors[3][i] = df[1][i * nquad0] * normal_3[0][i];
1705 }
1706
1707 for (int n = 1; n < ncoords; ++n)
1708 {
1709 for (int i = 0; i < nquad1; ++i)
1710 {
1711 d0factors[1][i] +=
1712 df[2 * n][(i + 1) * nquad0 - 1] * normal_1[n][i];
1713 d0factors[3][i] += df[2 * n][i * nquad0] * normal_3[n][i];
1714
1715 d1factors[1][i] +=
1716 df[2 * n + 1][(i + 1) * nquad0 - 1] * normal_1[n][i];
1717 d1factors[3][i] += df[2 * n + 1][i * nquad0] * normal_3[n][i];
1718 }
1719 }
1720 }
1721 else
1722 {
1723 // d xi_2/dx n_x
1724 for (int i = 0; i < nquad0; ++i)
1725 {
1726 d1factors[0][i] = df[1][0] * normal_0[0][i];
1727 d1factors[2][i] = df[1][0] * normal_2[0][i];
1728 }
1729
1730 // d xi_1/dx n_x
1731 for (int i = 0; i < nquad1; ++i)
1732 {
1733 d0factors[1][i] = df[0][0] * normal_1[0][i];
1734 d0factors[3][i] = df[0][0] * normal_3[0][i];
1735 }
1736
1737 for (int n = 1; n < ncoords; ++n)
1738 {
1739 // d xi_2/dy n_y
1740 // needs checking for 3D coords
1741 for (int i = 0; i < nquad0; ++i)
1742 {
1743 d1factors[0][i] += df[2 * n + 1][0] * normal_0[n][i];
1744 d1factors[2][i] += df[2 * n + 1][0] * normal_2[n][i];
1745 }
1746
1747 // d xi_1/dy n_y
1748 // needs checking for 3D coords
1749 for (int i = 0; i < nquad1; ++i)
1750 {
1751 d0factors[1][i] += df[2 * n][0] * normal_1[n][i];
1752 d0factors[3][i] += df[2 * n][0] * normal_3[n][i];
1753 }
1754 }
1755
1756 // d1factors
1757 // d xi_1/dx n_x
1758 for (int i = 0; i < nquad0; ++i)
1759 {
1760 d0factors[0][i] = df[0][0] * normal_0[0][i];
1761 d0factors[2][i] = df[0][0] * normal_2[0][i];
1762 }
1763
1764 // d xi_2/dx n_x
1765 for (int i = 0; i < nquad1; ++i)
1766 {
1767 d1factors[1][i] = df[1][0] * normal_1[0][i];
1768 d1factors[3][i] = df[1][0] * normal_3[0][i];
1769 }
1770
1771 for (int n = 1; n < ncoords; ++n)
1772 {
1773 // d xi_1/dy n_y
1774 // needs checking for 3D coords
1775 for (int i = 0; i < nquad0; ++i)
1776 {
1777 d0factors[0][i] += df[2 * n][0] * normal_0[n][i];
1778 d0factors[2][i] += df[2 * n][0] * normal_2[n][i];
1779 }
1780
1781 // d xi_2/dy n_y
1782 // needs checking for 3D coords
1783 for (int i = 0; i < nquad1; ++i)
1784 {
1785 d1factors[1][i] += df[2 * n + 1][0] * normal_1[n][i];
1786 d1factors[3][i] += df[2 * n + 1][0] * normal_3[n][i];
1787 }
1788 }
1789 }
1790}
const NormalVector & GetTraceNormal(const int id)
Definition: Expansion.cpp:251
References Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetNumPoints(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::LocalRegions::Expansion::GetTraceNormal(), and Nektar::LocalRegions::Expansion::m_metricinfo.
◆ v_NormVectorIProductWRTBase() [1/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 516 of file QuadExp.cpp.
519{
520 NormVectorIProductWRTBase(Fvec[0], Fvec[1], Fvec[2], outarray);
521}
void NormVectorIProductWRTBase(const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:622
References Nektar::StdRegions::StdExpansion::NormVectorIProductWRTBase().
◆ v_NormVectorIProductWRTBase() [2/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 488 of file QuadExp.cpp.
492{
494 Array<OneD, NekDouble> Fn(nq);
495
496 const Array<OneD, const Array<OneD, NekDouble>> &normals =
497 GetLeftAdjacentElementExp()->GetTraceNormal(
498 GetLeftAdjacentElementTrace());
499
501 {
502 Vmath::Vvtvvtp(nq, &normals[0][0], 1, &Fx[0], 1, &normals[1][0], 1,
503 &Fy[0], 1, &Fn[0], 1);
504 Vmath::Vvtvp(nq, &normals[2][0], 1, &Fz[0], 1, &Fn[0], 1, &Fn[0], 1);
505 }
506 else
507 {
508 Vmath::Svtsvtp(nq, normals[0][0], &Fx[0], 1, normals[1][0], &Fy[0], 1,
509 &Fn[0], 1);
510 Vmath::Svtvp(nq, normals[2][0], &Fz[0], 1, &Fn[0], 1, &Fn[0], 1);
511 }
512
513 IProductWRTBase(Fn, outarray);
514}
ExpansionSharedPtr GetLeftAdjacentElementExp() const
Definition: Expansion.h:441
int GetLeftAdjacentElementTrace() const
Definition: Expansion.h:454
void Svtsvtp(int n, const T alpha, const T *x, int incx, const T beta, const T *y, int incy, T *z, int incz)
Svtsvtp (scalar times vector plus scalar times vector):
Definition: Vmath.hpp:473
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.hpp:396
References Nektar::SpatialDomains::eDeformed, Nektar::LocalRegions::Expansion::GetLeftAdjacentElementExp(), Nektar::LocalRegions::Expansion::GetLeftAdjacentElementTrace(), Nektar::StdRegions::StdExpansion::IProductWRTBase(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Svtsvtp(), Vmath::Svtvp(), Vmath::Vvtvp(), and Vmath::Vvtvvtp().
◆ v_PhysDeriv() [1/2]
|
overrideprotectedvirtual |
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 92 of file QuadExp.cpp.
96{
99 int nqtot = nquad0 * nquad1;
100 const Array<TwoD, const NekDouble> &df =
102 Array<OneD, NekDouble> diff0(2 * nqtot);
103 Array<OneD, NekDouble> diff1(diff0 + nqtot);
104
105 StdQuadExp::v_PhysDeriv(inarray, diff0, diff1);
106
108 {
109 if (out_d0.size())
110 {
111 Vmath::Vmul(nqtot, df[0], 1, diff0, 1, out_d0, 1);
112 Vmath::Vvtvp(nqtot, df[1], 1, diff1, 1, out_d0, 1, out_d0, 1);
113 }
114
115 if (out_d1.size())
116 {
117 Vmath::Vmul(nqtot, df[2], 1, diff0, 1, out_d1, 1);
118 Vmath::Vvtvp(nqtot, df[3], 1, diff1, 1, out_d1, 1, out_d1, 1);
119 }
120
121 if (out_d2.size())
122 {
123 Vmath::Vmul(nqtot, df[4], 1, diff0, 1, out_d2, 1);
124 Vmath::Vvtvp(nqtot, df[5], 1, diff1, 1, out_d2, 1, out_d2, 1);
125 }
126 }
127 else // regular geometry
128 {
129 if (out_d0.size())
130 {
131 Vmath::Smul(nqtot, df[0][0], diff0, 1, out_d0, 1);
132 Blas::Daxpy(nqtot, df[1][0], diff1, 1, out_d0, 1);
133 }
134
135 if (out_d1.size())
136 {
137 Vmath::Smul(nqtot, df[2][0], diff0, 1, out_d1, 1);
138 Blas::Daxpy(nqtot, df[3][0], diff1, 1, out_d1, 1);
139 }
140
141 if (out_d2.size())
142 {
143 Vmath::Smul(nqtot, df[4][0], diff0, 1, out_d2, 1);
144 Blas::Daxpy(nqtot, df[5][0], diff1, 1, out_d2, 1);
145 }
146 }
147}
static void Daxpy(const int &n, const double &alpha, const double *x, const int &incx, const double *y, const int &incy)
BLAS level 1: y = alpha x plus y.
Definition: Blas.hpp:135
References Blas::Daxpy(), Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Smul(), Vmath::Vmul(), and Vmath::Vvtvp().
Referenced by v_PhysDeriv().
◆ v_PhysDeriv() [2/2]
|
overrideprotectedvirtual |
Calculate the derivative of the physical points in a given direction.
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 149 of file QuadExp.cpp.
152{
153 switch (dir)
154 {
155 case 0:
156 {
158 NullNekDouble1DArray);
159 }
160 break;
161 case 1:
162 {
164 NullNekDouble1DArray);
165 }
166 break;
167 case 2:
168 {
170 outarray);
171 }
172 break;
173 default:
174 {
176 }
177 break;
178 }
179}
void v_PhysDeriv(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2=NullNekDouble1DArray) override
Calculate the derivative of the physical points.
Definition: QuadExp.cpp:92
static Array< OneD, NekDouble > NullNekDouble1DArray
Definition: BasicUtils/SharedArray.hpp:919
References ASSERTL1, Nektar::NullNekDouble1DArray, and v_PhysDeriv().
◆ v_PhysDirectionalDeriv()
|
overrideprotectedvirtual |
Physical derivative along a direction vector.
D_v = d/dx_v^s + d/dx_v^r
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 181 of file QuadExp.cpp.
184{
187 int nqtot = nquad0 * nquad1;
188
189 const Array<TwoD, const NekDouble> &df =
191
192 Array<OneD, NekDouble> diff0(2 * nqtot);
193 Array<OneD, NekDouble> diff1(diff0 + nqtot);
194
195 StdQuadExp::v_PhysDeriv(inarray, diff0, diff1);
196
198 {
199 Array<OneD, Array<OneD, NekDouble>> tangmat(2);
200
201 // d/dx_v^s = v_x*ds/dx + v_y*ds/dy + v_z*dx/dz
202 for (int i = 0; i < 2; ++i)
203 {
204 tangmat[i] = Array<OneD, NekDouble>(nqtot, 0.0);
206 {
207 Vmath::Vvtvp(nqtot, &df[2 * k + i][0], 1, &direction[k * nqtot],
208 1, &tangmat[i][0], 1, &tangmat[i][0], 1);
209 }
210 }
211
212 /// D_v = d/dx_v^s + d/dx_v^r
213 if (out.size())
214 {
215 Vmath::Vmul(nqtot, &tangmat[0][0], 1, &diff0[0], 1, &out[0], 1);
216 Vmath::Vvtvp(nqtot, &tangmat[1][0], 1, &diff1[0], 1, &out[0], 1,
217 &out[0], 1);
218 }
219 }
220 else
221 {
223 "Wrong route");
224 }
225}
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().
◆ v_PhysEvalFirstDeriv()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 587 of file QuadExp.cpp.
591{
592 Array<OneD, NekDouble> Lcoord(2);
594 m_geom->GetLocCoords(coord, Lcoord);
595 return StdQuadExp::v_PhysEvalFirstDeriv(Lcoord, inarray, firstOrderDerivs);
596}
References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_PhysEvaluate()
|
overrideprotectedvirtual |
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
-
coords the coordinates of the single point
- Returns
- returns the value of the expansion at the single point
Reimplemented from Nektar::StdRegions::StdExpansion2D.
Definition at line 576 of file QuadExp.cpp.
578{
579 Array<OneD, NekDouble> Lcoord = Array<OneD, NekDouble>(2);
580
582 m_geom->GetLocCoords(coord, Lcoord);
583
584 return StdExpansion2D::v_PhysEvaluate(Lcoord, physvals);
585}
References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_ReduceOrderCoeffs()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1455 of file QuadExp.cpp.
1458{
1459 int n_coeffs = inarray.size();
1460
1461 Array<OneD, NekDouble> coeff(n_coeffs);
1462 Array<OneD, NekDouble> coeff_tmp(n_coeffs, 0.0);
1463 Array<OneD, NekDouble> tmp, tmp2;
1464
1467 int numMax = nmodes0;
1468
1469 Vmath::Vcopy(n_coeffs, inarray, 1, coeff_tmp, 1);
1470
1471 const LibUtilities::PointsKey Pkey0(nmodes0,
1473 const LibUtilities::PointsKey Pkey1(nmodes1,
1477 LibUtilities::BasisKey bortho0(LibUtilities::eOrtho_A, nmodes0, Pkey0);
1478 LibUtilities::BasisKey bortho1(LibUtilities::eOrtho_A, nmodes1, Pkey1);
1479
1480 LibUtilities::InterpCoeff2D(b0, b1, coeff_tmp, bortho0, bortho1, coeff);
1481
1482 Vmath::Zero(n_coeffs, coeff_tmp, 1);
1483
1484 int cnt = 0;
1485 for (int i = 0; i < numMin + 1; ++i)
1486 {
1487 Vmath::Vcopy(numMin, tmp = coeff + cnt, 1, tmp2 = coeff_tmp + cnt, 1);
1488
1489 cnt = i * numMax;
1490 }
1491
1492 LibUtilities::InterpCoeff2D(bortho0, bortho1, coeff_tmp, b0, b1, outarray);
1493}
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:67
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().
◆ v_StdPhysEvaluate()
|
overrideprotectedvirtual |
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 568 of file QuadExp.cpp.
571{
572 // Evaluate point in local (eta) coordinates.
573 return StdExpansion2D::v_PhysEvaluate(Lcoord, physvals);
574}
◆ v_SVVLaplacianFilter()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1590 of file QuadExp.cpp.
1592{
1594
1595 // Calculate sqrt of the Jacobian
1597 Array<OneD, NekDouble> sqrt_jac(nq);
1599 {
1600 Vmath::Vsqrt(nq, jac, 1, sqrt_jac, 1);
1601 }
1602 else
1603 {
1605 }
1606
1607 // Multiply array by sqrt(Jac)
1608 Vmath::Vmul(nq, sqrt_jac, 1, array, 1, array, 1);
1609
1610 // Apply std region filter
1611 StdQuadExp::v_SVVLaplacianFilter(array, mkey);
1612
1613 // Divide by sqrt(Jac)
1614 Vmath::Vdiv(nq, array, 1, sqrt_jac, 1, array, 1);
1615}
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.hpp:126
References Nektar::SpatialDomains::eDeformed, Vmath::Fill(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::LocalRegions::Expansion::m_metricinfo, tinysimd::sqrt(), Vmath::Vdiv(), Vmath::Vmul(), and Vmath::Vsqrt().
◆ v_WeakDerivMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1426 of file QuadExp.cpp.
1430{
1431 StdExpansion::WeakDerivMatrixOp_MatFree(i, inarray, outarray, mkey);
1432}
◆ v_WeakDirectionalDerivMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1434 of file QuadExp.cpp.
1437{
1438 StdExpansion::WeakDirectionalDerivMatrixOp_MatFree(inarray, outarray, mkey);
1439}
Member Data Documentation
◆ m_matrixManager
|
private |
Definition at line 244 of file QuadExp.h.
Referenced by GetEdgeInterpVals(), v_DropLocMatrix(), v_FwdTrans(), and v_GetLocMatrix().
◆ m_staticCondMatrixManager
|
private |
Definition at line 246 of file QuadExp.h.
Referenced by v_DropLocStaticCondMatrix(), v_FwdTransBndConstrained(), and v_GetLocStaticCondMatrix().
