#include <PyrExp.h>
Inheritance diagram for Nektar::LocalRegions::PyrExp:
Public Member Functions | |
| PyrExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::BasisKey &Bc, const SpatialDomains::PyrGeomSharedPtr &geom) | |
| Constructor using BasisKey class for quadrature points and order definition. More... | |
| PyrExp (const PyrExp &T) | |
| ~PyrExp () override=default | |
 Public Member Functions inherited from Nektar::StdRegions::StdPyrExp | |
| StdPyrExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::BasisKey &Bc) | |
| StdPyrExp (const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::BasisKey &Bc, NekDouble *coeffs, NekDouble *phys) | |
| StdPyrExp ()=default | |
| StdPyrExp (const StdPyrExp &T)=default | |
| ~StdPyrExp () override=default | |
| Public Member Functions inherited from Nektar::StdRegions::StdExpansion3D | |
| StdExpansion3D (int numcoeffs, const LibUtilities::BasisKey &Ba, const LibUtilities::BasisKey &Bb, const LibUtilities::BasisKey &Bc) | |
| StdExpansion3D ()=default | |
| StdExpansion3D (const StdExpansion3D &T)=default | |
| ~StdExpansion3D () override=default | |
| void | PhysTensorDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray_d1, Array< OneD, NekDouble > &outarray_d2, Array< OneD, NekDouble > &outarray_d3) |
| Calculate the 3D derivative in the local tensor/collapsed coordinate at the physical points. More... | |
| void | BwdTrans_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2) |
| void | IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2) |
| int | GetNedges () const |
| return the number of edges in 3D expansion More... | |
| int | GetEdgeNcoeffs (const int i) const |
| This function returns the number of expansion coefficients belonging to the i-th edge. More... | |
| void | GetEdgeInteriorToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards) |
| 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::Expansion3D | |
| Expansion3D (SpatialDomains::Geometry3DSharedPtr pGeom) | |
| ~Expansion3D () override=default | |
| void | SetTraceToGeomOrientation (Array< OneD, NekDouble > &inout) |
| Align trace orientation with the geometry orientation. More... | |
| void | SetFaceToGeomOrientation (const int face, Array< OneD, NekDouble > &inout) |
| Align face orientation with the geometry orientation. More... | |
| void | AddHDGHelmholtzFaceTerms (const NekDouble tau, const int edge, Array< OneD, NekDouble > &facePhys, const StdRegions::VarCoeffMap &dirForcing, Array< OneD, NekDouble > &outarray) |
| void | AddNormTraceInt (const int dir, Array< OneD, ExpansionSharedPtr > &FaceExp, Array< OneD, Array< OneD, NekDouble > > &faceCoeffs, Array< OneD, NekDouble > &outarray) |
| void | AddNormTraceInt (const int dir, Array< OneD, const NekDouble > &inarray, Array< OneD, ExpansionSharedPtr > &FaceExp, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs) |
| void | AddFaceBoundaryInt (const int face, ExpansionSharedPtr &FaceExp, Array< OneD, NekDouble > &facePhys, Array< OneD, NekDouble > &outarray, const StdRegions::VarCoeffMap &varcoeffs=StdRegions::NullVarCoeffMap) |
| SpatialDomains::Geometry3DSharedPtr | GetGeom3D () const |
| void | v_ReOrientTracePhysMap (const StdRegions::Orientation orient, Array< OneD, int > &idmap, const int nq0, const int nq1) override |
| void | v_NormVectorIProductWRTBase (const Array< OneD, const Array< OneD, NekDouble > > &Fvec, Array< OneD, NekDouble > &outarray) override |
| Array< OneD, unsigned int > | GetEdgeInverseBoundaryMap (int eid) |
| Array< OneD, unsigned int > | GetTraceInverseBoundaryMap (int fid, StdRegions::Orientation faceOrient=StdRegions::eNoOrientation, int P1=-1, int P2=-1) |
| void | GetInverseBoundaryMaps (Array< OneD, unsigned int > &vmap, Array< OneD, Array< OneD, unsigned int > > &emap, Array< OneD, Array< OneD, unsigned int > > &fmap) |
| DNekScalMatSharedPtr | CreateMatrix (const MatrixKey &mkey) |
| 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 |
| Integrate the physical point list inarray over pyramidic region and return the value. 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) override |
| Calculate the derivative of the physical points. More... | |
| void | v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Forward transform from physical quadrature space stored in inarray and evaluate the expansion coefficients and store in (this)->m_coeffs. More... | |
| 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*base2 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_IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Calculates the inner product \( I_{pqr} = (u,
\partial_{x_i} \phi_{pqr}) \). More... | |
| 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 |
| 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 |
| void | v_ExtractDataToCoeffs (const NekDouble *data, const std::vector< unsigned int > &nummodes, const int mode_offset, NekDouble *coeffs, std::vector< LibUtilities::BasisType > &fromType) override |
| NekDouble | v_StdPhysEvaluate (const Array< OneD, const NekDouble > &Lcoord, const Array< OneD, const NekDouble > &physvals) 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_GetTracePhysMap (const int face, Array< OneD, int > &outarray) override |
| void | v_ComputeTraceNormal (const int face) override |
| void | v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdRegions::StdMatrixKey &mkey) override |
| DNekMatSharedPtr | v_GenMatrix (const StdRegions::StdMatrixKey &mkey) override |
| DNekMatSharedPtr | v_CreateStdMatrix (const StdRegions::StdMatrixKey &mkey) override |
| DNekScalMatSharedPtr | v_GetLocMatrix (const MatrixKey &mkey) override |
| DNekScalBlkMatSharedPtr | v_GetLocStaticCondMatrix (const MatrixKey &mkey) override |
| void | v_DropLocMatrix (const MatrixKey &mkey) override |
| void | v_DropLocStaticCondMatrix (const MatrixKey &mkey) override |
| void | v_ComputeLaplacianMetric () override |
| void | v_NormalTraceDerivFactors (Array< OneD, Array< OneD, NekDouble > > &d0factors, Array< OneD, Array< OneD, NekDouble > > &d1factors, Array< OneD, Array< OneD, NekDouble > > &d2factors) override |
| : This method gets all of the factors which are required as part of the Gradient Jump Penalty stabilisation and involves the product of the normal and geometric factors along the element trace. More... | |
| Protected Member Functions inherited from Nektar::StdRegions::StdPyrExp | |
| void | v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2) 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) override |
| void | v_StdPhysDeriv (const int dir, 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_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Backward transformation is evaluated at the quadrature points. More... | |
| 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 > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2) override |
| void | v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Forward transform from physical quadrature space stored in inarray and evaluate the expansion coefficients and store in outarray. More... | |
| void | v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Inner product of inarray over region with respect to the expansion basis m_base[0]->GetBdata(),m_base[1]->GetBdata(), m_base[2]->GetBdata() and return in 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 > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2) 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_GetCoords (Array< OneD, NekDouble > &xi_x, Array< OneD, NekDouble > &xi_y, Array< OneD, NekDouble > &xi_z) override |
| void | v_FillMode (const int mode, Array< OneD, NekDouble > &outarray) override |
| void | v_GetTraceNumModes (const int fid, int &numModes0, int &numModes1, Orientation faceOrient=eDir1FwdDir1_Dir2FwdDir2) override |
| NekDouble | v_PhysEvaluateBasis (const Array< OneD, const NekDouble > &coords, int mode) final |
| NekDouble | v_PhysEvaluate (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) override |
| int | v_GetNverts () const override |
| int | v_GetNedges () const override |
| int | v_GetNtraces () const override |
| LibUtilities::ShapeType | v_DetShapeType () const override |
| int | v_NumBndryCoeffs () const override |
| int | v_NumDGBndryCoeffs () const override |
| int | v_GetTraceNcoeffs (const int i) const override |
| int | v_GetTraceIntNcoeffs (const int i) const override |
| int | v_GetTraceNumPoints (const int i) const override |
| int | v_GetEdgeNcoeffs (const int i) const override |
| int | v_CalcNumberOfCoefficients (const std::vector< unsigned int > &nummodes, int &modes_offset) override |
| const LibUtilities::BasisKey | v_GetTraceBasisKey (const int i, const int k) const override |
| int | v_GetVertexMap (int localVertexId, bool useCoeffPacking=false) override |
| void | v_GetInteriorMap (Array< OneD, unsigned int > &outarray) override |
| void | v_GetBoundaryMap (Array< OneD, unsigned int > &outarray) override |
| void | v_GetTraceCoeffMap (const unsigned int fid, Array< OneD, unsigned int > &maparray) override |
| void | v_GetElmtTraceToTraceMap (const unsigned int fid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation faceOrient, int P, int Q) override |
| void | v_GetEdgeInteriorToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation traceOrient=eDir1FwdDir1_Dir2FwdDir2) override |
| void | v_GetTraceInteriorToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, const Orientation traceOrient=eDir1FwdDir1_Dir2FwdDir2) override |
| DNekMatSharedPtr | v_GenMatrix (const StdMatrixKey &mkey) override |
| DNekMatSharedPtr | v_CreateStdMatrix (const StdMatrixKey &mkey) override |
| void | v_SVVLaplacianFilter (Array< OneD, NekDouble > &array, const StdMatrixKey &mkey) override |
| void | v_ReduceOrderCoeffs (int numMin, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override |
| Protected Member Functions inherited from Nektar::StdRegions::StdExpansion3D | |
| 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 > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2)=0 |
| virtual void | v_IProductWRTBase_SumFacKernel (const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2)=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 |
| NekDouble | v_Integral (const Array< OneD, const NekDouble > &inarray) override |
| Integrates the specified function over the domain. More... | |
| virtual int | v_GetNedges (void) const |
| virtual int | v_GetEdgeNcoeffs (const int i) const |
| NekDouble | BaryTensorDeriv (const Array< OneD, NekDouble > &coord, const Array< OneD, const NekDouble > &inarray, std::array< NekDouble, 3 > &firstOrderDerivs) |
| virtual void | v_GetEdgeInteriorToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient=eForwards) |
| void | v_GetTraceToElementMap (const int tid, Array< OneD, unsigned int > &maparray, Array< OneD, int > &signarray, Orientation traceOrient, int P, int Q) 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::Expansion3D | |
| void | v_DGDeriv (const int dir, const Array< OneD, const NekDouble > &incoeffs, Array< OneD, ExpansionSharedPtr > &FaceExp, Array< OneD, Array< OneD, NekDouble > > &faceCoeffs, Array< OneD, NekDouble > &out_d) override |
| Evaluate coefficients of weak deriviative in the direction dir given the input coefficicents incoeffs and the imposed boundary values in EdgeExp (which will have its phys space updated). More... | |
| DNekMatSharedPtr | v_GenMatrix (const StdRegions::StdMatrixKey &mkey) override |
| void | v_AddFaceNormBoundaryInt (const int face, const ExpansionSharedPtr &FaceExp, const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray) override |
| void | v_AddRobinMassMatrix (const int face, const Array< OneD, const NekDouble > &primCoeffs, DNekMatSharedPtr &inoutmat) override |
| StdRegions::Orientation | v_GetTraceOrient (int face) override |
| void | v_GetTracePhysVals (const int face, const StdRegions::StdExpansionSharedPtr &FaceExp, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, StdRegions::Orientation orient) override |
| Extract the physical values along face face from inarray into outarray following the local face orientation and point distribution defined by defined in FaceExp. More... | |
| void | v_GenTraceExp (const int traceid, ExpansionSharedPtr &exp) override |
| void | GetPhysFaceVarCoeffsFromElement (const int face, ExpansionSharedPtr &FaceExp, const Array< OneD, const NekDouble > &varcoeff, Array< OneD, NekDouble > &outarray) |
| DNekMatSharedPtr | v_BuildTransformationMatrix (const DNekScalMatSharedPtr &r_bnd, const StdRegions::MatrixType matrixType) override |
| DNekMatSharedPtr | v_BuildInverseTransformationMatrix (const DNekScalMatSharedPtr &transformationmatrix) override |
| Build inverse and inverse transposed transformation matrix: \(\mathbf{R^{-1}}\) and \(\mathbf{R^{-T}}\). More... | |
| DNekMatSharedPtr | v_BuildVertexMatrix (const DNekScalMatSharedPtr &r_bnd) 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 | v_LaplacianMatrixOp_MatFree_Kernel (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp) override |
Detailed Description
Constructor & Destructor Documentation
◆ PyrExp() [1/2]
| Nektar::LocalRegions::PyrExp::PyrExp | ( | const LibUtilities::BasisKey & | Ba, |
| const LibUtilities::BasisKey & | Bb, | ||
| const LibUtilities::BasisKey & | Bc, | ||
| const SpatialDomains::PyrGeomSharedPtr & | geom | ||
| ) |
Constructor using BasisKey class for quadrature points and order definition.
Definition at line 43 of file PyrExp.cpp.
48 Ba.GetNumModes(), Bb.GetNumModes(), Bc.GetNumModes()),
49 3, Ba, Bb, Bc),
51 Ba.GetNumModes(), Bb.GetNumModes(), Bc.GetNumModes()),
52 Ba, Bb, Bc),
54 m_matrixManager(
56 std::string("PyrExpMatrix")),
58 this, std::placeholders::_1),
59 std::string("PyrExpStaticCondMatrix"))
60{
61}
Expansion3D(SpatialDomains::Geometry3DSharedPtr pGeom)
Definition: Expansion3D.h:59
DNekScalMatSharedPtr CreateMatrix(const MatrixKey &mkey)
Definition: Expansion3D.cpp:407
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: PyrExp.h:174
LibUtilities::NekManager< MatrixKey, DNekScalBlkMat, MatrixKey::opLess > m_staticCondMatrixManager
Definition: PyrExp.h:176
StdExpansion3D()=default
StdPyrExp()=default
int getNumberOfCoefficients(int Na, int Nb, int Nc)
Definition: ShapeType.hpp:232
◆ PyrExp() [2/2]
| Nektar::LocalRegions::PyrExp::PyrExp | ( | const PyrExp & | T | ) |
Definition at line 63 of file PyrExp.cpp.
◆ ~PyrExp()
|
overridedefault |
Member Function Documentation
◆ v_AlignVectorToCollapsedDir()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 394 of file PyrExp.cpp.
397{
403 const int nqtot = nquad0 * nquad1 * nquad2;
404
408
409 Array<OneD, NekDouble> gfac0(nquad0);
410 Array<OneD, NekDouble> gfac1(nquad1);
411 Array<OneD, NekDouble> gfac2(nquad2);
412 Array<OneD, NekDouble> tmp5(nqtot);
413 Array<OneD, NekDouble> wsp(
414 std::max(nqtot, order0 * nquad2 * (nquad1 + order1)));
415
416 Array<OneD, NekDouble> tmp2 = outarray[0];
417 Array<OneD, NekDouble> tmp3 = outarray[1];
418 Array<OneD, NekDouble> tmp4 = outarray[2];
419
420 const Array<TwoD, const NekDouble> &df =
422
423 Array<OneD, NekDouble> tmp1;
424 tmp1 = inarray;
425
427 {
428 Vmath::Vmul(nqtot, &df[3 * dir][0], 1, tmp1.get(), 1, tmp2.get(), 1);
429 Vmath::Vmul(nqtot, &df[3 * dir + 1][0], 1, tmp1.get(), 1, tmp3.get(),
430 1);
431 Vmath::Vmul(nqtot, &df[3 * dir + 2][0], 1, tmp1.get(), 1, tmp4.get(),
432 1);
433 }
434 else
435 {
436 Vmath::Smul(nqtot, df[3 * dir][0], tmp1.get(), 1, tmp2.get(), 1);
437 Vmath::Smul(nqtot, df[3 * dir + 1][0], tmp1.get(), 1, tmp3.get(), 1);
438 Vmath::Smul(nqtot, df[3 * dir + 2][0], tmp1.get(), 1, tmp4.get(), 1);
439 }
440
441 // set up geometric factor: (1+z0)/2
442 for (int i = 0; i < nquad0; ++i)
443 {
444 gfac0[i] = 0.5 * (1 + z0[i]);
445 }
446
447 // set up geometric factor: (1+z1)/2
448 for (int i = 0; i < nquad1; ++i)
449 {
450 gfac1[i] = 0.5 * (1 + z1[i]);
451 }
452
453 // Set up geometric factor: 2/(1-z2)
454 for (int i = 0; i < nquad2; ++i)
455 {
456 gfac2[i] = 2.0 / (1 - z2[i]);
457 }
458
459 const int nq01 = nquad0 * nquad1;
460
461 for (int i = 0; i < nquad2; ++i)
462 {
463 Vmath::Smul(nq01, gfac2[i], &tmp2[0] + i * nq01, 1, &tmp2[0] + i * nq01,
464 1); // 2/(1-z2) for d/dxi_0
465 Vmath::Smul(nq01, gfac2[i], &tmp3[0] + i * nq01, 1, &tmp3[0] + i * nq01,
466 1); // 2/(1-z2) for d/dxi_1
467 Vmath::Smul(nq01, gfac2[i], &tmp4[0] + i * nq01, 1, &tmp5[0] + i * nq01,
468 1); // 2/(1-z2) for d/dxi_2
469 }
470
471 // (1+z0)/(1-z2) for d/d eta_0
472 for (int i = 0; i < nquad1 * nquad2; ++i)
473 {
474 Vmath::Vmul(nquad0, &gfac0[0], 1, &tmp5[0] + i * nquad0, 1,
475 &wsp[0] + i * nquad0, 1);
476 }
477
478 Vmath::Vadd(nqtot, &tmp2[0], 1, &wsp[0], 1, &tmp2[0], 1);
479
480 // (1+z1)/(1-z2) for d/d eta_1
481 for (int i = 0; i < nquad1 * nquad2; ++i)
482 {
483 Vmath::Smul(nquad0, gfac1[i % nquad1], &tmp5[0] + i * nquad0, 1,
484 &tmp5[0] + i * nquad0, 1);
485 }
486 Vmath::Vadd(nqtot, &tmp3[0], 1, &tmp5[0], 1, &tmp3[0], 1);
487}
SpatialDomains::GeomFactorsSharedPtr m_metricinfo
Definition: Expansion.h:274
const LibUtilities::PointsKeyVector GetPointsKeys() const
Definition: StdExpansion.h:1095
Array< OneD, LibUtilities::BasisSharedPtr > m_base
Definition: StdExpansion.h:1173
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 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
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 Nektar::SpatialDomains::eDeformed, Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Vmath::Smul(), Vmath::Vadd(), and Vmath::Vmul().
Referenced by v_IProductWRTDerivBase_SumFac().
◆ v_ComputeLaplacianMetric()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1095 of file PyrExp.cpp.
1096{
1098 {
1099 ComputeQuadratureMetric();
1100 }
1101
1102 int i, j;
1104 const unsigned int dim = 3;
1109
1110 for (unsigned int i = 0; i < dim; ++i)
1111 {
1112 for (unsigned int j = i; j < dim; ++j)
1113 {
1114 m_metrics[m[i][j]] = Array<OneD, NekDouble>(nqtot);
1115 }
1116 }
1117
1118 // Define shorthand synonyms for m_metrics storage
1119 Array<OneD, NekDouble> g0(m_metrics[m[0][0]]);
1120 Array<OneD, NekDouble> g1(m_metrics[m[1][1]]);
1121 Array<OneD, NekDouble> g2(m_metrics[m[2][2]]);
1122 Array<OneD, NekDouble> g3(m_metrics[m[0][1]]);
1123 Array<OneD, NekDouble> g4(m_metrics[m[0][2]]);
1124 Array<OneD, NekDouble> g5(m_metrics[m[1][2]]);
1125
1126 // Allocate temporary storage
1127 Array<OneD, NekDouble> alloc(9 * nqtot, 0.0);
1128 Array<OneD, NekDouble> h0(nqtot, alloc);
1129 Array<OneD, NekDouble> h1(nqtot, alloc + 1 * nqtot);
1130 Array<OneD, NekDouble> h2(nqtot, alloc + 2 * nqtot);
1131 Array<OneD, NekDouble> wsp1(nqtot, alloc + 3 * nqtot);
1132 Array<OneD, NekDouble> wsp2(nqtot, alloc + 4 * nqtot);
1133 Array<OneD, NekDouble> wsp3(nqtot, alloc + 5 * nqtot);
1134 Array<OneD, NekDouble> wsp4(nqtot, alloc + 6 * nqtot);
1135 Array<OneD, NekDouble> wsp5(nqtot, alloc + 7 * nqtot);
1136 Array<OneD, NekDouble> wsp6(nqtot, alloc + 8 * nqtot);
1137
1138 const Array<TwoD, const NekDouble> &df =
1146
1147 // Populate collapsed coordinate arrays h0, h1 and h2.
1148 for (j = 0; j < nquad2; ++j)
1149 {
1150 for (i = 0; i < nquad1; ++i)
1151 {
1152 Vmath::Fill(nquad0, 2.0 / (1.0 - z2[j]),
1153 &h0[0] + i * nquad0 + j * nquad0 * nquad1, 1);
1154 Vmath::Fill(nquad0, 1.0 / (1.0 - z2[j]),
1155 &h1[0] + i * nquad0 + j * nquad0 * nquad1, 1);
1156 Vmath::Fill(nquad0, (1.0 + z1[i]) / (1.0 - z2[j]),
1157 &h2[0] + i * nquad0 + j * nquad0 * nquad1, 1);
1158 }
1159 }
1160 for (i = 0; i < nquad0; i++)
1161 {
1162 Blas::Dscal(nquad1 * nquad2, 1 + z0[i], &h1[0] + i, nquad0);
1163 }
1164
1165 // Step 3. Construct combined metric terms for physical space to
1166 // collapsed coordinate system.
1167 // Order of construction optimised to minimise temporary storage
1169 {
1170 // f_{1k}
1171 Vmath::Vvtvvtp(nqtot, &df[0][0], 1, &h0[0], 1, &df[2][0], 1, &h1[0], 1,
1172 &wsp1[0], 1);
1173 Vmath::Vvtvvtp(nqtot, &df[3][0], 1, &h0[0], 1, &df[5][0], 1, &h1[0], 1,
1174 &wsp2[0], 1);
1175 Vmath::Vvtvvtp(nqtot, &df[6][0], 1, &h0[0], 1, &df[8][0], 1, &h1[0], 1,
1176 &wsp3[0], 1);
1177
1178 // g0
1179 Vmath::Vvtvvtp(nqtot, &wsp1[0], 1, &wsp1[0], 1, &wsp2[0], 1, &wsp2[0],
1180 1, &g0[0], 1);
1181 Vmath::Vvtvp(nqtot, &wsp3[0], 1, &wsp3[0], 1, &g0[0], 1, &g0[0], 1);
1182
1183 // g4
1184 Vmath::Vvtvvtp(nqtot, &df[2][0], 1, &wsp1[0], 1, &df[5][0], 1, &wsp2[0],
1185 1, &g4[0], 1);
1186 Vmath::Vvtvp(nqtot, &df[8][0], 1, &wsp3[0], 1, &g4[0], 1, &g4[0], 1);
1187
1188 // f_{2k}
1189 Vmath::Vvtvvtp(nqtot, &df[1][0], 1, &h0[0], 1, &df[2][0], 1, &h2[0], 1,
1190 &wsp4[0], 1);
1191 Vmath::Vvtvvtp(nqtot, &df[4][0], 1, &h0[0], 1, &df[5][0], 1, &h2[0], 1,
1192 &wsp5[0], 1);
1193 Vmath::Vvtvvtp(nqtot, &df[7][0], 1, &h0[0], 1, &df[8][0], 1, &h2[0], 1,
1194 &wsp6[0], 1);
1195
1196 // g1
1197 Vmath::Vvtvvtp(nqtot, &wsp4[0], 1, &wsp4[0], 1, &wsp5[0], 1, &wsp5[0],
1198 1, &g1[0], 1);
1199 Vmath::Vvtvp(nqtot, &wsp6[0], 1, &wsp6[0], 1, &g1[0], 1, &g1[0], 1);
1200
1201 // g3
1202 Vmath::Vvtvvtp(nqtot, &wsp1[0], 1, &wsp4[0], 1, &wsp2[0], 1, &wsp5[0],
1203 1, &g3[0], 1);
1204 Vmath::Vvtvp(nqtot, &wsp3[0], 1, &wsp6[0], 1, &g3[0], 1, &g3[0], 1);
1205
1206 // g5
1207 Vmath::Vvtvvtp(nqtot, &df[2][0], 1, &wsp4[0], 1, &df[5][0], 1, &wsp5[0],
1208 1, &g5[0], 1);
1209 Vmath::Vvtvp(nqtot, &df[8][0], 1, &wsp6[0], 1, &g5[0], 1, &g5[0], 1);
1210
1211 // g2
1212 Vmath::Vvtvvtp(nqtot, &df[2][0], 1, &df[2][0], 1, &df[5][0], 1,
1213 &df[5][0], 1, &g2[0], 1);
1214 Vmath::Vvtvp(nqtot, &df[8][0], 1, &df[8][0], 1, &g2[0], 1, &g2[0], 1);
1215 }
1216 else
1217 {
1218 // f_{1k}
1219 Vmath::Svtsvtp(nqtot, df[0][0], &h0[0], 1, df[2][0], &h1[0], 1,
1220 &wsp1[0], 1);
1221 Vmath::Svtsvtp(nqtot, df[3][0], &h0[0], 1, df[5][0], &h1[0], 1,
1222 &wsp2[0], 1);
1223 Vmath::Svtsvtp(nqtot, df[6][0], &h0[0], 1, df[8][0], &h1[0], 1,
1224 &wsp3[0], 1);
1225
1226 // g0
1227 Vmath::Vvtvvtp(nqtot, &wsp1[0], 1, &wsp1[0], 1, &wsp2[0], 1, &wsp2[0],
1228 1, &g0[0], 1);
1229 Vmath::Vvtvp(nqtot, &wsp3[0], 1, &wsp3[0], 1, &g0[0], 1, &g0[0], 1);
1230
1231 // g4
1232 Vmath::Svtsvtp(nqtot, df[2][0], &wsp1[0], 1, df[5][0], &wsp2[0], 1,
1233 &g4[0], 1);
1234 Vmath::Svtvp(nqtot, df[8][0], &wsp3[0], 1, &g4[0], 1, &g4[0], 1);
1235
1236 // f_{2k}
1237 Vmath::Svtsvtp(nqtot, df[1][0], &h0[0], 1, df[2][0], &h2[0], 1,
1238 &wsp4[0], 1);
1239 Vmath::Svtsvtp(nqtot, df[4][0], &h0[0], 1, df[5][0], &h2[0], 1,
1240 &wsp5[0], 1);
1241 Vmath::Svtsvtp(nqtot, df[7][0], &h0[0], 1, df[8][0], &h2[0], 1,
1242 &wsp6[0], 1);
1243
1244 // g1
1245 Vmath::Vvtvvtp(nqtot, &wsp4[0], 1, &wsp4[0], 1, &wsp5[0], 1, &wsp5[0],
1246 1, &g1[0], 1);
1247 Vmath::Vvtvp(nqtot, &wsp6[0], 1, &wsp6[0], 1, &g1[0], 1, &g1[0], 1);
1248
1249 // g3
1250 Vmath::Vvtvvtp(nqtot, &wsp1[0], 1, &wsp4[0], 1, &wsp2[0], 1, &wsp5[0],
1251 1, &g3[0], 1);
1252 Vmath::Vvtvp(nqtot, &wsp3[0], 1, &wsp6[0], 1, &g3[0], 1, &g3[0], 1);
1253
1254 // g5
1255 Vmath::Svtsvtp(nqtot, df[2][0], &wsp4[0], 1, df[5][0], &wsp5[0], 1,
1256 &g5[0], 1);
1257 Vmath::Svtvp(nqtot, df[8][0], &wsp6[0], 1, &g5[0], 1, &g5[0], 1);
1258
1259 // g2
1260 Vmath::Fill(nqtot,
1261 df[2][0] * df[2][0] + df[5][0] * df[5][0] +
1262 df[8][0] * df[8][0],
1263 &g2[0], 1);
1264 }
1265
1266 for (unsigned int i = 0; i < dim; ++i)
1267 {
1268 for (unsigned int j = i; j < dim; ++j)
1269 {
1271 }
1272 }
1273}
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
static void Dscal(const int &n, const double &alpha, double *x, const int &incx)
BLAS level 1: x = alpha x.
Definition: Blas.hpp:149
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
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 Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.hpp:54
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 Nektar::LocalRegions::Expansion::ComputeQuadratureMetric(), Blas::Dscal(), 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::StdRegions::StdExpansion::m_base, Nektar::LocalRegions::Expansion::m_metricinfo, Nektar::LocalRegions::Expansion::m_metrics, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), Vmath::Svtsvtp(), Vmath::Svtvp(), Vmath::Vvtvp(), and Vmath::Vvtvvtp().
◆ v_ComputeTraceNormal()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 729 of file PyrExp.cpp.
730{
732 GetGeom()->GetMetricInfo();
733
735 for (int i = 0; i < ptsKeys.size(); ++i)
736 {
737 // Need at least 2 points for computing normals
739 {
740 LibUtilities::PointsKey pKey(2, ptsKeys[i].GetPointsType());
741 ptsKeys[i] = pKey;
742 }
743 }
744
745 SpatialDomains::GeomType type = geomFactors->GetGtype();
746 const Array<TwoD, const NekDouble> &df =
747 geomFactors->GetDerivFactors(ptsKeys);
748 const Array<OneD, const NekDouble> &jac = geomFactors->GetJac(ptsKeys);
749
750 LibUtilities::BasisKey tobasis0 = GetTraceBasisKey(face, 0);
751 LibUtilities::BasisKey tobasis1 = GetTraceBasisKey(face, 1);
752
753 // Number of quadrature points in face expansion.
754 int nq_face = tobasis0.GetNumPoints() * tobasis1.GetNumPoints();
755
757 int i;
758
759 m_traceNormals[face] = Array<OneD, Array<OneD, NekDouble>>(vCoordDim);
760 Array<OneD, Array<OneD, NekDouble>> &normal = m_traceNormals[face];
761 for (i = 0; i < vCoordDim; ++i)
762 {
763 normal[i] = Array<OneD, NekDouble>(nq_face);
764 }
765
766 size_t nqb = nq_face;
767 size_t nbnd = face;
768 m_elmtBndNormDirElmtLen[nbnd] = Array<OneD, NekDouble>{nqb, 0.0};
769 Array<OneD, NekDouble> &length = m_elmtBndNormDirElmtLen[nbnd];
770
771 // Regular geometry case
773 type == SpatialDomains::eMovingRegular)
774 {
775 NekDouble fac;
776 // Set up normals
777 switch (face)
778 {
779 case 0:
780 {
781 for (i = 0; i < vCoordDim; ++i)
782 {
783 normal[i][0] = -df[3 * i + 2][0];
784 }
785 break;
786 }
787 case 1:
788 {
789 for (i = 0; i < vCoordDim; ++i)
790 {
791 normal[i][0] = -df[3 * i + 1][0];
792 }
793 break;
794 }
795 case 2:
796 {
797 for (i = 0; i < vCoordDim; ++i)
798 {
799 normal[i][0] = df[3 * i][0] + df[3 * i + 2][0];
800 }
801 break;
802 }
803 case 3:
804 {
805 for (i = 0; i < vCoordDim; ++i)
806 {
807 normal[i][0] = df[3 * i + 1][0] + df[3 * i + 2][0];
808 }
809 break;
810 }
811 case 4:
812 {
813 for (i = 0; i < vCoordDim; ++i)
814 {
815 normal[i][0] = -df[3 * i][0];
816 }
817 break;
818 }
819 default:
821 }
822
823 // Normalise resulting vector.
824 fac = 0.0;
825 for (i = 0; i < vCoordDim; ++i)
826 {
827 fac += normal[i][0] * normal[i][0];
828 }
829 fac = 1.0 / sqrt(fac);
830
831 Vmath::Fill(nqb, fac, length, 1);
832
833 for (i = 0; i < vCoordDim; ++i)
834 {
835 Vmath::Fill(nq_face, fac * normal[i][0], normal[i], 1);
836 }
837 }
838 else
839 {
840 // Set up deformed normals.
841 int j, k;
842
843 int nq0 = ptsKeys[0].GetNumPoints();
844 int nq1 = ptsKeys[1].GetNumPoints();
845 int nq2 = ptsKeys[2].GetNumPoints();
846 int nq01 = nq0 * nq1;
847 int nqtot;
848
849 // Determine number of quadrature points on the face.
850 if (face == 0)
851 {
852 nqtot = nq0 * nq1;
853 }
854 else if (face == 1 || face == 3)
855 {
856 nqtot = nq0 * nq2;
857 }
858 else
859 {
860 nqtot = nq1 * nq2;
861 }
862
863 LibUtilities::PointsKey points0;
864 LibUtilities::PointsKey points1;
865
866 Array<OneD, NekDouble> faceJac(nqtot);
867 Array<OneD, NekDouble> normals(vCoordDim * nqtot, 0.0);
868
869 // Extract Jacobian along face and recover local derivatives
870 // (dx/dr) for polynomial interpolation by multiplying m_gmat by
871 // jacobian
872 switch (face)
873 {
874 case 0:
875 {
876 for (j = 0; j < nq01; ++j)
877 {
878 normals[j] = -df[2][j] * jac[j];
879 normals[nqtot + j] = -df[5][j] * jac[j];
880 normals[2 * nqtot + j] = -df[8][j] * jac[j];
881 faceJac[j] = jac[j];
882 }
883
884 points0 = ptsKeys[0];
885 points1 = ptsKeys[1];
886 break;
887 }
888
889 case 1:
890 {
891 for (j = 0; j < nq0; ++j)
892 {
893 for (k = 0; k < nq2; ++k)
894 {
895 int tmp = j + nq01 * k;
896 normals[j + k * nq0] = -df[1][tmp] * jac[tmp];
897 normals[nqtot + j + k * nq0] = -df[4][tmp] * jac[tmp];
898 normals[2 * nqtot + j + k * nq0] =
899 -df[7][tmp] * jac[tmp];
900 faceJac[j + k * nq0] = jac[tmp];
901 }
902 }
903
904 points0 = ptsKeys[0];
905 points1 = ptsKeys[2];
906 break;
907 }
908
909 case 2:
910 {
911 for (j = 0; j < nq1; ++j)
912 {
913 for (k = 0; k < nq2; ++k)
914 {
915 int tmp = nq0 - 1 + nq0 * j + nq01 * k;
916 normals[j + k * nq1] =
917 (df[0][tmp] + df[2][tmp]) * jac[tmp];
918 normals[nqtot + j + k * nq1] =
919 (df[3][tmp] + df[5][tmp]) * jac[tmp];
920 normals[2 * nqtot + j + k * nq1] =
921 (df[6][tmp] + df[8][tmp]) * jac[tmp];
922 faceJac[j + k * nq1] = jac[tmp];
923 }
924 }
925
926 points0 = ptsKeys[1];
927 points1 = ptsKeys[2];
928 break;
929 }
930
931 case 3:
932 {
933 for (j = 0; j < nq0; ++j)
934 {
935 for (k = 0; k < nq2; ++k)
936 {
937 int tmp = nq0 * (nq1 - 1) + j + nq01 * k;
938 normals[j + k * nq0] =
939 (df[1][tmp] + df[2][tmp]) * jac[tmp];
940 normals[nqtot + j + k * nq0] =
941 (df[4][tmp] + df[5][tmp]) * jac[tmp];
942 normals[2 * nqtot + j + k * nq0] =
943 (df[7][tmp] + df[8][tmp]) * jac[tmp];
944 faceJac[j + k * nq0] = jac[tmp];
945 }
946 }
947
948 points0 = ptsKeys[0];
949 points1 = ptsKeys[2];
950 break;
951 }
952
953 case 4:
954 {
955 for (j = 0; j < nq1; ++j)
956 {
957 for (k = 0; k < nq2; ++k)
958 {
959 int tmp = j * nq0 + nq01 * k;
960 normals[j + k * nq1] = -df[0][tmp] * jac[tmp];
961 normals[nqtot + j + k * nq1] = -df[3][tmp] * jac[tmp];
962 normals[2 * nqtot + j + k * nq1] =
963 -df[6][tmp] * jac[tmp];
964 faceJac[j + k * nq1] = jac[tmp];
965 }
966 }
967
968 points0 = ptsKeys[1];
969 points1 = ptsKeys[2];
970 break;
971 }
972
973 default:
975 }
976
977 Array<OneD, NekDouble> work(nq_face, 0.0);
978 // Interpolate Jacobian and invert
979 LibUtilities::Interp2D(points0, points1, faceJac,
980 tobasis0.GetPointsKey(), tobasis1.GetPointsKey(),
981 work);
982 Vmath::Sdiv(nq_face, 1.0, &work[0], 1, &work[0], 1);
983
984 // Interpolate normal and multiply by inverse Jacobian.
985 for (i = 0; i < vCoordDim; ++i)
986 {
987 LibUtilities::Interp2D(points0, points1, &normals[i * nqtot],
988 tobasis0.GetPointsKey(),
989 tobasis1.GetPointsKey(), &normal[i][0]);
990 Vmath::Vmul(nq_face, work, 1, normal[i], 1, normal[i], 1);
991 }
992
993 // Normalise to obtain unit normals.
994 Vmath::Zero(nq_face, work, 1);
996 {
997 Vmath::Vvtvp(nq_face, normal[i], 1, normal[i], 1, work, 1, work, 1);
998 }
999
1000 Vmath::Vsqrt(nq_face, work, 1, work, 1);
1001 Vmath::Sdiv(nq_face, 1.0, work, 1, work, 1);
1002
1003 Vmath::Vcopy(nqb, work, 1, length, 1);
1004
1006 {
1007 Vmath::Vmul(nq_face, normal[i], 1, work, 1, normal[i], 1);
1008 }
1009 }
1010}
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
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 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
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 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 Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.hpp:825
References ASSERTL0, Nektar::SpatialDomains::eMovingRegular, Nektar::SpatialDomains::eRegular, Vmath::Fill(), Nektar::StdRegions::StdExpansion::GetCoordim(), Nektar::LocalRegions::Expansion::GetGeom(), Nektar::LibUtilities::BasisKey::GetNumPoints(), Nektar::StdRegions::StdExpansion::GetNumPoints(), Nektar::LibUtilities::BasisKey::GetPointsKey(), Nektar::StdRegions::StdExpansion::GetPointsKeys(), Nektar::StdRegions::StdExpansion::GetPointsType(), Nektar::StdRegions::StdExpansion::GetTraceBasisKey(), Nektar::LibUtilities::Interp2D(), Nektar::LocalRegions::Expansion::m_elmtBndNormDirElmtLen, Nektar::LocalRegions::Expansion::m_traceNormals, Vmath::Sdiv(), tinysimd::sqrt(), Vmath::Vcopy(), Vmath::Vmul(), Vmath::Vsqrt(), Vmath::Vvtvp(), and Vmath::Zero().
◆ v_CreateStdMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 1064 of file PyrExp.cpp.
1065{
1066 LibUtilities::BasisKey bkey0 = m_base[0]->GetBasisKey();
1067 LibUtilities::BasisKey bkey1 = m_base[1]->GetBasisKey();
1068 LibUtilities::BasisKey bkey2 = m_base[2]->GetBasisKey();
1069 StdRegions::StdPyrExpSharedPtr tmp =
1070 MemoryManager<StdPyrExp>::AllocateSharedPtr(bkey0, bkey1, bkey2);
1071
1072 return tmp->GetStdMatrix(mkey);
1073}
static std::shared_ptr< DataType > AllocateSharedPtr(const Args &...args)
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:166
std::shared_ptr< StdPyrExp > StdPyrExpSharedPtr
Definition: StdPyrExp.h:214
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_DropLocMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1080 of file PyrExp.cpp.
References m_matrixManager.
◆ v_DropLocStaticCondMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1090 of file PyrExp.cpp.
References m_staticCondMatrixManager.
◆ v_ExtractDataToCoeffs()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 541 of file PyrExp.cpp.
545{
546 int data_order0 = nummodes[mode_offset];
548 int data_order1 = nummodes[mode_offset + 1];
550 int fillorder1 = min(order1, data_order1);
551 int data_order2 = nummodes[mode_offset + 2];
553 int fillorder2 = min(order2, data_order2);
554
555 // Check if not same order or basis and if not make temp
556 // element to read in data
560 data_order1 != fillorder1 || data_order2 != fillorder2)
561 {
562 // Construct a pyr with the appropriate basis type at our
563 // quadrature points, and one more to do a forwards
564 // transform. We can then copy the output to coeffs.
565 StdRegions::StdPyrExp tmpPyr(
566 LibUtilities::BasisKey(fromType[0], data_order0,
567 m_base[0]->GetPointsKey()),
568 LibUtilities::BasisKey(fromType[1], data_order1,
569 m_base[1]->GetPointsKey()),
570 LibUtilities::BasisKey(fromType[2], data_order2,
571 m_base[2]->GetPointsKey()));
572
573 StdRegions::StdPyrExp tmpPyr2(m_base[0]->GetBasisKey(),
574 m_base[1]->GetBasisKey(),
575 m_base[2]->GetBasisKey());
576
577 Array<OneD, const NekDouble> tmpData(tmpPyr.GetNcoeffs(), data);
578 Array<OneD, NekDouble> tmpBwd(tmpPyr2.GetTotPoints());
579 Array<OneD, NekDouble> tmpOut(tmpPyr2.GetNcoeffs());
580
581 tmpPyr.BwdTrans(tmpData, tmpBwd);
582 tmpPyr2.FwdTrans(tmpBwd, tmpOut);
583 Vmath::Vcopy(tmpOut.size(), &tmpOut[0], 1, coeffs, 1);
584 }
585 else
586 {
588 }
589}
LibUtilities::BasisType GetBasisType(const int dir) const
This function returns the type of basis used in the dir direction.
Definition: StdExpansion.h:156
References Nektar::StdRegions::StdExpansion::BwdTrans(), Nektar::StdRegions::StdExpansion::FwdTrans(), Nektar::StdRegions::StdExpansion::GetBasisType(), Nektar::StdRegions::StdExpansion::GetNcoeffs(), Nektar::StdRegions::StdExpansion::GetTotPoints(), Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::m_ncoeffs, and Vmath::Vcopy().
◆ v_FwdTrans()
|
overrideprotectedvirtual |
Forward transform from physical quadrature space stored in inarray and evaluate the expansion coefficients and store in (this)->m_coeffs.
Inputs:
- inarray: array of physical quadrature points to be transformed
Outputs:
- (this)->_coeffs: updated array of expansion coefficients.
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 221 of file PyrExp.cpp.
223{
225 m_base[2]->Collocation())
226 {
228 }
229 else
230 {
231 v_IProductWRTBase(inarray, outarray);
232
233 // get Mass matrix inverse
236
237 // copy inarray in case inarray == outarray
240
241 out = (*matsys) * in;
242 }
243}
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*base2 and put into out...
Definition: PyrExp.cpp:276
int GetNcoeffs(void) const
This function returns the total number of coefficients used in the expansion.
Definition: StdExpansion.h:124
LibUtilities::ShapeType DetShapeType() const
This function returns the shape of the expansion domain.
Definition: StdExpansion.h:367
std::shared_ptr< DNekScalMat > DNekScalMatSharedPtr
Definition: NekMatrixFwd.hpp:66
References Nektar::StdRegions::StdExpansion::DetShapeType(), Nektar::StdRegions::eInvMass, Nektar::eWrapper, Nektar::StdRegions::StdExpansion::GetNcoeffs(), Nektar::StdRegions::StdExpansion::m_base, m_matrixManager, Nektar::StdRegions::StdExpansion::m_ncoeffs, v_IProductWRTBase(), and Vmath::Vcopy().
◆ v_GenMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 1043 of file PyrExp.cpp.
1044{
1045 DNekMatSharedPtr returnval;
1046
1047 switch (mkey.GetMatrixType())
1048 {
1055 returnval = Expansion3D::v_GenMatrix(mkey);
1056 break;
1057 default:
1058 returnval = StdPyrExp::v_GenMatrix(mkey);
1059 }
1060
1061 return returnval;
1062}
DNekMatSharedPtr v_GenMatrix(const StdRegions::StdMatrixKey &mkey) override
Definition: Expansion3D.cpp:876
References Nektar::StdRegions::eHybridDGHelmBndLam, Nektar::StdRegions::eHybridDGHelmholtz, Nektar::StdRegions::eHybridDGLamToQ0, Nektar::StdRegions::eHybridDGLamToQ1, Nektar::StdRegions::eHybridDGLamToQ2, Nektar::StdRegions::eHybridDGLamToU, Nektar::StdRegions::StdMatrixKey::GetMatrixType(), and Nektar::LocalRegions::Expansion3D::v_GenMatrix().
◆ v_GetCoord()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 517 of file PyrExp.cpp.
519{
520 int i;
521
522 ASSERTL1(Lcoords[0] <= -1.0 && Lcoords[0] >= 1.0 && Lcoords[1] <= -1.0 &&
523 Lcoords[1] >= 1.0 && Lcoords[2] <= -1.0 && Lcoords[2] >= 1.0,
524 "Local coordinates are not in region [-1,1]");
525
526 // m_geom->FillGeom(); // TODO: implement FillGeom()
527
529 {
530 coords[i] = m_geom->GetCoord(i, Lcoords);
531 }
532}
#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
References ASSERTL1, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_GetCoords()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 534 of file PyrExp.cpp.
537{
538 Expansion::v_GetCoords(coords_1, coords_2, coords_3);
539}
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 500 of file PyrExp.cpp.
501{
503 m_base[0]->GetPointsKey());
505 m_base[1]->GetPointsKey());
507 m_base[2]->GetPointsKey());
508
510 bkey2);
511}
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 1075 of file PyrExp.cpp.
References m_matrixManager.
◆ v_GetLocStaticCondMatrix()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1085 of file PyrExp.cpp.
References m_staticCondMatrixManager.
◆ v_GetStdExp()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 493 of file PyrExp.cpp.
References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Nektar::StdRegions::StdExpansion::m_base.
◆ v_GetTracePhysMap()
|
overrideprotectedvirtual |
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 631 of file PyrExp.cpp.
632{
636
637 int nq0 = 0;
638 int nq1 = 0;
639
640 switch (face)
641 {
642 case 0:
643 nq0 = nquad0;
644 nq1 = nquad1;
645 if (outarray.size() != nq0 * nq1)
646 {
647 outarray = Array<OneD, int>(nq0 * nq1);
648 }
649
650 // Directions A and B positive
651 for (int i = 0; i < nquad0 * nquad1; ++i)
652 {
653 outarray[i] = i;
654 }
655
656 break;
657 case 1:
658 nq0 = nquad0;
659 nq1 = nquad2;
660 if (outarray.size() != nq0 * nq1)
661 {
662 outarray = Array<OneD, int>(nq0 * nq1);
663 }
664
665 // Direction A and B positive
666 for (int k = 0; k < nquad2; k++)
667 {
668 for (int i = 0; i < nquad0; ++i)
669 {
670 outarray[k * nquad0 + i] = (nquad0 * nquad1 * k) + i;
671 }
672 }
673
674 break;
675 case 2:
676 nq0 = nquad1;
677 nq1 = nquad2;
678 if (outarray.size() != nq0 * nq1)
679 {
680 outarray = Array<OneD, int>(nq0 * nq1);
681 }
682
683 // Directions A and B positive
684 for (int j = 0; j < nquad1 * nquad2; ++j)
685 {
686 outarray[j] = nquad0 - 1 + j * nquad0;
687 }
688 break;
689 case 3:
690
691 nq0 = nquad0;
692 nq1 = nquad2;
693 if (outarray.size() != nq0 * nq1)
694 {
695 outarray = Array<OneD, int>(nq0 * nq1);
696 }
697
698 // Direction A and B positive
699 for (int k = 0; k < nquad2; k++)
700 {
701 for (int i = 0; i < nquad0; ++i)
702 {
703 outarray[k * nquad0 + i] =
704 nquad0 * (nquad1 - 1) + (nquad0 * nquad1 * k) + i;
705 }
706 }
707 break;
708 case 4:
709 nq0 = nquad1;
710 nq1 = nquad2;
711
712 if (outarray.size() != nq0 * nq1)
713 {
714 outarray = Array<OneD, int>(nq0 * nq1);
715 }
716
717 // Directions A and B positive
718 for (int j = 0; j < nquad1 * nquad2; ++j)
719 {
720 outarray[j] = j * nquad0;
721 }
722 break;
723 default:
725 break;
726 }
727}
References ASSERTL0, and Nektar::StdRegions::StdExpansion::m_base.
◆ v_Integral()
|
overrideprotectedvirtual |
Integrate the physical point list inarray over pyramidic region and return the value.
Inputs:
- inarray: definition of function to be returned at quadrature point of expansion.
Outputs:
- returns \(\int^1_{-1}\int^1_{-1}\int^1_{-1} u(\bar \eta_1,
\eta_2, \eta_3) J[i,j,k] d \bar \eta_1 d \eta_2 d \eta_3\)
\(= \sum_{i=0}^{Q_1 - 1} \sum_{j=0}^{Q_2 - 1} \sum_{k=0}^{Q_3 - 1} u(\bar \eta_{1i}^{0,0}, \eta_{2j}^{0,0},\eta_{3k}^{2,0})w_{i}^{0,0} w_{j}^{0,0} \hat w_{k}^{2,0} \)
where \(inarray[i,j, k] = u(\bar \eta_{1i},\eta_{2j}, \eta_{3k}) \),
\(\hat w_{k}^{2,0} = \frac {w^{2,0}} {2} \)
and \( J[i,j,k] \) is the Jacobian evaluated at the quadrature point.
Reimplemented from Nektar::StdRegions::StdExpansion3D.
Definition at line 94 of file PyrExp.cpp.
95{
100 Array<OneD, NekDouble> tmp(nquad0 * nquad1 * nquad2);
101
102 // multiply inarray with Jacobian
104 {
105 Vmath::Vmul(nquad0 * nquad1 * nquad2, &jac[0], 1,
106 (NekDouble *)&inarray[0], 1, &tmp[0], 1);
107 }
108 else
109 {
111 (NekDouble *)&inarray[0], 1, &tmp[0], 1);
112 }
113
114 // call StdPyrExp version;
115 return StdPyrExp::v_Integral(tmp);
116}
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*base2 and put into outarray:
\( \begin{array}{rcl} I_{pqr} = (\phi_{pqr}, u)_{\delta} & = &
\sum_{i=0}^{nq_0} \sum_{j=0}^{nq_1} \sum_{k=0}^{nq_2} \psi_{p}^{a}
(\bar \eta_{1i}) \psi_{q}^{a} (\eta_{2j}) \psi_{pqr}^{c}
(\eta_{3k}) w_i w_j w_k u(\bar \eta_{1,i} \eta_{2,j} \eta_{3,k})
J_{i,j,k}\\ & = & \sum_{i=0}^{nq_0} \psi_p^a(\bar \eta_{1,i})
\sum_{j=0}^{nq_1} \psi_{q}^a(\eta_{2,j}) \sum_{k=0}^{nq_2}
\psi_{pqr}^c u(\bar \eta_{1i},\eta_{2j},\eta_{3k}) J_{i,j,k}
\end{array} \)
where
\(\phi_{pqr} (\xi_1 , \xi_2 , \xi_3) = \psi_p^a (\bar \eta_1)
\psi_{q}^a (\eta_2) \psi_{pqr}^c (\eta_3) \)
which can be implemented as
\(f_{pqr} (\xi_{3k}) =
\sum_{k=0}^{nq_3} \psi_{pqr}^c u(\bar
\eta_{1i},\eta_{2j},\eta_{3k}) J_{i,j,k} = {\bf B_3 U} \)
\(
g_{pq} (\xi_{3k}) = \sum_{j=0}^{nq_1} \psi_{q}^a (\xi_{2j}) f_{pqr}
(\xi_{3k}) = {\bf B_2 F} \)
\( (\phi_{pqr}, u)_{\delta} =
\sum_{k=0}^{nq_0} \psi_{p}^a (\xi_{3k}) g_{pq} (\xi_{3k}) = {\bf
B_1 G} \)
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 276 of file PyrExp.cpp.
278{
279 v_IProductWRTBase_SumFac(inarray, outarray);
280}
void v_IProductWRTBase_SumFac(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool multiplybyweights=true) override
Definition: PyrExp.cpp:282
References v_IProductWRTBase_SumFac().
Referenced by v_FwdTrans().
◆ v_IProductWRTBase_SumFac()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 282 of file PyrExp.cpp.
285{
291
292 Array<OneD, NekDouble> wsp(order0 * nquad2 * (nquad1 + order1));
293
294 if (multiplybyweights)
295 {
296 Array<OneD, NekDouble> tmp(nquad0 * nquad1 * nquad2);
297
298 MultiplyByQuadratureMetric(inarray, tmp);
299
302 tmp, outarray, wsp, true, true, true);
303 }
304 else
305 {
308 inarray, outarray, wsp, true, true, true);
309 }
310}
void IProductWRTBase_SumFacKernel(const Array< OneD, const NekDouble > &base0, const Array< OneD, const NekDouble > &base1, const Array< OneD, const NekDouble > &base2, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, Array< OneD, NekDouble > &wsp, bool doCheckCollDir0, bool doCheckCollDir1, bool doCheckCollDir2)
Definition: StdExpansion3D.cpp:107
References Nektar::StdRegions::StdExpansion3D::IProductWRTBase_SumFacKernel(), Nektar::StdRegions::StdExpansion::m_base, and Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric().
Referenced by v_IProductWRTBase().
◆ v_IProductWRTDerivBase()
|
overrideprotectedvirtual |
Calculates the inner product \( I_{pqr} = (u, \partial_{x_i} \phi_{pqr}) \).
The derivative of the basis functions is performed using the chain rule in order to incorporate the geometric factors. Assuming that the basis functions are a tensor product \(\phi_{pqr}(\eta_1,\eta_2,\eta_3) = \phi_1(\eta_1)\phi_2(\eta_2)\phi_3(\eta_3)\), this yields the result
\[ I_{pqr} = \sum_{j=1}^3 \left(u, \frac{\partial u}{\partial \eta_j} \frac{\partial \eta_j}{\partial x_i}\right) \]
In the pyramid element, we must also incorporate a second set of geometric factors which incorporate the collapsed co-ordinate system, so that
\[ \frac{\partial\eta_j}{\partial x_i} = \sum_{k=1}^3 \frac{\partial\eta_j}{\partial\xi_k}\frac{\partial\xi_k}{\partial x_i} \]
These derivatives can be found on p152 of Sherwin & Karniadakis.
- Parameters
-
dir Direction in which to take the derivative. inarray The function \( u \). outarray Value of the inner product.
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 342 of file PyrExp.cpp.
345{
346 v_IProductWRTDerivBase_SumFac(dir, inarray, outarray);
347}
void v_IProductWRTDerivBase_SumFac(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray) override
Definition: PyrExp.cpp:349
References v_IProductWRTDerivBase_SumFac().
◆ v_IProductWRTDerivBase_SumFac()
|
overrideprotectedvirtual |
- Parameters
-
inarray Function evaluated at physical collocation points. outarray Inner product with respect to each basis function over the element.
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 349 of file PyrExp.cpp.
352{
358 const int nqtot = nquad0 * nquad1 * nquad2;
359
360 Array<OneD, NekDouble> tmp1(nqtot);
361 Array<OneD, NekDouble> tmp2(nqtot);
362 Array<OneD, NekDouble> tmp3(nqtot);
363 Array<OneD, NekDouble> tmp4(nqtot);
364 Array<OneD, NekDouble> tmp6(m_ncoeffs);
365 Array<OneD, NekDouble> wsp(
366 std::max(nqtot, order0 * nquad2 * (nquad1 + order1)));
367
368 MultiplyByQuadratureMetric(inarray, tmp1);
369
370 Array<OneD, Array<OneD, NekDouble>> tmp2D{3};
371 tmp2D[0] = tmp2;
372 tmp2D[1] = tmp3;
373 tmp2D[2] = tmp4;
374
375 PyrExp::v_AlignVectorToCollapsedDir(dir, tmp1, tmp2D);
376
378 m_base[2]->GetBdata(), tmp2, outarray, wsp,
379 false, true, true);
380
383 false, true);
384
386
389 true, false);
390
392}
void v_AlignVectorToCollapsedDir(const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray) override
Definition: PyrExp.cpp:394
References Nektar::StdRegions::StdExpansion3D::IProductWRTBase_SumFacKernel(), Nektar::StdRegions::StdExpansion::m_base, Nektar::StdRegions::StdExpansion::m_ncoeffs, Nektar::StdRegions::StdExpansion::MultiplyByQuadratureMetric(), v_AlignVectorToCollapsedDir(), and Vmath::Vadd().
Referenced by v_IProductWRTDerivBase().
◆ v_LaplacianMatrixOp_MatFree_Kernel()
|
overrideprivatevirtual |
Reimplemented from Nektar::StdRegions::StdExpansion.
Definition at line 1275 of file PyrExp.cpp.
1278{
1279 // This implementation is only valid when there are no coefficients
1280 // associated to the Laplacian operator
1282 {
1283 ComputeLaplacianMetric();
1284 }
1285
1289 int nqtot = nquad0 * nquad1 * nq2;
1290
1293
1294 const Array<OneD, const NekDouble> &base0 = m_base[0]->GetBdata();
1295 const Array<OneD, const NekDouble> &base1 = m_base[1]->GetBdata();
1296 const Array<OneD, const NekDouble> &base2 = m_base[2]->GetBdata();
1297 const Array<OneD, const NekDouble> &dbase0 = m_base[0]->GetDbdata();
1298 const Array<OneD, const NekDouble> &dbase1 = m_base[1]->GetDbdata();
1299 const Array<OneD, const NekDouble> &dbase2 = m_base[2]->GetDbdata();
1300 const Array<OneD, const NekDouble> &metric00 =
1301 m_metrics[eMetricLaplacian00];
1302 const Array<OneD, const NekDouble> &metric01 =
1303 m_metrics[eMetricLaplacian01];
1304 const Array<OneD, const NekDouble> &metric02 =
1305 m_metrics[eMetricLaplacian02];
1306 const Array<OneD, const NekDouble> &metric11 =
1307 m_metrics[eMetricLaplacian11];
1308 const Array<OneD, const NekDouble> &metric12 =
1309 m_metrics[eMetricLaplacian12];
1310 const Array<OneD, const NekDouble> &metric22 =
1311 m_metrics[eMetricLaplacian22];
1312
1313 // Allocate temporary storage
1314 Array<OneD, NekDouble> wsp0(2 * nqtot, wsp);
1315 Array<OneD, NekDouble> wsp1(nqtot, wsp + 1 * nqtot);
1316 Array<OneD, NekDouble> wsp2(nqtot, wsp + 2 * nqtot);
1317 Array<OneD, NekDouble> wsp3(nqtot, wsp + 3 * nqtot);
1318 Array<OneD, NekDouble> wsp4(nqtot, wsp + 4 * nqtot);
1319 Array<OneD, NekDouble> wsp5(nqtot, wsp + 5 * nqtot);
1320
1321 // LAPLACIAN MATRIX OPERATION
1322 // wsp1 = du_dxi1 = D_xi1 * inarray = D_xi1 * u
1323 // wsp2 = du_dxi2 = D_xi2 * inarray = D_xi2 * u
1324 // wsp2 = du_dxi3 = D_xi3 * inarray = D_xi3 * u
1325 StdExpansion3D::PhysTensorDeriv(inarray, wsp0, wsp1, wsp2);
1326
1327 // wsp0 = k = g0 * wsp1 + g1 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
1328 // wsp2 = l = g1 * wsp1 + g2 * wsp2 = g0 * du_dxi1 + g1 * du_dxi2
1329 // where g0, g1 and g2 are the metric terms set up in the GeomFactors class
1330 // especially for this purpose
1331 Vmath::Vvtvvtp(nqtot, &metric00[0], 1, &wsp0[0], 1, &metric01[0], 1,
1332 &wsp1[0], 1, &wsp3[0], 1);
1333 Vmath::Vvtvp(nqtot, &metric02[0], 1, &wsp2[0], 1, &wsp3[0], 1, &wsp3[0], 1);
1334 Vmath::Vvtvvtp(nqtot, &metric01[0], 1, &wsp0[0], 1, &metric11[0], 1,
1335 &wsp1[0], 1, &wsp4[0], 1);
1336 Vmath::Vvtvp(nqtot, &metric12[0], 1, &wsp2[0], 1, &wsp4[0], 1, &wsp4[0], 1);
1337 Vmath::Vvtvvtp(nqtot, &metric02[0], 1, &wsp0[0], 1, &metric12[0], 1,
1338 &wsp1[0], 1, &wsp5[0], 1);
1339 Vmath::Vvtvp(nqtot, &metric22[0], 1, &wsp2[0], 1, &wsp5[0], 1, &wsp5[0], 1);
1340
1341 // outarray = m = (D_xi1 * B)^T * k
1342 // wsp1 = n = (D_xi2 * B)^T * l
1343 IProductWRTBase_SumFacKernel(dbase0, base1, base2, wsp3, outarray, wsp0,
1344 false, true, true);
1345 IProductWRTBase_SumFacKernel(base0, dbase1, base2, wsp4, wsp2, wsp0, true,
1346 false, true);
1347 Vmath::Vadd(m_ncoeffs, wsp2.get(), 1, outarray.get(), 1, outarray.get(), 1);
1348 IProductWRTBase_SumFacKernel(base0, base1, dbase2, wsp5, wsp2, wsp0, true,
1349 true, false);
1350 Vmath::Vadd(m_ncoeffs, wsp2.get(), 1, outarray.get(), 1, outarray.get(), 1);
1351}
void ComputeLaplacianMetric()
Definition: Expansion.cpp:449
References ASSERTL1, Nektar::LocalRegions::Expansion::ComputeLaplacianMetric(), Nektar::LocalRegions::eMetricLaplacian00, Nektar::StdRegions::StdExpansion::m_base, and Nektar::LocalRegions::Expansion::m_metrics.
◆ v_NormalTraceDerivFactors()
|
overrideprotectedvirtual |
: This method gets all of the factors which are required as part of the Gradient Jump Penalty stabilisation and involves the product of the normal and geometric factors along the element trace.
Reimplemented from Nektar::LocalRegions::Expansion.
Definition at line 1358 of file PyrExp.cpp.
1362{
1366
1367 const Array<TwoD, const NekDouble> &df =
1369
1370 if (d0factors.size() != 5)
1371 {
1372 d0factors = Array<OneD, Array<OneD, NekDouble>>(5);
1373 d1factors = Array<OneD, Array<OneD, NekDouble>>(5);
1374 d2factors = Array<OneD, Array<OneD, NekDouble>>(5);
1375 }
1376
1377 if (d0factors[0].size() != nquad0 * nquad1)
1378 {
1379 d0factors[0] = Array<OneD, NekDouble>(nquad0 * nquad1);
1380 d1factors[0] = Array<OneD, NekDouble>(nquad0 * nquad1);
1381 d2factors[0] = Array<OneD, NekDouble>(nquad0 * nquad1);
1382 }
1383
1384 if (d0factors[1].size() != nquad0 * nquad2)
1385 {
1386 d0factors[1] = Array<OneD, NekDouble>(nquad0 * nquad2);
1387 d0factors[3] = Array<OneD, NekDouble>(nquad0 * nquad2);
1388 d1factors[1] = Array<OneD, NekDouble>(nquad0 * nquad2);
1389 d1factors[3] = Array<OneD, NekDouble>(nquad0 * nquad2);
1390 d2factors[1] = Array<OneD, NekDouble>(nquad0 * nquad2);
1391 d2factors[3] = Array<OneD, NekDouble>(nquad0 * nquad2);
1392 }
1393
1394 if (d0factors[2].size() != nquad1 * nquad2)
1395 {
1396 d0factors[2] = Array<OneD, NekDouble>(nquad1 * nquad2);
1397 d0factors[4] = Array<OneD, NekDouble>(nquad1 * nquad2);
1398 d1factors[2] = Array<OneD, NekDouble>(nquad1 * nquad2);
1399 d1factors[4] = Array<OneD, NekDouble>(nquad1 * nquad2);
1400 d2factors[2] = Array<OneD, NekDouble>(nquad1 * nquad2);
1401 d2factors[4] = Array<OneD, NekDouble>(nquad1 * nquad2);
1402 }
1403
1404 // Outwards normals
1405 const Array<OneD, const Array<OneD, NekDouble>> &normal_0 =
1406 GetTraceNormal(0);
1407 const Array<OneD, const Array<OneD, NekDouble>> &normal_1 =
1408 GetTraceNormal(1);
1409 const Array<OneD, const Array<OneD, NekDouble>> &normal_2 =
1410 GetTraceNormal(2);
1411 const Array<OneD, const Array<OneD, NekDouble>> &normal_3 =
1412 GetTraceNormal(3);
1413 const Array<OneD, const Array<OneD, NekDouble>> &normal_4 =
1414 GetTraceNormal(4);
1415
1416 int ncoords = normal_0.size();
1417
1418 // first gather together standard cartesian inner products
1420 {
1421 // face 0
1422 for (int i = 0; i < nquad0 * nquad1; ++i)
1423 {
1424 d0factors[0][i] = df[0][i] * normal_0[0][i];
1425 d1factors[0][i] = df[1][i] * normal_0[0][i];
1426 d2factors[0][i] = df[2][i] * normal_0[0][i];
1427 }
1428
1429 for (int n = 1; n < ncoords; ++n)
1430 {
1431 for (int i = 0; i < nquad0 * nquad1; ++i)
1432 {
1433 d0factors[0][i] += df[3 * n][i] * normal_0[n][i];
1434 d1factors[0][i] += df[3 * n + 1][i] * normal_0[n][i];
1435 d2factors[0][i] += df[3 * n + 2][i] * normal_0[n][i];
1436 }
1437 }
1438
1439 // faces 1 and 3
1440 for (int j = 0; j < nquad2; ++j)
1441 {
1442 for (int i = 0; i < nquad0; ++i)
1443 {
1444 d0factors[1][j * nquad0 + i] = df[0][j * nquad0 * nquad1 + i] *
1445 normal_1[0][j * nquad0 + i];
1446 d1factors[1][j * nquad0 + i] = df[1][j * nquad0 * nquad1 + i] *
1447 normal_1[0][j * nquad0 + i];
1448 d2factors[1][j * nquad0 + i] = df[2][j * nquad0 * nquad1 + i] *
1449 normal_1[0][j * nquad0 + i];
1450
1451 d0factors[3][j * nquad0 + i] =
1452 df[0][(j + 1) * nquad0 * nquad1 - nquad0 + i] *
1453 normal_3[0][j * nquad0 + i];
1454 d1factors[3][j * nquad0 + i] =
1455 df[1][(j + 1) * nquad0 * nquad1 - nquad0 + i] *
1456 normal_3[0][j * nquad0 + i];
1457 d2factors[3][j * nquad0 + i] =
1458 df[2][(j + 1) * nquad0 * nquad1 - nquad0 + i] *
1459 normal_3[0][j * nquad0 + i];
1460 }
1461 }
1462
1463 for (int n = 1; n < ncoords; ++n)
1464 {
1465 for (int j = 0; j < nquad2; ++j)
1466 {
1467 for (int i = 0; i < nquad0; ++i)
1468 {
1469 d0factors[1][j * nquad0 + i] +=
1470 df[3 * n][j * nquad0 * nquad1 + i] *
1471 normal_1[0][j * nquad0 + i];
1472 d1factors[1][j * nquad0 + i] +=
1473 df[3 * n + 1][j * nquad0 * nquad1 + i] *
1474 normal_1[0][j * nquad0 + i];
1475 d2factors[1][j * nquad0 + i] +=
1476 df[3 * n + 2][j * nquad0 * nquad1 + i] *
1477 normal_1[0][j * nquad0 + i];
1478
1479 d0factors[3][j * nquad0 + i] +=
1480 df[3 * n][(j + 1) * nquad0 * nquad1 - nquad0 + i] *
1481 normal_3[0][j * nquad0 + i];
1482 d1factors[3][j * nquad0 + i] +=
1483 df[3 * n + 1][(j + 1) * nquad0 * nquad1 - nquad0 + i] *
1484 normal_3[0][j * nquad0 + i];
1485 d2factors[3][j * nquad0 + i] +=
1486 df[3 * n + 2][(j + 1) * nquad0 * nquad1 - nquad0 + i] *
1487 normal_3[0][j * nquad0 + i];
1488 }
1489 }
1490 }
1491
1492 // faces 2 and 4
1493 for (int j = 0; j < nquad2; ++j)
1494 {
1495 for (int i = 0; i < nquad1; ++i)
1496 {
1497 d0factors[2][j * nquad1 + i] =
1498 df[0][j * nquad0 * nquad1 + (i + 1) * nquad0 - 1] *
1499 normal_2[0][j * nquad1 + i];
1500 d1factors[2][j * nquad1 + i] =
1501 df[1][j * nquad0 * nquad1 + (i + 1) * nquad0 - 1] *
1502 normal_2[0][j * nquad1 + i];
1503 d2factors[2][j * nquad1 + i] =
1504 df[2][j * nquad0 * nquad1 + (i + 1) * nquad0 - 1] *
1505 normal_2[0][j * nquad1 + i];
1506
1507 d0factors[4][j * nquad1 + i] =
1508 df[0][j * nquad0 * nquad1 + i * nquad0] *
1509 normal_4[0][j * nquad1 + i];
1510 d1factors[4][j * nquad1 + i] =
1511 df[1][j * nquad0 * nquad1 + i * nquad0] *
1512 normal_4[0][j * nquad1 + i];
1513 d2factors[4][j * nquad1 + i] =
1514 df[2][j * nquad0 * nquad1 + i * nquad0] *
1515 normal_4[0][j * nquad1 + i];
1516 }
1517 }
1518
1519 for (int n = 1; n < ncoords; ++n)
1520 {
1521 for (int j = 0; j < nquad2; ++j)
1522 {
1523 for (int i = 0; i < nquad1; ++i)
1524 {
1525 d0factors[2][j * nquad1 + i] +=
1526 df[3 * n][j * nquad0 * nquad1 + (i + 1) * nquad0 - 1] *
1527 normal_2[n][j * nquad1 + i];
1528 d1factors[2][j * nquad1 + i] +=
1529 df[3 * n + 1]
1530 [j * nquad0 * nquad1 + (i + 1) * nquad0 - 1] *
1531 normal_2[n][j * nquad1 + i];
1532 d2factors[2][j * nquad1 + i] +=
1533 df[3 * n + 2]
1534 [j * nquad0 * nquad1 + (i + 1) * nquad0 - 1] *
1535 normal_2[n][j * nquad1 + i];
1536
1537 d0factors[4][j * nquad1 + i] +=
1538 df[3 * n][i * nquad0 + j * nquad0 * nquad1] *
1539 normal_4[n][j * nquad1 + i];
1540 d1factors[4][j * nquad1 + i] +=
1541 df[3 * n + 1][i * nquad0 + j * nquad0 * nquad1] *
1542 normal_4[n][j * nquad1 + i];
1543 d2factors[4][j * nquad1 + i] +=
1544 df[3 * n + 2][i * nquad0 + j * nquad0 * nquad1] *
1545 normal_4[n][j * nquad1 + i];
1546 }
1547 }
1548 }
1549 }
1550 else
1551 {
1552 // Face 0
1553 for (int i = 0; i < nquad0 * nquad1; ++i)
1554 {
1555 d0factors[0][i] = df[0][0] * normal_0[0][i];
1556 d1factors[0][i] = df[1][0] * normal_0[0][i];
1557 d2factors[0][i] = df[2][0] * normal_0[0][i];
1558 }
1559
1560 for (int n = 1; n < ncoords; ++n)
1561 {
1562 for (int i = 0; i < nquad0 * nquad1; ++i)
1563 {
1564 d0factors[0][i] += df[3 * n][0] * normal_0[n][i];
1565 d1factors[0][i] += df[3 * n + 1][0] * normal_0[n][i];
1566 d2factors[0][i] += df[3 * n + 2][0] * normal_0[n][i];
1567 }
1568 }
1569
1570 // faces 1 and 3
1571 for (int i = 0; i < nquad0 * nquad2; ++i)
1572 {
1573 d0factors[1][i] = df[0][0] * normal_1[0][i];
1574 d0factors[3][i] = df[0][0] * normal_3[0][i];
1575
1576 d1factors[1][i] = df[1][0] * normal_1[0][i];
1577 d1factors[3][i] = df[1][0] * normal_3[0][i];
1578
1579 d2factors[1][i] = df[2][0] * normal_1[0][i];
1580 d2factors[3][i] = df[2][0] * normal_3[0][i];
1581 }
1582
1583 for (int n = 1; n < ncoords; ++n)
1584 {
1585 for (int i = 0; i < nquad0 * nquad2; ++i)
1586 {
1587 d0factors[1][i] += df[3 * n][0] * normal_1[n][i];
1588 d0factors[3][i] += df[3 * n][0] * normal_3[n][i];
1589
1590 d1factors[1][i] += df[3 * n + 1][0] * normal_1[n][i];
1591 d1factors[3][i] += df[3 * n + 1][0] * normal_3[n][i];
1592
1593 d2factors[1][i] += df[3 * n + 2][0] * normal_1[n][i];
1594 d2factors[3][i] += df[3 * n + 2][0] * normal_3[n][i];
1595 }
1596 }
1597
1598 // faces 2 and 4
1599 for (int i = 0; i < nquad1 * nquad2; ++i)
1600 {
1601 d0factors[2][i] = df[0][0] * normal_2[0][i];
1602 d0factors[4][i] = df[0][0] * normal_4[0][i];
1603
1604 d1factors[2][i] = df[1][0] * normal_2[0][i];
1605 d1factors[4][i] = df[1][0] * normal_4[0][i];
1606
1607 d2factors[2][i] = df[2][0] * normal_2[0][i];
1608 d2factors[4][i] = df[2][0] * normal_4[0][i];
1609 }
1610
1611 for (int n = 1; n < ncoords; ++n)
1612 {
1613 for (int i = 0; i < nquad1 * nquad2; ++i)
1614 {
1615 d0factors[2][i] += df[3 * n][0] * normal_2[n][i];
1616 d0factors[4][i] += df[3 * n][0] * normal_4[n][i];
1617
1618 d1factors[2][i] += df[3 * n + 1][0] * normal_2[n][i];
1619 d1factors[4][i] += df[3 * n + 1][0] * normal_4[n][i];
1620
1621 d2factors[2][i] += df[3 * n + 2][0] * normal_2[n][i];
1622 d2factors[4][i] += df[3 * n + 2][0] * normal_4[n][i];
1623 }
1624 }
1625 }
1626}
const NormalVector & GetTraceNormal(const int id)
Definition: Expansion.cpp:251
◆ v_PhysDeriv()
|
overrideprotectedvirtual |
Calculate the derivative of the physical points.
The derivative is evaluated at the nodal physical points. Derivatives with respect to the local Cartesian coordinates.
\(\begin{Bmatrix} \frac {\partial} {\partial \xi_1} \\ \frac {\partial} {\partial \xi_2} \\ \frac {\partial} {\partial \xi_3} \end{Bmatrix} = \begin{Bmatrix} \frac 2 {(1-\eta_3)} \frac \partial {\partial \bar \eta_1} \\ \frac {\partial} {\partial \xi_2} \ \ \frac {(1 + \bar \eta_1)} {(1 - \eta_3)} \frac \partial {\partial \bar \eta_1} + \frac {\partial} {\partial \eta_3} \end{Bmatrix}\)
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 122 of file PyrExp.cpp.
126{
130 Array<TwoD, const NekDouble> gmat =
132 Array<OneD, NekDouble> diff0(nquad0 * nquad1 * nquad2);
133 Array<OneD, NekDouble> diff1(nquad0 * nquad1 * nquad2);
134 Array<OneD, NekDouble> diff2(nquad0 * nquad1 * nquad2);
135
136 StdPyrExp::v_PhysDeriv(inarray, diff0, diff1, diff2);
137
139 {
140 if (out_d0.size())
141 {
142 Vmath::Vmul(nquad0 * nquad1 * nquad2, &gmat[0][0], 1, &diff0[0], 1,
143 &out_d0[0], 1);
144 Vmath::Vvtvp(nquad0 * nquad1 * nquad2, &gmat[1][0], 1, &diff1[0], 1,
145 &out_d0[0], 1, &out_d0[0], 1);
146 Vmath::Vvtvp(nquad0 * nquad1 * nquad2, &gmat[2][0], 1, &diff2[0], 1,
147 &out_d0[0], 1, &out_d0[0], 1);
148 }
149
150 if (out_d1.size())
151 {
152 Vmath::Vmul(nquad0 * nquad1 * nquad2, &gmat[3][0], 1, &diff0[0], 1,
153 &out_d1[0], 1);
154 Vmath::Vvtvp(nquad0 * nquad1 * nquad2, &gmat[4][0], 1, &diff1[0], 1,
155 &out_d1[0], 1, &out_d1[0], 1);
156 Vmath::Vvtvp(nquad0 * nquad1 * nquad2, &gmat[5][0], 1, &diff2[0], 1,
157 &out_d1[0], 1, &out_d1[0], 1);
158 }
159
160 if (out_d2.size())
161 {
162 Vmath::Vmul(nquad0 * nquad1 * nquad2, &gmat[6][0], 1, &diff0[0], 1,
163 &out_d2[0], 1);
164 Vmath::Vvtvp(nquad0 * nquad1 * nquad2, &gmat[7][0], 1, &diff1[0], 1,
165 &out_d2[0], 1, &out_d2[0], 1);
166 Vmath::Vvtvp(nquad0 * nquad1 * nquad2, &gmat[8][0], 1, &diff2[0], 1,
167 &out_d2[0], 1, &out_d2[0], 1);
168 }
169 }
170 else // regular geometry
171 {
172 if (out_d0.size())
173 {
174 Vmath::Smul(nquad0 * nquad1 * nquad2, gmat[0][0], &diff0[0], 1,
175 &out_d0[0], 1);
176 Blas::Daxpy(nquad0 * nquad1 * nquad2, gmat[1][0], &diff1[0], 1,
177 &out_d0[0], 1);
178 Blas::Daxpy(nquad0 * nquad1 * nquad2, gmat[2][0], &diff2[0], 1,
179 &out_d0[0], 1);
180 }
181
182 if (out_d1.size())
183 {
184 Vmath::Smul(nquad0 * nquad1 * nquad2, gmat[3][0], &diff0[0], 1,
185 &out_d1[0], 1);
186 Blas::Daxpy(nquad0 * nquad1 * nquad2, gmat[4][0], &diff1[0], 1,
187 &out_d1[0], 1);
188 Blas::Daxpy(nquad0 * nquad1 * nquad2, gmat[5][0], &diff2[0], 1,
189 &out_d1[0], 1);
190 }
191
192 if (out_d2.size())
193 {
194 Vmath::Smul(nquad0 * nquad1 * nquad2, gmat[6][0], &diff0[0], 1,
195 &out_d2[0], 1);
196 Blas::Daxpy(nquad0 * nquad1 * nquad2, gmat[7][0], &diff1[0], 1,
197 &out_d2[0], 1);
198 Blas::Daxpy(nquad0 * nquad1 * nquad2, gmat[8][0], &diff2[0], 1,
199 &out_d2[0], 1);
200 }
201 }
202}
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().
◆ v_PhysEvaluate() [1/2]
|
overrideprotectedvirtual |
This function evaluates the expansion at a single (arbitrary) point of the domain.
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 #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::StdExpansion3D.
Definition at line 604 of file PyrExp.cpp.
606{
607 Array<OneD, NekDouble> Lcoord(3);
608
610
611 // TODO: check GetLocCoords()
612 m_geom->GetLocCoords(coord, Lcoord);
613
614 return StdExpansion3D::v_PhysEvaluate(Lcoord, physvals);
615}
References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.
◆ v_PhysEvaluate() [2/2]
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 617 of file PyrExp.cpp.
620{
621 Array<OneD, NekDouble> Lcoord(3);
623 m_geom->GetLocCoords(coord, Lcoord);
624 return StdPyrExp::v_PhysEvaluate(Lcoord, inarray, firstOrderDerivs);
625}
References ASSERTL0, and Nektar::LocalRegions::Expansion::m_geom.
◆ 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 596 of file PyrExp.cpp.
599{
600 // Evaluate point in local coordinates.
601 return StdExpansion3D::v_PhysEvaluate(Lcoord, physvals);
602}
◆ v_SVVLaplacianFilter()
|
overrideprotectedvirtual |
Reimplemented from Nektar::StdRegions::StdPyrExp.
Definition at line 1012 of file PyrExp.cpp.
1014{
1016
1017 // Calculate sqrt of the Jacobian
1019 Array<OneD, NekDouble> sqrt_jac(nq);
1021 {
1022 Vmath::Vsqrt(nq, jac, 1, sqrt_jac, 1);
1023 }
1024 else
1025 {
1027 }
1028
1029 // Multiply array by sqrt(Jac)
1030 Vmath::Vmul(nq, sqrt_jac, 1, array, 1, array, 1);
1031
1032 // Apply std region filter
1033 StdPyrExp::v_SVVLaplacianFilter(array, mkey);
1034
1035 // Divide by sqrt(Jac)
1036 Vmath::Vdiv(nq, array, 1, sqrt_jac, 1, array, 1);
1037}
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().
Member Data Documentation
◆ m_matrixManager
|
private |
Definition at line 174 of file PyrExp.h.
Referenced by v_DropLocMatrix(), v_FwdTrans(), and v_GetLocMatrix().
◆ m_staticCondMatrixManager
|
private |
Definition at line 176 of file PyrExp.h.
Referenced by v_DropLocStaticCondMatrix(), and v_GetLocStaticCondMatrix().
