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

Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions. More...

#include <ExpList3DHomogeneous1D.h>

Inheritance diagram for Nektar::MultiRegions::ExpList3DHomogeneous1D:
Inheritance graph
[legend]
Collaboration diagram for Nektar::MultiRegions::ExpList3DHomogeneous1D:
Collaboration graph
[legend]

Public Member Functions

 ExpList3DHomogeneous1D ()
 Default constructor. More...
 
 ExpList3DHomogeneous1D (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::BasisKey &HomoBasis, const NekDouble lhom, const bool useFFT, const bool dealiasing)
 
 ExpList3DHomogeneous1D (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::BasisKey &HomoBasis, const NekDouble lhom, const bool useFFT, const bool dealiasing, const SpatialDomains::MeshGraphSharedPtr &graph2D, const std::string &var="DefaultVar", const Collections::ImplementationType ImpType=Collections::eNoImpType)
 Sets up a list of local expansions based on an input mesh. More...
 
 ExpList3DHomogeneous1D (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::BasisKey &HomoBasis, const NekDouble lhom, const bool useFFT, const bool dealiasing, const SpatialDomains::ExpansionMap &expansions, const Collections::ImplementationType ImpType=Collections::eNoImpType)
 Sets up a list of local expansions based on an mesh expansion. More...
 
 ExpList3DHomogeneous1D (const ExpList3DHomogeneous1D &In, const bool DeclarePlanesSetCoeffPhys=true)
 Copy constructor. More...
 
 ExpList3DHomogeneous1D (const ExpList3DHomogeneous1D &In, const std::vector< unsigned int > &eIDs, const bool DeclarePlanesSetCoeffPhys=true, const Collections::ImplementationType ImpType=Collections::eNoImpType)
 Constructor copying only elements defined in eIds. More...
 
virtual ~ExpList3DHomogeneous1D ()
 Destructor. More...
 
void GetCoords (Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coord_2=NullNekDouble1DArray)
 This function calculates the coordinates of all the elemental quadrature points $\boldsymbol{x}_i$. More...
 
void GetCoords (const int eid, Array< OneD, NekDouble > &xc0, Array< OneD, NekDouble > &xc1, Array< OneD, NekDouble > &xc2)
 
- Public Member Functions inherited from Nektar::MultiRegions::ExpListHomogeneous1D
 ExpListHomogeneous1D ()
 Default constructor. More...
 
 ExpListHomogeneous1D (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::BasisKey &HomoBasis, const NekDouble lz, const bool useFFT, const bool dealiasing)
 
 ExpListHomogeneous1D (const ExpListHomogeneous1D &In)
 Copy constructor. More...
 
 ExpListHomogeneous1D (const ExpListHomogeneous1D &In, const std::vector< unsigned int > &eIDs)
 
virtual ~ExpListHomogeneous1D ()
 Destructor. More...
 
void Homogeneous1DTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool IsForwards, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
void HomogeneousFwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
void HomogeneousBwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
LibUtilities::BasisSharedPtr GetHomogeneousBasis (void)
 
void PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2)
 
void PhysDeriv (Direction edir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
ExpListSharedPtr & GetPlane (int n)
 
- Public Member Functions inherited from Nektar::MultiRegions::ExpList
 ExpList ()
 The default constructor. More...
 
 ExpList (const LibUtilities::SessionReaderSharedPtr &pSession)
 The default constructor. More...
 
 ExpList (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 The default constructor. More...
 
 ExpList (const ExpList &in, const std::vector< unsigned int > &eIDs, const bool DeclareCoeffPhysArrays=true)
 Constructor copying only elements defined in eIds. More...
 
 ExpList (const ExpList &in, const bool DeclareCoeffPhysArrays=true)
 The copy constructor. More...
 
virtual ~ExpList ()
 The default destructor. More...
 
int GetNcoeffs (void) const
 Returns the total number of local degrees of freedom $N_{\mathrm{eof}}=\sum_{e=1}^{{N_{\mathrm{el}}}}N^{e}_m$. More...
 
int GetNcoeffs (const int eid) const
 Returns the total number of local degrees of freedom for element eid. More...
 
ExpansionType GetExpType (void)
 Returns the type of the expansion. More...
 
void SetExpType (ExpansionType Type)
 Returns the type of the expansion. More...
 
int EvalBasisNumModesMax (void) const
 Evaulates the maximum number of modes in the elemental basis order over all elements. More...
 
const Array< OneD, int > EvalBasisNumModesMaxPerExp (void) const
 Returns the vector of the number of modes in the elemental basis order over all elements. More...
 
int GetTotPoints (void) const
 Returns the total number of quadrature points m_npoints $=Q_{\mathrm{tot}}$. More...
 
int GetTotPoints (const int eid) const
 Returns the total number of quadrature points for eid's element $=Q_{\mathrm{tot}}$. More...
 
int GetNpoints (void) const
 Returns the total number of quadrature points m_npoints $=Q_{\mathrm{tot}}$. More...
 
int Get1DScaledTotPoints (const NekDouble scale) const
 Returns the total number of qudature points scaled by the factor scale on each 1D direction. More...
 
void SetWaveSpace (const bool wavespace)
 Sets the wave space to the one of the possible configuration true or false. More...
 
void SetModifiedBasis (const bool modbasis)
 Set Modified Basis for the stability analysis. More...
 
void SetPhys (int i, NekDouble val)
 Set the i th value of m_phys to value val. More...
 
bool GetWaveSpace (void) const
 This function returns the third direction expansion condition, which can be in wave space (coefficient) or not It is stored in the variable m_WaveSpace. More...
 
void SetPhys (const Array< OneD, const NekDouble > &inarray)
 Fills the array m_phys. More...
 
void SetPhysArray (Array< OneD, NekDouble > &inarray)
 Sets the array m_phys. More...
 
void SetPhysState (const bool physState)
 This function manually sets whether the array of physical values $\boldsymbol{u}_l$ (implemented as m_phys) is filled or not. More...
 
bool GetPhysState (void) const
 This function indicates whether the array of physical values $\boldsymbol{u}_l$ (implemented as m_phys) is filled or not. More...
 
NekDouble PhysIntegral (void)
 This function integrates a function $f(\boldsymbol{x})$ over the domain consisting of all the elements of the expansion. More...
 
NekDouble PhysIntegral (const Array< OneD, const NekDouble > &inarray)
 This function integrates a function $f(\boldsymbol{x})$ over the domain consisting of all the elements of the expansion. More...
 
void IProductWRTBase_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function calculates the inner product of a function $f(\boldsymbol{x})$ with respect to all {local} expansion modes $\phi_n^e(\boldsymbol{x})$. More...
 
void IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function calculates the inner product of a function $f(\boldsymbol{x})$ with respect to the derivative (in direction. More...
 
void IProductWRTDerivBase (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, NekDouble > &outarray)
 This function calculates the inner product of a function $f(\boldsymbol{x})$ with respect to the derivative (in direction. More...
 
void FwdTrans_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function elementally evaluates the forward transformation of a function $u(\boldsymbol{x})$ onto the global spectral/hp expansion. More...
 
void FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void MultiplyByElmtInvMass (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function elementally mulplies the coefficient space of Sin my the elemental inverse of the mass matrix. More...
 
void MultiplyByInvMassMatrix (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void SmoothField (Array< OneD, NekDouble > &field)
 Smooth a field across elements. More...
 
void HelmSolve (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const FlagList &flags, const StdRegions::ConstFactorMap &factors, const StdRegions::VarCoeffMap &varcoeff=StdRegions::NullVarCoeffMap, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray, const bool PhysSpaceForcing=true)
 Solve helmholtz problem. More...
 
void LinearAdvectionDiffusionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 Solve Advection Diffusion Reaction. More...
 
void LinearAdvectionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 Solve Advection Diffusion Reaction. More...
 
void FwdTrans_BndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void BwdTrans_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function elementally evaluates the backward transformation of the global spectral/hp element expansion. More...
 
void BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void GetCoords (Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coord_2=NullNekDouble1DArray)
 This function calculates the coordinates of all the elemental quadrature points $\boldsymbol{x}_i$. More...
 
void HomogeneousFwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
void HomogeneousBwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
void DealiasedProd (const Array< OneD, NekDouble > &inarray1, const Array< OneD, NekDouble > &inarray2, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void DealiasedDotProd (const Array< OneD, Array< OneD, NekDouble > > &inarray1, const Array< OneD, Array< OneD, NekDouble > > &inarray2, Array< OneD, Array< OneD, NekDouble > > &outarray, CoeffState coeffstate=eLocal)
 
void GetBCValues (Array< OneD, NekDouble > &BndVals, const Array< OneD, NekDouble > &TotField, int BndID)
 
void NormVectorIProductWRTBase (Array< OneD, const NekDouble > &V1, Array< OneD, const NekDouble > &V2, Array< OneD, NekDouble > &outarray, int BndID)
 
void NormVectorIProductWRTBase (Array< OneD, Array< OneD, NekDouble > > &V, Array< OneD, NekDouble > &outarray)
 
void ApplyGeomInfo ()
 Apply geometry information to each expansion. More...
 
void Reset ()
 Reset geometry information and reset matrices. More...
 
void WriteTecplotHeader (std::ostream &outfile, std::string var="")
 
void WriteTecplotZone (std::ostream &outfile, int expansion=-1)
 
void WriteTecplotField (std::ostream &outfile, int expansion=-1)
 
void WriteTecplotConnectivity (std::ostream &outfile, int expansion=-1)
 
void WriteVtkHeader (std::ostream &outfile)
 
void WriteVtkFooter (std::ostream &outfile)
 
void WriteVtkPieceHeader (std::ostream &outfile, int expansion, int istrip=0)
 
void WriteVtkPieceFooter (std::ostream &outfile, int expansion)
 
void WriteVtkPieceData (std::ostream &outfile, int expansion, std::string var="v")
 
int GetCoordim (int eid)
 This function returns the dimension of the coordinates of the element eid. More...
 
void SetCoeff (int i, NekDouble val)
 Set the i th coefficiient in m_coeffs to value val. More...
 
void SetCoeffs (int i, NekDouble val)
 Set the i th coefficiient in m_coeffs to value val. More...
 
void SetCoeffsArray (Array< OneD, NekDouble > &inarray)
 Set the m_coeffs array to inarray. More...
 
const Array< OneD, const
NekDouble > & 		GetCoeffs () const
 This function returns (a reference to) the array $\boldsymbol{\hat{u}}_l$ (implemented as m_coeffs) containing all local expansion coefficients. More...
 
void ImposeDirichletConditions (Array< OneD, NekDouble > &outarray)
 Impose Dirichlet Boundary Conditions onto Array. More...
 
void FillBndCondFromField (void)
 Fill Bnd Condition expansion from the values stored in expansion. More...
 
void FillBndCondFromField (const int nreg)
 Fill Bnd Condition expansion in nreg from the values stored in expansion. More...
 
void LocalToGlobal (bool useComm=true)
 Gathers the global coefficients $\boldsymbol{\hat{u}}_g$ from the local coefficients $\boldsymbol{\hat{u}}_l$. More...
 
void LocalToGlobal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool useComm=true)
 
void GlobalToLocal (void)
 Scatters from the global coefficients $\boldsymbol{\hat{u}}_g$ to the local coefficients $\boldsymbol{\hat{u}}_l$. More...
 
void GlobalToLocal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
NekDouble GetCoeff (int i)
 Get the i th value (coefficient) of m_coeffs. More...
 
NekDouble GetCoeffs (int i)
 Get the i th value (coefficient) of m_coeffs. More...
 
const Array< OneD, const
NekDouble > & 		GetPhys () const
 This function returns (a reference to) the array $\boldsymbol{u}_l$ (implemented as m_phys) containing the function $u^{\delta}(\boldsymbol{x})$ evaluated at the quadrature points. More...
 
NekDouble Linf (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 This function calculates the $L_\infty$ error of the global spectral/hp element approximation. More...
 
NekDouble L2 (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 This function calculates the $L_2$ error with respect to soln of the global spectral/hp element approximation. More...
 
NekDouble H1 (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 Calculates the $H^1$ error of the global spectral/hp element approximation. More...
 
NekDouble Integral (const Array< OneD, const NekDouble > &inarray)
 
Array< OneD, const NekDouble > HomogeneousEnergy (void)
 This function calculates the energy associated with each one of the modesof a 3D homogeneous nD expansion. More...
 
void SetHomo1DSpecVanVisc (Array< OneD, NekDouble > visc)
 This function sets the Spectral Vanishing Viscosity in homogeneous1D expansion. More...
 
Array< OneD, const unsigned int > GetZIDs (void)
 This function returns a vector containing the wave numbers in z-direction associated with the 3D homogenous expansion. Required if a parellelisation is applied in the Fourier direction. More...
 
LibUtilities::TranspositionSharedPtr GetTransposition (void)
 This function returns the transposition class associaed with the homogeneous expansion. More...
 
NekDouble GetHomoLen (void)
 This function returns the Width of homogeneous direction associaed with the homogeneous expansion. More...
 
Array< OneD, const unsigned int > GetYIDs (void)
 This function returns a vector containing the wave numbers in y-direction associated with the 3D homogenous expansion. Required if a parellelisation is applied in the Fourier direction. More...
 
void PhysInterp1DScaled (const NekDouble scale, const Array< OneD, NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function interpolates the physical space points in inarray to outarray using the same points defined in the expansion but where the number of points are rescaled by 1DScale. More...
 
void PhysGalerkinProjection1DScaled (const NekDouble scale, const Array< OneD, NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function Galerkin projects the physical space points in inarray to outarray where inarray is assumed to be defined in the expansion but where the number of points are rescaled by 1DScale. More...
 
int GetExpSize (void)
 This function returns the number of elements in the expansion. More...
 
int GetNumElmts (void)
 This function returns the number of elements in the expansion which may be different for a homogeoenous extended expansionp. More...
 
const boost::shared_ptr
< LocalRegions::ExpansionVector > 		GetExp () const
 This function returns the vector of elements in the expansion. More...
 
LocalRegions::ExpansionSharedPtr & GetExp (int n) const
 This function returns (a shared pointer to) the local elemental expansion of the $n^{\mathrm{th}}$ element. More...
 
LocalRegions::ExpansionSharedPtr & GetExp (const Array< OneD, const NekDouble > &gloCoord)
 This function returns (a shared pointer to) the local elemental expansion containing the arbitrary point given by gloCoord. More...
 
int GetExpIndex (const Array< OneD, const NekDouble > &gloCoord, NekDouble tol=0.0, bool returnNearestElmt=false)
 
int GetExpIndex (const Array< OneD, const NekDouble > &gloCoords, Array< OneD, NekDouble > &locCoords, NekDouble tol=0.0, bool returnNearestElmt=false)
 
int GetCoeff_Offset (int n) const
 Get the start offset position for a global list of m_coeffs correspoinding to element n. More...
 
int GetPhys_Offset (int n) const
 Get the start offset position for a global list of m_phys correspoinding to element n. More...
 
Array< OneD, NekDouble > & UpdateCoeffs ()
 This function returns (a reference to) the array $\boldsymbol{\hat{u}}_l$ (implemented as m_coeffs) containing all local expansion coefficients. More...
 
Array< OneD, NekDouble > & UpdatePhys ()
 This function returns (a reference to) the array $\boldsymbol{u}_l$ (implemented as m_phys) containing the function $u^{\delta}(\boldsymbol{x})$ evaluated at the quadrature points. More...
 
void PhysDeriv (Direction edir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
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)
 This function discretely evaluates the derivative of a function $f(\boldsymbol{x})$ on the domain consisting of all elements of the expansion. More...
 
void PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
void CurlCurl (Array< OneD, Array< OneD, NekDouble > > &Vel, Array< OneD, Array< OneD, NekDouble > > &Q)
 
const Array< OneD, const
boost::shared_ptr< ExpList > > & 		GetBndCondExpansions ()
 
boost::shared_ptr< ExpList > & UpdateBndCondExpansion (int i)
 
void Upwind (const Array< OneD, const Array< OneD, NekDouble > > &Vec, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &Upwind)
 
void Upwind (const Array< OneD, const NekDouble > &Vn, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &Upwind)
 
boost::shared_ptr< ExpList > & GetTrace ()
 
boost::shared_ptr
< AssemblyMapDG > & 		GetTraceMap (void)
 
const Array< OneD, const int > & GetTraceBndMap (void)
 
void GetNormals (Array< OneD, Array< OneD, NekDouble > > &normals)
 
void AddTraceIntegral (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
 
void AddTraceIntegral (const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 
void AddFwdBwdTraceIntegral (const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &outarray)
 
void GetFwdBwdTracePhys (Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
 
void GetFwdBwdTracePhys (const Array< OneD, const NekDouble > &field, Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
 
const std::vector< bool > & GetLeftAdjacentFaces (void) const
 
void ExtractTracePhys (Array< OneD, NekDouble > &outarray)
 
void ExtractTracePhys (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
const Array< OneD, const
SpatialDomains::BoundaryConditionShPtr > & 		GetBndConditions ()
 
Array< OneD,
SpatialDomains::BoundaryConditionShPtr > & 		UpdateBndConditions ()
 
void EvaluateBoundaryConditions (const NekDouble time=0.0, const std::string varName="", const NekDouble=NekConstants::kNekUnsetDouble, const NekDouble=NekConstants::kNekUnsetDouble)
 
void GeneralMatrixOp (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 This function calculates the result of the multiplication of a matrix of type specified by mkey with a vector given by inarray. More...
 
void GeneralMatrixOp_IterPerExp (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void SetUpPhysNormals ()
 
void GetBoundaryToElmtMap (Array< OneD, int > &ElmtID, Array< OneD, int > &EdgeID)
 
void GetBndElmtExpansion (int i, boost::shared_ptr< ExpList > &result, const bool DeclareCoeffPhysArrays=true)
 
void ExtractElmtToBndPhys (int i, Array< OneD, NekDouble > &elmt, Array< OneD, NekDouble > &boundary)
 
void ExtractPhysToBndElmt (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bndElmt)
 
void ExtractPhysToBnd (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bnd)
 
void GetBoundaryNormals (int i, Array< OneD, Array< OneD, NekDouble > > &normals)
 
void GeneralGetFieldDefinitions (std::vector< LibUtilities::FieldDefinitionsSharedPtr > &fielddef, int NumHomoDir=0, Array< OneD, LibUtilities::BasisSharedPtr > &HomoBasis=LibUtilities::NullBasisSharedPtr1DArray, std::vector< NekDouble > &HomoLen=LibUtilities::NullNekDoubleVector, bool homoStrips=false, std::vector< unsigned int > &HomoSIDs=LibUtilities::NullUnsignedIntVector, std::vector< unsigned int > &HomoZIDs=LibUtilities::NullUnsignedIntVector, std::vector< unsigned int > &HomoYIDs=LibUtilities::NullUnsignedIntVector)
 
const
NekOptimize::GlobalOptParamSharedPtr & 		GetGlobalOptParam (void)
 
std::map< int,
RobinBCInfoSharedPtr > 		GetRobinBCInfo ()
 
void GetPeriodicEntities (PeriodicMap &periodicVerts, PeriodicMap &periodicEdges, PeriodicMap &periodicFaces=NullPeriodicMap)
 
std::vector
< LibUtilities::FieldDefinitionsSharedPtr > 		GetFieldDefinitions ()
 
void GetFieldDefinitions (std::vector< LibUtilities::FieldDefinitionsSharedPtr > &fielddef)
 
void AppendFieldData (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata)
 Append the element data listed in elements fielddef->m_ElementIDs onto fielddata. More...
 
void AppendFieldData (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, Array< OneD, NekDouble > &coeffs)
 Append the data in coeffs listed in elements fielddef->m_ElementIDs onto fielddata. More...
 
void ExtractElmtDataToCoeffs (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, std::string &field, Array< OneD, NekDouble > &coeffs)
 Extract the data in fielddata into the coeffs using the basic ExpList Elemental expansions rather than planes in homogeneous case. More...
 
void ExtractCoeffsToCoeffs (const boost::shared_ptr< ExpList > &fromExpList, const Array< OneD, const NekDouble > &fromCoeffs, Array< OneD, NekDouble > &toCoeffs)
 Extract the data from fromField using fromExpList the coeffs using the basic ExpList Elemental expansions rather than planes in homogeneous case. More...
 
void ExtractDataToCoeffs (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, std::string &field, Array< OneD, NekDouble > &coeffs)
 Extract the data in fielddata into the coeffs. More...
 
boost::shared_ptr< ExpList > GetSharedThisPtr ()
 Returns a shared pointer to the current object. More...
 
boost::shared_ptr
< LibUtilities::SessionReader > 		GetSession () const
 Returns the session object. More...
 
boost::shared_ptr
< LibUtilities::Comm > 		GetComm ()
 Returns the comm object. More...
 
SpatialDomains::MeshGraphSharedPtr GetGraph ()
 
LibUtilities::BasisSharedPtr GetHomogeneousBasis (void)
 
boost::shared_ptr< ExpList > & GetPlane (int n)
 
void CreateCollections (Collections::ImplementationType ImpType=Collections::eNoImpType)
 Construct collections of elements containing a single element type and polynomial order from the list of expansions. More...
 
void ClearGlobalLinSysManager (void)
 

Protected Member Functions

void SetCoeffPhys (void)
 Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m_coeff and m_phys. More...
 
virtual void v_GetCoords (Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1, Array< OneD, NekDouble > &coord_2)
 
virtual void v_WriteTecplotConnectivity (std::ostream &outfile, int expansion)
 
virtual void v_WriteVtkPieceHeader (std::ostream &outfile, int expansion, int istrip)
 
virtual NekDouble v_L2 (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 
virtual Array< OneD, const
NekDouble > 		v_HomogeneousEnergy (void)
 
virtual void v_GetPeriodicEntities (PeriodicMap &periodicVerts, PeriodicMap &periodicEdges, PeriodicMap &periodicFaces)
 
- Protected Member Functions inherited from Nektar::MultiRegions::ExpListHomogeneous1D
DNekBlkMatSharedPtr GenHomogeneous1DBlockMatrix (Homogeneous1DMatType mattype, CoeffState coeffstate=eLocal) const
 
DNekBlkMatSharedPtr GetHomogeneous1DBlockMatrix (Homogeneous1DMatType mattype, CoeffState coeffstate=eLocal) const
 
NekDouble GetSpecVanVisc (const int k)
 
virtual void v_SetHomo1DSpecVanVisc (Array< OneD, NekDouble > visc)
 
virtual int v_GetNumElmts (void)
 
virtual
LibUtilities::BasisSharedPtr 		v_GetHomogeneousBasis (void)
 
virtual void v_FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate)
 
virtual void v_FwdTrans_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate)
 
virtual void v_BwdTrans_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate)
 
virtual void v_IProductWRTBase_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual std::vector
< LibUtilities::FieldDefinitionsSharedPtr > 		v_GetFieldDefinitions (void)
 
virtual void v_GetFieldDefinitions (std::vector< LibUtilities::FieldDefinitionsSharedPtr > &fielddef)
 
virtual void v_AppendFieldData (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata)
 
virtual void v_AppendFieldData (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, Array< OneD, NekDouble > &coeffs)
 
virtual void v_ExtractDataToCoeffs (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, std::string &field, Array< OneD, NekDouble > &coeffs)
 Extract data from raw field data into expansion list. More...
 
virtual void v_ExtractCoeffsToCoeffs (const boost::shared_ptr< ExpList > &fromExpList, const Array< OneD, const NekDouble > &fromCoeffs, Array< OneD, NekDouble > &toCoeffs)
 
virtual void v_WriteVtkPieceData (std::ostream &outfile, int expansion, std::string var)
 
virtual void v_PhysInterp1DScaled (const NekDouble scale, const Array< OneD, NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_PhysGalerkinProjection1DScaled (const NekDouble scale, const Array< OneD, NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_HomogeneousFwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
virtual void v_HomogeneousBwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
virtual void v_DealiasedProd (const Array< OneD, NekDouble > &inarray1, const Array< OneD, NekDouble > &inarray2, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
virtual void v_DealiasedDotProd (const Array< OneD, Array< OneD, NekDouble > > &inarray1, const Array< OneD, Array< OneD, NekDouble > > &inarray2, Array< OneD, Array< OneD, NekDouble > > &outarray, CoeffState coeffstate=eLocal)
 
virtual void v_PhysDeriv (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2)
 
virtual void v_PhysDeriv (Direction edir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
virtual
LibUtilities::TranspositionSharedPtr 		v_GetTransposition (void)
 
virtual Array< OneD, const
unsigned int > 		v_GetZIDs (void)
 
virtual ExpListSharedPtr & v_GetPlane (int n)
 
virtual NekDouble v_GetHomoLen (void)
 
- Protected Member Functions inherited from Nektar::MultiRegions::ExpList
void SetCoeffPhysOffsets ()
 Definition of the total number of degrees of freedom and quadrature points and offsets to access data. More...
 
boost::shared_ptr< DNekMat > GenGlobalMatrixFull (const GlobalLinSysKey &mkey, const boost::shared_ptr< AssemblyMapCG > &locToGloMap)
 
const DNekScalBlkMatSharedPtr GenBlockMatrix (const GlobalMatrixKey &gkey)
 This function assembles the block diagonal matrix of local matrices of the type mtype. More...
 
const DNekScalBlkMatSharedPtr & GetBlockMatrix (const GlobalMatrixKey &gkey)
 
void MultiplyByBlockMatrix (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
boost::shared_ptr< GlobalMatrix > GenGlobalMatrix (const GlobalMatrixKey &mkey, const boost::shared_ptr< AssemblyMapCG > &locToGloMap)
 Generates a global matrix from the given key and map. More...
 
void GlobalEigenSystem (const boost::shared_ptr< DNekMat > &Gmat, Array< OneD, NekDouble > &EigValsReal, Array< OneD, NekDouble > &EigValsImag, Array< OneD, NekDouble > &EigVecs=NullNekDouble1DArray)
 
boost::shared_ptr< GlobalLinSys > GenGlobalLinSys (const GlobalLinSysKey &mkey, const boost::shared_ptr< AssemblyMapCG > &locToGloMap)
 This operation constructs the global linear system of type mkey. More...
 
boost::shared_ptr< GlobalLinSys > GenGlobalBndLinSys (const GlobalLinSysKey &mkey, const AssemblyMapSharedPtr &locToGloMap)
 Generate a GlobalLinSys from information provided by the key "mkey" and the mapping provided in LocToGloBaseMap. More...
 
void ReadGlobalOptimizationParameters ()
 
virtual const Array< OneD,
const boost::shared_ptr
< ExpList > > & 		v_GetBndCondExpansions (void)
 
virtual boost::shared_ptr
< ExpList > & 		v_UpdateBndCondExpansion (int i)
 
virtual void v_Upwind (const Array< OneD, const Array< OneD, NekDouble > > &Vec, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &Upwind)
 
virtual void v_Upwind (const Array< OneD, const NekDouble > &Vn, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &Upwind)
 
virtual boost::shared_ptr
< ExpList > & 		v_GetTrace ()
 
virtual boost::shared_ptr
< AssemblyMapDG > & 		v_GetTraceMap ()
 
virtual const Array< OneD,
const int > & 		v_GetTraceBndMap ()
 
virtual void v_GetNormals (Array< OneD, Array< OneD, NekDouble > > &normals)
 
virtual void v_AddTraceIntegral (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
 
virtual void v_AddTraceIntegral (const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 
virtual void v_AddFwdBwdTraceIntegral (const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &outarray)
 
virtual void v_GetFwdBwdTracePhys (Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
 
virtual void v_GetFwdBwdTracePhys (const Array< OneD, const NekDouble > &field, Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
 
virtual const std::vector< bool > & v_GetLeftAdjacentFaces (void) const
 
virtual void v_ExtractTracePhys (Array< OneD, NekDouble > &outarray)
 
virtual void v_ExtractTracePhys (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_MultiplyByInvMassMatrix (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate)
 
virtual void v_HelmSolve (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const FlagList &flags, const StdRegions::ConstFactorMap &factors, const StdRegions::VarCoeffMap &varcoeff, const Array< OneD, const NekDouble > &dirForcing, const bool PhysSpaceForcing)
 
virtual void v_LinearAdvectionDiffusionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 
virtual void v_LinearAdvectionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 
virtual void v_ImposeDirichletConditions (Array< OneD, NekDouble > &outarray)
 
virtual void v_FillBndCondFromField ()
 
virtual void v_FillBndCondFromField (const int nreg)
 
virtual void v_Reset ()
 Reset geometry information, metrics, matrix managers and geometry information. More...
 
virtual void v_LocalToGlobal (bool UseComm)
 
virtual void v_LocalToGlobal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool UseComm)
 
virtual void v_GlobalToLocal (void)
 
virtual void v_GlobalToLocal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_SmoothField (Array< OneD, NekDouble > &field)
 
virtual void v_GeneralMatrixOp (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate)
 
virtual void v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
virtual void v_CurlCurl (Array< OneD, Array< OneD, NekDouble > > &Vel, Array< OneD, Array< OneD, NekDouble > > &Q)
 
virtual void v_GetBCValues (Array< OneD, NekDouble > &BndVals, const Array< OneD, NekDouble > &TotField, int BndID)
 
virtual void v_NormVectorIProductWRTBase (Array< OneD, const NekDouble > &V1, Array< OneD, const NekDouble > &V2, Array< OneD, NekDouble > &outarray, int BndID)
 
virtual void v_NormVectorIProductWRTBase (Array< OneD, Array< OneD, NekDouble > > &V, Array< OneD, NekDouble > &outarray)
 
virtual void v_SetUpPhysNormals ()
 
virtual void v_GetBoundaryToElmtMap (Array< OneD, int > &ElmtID, Array< OneD, int > &EdgeID)
 
virtual void v_GetBndElmtExpansion (int i, boost::shared_ptr< ExpList > &result, const bool DeclareCoeffPhysArrays)
 
virtual void v_ExtractElmtToBndPhys (int i, Array< OneD, NekDouble > &elmt, Array< OneD, NekDouble > &boundary)
 
virtual void v_ExtractPhysToBndElmt (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bndElmt)
 
virtual void v_ExtractPhysToBnd (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bnd)
 
virtual void v_GetBoundaryNormals (int i, Array< OneD, Array< OneD, NekDouble > > &normals)
 
virtual void v_ReadGlobalOptimizationParameters ()
 
virtual void v_WriteTecplotHeader (std::ostream &outfile, std::string var="")
 
virtual void v_WriteTecplotZone (std::ostream &outfile, int expansion)
 
virtual void v_WriteTecplotField (std::ostream &outfile, int expansion)
 
virtual NekDouble v_Integral (const Array< OneD, const NekDouble > &inarray)
 
virtual Array< OneD, const
unsigned int > 		v_GetYIDs (void)
 
virtual void v_ClearGlobalLinSysManager (void)
 
void ExtractFileBCs (const std::string &fileName, LibUtilities::CommSharedPtr comm, const std::string &varName, const boost::shared_ptr< ExpList > locExpList)
 

Private Member Functions

void GenExpList3DHomogeneous1D (const SpatialDomains::ExpansionMap &expansions, const Collections::ImplementationType ImpType)
 

Additional Inherited Members

- Public Attributes inherited from Nektar::MultiRegions::ExpListHomogeneous1D
LibUtilities::TranspositionSharedPtr m_transposition
 
LibUtilities::CommSharedPtr m_StripZcomm
 
- Public Attributes inherited from Nektar::MultiRegions::ExpList
ExpansionType m_expType
 
- Static Protected Member Functions inherited from Nektar::MultiRegions::ExpList
static
SpatialDomains::BoundaryConditionShPtr 		GetBoundaryCondition (const SpatialDomains::BoundaryConditionCollection &collection, unsigned int index, const std::string &variable)
 
- Protected Attributes inherited from Nektar::MultiRegions::ExpListHomogeneous1D
bool m_useFFT
 FFT variables. More...
 
LibUtilities::NektarFFTSharedPtr m_FFT
 
LibUtilities::NektarFFTSharedPtr m_FFT_deal
 
Array< OneD, NekDouble > m_tmpIN
 
Array< OneD, NekDouble > m_tmpOUT
 
LibUtilities::BasisSharedPtr m_homogeneousBasis
 Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m_coeff and m_phys. More...
 
NekDouble m_lhom
 Width of homogeneous direction. More...
 
Homo1DBlockMatrixMapShPtr m_homogeneous1DBlockMat
 
Array< OneD, ExpListSharedPtr > m_planes
 
boost::unordered_map< int, int > m_zIdToPlane
 
- Protected Attributes inherited from Nektar::MultiRegions::ExpList
LibUtilities::CommSharedPtr m_comm
 Communicator. More...
 
LibUtilities::SessionReaderSharedPtr m_session
 Session. More...
 
SpatialDomains::MeshGraphSharedPtr m_graph
 Mesh associated with this expansion list. More...
 
int m_ncoeffs
 The total number of local degrees of freedom. m_ncoeffs $=N_{\mathrm{eof}}=\sum_{e=1}^{{N_{\mathrm{el}}}}N^{e}_l$. More...
 
int m_npoints
 
Array< OneD, NekDouble > m_coeffs
 Concatenation of all local expansion coefficients. More...
 
Array< OneD, NekDouble > m_phys
 The global expansion evaluated at the quadrature points. More...
 
bool m_physState
 The state of the array m_phys. More...
 
boost::shared_ptr
< LocalRegions::ExpansionVector > 		m_exp
 The list of local expansions. More...
 
Collections::CollectionVector m_collections
 
std::vector< int > m_coll_coeff_offset
 Offset of elemental data into the array m_coeffs. More...
 
std::vector< int > m_coll_phys_offset
 Offset of elemental data into the array m_phys. More...
 
Array< OneD, int > m_coeff_offset
 Offset of elemental data into the array m_coeffs. More...
 
Array< OneD, int > m_phys_offset
 Offset of elemental data into the array m_phys. More...
 
NekOptimize::GlobalOptParamSharedPtr m_globalOptParam
 
BlockMatrixMapShPtr m_blockMat
 
bool m_WaveSpace
 
boost::unordered_map< int, int > m_elmtToExpId
 Mapping from geometry ID of element to index inside m_exp. More...
 

Detailed Description

Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions.

Definition at line 61 of file ExpList3DHomogeneous1D.h.

Constructor & Destructor Documentation

Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D ( )

Default constructor.

Definition at line 47 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::e3DH1D, and Nektar::MultiRegions::ExpList::SetExpType().

47  :
48  ExpListHomogeneous1D()
49  {
50  SetExpType(e3DH1D);
51  }
Nektar::MultiRegions::ExpList::SetExpType
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
Definition: ExpList.cpp:277
Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D ( const LibUtilities::SessionReaderSharedPtr &  pSession,
const LibUtilities::BasisKey &  HomoBasis,
const NekDouble  lhom,
const bool  useFFT,
const bool  dealiasing 
)

Definition at line 53 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::e3DH1D, and Nektar::MultiRegions::ExpList::SetExpType().

58  :
59  ExpListHomogeneous1D(pSession,HomoBasis,lhom,useFFT,dealiasing)
60  {
61  SetExpType(e3DH1D);
62  }
Nektar::MultiRegions::ExpList::SetExpType
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
Definition: ExpList.cpp:277
Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D ( const LibUtilities::SessionReaderSharedPtr &  pSession,
const LibUtilities::BasisKey &  HomoBasis,
const NekDouble  lhom,
const bool  useFFT,
const bool  dealiasing,
const SpatialDomains::MeshGraphSharedPtr &  graph2D,
const std::string &  var = "DefaultVar",
const Collections::ImplementationType  ImpType = Collections::eNoImpType 
)

Sets up a list of local expansions based on an input mesh.

Definition at line 65 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::e3DH1D, GenExpList3DHomogeneous1D(), and Nektar::MultiRegions::ExpList::SetExpType().

73  :
74  ExpListHomogeneous1D(pSession,HomoBasis,lhom,useFFT,dealiasing)
75  {
76  SetExpType(e3DH1D);
77 
78  GenExpList3DHomogeneous1D(graph2D->GetExpansions(var),ImpType);
79  }
Nektar::MultiRegions::ExpList3DHomogeneous1D::GenExpList3DHomogeneous1D
void GenExpList3DHomogeneous1D(const SpatialDomains::ExpansionMap &expansions, const Collections::ImplementationType ImpType)
Definition: ExpList3DHomogeneous1D.cpp:97
Nektar::MultiRegions::ExpList::SetExpType
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
Definition: ExpList.cpp:277
Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D ( const LibUtilities::SessionReaderSharedPtr &  pSession,
const LibUtilities::BasisKey &  HomoBasis,
const NekDouble  lhom,
const bool  useFFT,
const bool  dealiasing,
const SpatialDomains::ExpansionMap &  expansions,
const Collections::ImplementationType  ImpType = Collections::eNoImpType 
)

Sets up a list of local expansions based on an mesh expansion.

Definition at line 82 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::e3DH1D, GenExpList3DHomogeneous1D(), and Nektar::MultiRegions::ExpList::SetExpType().

89  :
90  ExpListHomogeneous1D(pSession,HomoBasis,lhom,useFFT,dealiasing)
91  {
92  SetExpType(e3DH1D);
93 
94  GenExpList3DHomogeneous1D(expansions,ImpType);
95  }
Nektar::MultiRegions::ExpList3DHomogeneous1D::GenExpList3DHomogeneous1D
void GenExpList3DHomogeneous1D(const SpatialDomains::ExpansionMap &expansions, const Collections::ImplementationType ImpType)
Definition: ExpList3DHomogeneous1D.cpp:97
Nektar::MultiRegions::ExpList::SetExpType
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
Definition: ExpList.cpp:277
Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D ( const ExpList3DHomogeneous1D &  In,
const bool  DeclarePlanesSetCoeffPhys = true 
)

Copy constructor.

Parameters
InExpList3DHomogeneous1D object to copy.

Definition at line 137 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::e3DH1D, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, SetCoeffPhys(), and Nektar::MultiRegions::ExpList::SetExpType().

137  :
138  ExpListHomogeneous1D(In)
139  {
140  SetExpType(e3DH1D);
141 
142  if(DeclarePlanesSetCoeffPhys)
143  {
144  bool False = false;
145  ExpList2DSharedPtr zero_plane = boost::dynamic_pointer_cast<ExpList2D> (In.m_planes[0]);
146 
147  for(int n = 0; n < m_planes.num_elements(); ++n)
148  {
149  m_planes[n] = MemoryManager<ExpList2D>::AllocateSharedPtr(*zero_plane,False);
150  }
151 
152  SetCoeffPhys();
153  }
154  }
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::MultiRegions::ExpList3DHomogeneous1D::SetCoeffPhys
void SetCoeffPhys(void)
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpList3DHomogeneous1D.cpp:202
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Nektar::MultiRegions::ExpList2DSharedPtr
boost::shared_ptr< ExpList2D > ExpList2DSharedPtr
Shared pointer to an ExpList2D object.
Definition: ExpList2D.h:49
Nektar::MultiRegions::ExpList::SetExpType
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
Definition: ExpList.cpp:277
Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D ( const ExpList3DHomogeneous1D &  In,
const std::vector< unsigned int > &  eIDs,
const bool  DeclarePlanesSetCoeffPhys = true,
const Collections::ImplementationType  ImpType = Collections::eNoImpType 
)

Constructor copying only elements defined in eIds.

Parameters
InExpList3DHomogeneous1D object to copy.
eIDsId of elements that should be copied.

Definition at line 160 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::e3DH1D, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, SetCoeffPhys(), and Nektar::MultiRegions::ExpList::SetExpType().

164  :
165  ExpListHomogeneous1D(In, eIDs)
166  {
167  SetExpType(e3DH1D);
168 
169  if(DeclarePlanesSetCoeffPhys)
170  {
171  bool False = false;
172  std::vector<unsigned int> eIDsPlane;
173  int nel = eIDs.size()/m_planes.num_elements();
174 
175  for (int i = 0; i < nel; ++i)
176  {
177  eIDsPlane.push_back(eIDs[i]);
178  }
179 
180  ExpList2DSharedPtr zero_plane_old =
181  boost::dynamic_pointer_cast<ExpList2D> (In.m_planes[0]);
182 
183  ExpList2DSharedPtr zero_plane =
184  MemoryManager<ExpList2D>::AllocateSharedPtr(*(zero_plane_old), eIDsPlane, ImpType);
185 
186  for(int n = 0; n < m_planes.num_elements(); ++n)
187  {
188  m_planes[n] = MemoryManager<ExpList2D>::AllocateSharedPtr(*zero_plane,False);
189  }
190 
191  SetCoeffPhys();
192  }
193  }
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::MultiRegions::ExpList3DHomogeneous1D::SetCoeffPhys
void SetCoeffPhys(void)
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpList3DHomogeneous1D.cpp:202
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Nektar::MultiRegions::ExpList2DSharedPtr
boost::shared_ptr< ExpList2D > ExpList2DSharedPtr
Shared pointer to an ExpList2D object.
Definition: ExpList2D.h:49
Nektar::MultiRegions::ExpList::SetExpType
void SetExpType(ExpansionType Type)
Returns the type of the expansion.
Definition: ExpList.cpp:277
Nektar::MultiRegions::ExpList3DHomogeneous1D::~ExpList3DHomogeneous1D ( )
virtual

Destructor.

Destructor

Definition at line 198 of file ExpList3DHomogeneous1D.cpp.

199  {
200  }

Member Function Documentation

void Nektar::MultiRegions::ExpList3DHomogeneous1D::GenExpList3DHomogeneous1D ( const SpatialDomains::ExpansionMap &  expansions,
const Collections::ImplementationType  ImpType 
)
private

Definition at line 97 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::ExpList::GetExp(), Nektar::MultiRegions::ExpList::GetExpSize(), Nektar::MultiRegions::ExpList::m_exp, Nektar::MultiRegions::ExpList::m_globalOptParam, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, Nektar::MultiRegions::ExpList::m_session, and SetCoeffPhys().

Referenced by ExpList3DHomogeneous1D().

100  {
101  int n,j,nel;
102  bool False = false;
103  ExpList2DSharedPtr plane_zero;
104 
105  // note that nzplanes can be larger than nzmodes
106  m_planes[0] = plane_zero = MemoryManager<ExpList2D>::AllocateSharedPtr(m_session, expansions, False,ImpType);
107 
108  m_exp = MemoryManager<LocalRegions::ExpansionVector>::AllocateSharedPtr();
109  nel = m_planes[0]->GetExpSize();
110 
111  for(j = 0; j < nel; ++j)
112  {
113  (*m_exp).push_back(m_planes[0]->GetExp(j));
114  }
115 
116  for(n = 1; n < m_planes.num_elements(); ++n)
117  {
118  m_planes[n] = MemoryManager<ExpList2D>::AllocateSharedPtr(*plane_zero,False);
119  for(j = 0; j < nel; ++j)
120  {
121  (*m_exp).push_back((*m_exp)[j]);
122  }
123  }
124 
125  // Setup Default optimisation information.
126  nel = GetExpSize();
127  m_globalOptParam = MemoryManager<NekOptimize::GlobalOptParam>
128  ::AllocateSharedPtr(nel);
129 
130  SetCoeffPhys();
131  }
Nektar::MemoryManager::AllocateSharedPtr
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Definition: NekMemoryManager.hpp:235
Nektar::MultiRegions::ExpList::m_globalOptParam
NekOptimize::GlobalOptParamSharedPtr m_globalOptParam
Definition: ExpList.h:1052
Nektar::MultiRegions::ExpList::GetExp
const boost::shared_ptr< LocalRegions::ExpansionVector > GetExp() const
This function returns the vector of elements in the expansion.
Definition: ExpList.h:2067
Nektar::MultiRegions::ExpList::GetExpSize
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2046
Nektar::MultiRegions::ExpList3DHomogeneous1D::SetCoeffPhys
void SetCoeffPhys(void)
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpList3DHomogeneous1D.cpp:202
Nektar::MultiRegions::ExpList::m_exp
boost::shared_ptr< LocalRegions::ExpansionVector > m_exp
The list of local expansions.
Definition: ExpList.h:1036
Nektar::MultiRegions::ExpList::m_session
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:969
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Nektar::MultiRegions::ExpList2DSharedPtr
boost::shared_ptr< ExpList2D > ExpList2DSharedPtr
Shared pointer to an ExpList2D object.
Definition: ExpList2D.h:49
void Nektar::MultiRegions::ExpList3DHomogeneous1D::GetCoords ( Array< OneD, NekDouble > &  coord_0,
Array< OneD, NekDouble > &  coord_1 = NullNekDouble1DArray,
Array< OneD, NekDouble > &  coord_2 = NullNekDouble1DArray 
)
inline

This function calculates the coordinates of all the elemental quadrature points $\boldsymbol{x}_i$.

Definition at line 164 of file ExpList3DHomogeneous1D.h.

References v_GetCoords().

Referenced by v_WriteVtkPieceHeader().

168  {
169  v_GetCoords(coord_0,coord_1,coord_2);
170  }
Nektar::MultiRegions::ExpList3DHomogeneous1D::v_GetCoords
virtual void v_GetCoords(Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1, Array< OneD, NekDouble > &coord_2)
Definition: ExpList3DHomogeneous1D.cpp:288
void Nektar::MultiRegions::ExpList3DHomogeneous1D::GetCoords ( const int  eid,
Array< OneD, NekDouble > &  xc0,
Array< OneD, NekDouble > &  xc1,
Array< OneD, NekDouble > &  xc2 
)

Definition at line 235 of file ExpList3DHomogeneous1D.cpp.

References Vmath::Fill(), Nektar::MultiRegions::ExpList::GetTotPoints(), Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis, Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition, Vmath::Sadd(), Vmath::Smul(), and Vmath::Vcopy().

239  {
240  int n;
241  Array<OneD, NekDouble> tmp_xc;
242  int nzplanes = m_planes.num_elements();
243  int npoints = GetTotPoints(eid);
244 
245  (*m_exp)[eid]->GetCoords(xc0,xc1);
246 
247  // Fill z-direction
248  Array<OneD, const NekDouble> pts = m_homogeneousBasis->GetZ();
249  Array<OneD, NekDouble> local_pts(m_planes.num_elements());
250 
251  for(n = 0; n < m_planes.num_elements(); n++)
252  {
253  local_pts[n] = pts[m_transposition->GetPlaneID(n)];
254  }
255 
256  Array<OneD, NekDouble> z(nzplanes);
257 
258  Vmath::Smul(nzplanes,m_lhom/2.0,local_pts,1,z,1);
259  Vmath::Sadd(nzplanes,m_lhom/2.0,z,1,z,1);
260 
261  for(n = 0; n < nzplanes; ++n)
262  {
263  Vmath::Fill(npoints,z[n],tmp_xc = xc2 + npoints*n,1);
264  if(n)
265  {
266  Vmath::Vcopy(npoints,xc0,1,tmp_xc = xc0+npoints*n,1);
267  Vmath::Vcopy(npoints,xc1,1,tmp_xc = xc1+npoints*n,1);
268  }
269  }
270  }
Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition
LibUtilities::TranspositionSharedPtr m_transposition
Definition: ExpListHomogeneous1D.h:134
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:46
Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom
NekDouble m_lhom
Width of homogeneous direction.
Definition: ExpListHomogeneous1D.h:151
Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpListHomogeneous1D.h:150
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:213
Nektar::MultiRegions::ExpList::GetTotPoints
int GetTotPoints(void) const
Returns the total number of quadrature points m_npoints .
Definition: ExpList.h:1535
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Vmath::Sadd
void Sadd(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Add vector y = alpha + x.
Definition: Vmath.cpp:315
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1061
void Nektar::MultiRegions::ExpList3DHomogeneous1D::SetCoeffPhys ( void  )
protected

Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m_coeff and m_phys.

Definition at line 202 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::ExpList::m_coeff_offset, Nektar::MultiRegions::ExpList::m_coeffs, Nektar::MultiRegions::ExpList::m_ncoeffs, Nektar::MultiRegions::ExpList::m_npoints, Nektar::MultiRegions::ExpList::m_phys, Nektar::MultiRegions::ExpList::m_phys_offset, and Nektar::MultiRegions::ExpListHomogeneous1D::m_planes.

Referenced by Nektar::MultiRegions::ContField3DHomogeneous1D::ContField3DHomogeneous1D(), Nektar::MultiRegions::DisContField3DHomogeneous1D::DisContField3DHomogeneous1D(), ExpList3DHomogeneous1D(), and GenExpList3DHomogeneous1D().

203  {
204  int i,n,cnt;
205  int ncoeffs_per_plane = m_planes[0]->GetNcoeffs();
206  int npoints_per_plane = m_planes[0]->GetTotPoints();
207 
208  int nzplanes = m_planes.num_elements();
209 
210  // Set total coefficients and points
211  m_ncoeffs = ncoeffs_per_plane*nzplanes;
212  m_npoints = npoints_per_plane*nzplanes;
213 
214  m_coeffs = Array<OneD, NekDouble> (m_ncoeffs,0.0);
215  m_phys = Array<OneD, NekDouble> (m_npoints,0.0);
216 
217  int nel = m_planes[0]->GetExpSize();
218  m_coeff_offset = Array<OneD,int>(nel*nzplanes);
219  m_phys_offset = Array<OneD,int>(nel*nzplanes);
220  Array<OneD, NekDouble> tmparray;
221 
222  for(cnt = n = 0; n < nzplanes; ++n)
223  {
224  m_planes[n]->SetCoeffsArray(tmparray= m_coeffs + ncoeffs_per_plane*n);
225  m_planes[n]->SetPhysArray(tmparray = m_phys + npoints_per_plane*n);
226 
227  for(i = 0; i < nel; ++i)
228  {
229  m_coeff_offset[cnt] = m_planes[n]->GetCoeff_Offset(i) + n*ncoeffs_per_plane;
230  m_phys_offset[cnt++] = m_planes[n]->GetPhys_Offset(i) + n*npoints_per_plane;
231  }
232  }
233  }
Nektar::MultiRegions::ExpList::m_phys
Array< OneD, NekDouble > m_phys
The global expansion evaluated at the quadrature points.
Definition: ExpList.h:1015
Nektar::MultiRegions::ExpList::m_coeffs
Array< OneD, NekDouble > m_coeffs
Concatenation of all local expansion coefficients.
Definition: ExpList.h:998
Nektar::MultiRegions::ExpList::m_coeff_offset
Array< OneD, int > m_coeff_offset
Offset of elemental data into the array m_coeffs.
Definition: ExpList.h:1047
Nektar::MultiRegions::ExpList::m_phys_offset
Array< OneD, int > m_phys_offset
Offset of elemental data into the array m_phys.
Definition: ExpList.h:1050
Nektar::MultiRegions::ExpList::m_ncoeffs
int m_ncoeffs
The total number of local degrees of freedom. m_ncoeffs .
Definition: ExpList.h:976
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
void Nektar::MultiRegions::ExpList3DHomogeneous1D::v_GetCoords ( Array< OneD, NekDouble > &  xc0,
Array< OneD, NekDouble > &  xc1,
Array< OneD, NekDouble > &  xc2 
)
protectedvirtual

The operation calls the 2D plane coordinates through the function ExpList::GetCoords and then evaluated the third coordinate using the member m_lhom

Parameters
coord_0After calculation, the $x_1$ coordinate will be stored in this array.
coord_1After calculation, the $x_2$ coordinate will be stored in this array.
coord_2After calculation, the $x_3$ coordinate will be stored in this array. This coordinate is evaluated using the predefined value m_lhom

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 288 of file ExpList3DHomogeneous1D.cpp.

References Vmath::Fill(), Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis, Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition, Vmath::Sadd(), Vmath::Smul(), and Vmath::Vcopy().

Referenced by GetCoords().

291  {
292  int n;
293  Array<OneD, NekDouble> tmp_xc;
294  int nzplanes = m_planes.num_elements();
295  int npoints = m_planes[0]->GetTotPoints();
296 
297  m_planes[0]->GetCoords(xc0,xc1);
298 
299  // Fill z-direction
300  Array<OneD, const NekDouble> pts = m_homogeneousBasis->GetZ();
301 
302  Array<OneD, NekDouble> local_pts(m_planes.num_elements());
303 
304  for(n = 0; n < m_planes.num_elements(); n++)
305  {
306  local_pts[n] = pts[m_transposition->GetPlaneID(n)];
307  }
308 
309  Array<OneD, NekDouble> z(nzplanes);
310 
311  Vmath::Smul(nzplanes,m_lhom/2.0,local_pts,1,z,1);
312  Vmath::Sadd(nzplanes,m_lhom/2.0,z,1,z,1);
313 
314  for(n = 0; n < nzplanes; ++n)
315  {
316  Vmath::Fill(npoints,z[n],tmp_xc = xc2 + npoints*n,1);
317  if(n)
318  {
319  Vmath::Vcopy(npoints,xc0,1,tmp_xc = xc0+npoints*n,1);
320  Vmath::Vcopy(npoints,xc1,1,tmp_xc = xc1+npoints*n,1);
321  }
322  }
323  }
Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition
LibUtilities::TranspositionSharedPtr m_transposition
Definition: ExpListHomogeneous1D.h:134
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:46
Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom
NekDouble m_lhom
Width of homogeneous direction.
Definition: ExpListHomogeneous1D.h:151
Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpListHomogeneous1D.h:150
Vmath::Smul
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:213
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Vmath::Sadd
void Sadd(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Add vector y = alpha + x.
Definition: Vmath.cpp:315
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1061
virtual void Nektar::MultiRegions::ExpList3DHomogeneous1D::v_GetPeriodicEntities ( PeriodicMap &  periodicVerts,
PeriodicMap &  periodicEdges,
PeriodicMap &  periodicFaces 
)
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 149 of file ExpList3DHomogeneous1D.h.

References Nektar::MultiRegions::ExpListHomogeneous1D::m_planes.

154  {
155  m_planes[0]->GetPeriodicEntities(periodicVerts,periodicEdges);
156  }
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Array< OneD, const NekDouble > Nektar::MultiRegions::ExpList3DHomogeneous1D::v_HomogeneousEnergy ( void  )
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 519 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::ExpList::GetCoeff_Offset(), Nektar::MultiRegions::ExpList::GetCoeffs(), Nektar::MultiRegions::ExpList::GetExp(), Nektar::MultiRegions::ExpList::GetTotPoints(), Nektar::MultiRegions::ExpList::m_comm, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, and Nektar::LibUtilities::ReduceSum.

520  {
521  Array<OneD, NekDouble> energy(m_planes.num_elements()/2);
522  NekDouble area = 0.0;
523 
524  // Calculate total area of elements.
525  for (int n = 0; n < m_planes[0]->GetExpSize(); ++n)
526  {
527  Array<OneD, NekDouble> inarray(m_planes[0]->GetExp(n)->GetTotPoints(), 1.0);
528  area += m_planes[0]->GetExp(n)->Integral(inarray);
529  }
530 
531  m_comm->GetRowComm()->AllReduce(area, LibUtilities::ReduceSum);
532 
533  // Calculate L2 norm of real/imaginary planes.
534  for (int n = 0; n < m_planes.num_elements(); n += 2)
535  {
536  NekDouble err;
537 
538  energy[n/2] = 0;
539 
540  for(int i = 0; i < m_planes[n]->GetExpSize(); ++i)
541  {
542  StdRegions::StdExpansionSharedPtr exp = m_planes[n]->GetExp(i);
543  Array<OneD, NekDouble> phys(exp->GetTotPoints());
544  exp->BwdTrans(m_planes[n]->GetCoeffs()+m_planes[n]->GetCoeff_Offset(i),
545  phys);
546  err = exp->L2(phys);
547  energy[n/2] += err*err;
548 
549  exp = m_planes[n+1]->GetExp(i);
550  exp->BwdTrans(m_planes[n+1]->GetCoeffs()+m_planes[n+1]->GetCoeff_Offset(i),
551  phys);
552  err = exp->L2(phys);
553  energy[n/2] += err*err;
554  }
555 
556  m_comm->GetRowComm()->AllReduce(energy[n/2], LibUtilities::ReduceSum);
557  energy[n/2] /= 2.0*area;
558  }
559 
560  return energy;
561  }
Nektar::MultiRegions::ExpList::GetCoeffs
const Array< OneD, const NekDouble > & GetCoeffs() const
This function returns (a reference to) the array (implemented as m_coeffs) containing all local expa...
Definition: ExpList.h:1938
Nektar::MultiRegions::ExpList::GetCoeff_Offset
int GetCoeff_Offset(int n) const
Get the start offset position for a global list of m_coeffs correspoinding to element n...
Definition: ExpList.h:2076
Nektar::MultiRegions::ExpList::GetExp
const boost::shared_ptr< LocalRegions::ExpansionVector > GetExp() const
This function returns the vector of elements in the expansion.
Definition: ExpList.h:2067
Nektar::MultiRegions::ExpList::GetTotPoints
int GetTotPoints(void) const
Returns the total number of quadrature points m_npoints .
Definition: ExpList.h:1535
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::MultiRegions::ExpList::m_comm
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:966
Nektar::StdRegions::StdExpansionSharedPtr
boost::shared_ptr< StdExpansion > StdExpansionSharedPtr
Definition: StdExpansion.h:1878
NekDouble Nektar::MultiRegions::ExpList3DHomogeneous1D::v_L2 ( const Array< OneD, const NekDouble > &  inarray,
const Array< OneD, const NekDouble > &  soln = NullNekDouble1DArray 
)
protectedvirtual

Given a spectral/hp approximation $u^{\delta}(\boldsymbol{x})$ evaluated at the quadrature points (which should be contained in m_phys), this function calculates the $L_2$ error of this approximation with respect to an exact solution. The local distribution of the quadrature points allows an elemental evaluation of this operation through the functions StdRegions::StdExpansion::L2.

The exact solution, also evaluated at the quadrature points, should be contained in the variable m_phys of the ExpList object Sol.

Parameters
SolAn ExpList, containing the discrete evaluation of the exact solution at the quadrature points in its array m_phys.
Returns
The $L_2$ error of the approximation.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 482 of file ExpList3DHomogeneous1D.cpp.

References Nektar::MultiRegions::ExpList::m_comm, Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis, Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition, Nektar::NullNekDouble1DArray, and Nektar::LibUtilities::ReduceSum.

485  {
486  int cnt = 0;
487  NekDouble errL2, err = 0.0;
488  Array<OneD, const NekDouble> w = m_homogeneousBasis->GetW();
489  Array<OneD, NekDouble> local_w(m_planes.num_elements());
490 
491  for(int n = 0; n < m_planes.num_elements(); n++)
492  {
493  local_w[n] = w[m_transposition->GetPlaneID(n)];
494  }
495 
496  if (soln == NullNekDouble1DArray)
497  {
498  for(int n = 0; n < m_planes.num_elements(); ++n)
499  {
500  errL2 = m_planes[n]->L2(inarray + cnt);
501  cnt += m_planes[n]->GetTotPoints();
502  err += errL2*errL2*local_w[n]*m_lhom*0.5;
503  }
504  }
505  else
506  {
507  for(int n = 0; n < m_planes.num_elements(); ++n)
508  {
509  errL2 = m_planes[n]->L2(inarray + cnt, soln + cnt);
510  cnt += m_planes[n]->GetTotPoints();
511  err += errL2*errL2*local_w[n]*m_lhom*0.5;
512  }
513  }
514  m_comm->GetColumnComm()->AllReduce(err, LibUtilities::ReduceSum);
515 
516  return sqrt(err);
517  }
Nektar::NullNekDouble1DArray
static Array< OneD, NekDouble > NullNekDouble1DArray
Definition: SharedArray.hpp:785
Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition
LibUtilities::TranspositionSharedPtr m_transposition
Definition: ExpListHomogeneous1D.h:134
Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom
NekDouble m_lhom
Width of homogeneous direction.
Definition: ExpListHomogeneous1D.h:151
Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpListHomogeneous1D.h:150
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::MultiRegions::ExpList::m_comm
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:966
void Nektar::MultiRegions::ExpList3DHomogeneous1D::v_WriteTecplotConnectivity ( std::ostream &  outfile,
int  expansion 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 325 of file ExpList3DHomogeneous1D.cpp.

References ASSERTL0, and Nektar::MultiRegions::ExpListHomogeneous1D::m_planes.

326  {
327  ASSERTL0(expansion == -1, "Multi-zone output not supported for homogeneous expansions.");
328 
329  const int nPtsPlane = m_planes[0]->GetNpoints();
330  const int nElmt = m_planes[0]->GetExpSize();
331  const int nPlanes = m_planes.num_elements();
332 
333  int cnt = 0;
334  int cnt2 = 0;
335  for (int i = 0; i < nElmt; ++i)
336  {
337  const int np0 = (*m_exp)[i]->GetNumPoints(0);
338  const int np1 = (*m_exp)[i]->GetNumPoints(1);
339 
340  for (int n = 1; n < nPlanes; ++n)
341  {
342  const int o1 = (n-1) * nPtsPlane;
343  const int o2 = n * nPtsPlane;
344  for (int j = 1; j < np1; ++j)
345  {
346  for(int k = 1; k < np0; ++k)
347  {
348  outfile << cnt + (j-1)*np0 + (k-1) + o1 + 1 << " ";
349  outfile << cnt + (j-1)*np0 + (k-1) + o2 + 1 << " ";
350  outfile << cnt + (j-1)*np0 + k + o2 + 1 << " ";
351  outfile << cnt + (j-1)*np0 + k + o1 + 1 << " ";
352  outfile << cnt + j *np0 + (k-1) + o1 + 1 << " ";
353  outfile << cnt + j *np0 + (k-1) + o2 + 1 << " ";
354  outfile << cnt + j *np0 + k + o2 + 1 << " ";
355  outfile << cnt + j *np0 + k + o1 + 1 << endl;
356  cnt2++;
357  }
358  }
359  }
360 
361  cnt += np0*np1;
362  }
363  }
ASSERTL0
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
void Nektar::MultiRegions::ExpList3DHomogeneous1D::v_WriteVtkPieceHeader ( std::ostream &  outfile,
int  expansion,
int  istrip 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 366 of file ExpList3DHomogeneous1D.cpp.

References Nektar::LibUtilities::eFourier, Nektar::LibUtilities::eFourierEvenlySpaced, Vmath::Fill(), GetCoords(), Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis, Nektar::MultiRegions::ExpListHomogeneous1D::m_planes, Nektar::MultiRegions::ExpList::m_session, and Vmath::Vcopy().

370  {
371  // If there is only one plane (e.g. HalfMode), we write a 2D plane.
372  if (m_planes.num_elements() == 1)
373  {
374  m_planes[0]->WriteVtkPieceHeader(outfile, expansion);
375  return;
376  }
377 
378  // If we are using Fourier points, output extra plane to fill domain
379  int outputExtraPlane = 0;
380  if ( m_homogeneousBasis->GetBasisType() == LibUtilities::eFourier
381  && m_homogeneousBasis->GetPointsType() ==
382  LibUtilities::eFourierEvenlySpaced)
383  {
384  outputExtraPlane = 1;
385  }
386 
387  int i,j,k;
388  int nq0 = (*m_exp)[expansion]->GetNumPoints(0);
389  int nq1 = (*m_exp)[expansion]->GetNumPoints(1);
390  int nq2 = m_planes.num_elements() + outputExtraPlane;
391  int ntot = nq0*nq1*nq2;
392  int ntotminus = (nq0-1)*(nq1-1)*(nq2-1);
393 
394  Array<OneD,NekDouble> coords[3];
395  coords[0] = Array<OneD,NekDouble>(ntot);
396  coords[1] = Array<OneD,NekDouble>(ntot);
397  coords[2] = Array<OneD,NekDouble>(ntot);
398  GetCoords(expansion,coords[0],coords[1],coords[2]);
399 
400  if (outputExtraPlane)
401  {
402  // Copy coords[0] and coords[1] to extra plane
403  Array<OneD,NekDouble> tmp;
404  Vmath::Vcopy (nq0*nq1, coords[0], 1,
405  tmp = coords[0] + (nq2-1)*nq0*nq1, 1);
406  Vmath::Vcopy (nq0*nq1, coords[1], 1,
407  tmp = coords[1] + (nq2-1)*nq0*nq1, 1);
408  // Fill coords[2] for extra plane
409  NekDouble z = coords[2][nq0*nq1*m_planes.num_elements()-1] +
410  (coords[2][nq0*nq1] - coords[2][0]);
411  Vmath::Fill(nq0*nq1, z, tmp = coords[2] + (nq2-1)*nq0*nq1, 1);
412  }
413 
414  NekDouble DistStrip;
415  m_session->LoadParameter("DistStrip", DistStrip, 0);
416  // Reset the z-coords for homostrips
417  for(int i = 0; i < ntot; i++)
418  {
419  coords[2][i] += istrip*DistStrip;
420  }
421 
422  outfile << " <Piece NumberOfPoints=\""
423  << ntot << "\" NumberOfCells=\""
424  << ntotminus << "\">" << endl;
425  outfile << " <Points>" << endl;
426  outfile << " <DataArray type=\"Float64\" "
427  << "NumberOfComponents=\"3\" format=\"ascii\">" << endl;
428  outfile << " ";
429  for (i = 0; i < ntot; ++i)
430  {
431  for (j = 0; j < 3; ++j)
432  {
433  outfile << coords[j][i] << " ";
434  }
435  outfile << endl;
436  }
437  outfile << endl;
438  outfile << " </DataArray>" << endl;
439  outfile << " </Points>" << endl;
440  outfile << " <Cells>" << endl;
441  outfile << " <DataArray type=\"Int32\" "
442  << "Name=\"connectivity\" format=\"ascii\">" << endl;
443  for (i = 0; i < nq0-1; ++i)
444  {
445  for (j = 0; j < nq1-1; ++j)
446  {
447  for (k = 0; k < nq2-1; ++k)
448  {
449  outfile << k*nq0*nq1 + j*nq0 + i << " "
450  << k*nq0*nq1 + j*nq0 + i + 1 << " "
451  << k*nq0*nq1 + (j+1)*nq0 + i + 1 << " "
452  << k*nq0*nq1 + (j+1)*nq0 + i << " "
453  << (k+1)*nq0*nq1 + j*nq0 + i << " "
454  << (k+1)*nq0*nq1 + j*nq0 + i + 1 << " "
455  << (k+1)*nq0*nq1 + (j+1)*nq0 + i + 1 << " "
456  << (k+1)*nq0*nq1 + (j+1)*nq0 + i << endl;
457  }
458  }
459  }
460  outfile << endl;
461  outfile << " </DataArray>" << endl;
462  outfile << " <DataArray type=\"Int32\" "
463  << "Name=\"offsets\" format=\"ascii\">" << endl;
464  for (i = 0; i < ntotminus; ++i)
465  {
466  outfile << i*8+8 << " ";
467  }
468  outfile << endl;
469  outfile << " </DataArray>" << endl;
470  outfile << " <DataArray type=\"UInt8\" "
471  << "Name=\"types\" format=\"ascii\">" << endl;
472  for (i = 0; i < ntotminus; ++i)
473  {
474  outfile << "12 ";
475  }
476  outfile << endl;
477  outfile << " </DataArray>" << endl;
478  outfile << " </Cells>" << endl;
479  outfile << " <PointData>" << endl;
480  }
Vmath::Fill
void Fill(int n, const T alpha, T *x, const int incx)
Fill a vector with a constant value.
Definition: Vmath.cpp:46
Nektar::LibUtilities::eFourier
Fourier Expansion .
Definition: BasisType.h:52
Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Definition: ExpListHomogeneous1D.h:150
Nektar::LibUtilities::eFourierEvenlySpaced
1D Evenly-spaced points using Fourier Fit
Definition: PointsType.h:66
Nektar::MultiRegions::ExpList::m_session
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:969
Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
Array< OneD, ExpListSharedPtr > m_planes
Definition: ExpListHomogeneous1D.h:153
Nektar::NekDouble
double NekDouble
Definition: NektarUnivTypeDefs.hpp:44
Nektar::MultiRegions::ExpList3DHomogeneous1D::GetCoords
void GetCoords(Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coord_2=NullNekDouble1DArray)
This function calculates the coordinates of all the elemental quadrature points . ...
Definition: ExpList3DHomogeneous1D.h:164
Vmath::Vcopy
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1061
Generated on Sat Nov 25 2017 16:19:01 for Nektar++ by   doxygen 1.8.8