#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) 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 | 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... | |
| 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_PhysEvaluate (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) 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_PhysEvaluate (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_PhysEvaluate (const Array< OneD, DNekMatSharedPtr > &I, const Array< OneD, const NekDouble > &physvals) override |
| NekDouble | v_PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) 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 |
| 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 void | v_SetCoeffsToOrientation (Array< OneD, NekDouble > &coeffs, StdRegions::Orientation dir) |
| virtual NekDouble | v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals) |
| 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 657 of file QuadExp.cpp.
660{
661 int i;
664
667
670
672
673 switch (edge)
674 {
675 case 0:
676 {
677 for (i = 0; i < nq0; i++)
678 {
679 outarray[i] =
680 Blas::Ddot(nq1, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
681 1, &inarray[i], nq0);
682 }
683 break;
684 }
685 case 1:
686 {
687 for (i = 0; i < nq1; i++)
688 {
689 outarray[i] =
690 Blas::Ddot(nq0, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
691 1, &inarray[i * nq0], 1);
692 }
693 break;
694 }
695 case 2:
696 {
697 for (i = 0; i < nq0; i++)
698 {
699 outarray[i] =
700 Blas::Ddot(nq1, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
701 1, &inarray[i], nq0);
702 }
703 break;
704 }
705 case 3:
706 {
707 for (i = 0; i < nq1; i++)
708 {
709 outarray[i] =
710 Blas::Ddot(nq0, mat_gauss->GetOwnedMatrix()->GetPtr().get(),
711 1, &inarray[i * nq0], 1);
712 }
713 break;
714 }
715 default:
717 break;
718 }
719}
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:367
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1173
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:402
StdRegions::ConstFactorMap factors
Definition: TestVarcoeffHashing.cpp:51
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.get(), 1, tmp1.get(), 1);
476 Vmath::Vmul(nqtot, &df[2 * dir + 1][0], 1, inarray.get(), 1, tmp2.get(),
477 1);
478 }
479 else
480 {
481 Vmath::Smul(nqtot, df[2 * dir][0], inarray.get(), 1, tmp1.get(), 1);
482 Vmath::Smul(nqtot, df[2 * dir + 1][0], inarray.get(), 1, tmp2.get(), 1);
483 }
484}
#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:1095
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 1549 of file QuadExp.cpp.
1550{
1552 {
1553 ComputeQuadratureMetric();
1554 }
1555
1558 const unsigned int dim = 2;
1563
1564 const Array<TwoD, const NekDouble> gmat =
1566 for (unsigned int i = 0; i < dim; ++i)
1567 {
1568 for (unsigned int j = i; j < dim; ++j)
1569 {
1570 m_metrics[m[i][j]] = Array<OneD, NekDouble>(nqtot);
1572 {
1573 Vmath::Vcopy(nqtot, &gmat[i * dim + j][0], 1,
1574 &m_metrics[m[i][j]][0], 1);
1575 }
1576 else
1577 {
1579 1);
1580 }
1582 }
1583 }
1584}
void ComputeQuadratureMetric()
Definition: Expansion.cpp:454
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:723
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 977 of file QuadExp.cpp.
978{
979 int i;
981 GetGeom()->GetMetricInfo();
982 SpatialDomains::GeomType type = geomFactors->GetGtype();
983
985 for (i = 0; i < ptsKeys.size(); ++i)
986 {
987 // Need at least 2 points for computing normals
989 {
990 LibUtilities::PointsKey pKey(2, ptsKeys[i].GetPointsType());
991 ptsKeys[i] = pKey;
992 }
993 }
994
995 const Array<TwoD, const NekDouble> &df =
996 geomFactors->GetDerivFactors(ptsKeys);
997 const Array<OneD, const NekDouble> &jac = geomFactors->GetJac(ptsKeys);
998
999 // The points of normals should follow trace basis, not local basis.
1000 LibUtilities::BasisKey tobasis = GetTraceBasisKey(edge);
1001
1002 int nqe = tobasis.GetNumPoints();
1004
1005 m_traceNormals[edge] = Array<OneD, Array<OneD, NekDouble>>(vCoordDim);
1006 Array<OneD, Array<OneD, NekDouble>> &normal = m_traceNormals[edge];
1007 for (i = 0; i < vCoordDim; ++i)
1008 {
1009 normal[i] = Array<OneD, NekDouble>(nqe);
1010 }
1011
1012 size_t nqb = nqe;
1013 size_t nbnd = edge;
1014 m_elmtBndNormDirElmtLen[nbnd] = Array<OneD, NekDouble>{nqb, 0.0};
1015 Array<OneD, NekDouble> &length = m_elmtBndNormDirElmtLen[nbnd];
1016
1017 // Regular geometry case
1019 (type == SpatialDomains::eMovingRegular))
1020 {
1021 NekDouble fac;
1022 // Set up normals
1023 switch (edge)
1024 {
1025 case 0:
1026 for (i = 0; i < vCoordDim; ++i)
1027 {
1028 Vmath::Fill(nqe, -df[2 * i + 1][0], normal[i], 1);
1029 }
1030 break;
1031 case 1:
1032 for (i = 0; i < vCoordDim; ++i)
1033 {
1034 Vmath::Fill(nqe, df[2 * i][0], normal[i], 1);
1035 }
1036 break;
1037 case 2:
1038 for (i = 0; i < vCoordDim; ++i)
1039 {
1040 Vmath::Fill(nqe, df[2 * i + 1][0], normal[i], 1);
1041 }
1042 break;
1043 case 3:
1044 for (i = 0; i < vCoordDim; ++i)
1045 {
1046 Vmath::Fill(nqe, -df[2 * i][0], normal[i], 1);
1047 }
1048 break;
1049 default:
1051 }
1052
1053 // normalise
1054 fac = 0.0;
1055 for (i = 0; i < vCoordDim; ++i)
1056 {
1057 fac += normal[i][0] * normal[i][0];
1058 }
1059 fac = 1.0 / sqrt(fac);
1060
1061 Vmath::Fill(nqb, fac, length, 1);
1062
1063 for (i = 0; i < vCoordDim; ++i)
1064 {
1065 Vmath::Smul(nqe, fac, normal[i], 1, normal[i], 1);
1066 }
1067 }
1068 else // Set up deformed normals
1069 {
1070 int j;
1071
1072 int nquad0 = ptsKeys[0].GetNumPoints();
1073 int nquad1 = ptsKeys[1].GetNumPoints();
1074
1075 LibUtilities::PointsKey from_key;
1076
1077 Array<OneD, NekDouble> normals(vCoordDim * max(nquad0, nquad1), 0.0);
1078 Array<OneD, NekDouble> edgejac(vCoordDim * max(nquad0, nquad1), 0.0);
1079
1080 // Extract Jacobian along edges and recover local
1081 // derivates (dx/dr) for polynomial interpolation by
1082 // multiplying m_gmat by jacobian
1083
1084 // Implementation for all the basis except Gauss points
1087 {
1088 switch (edge)
1089 {
1090 case 0:
1091 for (j = 0; j < nquad0; ++j)
1092 {
1093 edgejac[j] = jac[j];
1094 for (i = 0; i < vCoordDim; ++i)
1095 {
1096 normals[i * nquad0 + j] =
1097 -df[2 * i + 1][j] * edgejac[j];
1098 }
1099 }
1100 from_key = ptsKeys[0];
1101 break;
1102 case 1:
1103 for (j = 0; j < nquad1; ++j)
1104 {
1105 edgejac[j] = jac[nquad0 * j + nquad0 - 1];
1106 for (i = 0; i < vCoordDim; ++i)
1107 {
1108 normals[i * nquad1 + j] =
1109 df[2 * i][nquad0 * j + nquad0 - 1] * edgejac[j];
1110 }
1111 }
1112 from_key = ptsKeys[1];
1113 break;
1114 case 2:
1115 for (j = 0; j < nquad0; ++j)
1116 {
1117 edgejac[j] = jac[nquad0 * (nquad1 - 1) + j];
1118 for (i = 0; i < vCoordDim; ++i)
1119 {
1120 normals[i * nquad0 + j] =
1121 (df[2 * i + 1][nquad0 * (nquad1 - 1) + j]) *
1122 edgejac[j];
1123 }
1124 }
1125 from_key = ptsKeys[0];
1126 break;
1127 case 3:
1128 for (j = 0; j < nquad1; ++j)
1129 {
1130 edgejac[j] = jac[nquad0 * j];
1131 for (i = 0; i < vCoordDim; ++i)
1132 {
1133 normals[i * nquad1 + j] =
1134 -df[2 * i][nquad0 * j] * edgejac[j];
1135 }
1136 }
1137 from_key = ptsKeys[1];
1138 break;
1139 default:
1141 }
1142 }
1143 else
1144 {
1145 int nqtot = nquad0 * nquad1;
1146 Array<OneD, NekDouble> tmp_gmat(nqtot, 0.0);
1147 Array<OneD, NekDouble> tmp_gmat_edge(nqe, 0.0);
1148
1149 switch (edge)
1150 {
1151 case 0:
1152 for (j = 0; j < nquad0; ++j)
1153 {
1154 for (i = 0; i < vCoordDim; ++i)
1155 {
1156 Vmath::Vmul(nqtot, &(df[2 * i + 1][0]), 1, &jac[0],
1157 1, &(tmp_gmat[0]), 1);
1158 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1159 tmp_gmat_edge);
1160 normals[i * nquad0 + j] = -tmp_gmat_edge[j];
1161 }
1162 }
1163 from_key = ptsKeys[0];
1164 break;
1165 case 1:
1166 for (j = 0; j < nquad1; ++j)
1167 {
1168 for (i = 0; i < vCoordDim; ++i)
1169 {
1170 Vmath::Vmul(nqtot, &(df[2 * i][0]), 1, &jac[0], 1,
1171 &(tmp_gmat[0]), 1);
1172 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1173 tmp_gmat_edge);
1174 normals[i * nquad1 + j] = tmp_gmat_edge[j];
1175 }
1176 }
1177 from_key = ptsKeys[1];
1178 break;
1179 case 2:
1180 for (j = 0; j < nquad0; ++j)
1181 {
1182 for (i = 0; i < vCoordDim; ++i)
1183 {
1184 Vmath::Vmul(nqtot, &(df[2 * i + 1][0]), 1, &jac[0],
1185 1, &(tmp_gmat[0]), 1);
1186 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1187 tmp_gmat_edge);
1188 normals[i * nquad0 + j] = tmp_gmat_edge[j];
1189 }
1190 }
1191 from_key = ptsKeys[0];
1192 break;
1193 case 3:
1194 for (j = 0; j < nquad1; ++j)
1195 {
1196 for (i = 0; i < vCoordDim; ++i)
1197 {
1198 Vmath::Vmul(nqtot, &(df[2 * i][0]), 1, &jac[0], 1,
1199 &(tmp_gmat[0]), 1);
1200 QuadExp::GetEdgeInterpVals(edge, tmp_gmat,
1201 tmp_gmat_edge);
1202 normals[i * nquad1 + j] = -tmp_gmat_edge[j];
1203 }
1204 }
1205 from_key = ptsKeys[1];
1206 break;
1207 default:
1209 }
1210 }
1211
1212 int nq = from_key.GetNumPoints();
1213 Array<OneD, NekDouble> work(nqe, 0.0);
1214
1215 // interpolate Jacobian and invert
1216 LibUtilities::Interp1D(from_key, jac, tobasis.GetPointsKey(), work);
1217 Vmath::Sdiv(nqe, 1.0, &work[0], 1, &work[0], 1);
1218
1219 // interpolate
1221 {
1222 LibUtilities::Interp1D(from_key, &normals[i * nq],
1223 tobasis.GetPointsKey(), &normal[i][0]);
1224 Vmath::Vmul(nqe, work, 1, normal[i], 1, normal[i], 1);
1225 }
1226
1227 // normalise normal vectors
1228 Vmath::Zero(nqe, work, 1);
1230 {
1231 Vmath::Vvtvp(nqe, normal[i], 1, normal[i], 1, work, 1, work, 1);
1232 }
1233
1234 Vmath::Vsqrt(nqe, work, 1, work, 1);
1235 Vmath::Sdiv(nqe, 1.0, work, 1, work, 1);
1236
1237 Vmath::Vcopy(nqb, work, 1, length, 1);
1238
1240 {
1241 Vmath::Vmul(nqe, normal[i], 1, work, 1, normal[i], 1);
1242 }
1243 }
1245 {
1246 for (i = 0; i < vCoordDim; ++i)
1247 {
1249 {
1250 Vmath::Reverse(nqe, normal[i], 1, normal[i], 1);
1251 }
1252 }
1253 }
1254}
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:657
const LibUtilities::BasisKey GetTraceBasisKey(const int i, int k=-1) const
This function returns the basis key belonging to the i-th trace.
Definition: StdExpansion.h:299
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
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 1371 of file QuadExp.cpp.
1373{
1374 LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
1375 LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
1376 StdRegions::StdQuadExpSharedPtr tmp =
1377 MemoryManager<StdQuadExp>::AllocateSharedPtr(bkey0, bkey1);
1378 return tmp->GetStdMatrix(mkey);
1379}
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:224
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_DropLocMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1386 of file QuadExp.cpp.
References m_matrixManager.
◆ v_DropLocStaticCondMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1396 of file QuadExp.cpp.
References m_staticCondMatrixManager.
◆ v_ExtractDataToCoeffs()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1256 of file QuadExp.cpp.
1260{
1261 int data_order0 = nummodes[mode_offset];
1263
1264 int data_order1 = nummodes[mode_offset + 1];
1266 int fillorder1 = min(order1, data_order1);
1267
1268 // Check if same basis
1271 {
1272 // Construct a quad with the appropriate basis type at our
1273 // quadrature points, and one more to do a forwards
1274 // transform. We can then copy the output to coeffs.
1275 StdRegions::StdQuadExp tmpQuad(
1276 LibUtilities::BasisKey(fromType[0], data_order0,
1277 m_base[0]->GetPointsKey()),
1278 LibUtilities::BasisKey(fromType[1], data_order1,
1279 m_base[1]->GetPointsKey()));
1280 StdRegions::StdQuadExp tmpQuad2(m_base[0]->GetBasisKey(),
1281 m_base[1]->GetBasisKey());
1282
1283 Array<OneD, const NekDouble> tmpData(tmpQuad.GetNcoeffs(), data);
1284 Array<OneD, NekDouble> tmpBwd(tmpQuad2.GetTotPoints());
1285 Array<OneD, NekDouble> tmpOut(tmpQuad2.GetNcoeffs());
1286
1287 tmpQuad.BwdTrans(tmpData, tmpBwd);
1288 tmpQuad2.FwdTrans(tmpBwd, tmpOut);
1289 Vmath::Vcopy(tmpOut.size(), &tmpOut[0], 1, coeffs, 1);
1290
1291 return;
1292 }
1293
1295 {
1297 {
1298 int i;
1299 int cnt = 0;
1300 int cnt1 = 0;
1301
1303 "Extraction routine not set up for this basis");
1304
1306 for (i = 0; i < fillorder0; ++i)
1307 {
1308 Vmath::Vcopy(fillorder1, data + cnt, 1, coeffs + cnt1, 1);
1309 cnt += data_order1;
1310 cnt1 += order1;
1311 }
1312 }
1313 break;
1315 {
1316 LibUtilities::PointsKey p0(nummodes[0],
1318 LibUtilities::PointsKey p1(nummodes[1],
1320 LibUtilities::PointsKey t0(m_base[0]->GetNumModes(),
1322 LibUtilities::PointsKey t1(m_base[1]->GetNumModes(),
1324 LibUtilities::Interp2D(p0, p1, data, t0, t1, coeffs);
1325 }
1326 break;
1328 {
1329 // Assume that input is also Gll_Lagrange but no way to check;
1330 LibUtilities::PointsKey p0(nummodes[0],
1332 LibUtilities::PointsKey p1(nummodes[1],
1334 LibUtilities::PointsKey t0(m_base[0]->GetNumModes(),
1336 LibUtilities::PointsKey t1(m_base[1]->GetNumModes(),
1338 LibUtilities::Interp2D(p0, p1, data, t0, t1, coeffs);
1339 }
1340 break;
1341 default:
1343 }
1344}
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:528
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.get(), outarray.get() + 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]).get(),
299 (physEdge[i]).get() + 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.get(), 1, 0.0, result.get(), 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:322
void MassMatrixOp(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:752
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:684
void GetInteriorMap(Array< OneD, unsigned int > &outarray)
Definition: StdExpansion.h:674
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:248
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 1351 of file QuadExp.cpp.
1352{
1353 DNekMatSharedPtr returnval;
1354 switch (mkey.GetMatrixType())
1355 {
1363 returnval = Expansion2D::v_GenMatrix(mkey);
1364 break;
1365 default:
1366 returnval = StdQuadExp::v_GenMatrix(mkey);
1367 }
1368 return returnval;
1369}
DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey) override
Definition: Expansion2D.cpp:1250
@ 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 545 of file QuadExp.cpp.
547{
548 int i;
549
550 ASSERTL1(Lcoords[0] >= -1.0 && Lcoords[1] <= 1.0 && Lcoords[1] >= -1.0 &&
551 Lcoords[1] <= 1.0,
552 "Local coordinates are not in region [-1,1]");
553
554 m_geom->FillGeom();
556 {
557 coords[i] = m_geom->GetCoord(i, Lcoords);
558 }
559}
References ASSERTL1, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_GetCoords()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 538 of file QuadExp.cpp.
541{
542 Expansion::v_GetCoords(coords_0, coords_1, coords_2);
543}
void v_GetCoords(Array< OneD, NekDouble > &coords_1, Array< OneD, NekDouble > &coords_2, Array< OneD, NekDouble > &coords_3) override
Definition: Expansion.cpp:530
References Nektar::LocalRegions::Expansion::v_GetCoords().
◆ v_GetLinStdExp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 527 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 1381 of file QuadExp.cpp.
References m_matrixManager.
◆ v_GetLocStaticCondMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1391 of file QuadExp.cpp.
References m_staticCondMatrixManager.
◆ v_GetStdExp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 521 of file QuadExp.cpp.
522{
525}
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_GetTraceOrient()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1346 of file QuadExp.cpp.
References Nektar::LocalRegions::Expansion::m_geom.
◆ v_GetTracePhysMap()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 721 of file QuadExp.cpp.
722{
725
726 // Get points in Cartesian orientation
727 switch (edge)
728 {
729 case 0:
730 outarray = Array<OneD, int>(nquad0);
731 for (int i = 0; i < nquad0; ++i)
732 {
733 outarray[i] = i;
734 }
735 break;
736 case 1:
737 outarray = Array<OneD, int>(nquad1);
738 for (int i = 0; i < nquad1; ++i)
739 {
740 outarray[i] = (nquad0 - 1) + i * nquad0;
741 }
742 break;
743 case 2:
744 outarray = Array<OneD, int>(nquad0);
745 for (int i = 0; i < nquad0; ++i)
746 {
747 outarray[i] = i + nquad0 * (nquad1 - 1);
748 }
749 break;
750 case 3:
751 outarray = Array<OneD, int>(nquad1);
752 for (int i = 0; i < nquad1; ++i)
753 {
754 outarray[i] = i * nquad0;
755 }
756 break;
757 default:
759 break;
760 }
761}
References ASSERTL0, and Nektar::StdRegions::StdExpansion::m_base.
◆ v_GetTracePhysVals()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 595 of file QuadExp.cpp.
599{
602
603 // Implementation for all the basis except Gauss points
606 {
607 switch (edge)
608 {
609 case 0:
610 Vmath::Vcopy(nquad0, &(inarray[0]), 1, &(outarray[0]), 1);
611 break;
612 case 1:
613 Vmath::Vcopy(nquad1, &(inarray[0]) + (nquad0 - 1), nquad0,
614 &(outarray[0]), 1);
615 break;
616 case 2:
617 Vmath::Vcopy(nquad0, &(inarray[0]) + nquad0 * (nquad1 - 1), 1,
618 &(outarray[0]), 1);
619 break;
620 case 3:
621 Vmath::Vcopy(nquad1, &(inarray[0]), nquad0, &(outarray[0]), 1);
622 break;
623 default:
625 break;
626 }
627 }
628 else
629 {
630 QuadExp::GetEdgeInterpVals(edge, inarray, outarray);
631 }
632
633 // Interpolate if required
635 EdgeExp->GetBasis(0)->GetPointsKey())
636 {
637 Array<OneD, NekDouble> outtmp(max(nquad0, nquad1));
638
639 outtmp = outarray;
640
642 EdgeExp->GetBasis(0)->GetPointsKey(), outarray);
643 }
644
646 {
647 orient = GetTraceOrient(edge);
648 }
649 // Reverse data if necessary
651 {
652 Vmath::Reverse(EdgeExp->GetNumPoints(0), &outarray[0], 1, &outarray[0],
653 1);
654 }
655}
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 763 of file QuadExp.cpp.
765{
766 int i;
769
772 const Array<TwoD, const NekDouble> &df =
773 m_metricinfo->GetDerivFactors(ptsKeys);
774
775 Array<OneD, NekDouble> j(max(nquad0, nquad1), 0.0);
776 Array<OneD, NekDouble> g0(max(nquad0, nquad1), 0.0);
777 Array<OneD, NekDouble> g1(max(nquad0, nquad1), 0.0);
778 Array<OneD, NekDouble> g2(max(nquad0, nquad1), 0.0);
779 Array<OneD, NekDouble> g3(max(nquad0, nquad1), 0.0);
780
782 {
783 // Implementation for all the basis except Gauss points
786 {
787 switch (edge)
788 {
789 case 0:
790 Vmath::Vcopy(nquad0, &(df[1][0]), 1, &(g1[0]), 1);
791 Vmath::Vcopy(nquad0, &(df[3][0]), 1, &(g3[0]), 1);
792 Vmath::Vcopy(nquad0, &(jac[0]), 1, &(j[0]), 1);
793
794 for (i = 0; i < nquad0; ++i)
795 {
796 outarray[i] =
797 j[i] * sqrt(g1[i] * g1[i] + g3[i] * g3[i]);
798 }
799 break;
800 case 1:
801 Vmath::Vcopy(nquad1, &(df[0][0]) + (nquad0 - 1), nquad0,
802 &(g0[0]), 1);
803
804 Vmath::Vcopy(nquad1, &(df[2][0]) + (nquad0 - 1), nquad0,
805 &(g2[0]), 1);
806
807 Vmath::Vcopy(nquad1, &(jac[0]) + (nquad0 - 1), nquad0,
808 &(j[0]), 1);
809
810 for (i = 0; i < nquad1; ++i)
811 {
812 outarray[i] =
813 j[i] * sqrt(g0[i] * g0[i] + g2[i] * g2[i]);
814 }
815 break;
816 case 2:
817
818 Vmath::Vcopy(nquad0, &(df[1][0]) + (nquad0 * (nquad1 - 1)),
819 1, &(g1[0]), 1);
820
821 Vmath::Vcopy(nquad0, &(df[3][0]) + (nquad0 * (nquad1 - 1)),
822 1, &(g3[0]), 1);
823
824 Vmath::Vcopy(nquad0, &(jac[0]) + (nquad0 * (nquad1 - 1)), 1,
825 &(j[0]), 1);
826
827 for (i = 0; i < nquad0; ++i)
828 {
829 outarray[i] =
830 j[i] * sqrt(g1[i] * g1[i] + g3[i] * g3[i]);
831 }
832 break;
833 case 3:
834
835 Vmath::Vcopy(nquad1, &(df[0][0]), nquad0, &(g0[0]), 1);
836 Vmath::Vcopy(nquad1, &(df[2][0]), nquad0, &(g2[0]), 1);
837 Vmath::Vcopy(nquad1, &(jac[0]), nquad0, &(j[0]), 1);
838
839 for (i = 0; i < nquad1; ++i)
840 {
841 outarray[i] =
842 j[i] * sqrt(g0[i] * g0[i] + g2[i] * g2[i]);
843 }
844 break;
845 default:
847 break;
848 }
849 }
850 else
851 {
852 int nqtot = nquad0 * nquad1;
853 Array<OneD, NekDouble> tmp_gmat0(nqtot, 0.0);
854 Array<OneD, NekDouble> tmp_gmat1(nqtot, 0.0);
855 Array<OneD, NekDouble> tmp_gmat2(nqtot, 0.0);
856 Array<OneD, NekDouble> tmp_gmat3(nqtot, 0.0);
857 Array<OneD, NekDouble> g0_edge(max(nquad0, nquad1), 0.0);
858 Array<OneD, NekDouble> g1_edge(max(nquad0, nquad1), 0.0);
859 Array<OneD, NekDouble> g2_edge(max(nquad0, nquad1), 0.0);
860 Array<OneD, NekDouble> g3_edge(max(nquad0, nquad1), 0.0);
861 Array<OneD, NekDouble> jac_edge(max(nquad0, nquad1), 0.0);
862
863 switch (edge)
864 {
865 case 0:
866 Vmath::Vmul(nqtot, &(df[1][0]), 1, &jac[0], 1,
867 &(tmp_gmat1[0]), 1);
868 Vmath::Vmul(nqtot, &(df[3][0]), 1, &jac[0], 1,
869 &(tmp_gmat3[0]), 1);
870 QuadExp::GetEdgeInterpVals(edge, tmp_gmat1, g1_edge);
871 QuadExp::GetEdgeInterpVals(edge, tmp_gmat3, g3_edge);
872
873 for (i = 0; i < nquad0; ++i)
874 {
875 outarray[i] = sqrt(g1_edge[i] * g1_edge[i] +
876 g3_edge[i] * g3_edge[i]);
877 }
878 break;
879
880 case 1:
881 Vmath::Vmul(nqtot, &(df[0][0]), 1, &jac[0], 1,
882 &(tmp_gmat0[0]), 1);
883 Vmath::Vmul(nqtot, &(df[2][0]), 1, &jac[0], 1,
884 &(tmp_gmat2[0]), 1);
885 QuadExp::GetEdgeInterpVals(edge, tmp_gmat0, g0_edge);
886 QuadExp::GetEdgeInterpVals(edge, tmp_gmat2, g2_edge);
887
888 for (i = 0; i < nquad1; ++i)
889 {
890 outarray[i] = sqrt(g0_edge[i] * g0_edge[i] +
891 g2_edge[i] * g2_edge[i]);
892 }
893
894 break;
895 case 2:
896
897 Vmath::Vmul(nqtot, &(df[1][0]), 1, &jac[0], 1,
898 &(tmp_gmat1[0]), 1);
899 Vmath::Vmul(nqtot, &(df[3][0]), 1, &jac[0], 1,
900 &(tmp_gmat3[0]), 1);
901 QuadExp::GetEdgeInterpVals(edge, tmp_gmat1, g1_edge);
902 QuadExp::GetEdgeInterpVals(edge, tmp_gmat3, g3_edge);
903
904 for (i = 0; i < nquad0; ++i)
905 {
906 outarray[i] = sqrt(g1_edge[i] * g1_edge[i] +
907 g3_edge[i] * g3_edge[i]);
908 }
909
910 Vmath::Reverse(nquad0, &outarray[0], 1, &outarray[0], 1);
911
912 break;
913 case 3:
914 Vmath::Vmul(nqtot, &(df[0][0]), 1, &jac[0], 1,
915 &(tmp_gmat0[0]), 1);
916 Vmath::Vmul(nqtot, &(df[2][0]), 1, &jac[0], 1,
917 &(tmp_gmat2[0]), 1);
918 QuadExp::GetEdgeInterpVals(edge, tmp_gmat0, g0_edge);
919 QuadExp::GetEdgeInterpVals(edge, tmp_gmat2, g2_edge);
920
921 for (i = 0; i < nquad1; ++i)
922 {
923 outarray[i] = sqrt(g0_edge[i] * g0_edge[i] +
924 g2_edge[i] * g2_edge[i]);
925 }
926
927 Vmath::Reverse(nquad1, &outarray[0], 1, &outarray[0], 1);
928
929 break;
930 default:
932 break;
933 }
934 }
935 }
936 else
937 {
938
939 switch (edge)
940 {
941 case 0:
942
943 for (i = 0; i < nquad0; ++i)
944 {
945 outarray[i] = jac[0] * sqrt(df[1][0] * df[1][0] +
946 df[3][0] * df[3][0]);
947 }
948 break;
949 case 1:
950 for (i = 0; i < nquad1; ++i)
951 {
952 outarray[i] = jac[0] * sqrt(df[0][0] * df[0][0] +
953 df[2][0] * df[2][0]);
954 }
955 break;
956 case 2:
957 for (i = 0; i < nquad0; ++i)
958 {
959 outarray[i] = jac[0] * sqrt(df[1][0] * df[1][0] +
960 df[3][0] * df[3][0]);
961 }
962 break;
963 case 3:
964 for (i = 0; i < nquad1; ++i)
965 {
966 outarray[i] = jac[0] * sqrt(df[0][0] * df[0][0] +
967 df[2][0] * df[2][0]);
968 }
969 break;
970 default:
972 break;
973 }
974 }
975}
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 1445 of file QuadExp.cpp.
1448{
1449 QuadExp::HelmholtzMatrixOp_MatFree(inarray, outarray, mkey);
1450}
void HelmholtzMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:1284
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:1158
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:1209
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:198
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 1408 of file QuadExp.cpp.
1411{
1412 QuadExp::LaplacianMatrixOp_MatFree(inarray, outarray, mkey);
1413}
void LaplacianMatrixOp_MatFree(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const StdMatrixKey &mkey)
Definition: StdExpansion.h:1241
References Nektar::StdRegions::StdExpansion::LaplacianMatrixOp_MatFree().
◆ v_LaplacianMatrixOp() [2/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1415 of file QuadExp.cpp.
1419{
1420 StdExpansion::LaplacianMatrixOp_MatFree(k1, k2, inarray, outarray, mkey);
1421}
◆ v_LaplacianMatrixOp_MatFree_Kernel()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1492 of file QuadExp.cpp.
1495{
1497 {
1498 ComputeLaplacianMetric();
1499 }
1500
1503 int nqtot = nquad0 * nquad1;
1506 int wspsize =
1507 max(max(max(nqtot, m_ncoeffs), nquad1 * nmodes0), nquad0 * nmodes1);
1508
1510
1515 const Array<OneD, const NekDouble> &metric00 =
1517 const Array<OneD, const NekDouble> &metric01 =
1519 const Array<OneD, const NekDouble> &metric11 =
1521
1522 // Allocate temporary storage
1523 Array<OneD, NekDouble> wsp0(wsp);
1524 Array<OneD, NekDouble> wsp1(wsp + wspsize);
1525 Array<OneD, NekDouble> wsp2(wsp + 2 * wspsize);
1526
1527 StdExpansion2D::PhysTensorDeriv(inarray, wsp1, wsp2);
1528
1529 // wsp0 = k = g0 * wsp1 + g1 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
1530 // wsp2 = l = g1 * wsp1 + g2 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
1531 // where g0, g1 and g2 are the metric terms set up in the GeomFactors class
1532 // especially for this purpose
1533 Vmath::Vvtvvtp(nqtot, &metric00[0], 1, &wsp1[0], 1, &metric01[0], 1,
1534 &wsp2[0], 1, &wsp0[0], 1);
1535 Vmath::Vvtvvtp(nqtot, &metric01[0], 1, &wsp1[0], 1, &metric11[0], 1,
1536 &wsp2[0], 1, &wsp2[0], 1);
1537
1538 // outarray = m = (D_xi1 * B)^T * k
1539 // wsp1 = n = (D_xi2 * B)^T * l
1541 true);
1543
1544 // outarray = outarray + wsp1
1545 // = L * u_hat
1547}
void ComputeLaplacianMetric()
Definition: Expansion.cpp:449
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 1438 of file QuadExp.cpp.
1441{
1442 StdExpansion::MassLevelCurvatureMatrixOp_MatFree(inarray, outarray, mkey);
1443}
◆ v_MassMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1401 of file QuadExp.cpp.
1404{
1405 StdExpansion::MassMatrixOp_MatFree(inarray, outarray, mkey);
1406}
◆ 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 1618 of file QuadExp.cpp.
1622{
1625
1626 const Array<TwoD, const NekDouble> &df =
1628
1629 if (d0factors.size() != 4)
1630 {
1631 d0factors = Array<OneD, Array<OneD, NekDouble>>(4);
1632 d1factors = Array<OneD, Array<OneD, NekDouble>>(4);
1633 }
1634
1635 if (d0factors[0].size() != nquad0)
1636 {
1637 d0factors[0] = Array<OneD, NekDouble>(nquad0);
1638 d0factors[2] = Array<OneD, NekDouble>(nquad0);
1639 d1factors[0] = Array<OneD, NekDouble>(nquad0);
1640 d1factors[2] = Array<OneD, NekDouble>(nquad0);
1641 }
1642
1643 if (d0factors[1].size() != nquad1)
1644 {
1645 d0factors[1] = Array<OneD, NekDouble>(nquad1);
1646 d0factors[3] = Array<OneD, NekDouble>(nquad1);
1647 d1factors[1] = Array<OneD, NekDouble>(nquad1);
1648 d1factors[3] = Array<OneD, NekDouble>(nquad1);
1649 }
1650
1651 // Outwards normals
1652 const Array<OneD, const Array<OneD, NekDouble>> &normal_0 =
1653 GetTraceNormal(0);
1654 const Array<OneD, const Array<OneD, NekDouble>> &normal_1 =
1655 GetTraceNormal(1);
1656 const Array<OneD, const Array<OneD, NekDouble>> &normal_2 =
1657 GetTraceNormal(2);
1658 const Array<OneD, const Array<OneD, NekDouble>> &normal_3 =
1659 GetTraceNormal(3);
1660
1661 int ncoords = normal_0.size();
1662
1664 {
1665 // needs checking for 3D coords
1666
1667 // factors 0 and 2
1668 for (int i = 0; i < nquad0; ++i)
1669 {
1670 d0factors[0][i] = df[0][i] * normal_0[0][i];
1671 d0factors[2][i] = df[0][nquad0 * (nquad1 - 1) + i] * normal_2[0][i];
1672
1673 d1factors[0][i] = df[1][i] * normal_0[0][i];
1674 d1factors[2][i] = df[1][nquad0 * (nquad1 - 1) + i] * normal_2[0][i];
1675 }
1676
1677 for (int n = 1; n < ncoords; ++n)
1678 {
1679 // d xi_1/dy n_y
1680 // needs checking for 3D coords
1681 for (int i = 0; i < nquad0; ++i)
1682 {
1683 d0factors[0][i] += df[2 * n][i] * normal_0[n][i];
1684 d0factors[2][i] +=
1685 df[2 * n][nquad0 * (nquad1 - 1) + i] * normal_2[n][i];
1686
1687 d1factors[0][i] += df[2 * n + 1][i] * normal_0[n][i];
1688 d1factors[2][i] +=
1689 df[2 * n + 1][nquad0 * (nquad1 - 1) + i] * normal_2[n][i];
1690 }
1691 }
1692
1693 // faces 1 and 3
1694 for (int i = 0; i < nquad1; ++i)
1695 {
1696 d0factors[1][i] = df[0][(i + 1) * nquad0 - 1] * normal_1[0][i];
1697 d0factors[3][i] = df[0][i * nquad0] * normal_3[0][i];
1698
1699 d1factors[1][i] = df[1][(i + 1) * nquad0 - 1] * normal_1[0][i];
1700 d1factors[3][i] = df[1][i * nquad0] * normal_3[0][i];
1701 }
1702
1703 for (int n = 1; n < ncoords; ++n)
1704 {
1705 for (int i = 0; i < nquad1; ++i)
1706 {
1707 d0factors[1][i] +=
1708 df[2 * n][(i + 1) * nquad0 - 1] * normal_1[n][i];
1709 d0factors[3][i] += df[2 * n][i * nquad0] * normal_3[n][i];
1710
1711 d1factors[1][i] +=
1712 df[2 * n + 1][(i + 1) * nquad0 - 1] * normal_1[n][i];
1713 d1factors[3][i] += df[2 * n + 1][i * nquad0] * normal_3[n][i];
1714 }
1715 }
1716 }
1717 else
1718 {
1719 // d xi_2/dx n_x
1720 for (int i = 0; i < nquad0; ++i)
1721 {
1722 d1factors[0][i] = df[1][0] * normal_0[0][i];
1723 d1factors[2][i] = df[1][0] * normal_2[0][i];
1724 }
1725
1726 // d xi_1/dx n_x
1727 for (int i = 0; i < nquad1; ++i)
1728 {
1729 d0factors[1][i] = df[0][0] * normal_1[0][i];
1730 d0factors[3][i] = df[0][0] * normal_3[0][i];
1731 }
1732
1733 for (int n = 1; n < ncoords; ++n)
1734 {
1735 // d xi_2/dy n_y
1736 // needs checking for 3D coords
1737 for (int i = 0; i < nquad0; ++i)
1738 {
1739 d1factors[0][i] += df[2 * n + 1][0] * normal_0[n][i];
1740 d1factors[2][i] += df[2 * n + 1][0] * normal_2[n][i];
1741 }
1742
1743 // d xi_1/dy n_y
1744 // needs checking for 3D coords
1745 for (int i = 0; i < nquad1; ++i)
1746 {
1747 d0factors[1][i] += df[2 * n][0] * normal_1[n][i];
1748 d0factors[3][i] += df[2 * n][0] * normal_3[n][i];
1749 }
1750 }
1751
1752 // d1factors
1753 // d xi_1/dx n_x
1754 for (int i = 0; i < nquad0; ++i)
1755 {
1756 d0factors[0][i] = df[0][0] * normal_0[0][i];
1757 d0factors[2][i] = df[0][0] * normal_2[0][i];
1758 }
1759
1760 // d xi_2/dx n_x
1761 for (int i = 0; i < nquad1; ++i)
1762 {
1763 d1factors[1][i] = df[1][0] * normal_1[0][i];
1764 d1factors[3][i] = df[1][0] * normal_3[0][i];
1765 }
1766
1767 for (int n = 1; n < ncoords; ++n)
1768 {
1769 // d xi_1/dy n_y
1770 // needs checking for 3D coords
1771 for (int i = 0; i < nquad0; ++i)
1772 {
1773 d0factors[0][i] += df[2 * n][0] * normal_0[n][i];
1774 d0factors[2][i] += df[2 * n][0] * normal_2[n][i];
1775 }
1776
1777 // d xi_2/dy n_y
1778 // needs checking for 3D coords
1779 for (int i = 0; i < nquad1; ++i)
1780 {
1781 d1factors[1][i] += df[2 * n + 1][0] * normal_1[n][i];
1782 d1factors[3][i] += df[2 * n + 1][0] * normal_3[n][i];
1783 }
1784 }
1785 }
1786}
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 514 of file QuadExp.cpp.
517{
518 NormVectorIProductWRTBase(Fvec[0], Fvec[1], Fvec[2], outarray);
519}
void NormVectorIProductWRTBase(const Array< OneD, const NekDouble > &Fx, Array< OneD, NekDouble > &outarray)
Definition: StdExpansion.h:613
References Nektar::StdRegions::StdExpansion::NormVectorIProductWRTBase().
◆ v_NormVectorIProductWRTBase() [2/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 486 of file QuadExp.cpp.
490{
492 Array<OneD, NekDouble> Fn(nq);
493
494 const Array<OneD, const Array<OneD, NekDouble>> &normals =
495 GetLeftAdjacentElementExp()->GetTraceNormal(
496 GetLeftAdjacentElementTrace());
497
499 {
500 Vmath::Vvtvvtp(nq, &normals[0][0], 1, &Fx[0], 1, &normals[1][0], 1,
501 &Fy[0], 1, &Fn[0], 1);
502 Vmath::Vvtvp(nq, &normals[2][0], 1, &Fz[0], 1, &Fn[0], 1, &Fn[0], 1);
503 }
504 else
505 {
506 Vmath::Svtsvtp(nq, normals[0][0], &Fx[0], 1, normals[1][0], &Fy[0], 1,
507 &Fn[0], 1);
508 Vmath::Svtvp(nq, normals[2][0], &Fz[0], 1, &Fn[0], 1, &Fn[0], 1);
509 }
510
511 IProductWRTBase(Fn, outarray);
512}
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: SharedArray.hpp:888
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_PhysEvaluate() [1/2]
|
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 574 of file QuadExp.cpp.
576{
577 Array<OneD, NekDouble> Lcoord = Array<OneD, NekDouble>(2);
578
580 m_geom->GetLocCoords(coord, Lcoord);
581
582 return StdExpansion2D::v_PhysEvaluate(Lcoord, physvals);
583}
References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_PhysEvaluate() [2/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 585 of file QuadExp.cpp.
588{
589 Array<OneD, NekDouble> Lcoord(2);
591 m_geom->GetLocCoords(coord, Lcoord);
592 return StdQuadExp::v_PhysEvaluate(Lcoord, inarray, firstOrderDerivs);
593}
References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_ReduceOrderCoeffs()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1452 of file QuadExp.cpp.
1455{
1456 int n_coeffs = inarray.size();
1457
1458 Array<OneD, NekDouble> coeff(n_coeffs);
1459 Array<OneD, NekDouble> coeff_tmp(n_coeffs, 0.0);
1460 Array<OneD, NekDouble> tmp, tmp2;
1461
1464 int numMax = nmodes0;
1465
1466 Vmath::Vcopy(n_coeffs, inarray, 1, coeff_tmp, 1);
1467
1468 const LibUtilities::PointsKey Pkey0(nmodes0,
1470 const LibUtilities::PointsKey Pkey1(nmodes1,
1474 LibUtilities::BasisKey bortho0(LibUtilities::eOrtho_A, nmodes0, Pkey0);
1475 LibUtilities::BasisKey bortho1(LibUtilities::eOrtho_A, nmodes1, Pkey1);
1476
1477 LibUtilities::InterpCoeff2D(b0, b1, coeff_tmp, bortho0, bortho1, coeff);
1478
1479 Vmath::Zero(n_coeffs, coeff_tmp, 1);
1480
1481 int cnt = 0;
1482 for (int i = 0; i < numMin + 1; ++i)
1483 {
1484 Vmath::Vcopy(numMin, tmp = coeff + cnt, 1, tmp2 = coeff_tmp + cnt, 1);
1485
1486 cnt = i * numMax;
1487 }
1488
1489 LibUtilities::InterpCoeff2D(bortho0, bortho1, coeff_tmp, b0, b1, outarray);
1490}
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 566 of file QuadExp.cpp.
569{
570 // Evaluate point in local (eta) coordinates.
571 return StdExpansion2D::v_PhysEvaluate(Lcoord, physvals);
572}
◆ v_SVVLaplacianFilter()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdQuadExp.
Definition at line 1586 of file QuadExp.cpp.
1588{
1590
1591 // Calculate sqrt of the Jacobian
1593 Array<OneD, NekDouble> sqrt_jac(nq);
1595 {
1596 Vmath::Vsqrt(nq, jac, 1, sqrt_jac, 1);
1597 }
1598 else
1599 {
1601 }
1602
1603 // Multiply array by sqrt(Jac)
1604 Vmath::Vmul(nq, sqrt_jac, 1, array, 1, array, 1);
1605
1606 // Apply std region filter
1607 StdQuadExp::v_SVVLaplacianFilter(array, mkey);
1608
1609 // Divide by sqrt(Jac)
1610 Vmath::Vdiv(nq, array, 1, sqrt_jac, 1, array, 1);
1611}
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 1423 of file QuadExp.cpp.
1427{
1428 StdExpansion::WeakDerivMatrixOp_MatFree(i, inarray, outarray, mkey);
1429}
◆ v_WeakDirectionalDerivMatrixOp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1431 of file QuadExp.cpp.
1434{
1435 StdExpansion::WeakDirectionalDerivMatrixOp_MatFree(inarray, outarray, mkey);
1436}
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().
