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

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

#include <ExpListHomogeneous1D.h>

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

Public Member Functions

 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)
 
void DealiasedProd (const Array< OneD, NekDouble > &inarray1, const Array< OneD, NekDouble > &inarray2, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
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)
 
ExpListSharedPtrGetPlane (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)
 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 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 LocalToGlobal (void)
 Put the coefficients into global ordering using m_coeffs. More...
 
void GlobalToLocal (void)
 Put the coefficients into local ordering and place in m_coeffs. More...
 
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 NekDoubleHomogeneousEnergy (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::ExpansionSharedPtrGetExp (int n) const
 This function returns (a shared pointer to) the local elemental expansion of the $n^{\mathrm{th}}$ element. More...
 
LocalRegions::ExpansionSharedPtrGetExp (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...
 
int GetOffset_Elmt_Id (int n) const
 Get the element id associated with the n th consecutive block of data in m_phys and m_coeffs. 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)
 
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 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)
 
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 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)
 
map< int, RobinBCInfoSharedPtrGetRobinBCInfo ()
 
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< ExpListGetSharedThisPtr ()
 Returns a shared pointer to the current object. More...
 
boost::shared_ptr
< LibUtilities::SessionReader
GetSession ()
 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)
 

Public Attributes

LibUtilities::TranspositionSharedPtr m_transposition
 
LibUtilities::CommSharedPtr m_StripZcomm
 
- Public Attributes inherited from Nektar::MultiRegions::ExpList
ExpansionType m_expType
 

Protected Member Functions

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_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 ExpListSharedPtrv_GetPlane (int n)
 
virtual NekDouble v_GetHomoLen (void)
 
- Protected Member Functions inherited from Nektar::MultiRegions::ExpList
boost::shared_ptr< DNekMatGenGlobalMatrixFull (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 DNekScalBlkMatSharedPtrGetBlockMatrix (const GlobalMatrixKey &gkey)
 
void MultiplyByBlockMatrix (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
boost::shared_ptr< GlobalMatrixGenGlobalMatrix (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< GlobalLinSysGenGlobalLinSys (const GlobalLinSysKey &mkey, const boost::shared_ptr< AssemblyMapCG > &locToGloMap)
 This operation constructs the global linear system of type mkey. More...
 
boost::shared_ptr< GlobalLinSysGenGlobalBndLinSys (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 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)
 
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_Reset ()
 Reset geometry information, metrics, matrix managers and geometry information. More...
 
virtual void v_LocalToGlobal (void)
 
virtual void v_GlobalToLocal (void)
 
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_GetCoords (Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1, Array< OneD, NekDouble > &coord_2=NullNekDouble1DArray)
 
virtual void v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
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)
 
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_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 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 > &phys, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 
virtual NekDouble v_Integral (const Array< OneD, const NekDouble > &inarray)
 
virtual Array< OneD, const
NekDouble
v_HomogeneousEnergy (void)
 
virtual Array< OneD, const
unsigned int > 
v_GetYIDs (void)
 
virtual void v_ClearGlobalLinSysManager (void)
 
void ExtractFileBCs (const std::string &fileName, const std::string &varName, const boost::shared_ptr< ExpList > locExpList)
 

Protected Attributes

bool m_useFFT
 FFT variables. More...
 
LibUtilities::NektarFFTSharedPtr m_FFT
 
LibUtilities::NektarFFTSharedPtr m_FFT_deal
 
Array< OneD, NekDoublem_tmpIN
 
Array< OneD, NekDoublem_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, ExpListSharedPtrm_planes
 
- 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, NekDoublem_coeffs
 Concatenation of all local expansion coefficients. More...
 
Array< OneD, NekDoublem_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...
 
Array< OneD, int > m_offset_elmt_id
 Array containing the element id m_offset_elmt_id[n] that the n^th consecutive block of data in m_coeffs and m_phys is associated, i.e. for an array of constant expansion size and single shape elements m_phys[n*m_npoints] is the data related to m_exp[m_offset_elmt_id[n]];. More...
 
NekOptimize::GlobalOptParamSharedPtr m_globalOptParam
 
BlockMatrixMapShPtr m_blockMat
 
bool m_WaveSpace
 

Private Attributes

bool m_dealiasing
 
int m_padsize
 
Array< OneD, NekDoublem_specVanVisc
 Spectral vanishing Viscosity coefficient for stabilisation. More...
 

Additional Inherited Members

- Static Protected Member Functions inherited from Nektar::MultiRegions::ExpList
static
SpatialDomains::BoundaryConditionShPtr 
GetBoundaryCondition (const SpatialDomains::BoundaryConditionCollection &collection, unsigned int index, const std::string &variable)
 

Detailed Description

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

Definition at line 77 of file ExpListHomogeneous1D.h.

Constructor & Destructor Documentation

Nektar::MultiRegions::ExpListHomogeneous1D::ExpListHomogeneous1D ( )

Default constructor.

Definition at line 48 of file ExpListHomogeneous1D.cpp.

48  :
49  ExpList(),
51  m_lhom(1),
53  {
54  }
ExpList()
The default constructor.
Definition: ExpList.cpp:93
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
static BasisSharedPtr NullBasisSharedPtr
Definition: Basis.h:358
NekDouble m_lhom
Width of homogeneous direction.
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Nektar::MultiRegions::ExpListHomogeneous1D::ExpListHomogeneous1D ( const LibUtilities::SessionReaderSharedPtr pSession,
const LibUtilities::BasisKey HomoBasis,
const NekDouble  lz,
const bool  useFFT,
const bool  dealiasing 
)

Definition at line 56 of file ExpListHomogeneous1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, ASSERTL2, Nektar::LibUtilities::BasisManager(), Nektar::LibUtilities::NekFactory< tKey, tBase, >::CreateInstance(), Nektar::LibUtilities::GetNektarFFTFactory(), Nektar::MultiRegions::ExpList::m_comm, m_dealiasing, m_FFT, m_FFT_deal, m_homogeneousBasis, m_padsize, m_planes, Nektar::MultiRegions::ExpList::m_session, m_StripZcomm, m_transposition, m_useFFT, and Nektar::LibUtilities::NullBasisKey().

57  :
58  ExpList(pSession),
59  m_useFFT(useFFT),
60  m_lhom(lhom),
62  m_dealiasing(dealiasing)
63  {
64  ASSERTL2(HomoBasis != LibUtilities::NullBasisKey,"Homogeneous Basis is a null basis");
65 
67 
68  m_StripZcomm = m_session->DefinesSolverInfo("HomoStrip") ?
69  m_comm->GetColumnComm()->GetColumnComm() :
70  m_comm->GetColumnComm();
71 
74 
76  m_homogeneousBasis->GetNumPoints() /
77  m_StripZcomm->GetSize());
78 
79  if(m_useFFT)
80  {
82  "NekFFTW", m_homogeneousBasis->GetNumPoints());
83  }
84 
85  if(m_dealiasing)
86  {
87  if(m_useFFT)
88  {
89  NekDouble size = 1.5*m_homogeneousBasis->GetNumPoints();
90  m_padsize = int(size);
92  .CreateInstance("NekFFTW", m_padsize);
93  }
94  else
95  {
96  ASSERTL0(false, "Dealiasing available just in combination "
97  "with FFTW");
98  }
99  }
100  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
ExpList()
The default constructor.
Definition: ExpList.cpp:93
tBaseSharedPtr CreateInstance(tKey idKey BOOST_PP_COMMA_IF(MAX_PARAM) BOOST_PP_ENUM_BINARY_PARAMS(MAX_PARAM, tParam, x))
Create an instance of the class referred to by idKey.
Definition: NekFactory.hpp:162
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
LibUtilities::TranspositionSharedPtr m_transposition
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
LibUtilities::NektarFFTSharedPtr m_FFT_deal
NekDouble m_lhom
Width of homogeneous direction.
NektarFFTFactory & GetNektarFFTFactory()
Definition: NektarFFT.cpp:69
BasisManagerT & BasisManager(void)
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
double NekDouble
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:907
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Definition: ErrorUtil.hpp:213
static const BasisKey NullBasisKey(eNoBasisType, 0, NullPointsKey)
Defines a null basis with no type or points.
LibUtilities::NektarFFTSharedPtr m_FFT
Nektar::MultiRegions::ExpListHomogeneous1D::ExpListHomogeneous1D ( const ExpListHomogeneous1D In)

Copy constructor.

Parameters
InExpListHomogeneous1D object to copy.

Definition at line 106 of file ExpListHomogeneous1D.cpp.

References m_planes.

106  :
107  ExpList(In,false),
108  m_transposition(In.m_transposition),
109  m_StripZcomm(In.m_StripZcomm),
110  m_useFFT(In.m_useFFT),
111  m_FFT(In.m_FFT),
112  m_tmpIN(In.m_tmpIN),
113  m_tmpOUT(In.m_tmpOUT),
114  m_homogeneousBasis(In.m_homogeneousBasis),
115  m_lhom(In.m_lhom),
116  m_homogeneous1DBlockMat(In.m_homogeneous1DBlockMat),
117  m_dealiasing(In.m_dealiasing),
118  m_padsize(In.m_padsize)
119  {
120  m_planes = Array<OneD, ExpListSharedPtr>(In.m_planes.num_elements());
121  }
ExpList()
The default constructor.
Definition: ExpList.cpp:93
LibUtilities::TranspositionSharedPtr m_transposition
NekDouble m_lhom
Width of homogeneous direction.
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Array< OneD, ExpListSharedPtr > m_planes
LibUtilities::NektarFFTSharedPtr m_FFT
Nektar::MultiRegions::ExpListHomogeneous1D::ExpListHomogeneous1D ( const ExpListHomogeneous1D In,
const std::vector< unsigned int > &  eIDs 
)

Definition at line 123 of file ExpListHomogeneous1D.cpp.

References m_planes.

124  :
125  ExpList(In,eIDs,false),
126  m_transposition(In.m_transposition),
127  m_useFFT(In.m_useFFT),
128  m_FFT(In.m_FFT),
129  m_tmpIN(In.m_tmpIN),
130  m_tmpOUT(In.m_tmpOUT),
131  m_homogeneousBasis(In.m_homogeneousBasis),
132  m_lhom(In.m_lhom),
134  m_dealiasing(In.m_dealiasing),
135  m_padsize(In.m_padsize)
136  {
137  m_planes = Array<OneD, ExpListSharedPtr>(In.m_planes.num_elements());
138  }
ExpList()
The default constructor.
Definition: ExpList.cpp:93
LibUtilities::TranspositionSharedPtr m_transposition
General purpose memory allocation routines with the ability to allocate from thread specific memory p...
NekDouble m_lhom
Width of homogeneous direction.
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Array< OneD, ExpListSharedPtr > m_planes
LibUtilities::NektarFFTSharedPtr m_FFT
Nektar::MultiRegions::ExpListHomogeneous1D::~ExpListHomogeneous1D ( )
virtual

Destructor.

Destructor

Definition at line 143 of file ExpListHomogeneous1D.cpp.

144  {
145  }

Member Function Documentation

void Nektar::MultiRegions::ExpListHomogeneous1D::DealiasedProd ( const Array< OneD, NekDouble > &  inarray1,
const Array< OneD, NekDouble > &  inarray2,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate = eLocal 
)
inline

Definition at line 300 of file ExpListHomogeneous1D.h.

References v_DealiasedProd().

304  {
305  v_DealiasedProd(inarray1,inarray2,outarray,coeffstate);
306  }
virtual void v_DealiasedProd(const Array< OneD, NekDouble > &inarray1, const Array< OneD, NekDouble > &inarray2, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
DNekBlkMatSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::GenHomogeneous1DBlockMatrix ( Homogeneous1DMatType  mattype,
CoeffState  coeffstate = eLocal 
) const
protected

Definition at line 524 of file ExpListHomogeneous1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::eBackwardsCoeffSpace1D, Nektar::StdRegions::eBwdTrans, Nektar::eDIAGONAL, Nektar::MultiRegions::eForwardsCoeffSpace1D, Nektar::MultiRegions::eForwardsPhysSpace1D, Nektar::LibUtilities::eFourierHalfModeIm, Nektar::LibUtilities::eFourierHalfModeRe, Nektar::StdRegions::eFwdTrans, Nektar::MultiRegions::ExpList::m_comm, m_homogeneousBasis, and m_planes.

Referenced by GetHomogeneous1DBlockMatrix().

525  {
526  DNekMatSharedPtr loc_mat;
527  DNekBlkMatSharedPtr BlkMatrix;
528  int n_exp = 0;
529  int num_trans_per_proc = 0;
530 
531  if((mattype == eForwardsCoeffSpace1D)
532  ||(mattype == eBackwardsCoeffSpace1D)) // will operate on m_coeffs
533  {
534  n_exp = m_planes[0]->GetNcoeffs();
535  }
536  else
537  {
538  n_exp = m_planes[0]->GetTotPoints(); // will operatore on m_phys
539  }
540 
541  num_trans_per_proc = n_exp/m_comm->GetColumnComm()->GetSize() + (n_exp%m_comm->GetColumnComm()->GetSize() > 0);
542 
543  Array<OneD,unsigned int> nrows(num_trans_per_proc);
544  Array<OneD,unsigned int> ncols(num_trans_per_proc);
545 
546  if((mattype == eForwardsCoeffSpace1D)||(mattype == eForwardsPhysSpace1D))
547  {
548  nrows = Array<OneD, unsigned int>(num_trans_per_proc,m_homogeneousBasis->GetNumModes());
549  ncols = Array<OneD, unsigned int>(num_trans_per_proc,m_homogeneousBasis->GetNumPoints());
550  }
551  else
552  {
553  nrows = Array<OneD, unsigned int>(num_trans_per_proc,m_homogeneousBasis->GetNumPoints());
554  ncols = Array<OneD, unsigned int>(num_trans_per_proc,m_homogeneousBasis->GetNumModes());
555  }
556 
557  MatrixStorage blkmatStorage = eDIAGONAL;
558  BlkMatrix = MemoryManager<DNekBlkMat>
559  ::AllocateSharedPtr(nrows,ncols,blkmatStorage);
560 
561  //Half Mode
563  {
564  StdRegions::StdPointExp StdPoint(m_homogeneousBasis->GetBasisKey());
565 
566  if((mattype == eForwardsCoeffSpace1D)||(mattype == eForwardsPhysSpace1D))
567  {
569  StdPoint.DetShapeType(),
570  StdPoint);
571 
572  loc_mat = StdPoint.GetStdMatrix(matkey);
573  }
574  else
575  {
577  StdPoint.DetShapeType(),
578  StdPoint);
579 
580  loc_mat = StdPoint.GetStdMatrix(matkey);
581  }
582  }
583  //other cases
584  else
585  {
586  StdRegions::StdSegExp StdSeg(m_homogeneousBasis->GetBasisKey());
587 
588  if((mattype == eForwardsCoeffSpace1D)||(mattype == eForwardsPhysSpace1D))
589  {
591  StdSeg.DetShapeType(),
592  StdSeg);
593 
594  loc_mat = StdSeg.GetStdMatrix(matkey);
595  }
596  else
597  {
599  StdSeg.DetShapeType(),
600  StdSeg);
601 
602  loc_mat = StdSeg.GetStdMatrix(matkey);
603  }
604 
605  }
606 
607  // set up array of block matrices.
608  for(int i = 0; i < num_trans_per_proc; ++i)
609  {
610  BlkMatrix->SetBlock(i,i,loc_mat);
611  }
612 
613  return BlkMatrix;
614  }
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
Class representing a segment element in reference space.
Definition: StdSegExp.h:54
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Fourier Modified expansions with just the real part of the first mode .
Definition: BasisType.h:59
Array< OneD, ExpListSharedPtr > m_planes
Fourier Modified expansions with just the imaginary part of the first mode .
Definition: BasisType.h:60
boost::shared_ptr< DNekBlkMat > DNekBlkMatSharedPtr
Definition: NekTypeDefs.hpp:72
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:907
DNekBlkMatSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::GetHomogeneous1DBlockMatrix ( Homogeneous1DMatType  mattype,
CoeffState  coeffstate = eLocal 
) const
protected

Definition at line 508 of file ExpListHomogeneous1D.cpp.

References GenHomogeneous1DBlockMatrix(), Nektar::iterator, and m_homogeneous1DBlockMat.

Referenced by Homogeneous1DTrans().

509  {
510  Homo1DBlockMatrixMap::iterator matrixIter = m_homogeneous1DBlockMat->find(mattype);
511 
512  if(matrixIter == m_homogeneous1DBlockMat->end())
513  {
514  return ((*m_homogeneous1DBlockMat)[mattype] =
515  GenHomogeneous1DBlockMatrix(mattype,coeffstate));
516  }
517  else
518  {
519  return matrixIter->second;
520  }
521  }
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
DNekBlkMatSharedPtr GenHomogeneous1DBlockMatrix(Homogeneous1DMatType mattype, CoeffState coeffstate=eLocal) const
LibUtilities::BasisSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::GetHomogeneousBasis ( void  )
inline

Definition at line 120 of file ExpListHomogeneous1D.h.

References m_homogeneousBasis.

Referenced by v_GetHomogeneousBasis().

121  {
122  return m_homogeneousBasis;
123  }
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
ExpListSharedPtr& Nektar::MultiRegions::ExpListHomogeneous1D::GetPlane ( int  n)
inline

Definition at line 134 of file ExpListHomogeneous1D.h.

References m_planes.

Referenced by v_GetPlane().

135  {
136  return m_planes[n];
137  }
Array< OneD, ExpListSharedPtr > m_planes
NekDouble Nektar::MultiRegions::ExpListHomogeneous1D::GetSpecVanVisc ( const int  k)
inlineprotected

Definition at line 165 of file ExpListHomogeneous1D.h.

References m_specVanVisc.

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

166  {
167  NekDouble returnval = 0.0;
168 
169  if(m_specVanVisc.num_elements())
170  {
171  returnval = m_specVanVisc[k];
172  }
173 
174  return returnval;
175  }
Array< OneD, NekDouble > m_specVanVisc
Spectral vanishing Viscosity coefficient for stabilisation.
double NekDouble
void Nektar::MultiRegions::ExpListHomogeneous1D::Homogeneous1DTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
bool  IsForwards,
CoeffState  coeffstate = eLocal,
bool  Shuff = true,
bool  UnShuff = true 
)

Homogeneous transform Bwd/Fwd (MVM and FFT)

Definition at line 375 of file ExpListHomogeneous1D.cpp.

References Nektar::MultiRegions::eBackwardsCoeffSpace1D, Nektar::MultiRegions::eBackwardsPhysSpace1D, Nektar::MultiRegions::eForwardsCoeffSpace1D, Nektar::MultiRegions::eForwardsPhysSpace1D, Nektar::eWrapper, Nektar::LibUtilities::eXYtoZ, Nektar::LibUtilities::eZtoXY, GetHomogeneous1DBlockMatrix(), Nektar::MultiRegions::ExpList::m_comm, m_FFT, m_homogeneousBasis, Nektar::MultiRegions::ExpList::m_npoints, m_planes, Nektar::MultiRegions::ExpList::m_session, m_StripZcomm, m_tmpIN, m_tmpOUT, m_transposition, m_useFFT, and Vmath::Vcopy().

Referenced by v_HomogeneousBwdTrans(), and v_HomogeneousFwdTrans().

380  {
381  int num_dofs;
382 
383  if(IsForwards)
384  {
385  num_dofs = inarray.num_elements();
386  }
387  else
388  {
389  num_dofs = outarray.num_elements();
390  }
391 
392  if(m_useFFT)
393  {
394  int num_points_per_plane = num_dofs/m_planes.num_elements();
395  int num_dfts_per_proc;
396  if(!m_session->DefinesSolverInfo("HomoStrip"))
397  {
398  int nP = m_comm->GetColumnComm()->GetSize();
399  num_dfts_per_proc = num_points_per_plane / nP
400  + (num_points_per_plane % nP > 0);
401  }
402  else
403  {
404  int nP = m_StripZcomm->GetSize();
405  num_dfts_per_proc = num_points_per_plane / nP
406  + (num_points_per_plane % nP > 0);
407  }
408 
409  Array<OneD, NekDouble> fft_in (num_dfts_per_proc*m_homogeneousBasis->GetNumPoints(),0.0);
410  Array<OneD, NekDouble> fft_out(num_dfts_per_proc*m_homogeneousBasis->GetNumPoints(),0.0);
411 
412  if(Shuff)
413  {
414  m_transposition->Transpose(inarray,fft_in,false,LibUtilities::eXYtoZ);
415  }
416  else
417  {
418  Vmath::Vcopy(num_dfts_per_proc*m_homogeneousBasis->GetNumPoints(),
419  inarray,1,fft_in,1);
420  }
421 
422  if(IsForwards)
423  {
424  for(int i = 0 ; i < num_dfts_per_proc ; i++)
425  {
426  m_FFT->FFTFwdTrans(m_tmpIN = fft_in + i*m_homogeneousBasis->GetNumPoints(), m_tmpOUT = fft_out + i*m_homogeneousBasis->GetNumPoints());
427  }
428  }
429  else
430  {
431  for(int i = 0 ; i < num_dfts_per_proc ; i++)
432  {
433  m_FFT->FFTBwdTrans(m_tmpIN = fft_in + i*m_homogeneousBasis->GetNumPoints(), m_tmpOUT = fft_out + i*m_homogeneousBasis->GetNumPoints());
434  }
435  }
436 
437  if(UnShuff)
438  {
439  m_transposition->Transpose(fft_out,outarray,false,LibUtilities::eZtoXY);
440  }
441  else
442  {
443  Vmath::Vcopy(num_dfts_per_proc*m_homogeneousBasis->GetNumPoints(),
444  fft_out,1,outarray,1);
445  }
446  }
447  else
448  {
449  DNekBlkMatSharedPtr blkmat;
450 
451  if(num_dofs == m_npoints) //transform phys space
452  {
453  if(IsForwards)
454  {
456  }
457  else
458  {
460  }
461  }
462  else
463  {
464  if(IsForwards)
465  {
467  }
468  else
469  {
471  }
472  }
473 
474  int nrows = blkmat->GetRows();
475  int ncols = blkmat->GetColumns();
476 
477  Array<OneD, NekDouble> sortedinarray(ncols,0.0);
478  Array<OneD, NekDouble> sortedoutarray(nrows,0.0);
479 
480  if(Shuff)
481  {
482  m_transposition->Transpose(inarray,sortedinarray,!IsForwards,LibUtilities::eXYtoZ);
483  }
484  else
485  {
486  Vmath::Vcopy(ncols,inarray,1,sortedinarray,1);
487  }
488 
489  // Create NekVectors from the given data arrays
490  NekVector<NekDouble> in (ncols,sortedinarray,eWrapper);
491  NekVector<NekDouble> out(nrows,sortedoutarray,eWrapper);
492 
493  // Perform matrix-vector multiply.
494  out = (*blkmat)*in;
495 
496  if(UnShuff)
497  {
498  m_transposition->Transpose(sortedoutarray,outarray,IsForwards,LibUtilities::eZtoXY);
499  }
500  else
501  {
502  Vmath::Vcopy(nrows,sortedinarray,1,outarray,1);
503  }
504 
505  }
506  }
DNekBlkMatSharedPtr GetHomogeneous1DBlockMatrix(Homogeneous1DMatType mattype, CoeffState coeffstate=eLocal) const
LibUtilities::TranspositionSharedPtr m_transposition
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
boost::shared_ptr< DNekBlkMat > DNekBlkMatSharedPtr
Definition: NekTypeDefs.hpp:72
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:907
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
LibUtilities::NektarFFTSharedPtr m_FFT
void Nektar::MultiRegions::ExpListHomogeneous1D::HomogeneousBwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate = eLocal,
bool  Shuff = true,
bool  UnShuff = true 
)
inline

Definition at line 291 of file ExpListHomogeneous1D.h.

References v_HomogeneousBwdTrans().

Referenced by v_BwdTrans(), v_BwdTrans_IterPerExp(), v_DealiasedProd(), and v_PhysDeriv().

296  {
297  v_HomogeneousBwdTrans(inarray,outarray,coeffstate,Shuff,UnShuff);
298  }
virtual void v_HomogeneousBwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
void Nektar::MultiRegions::ExpListHomogeneous1D::HomogeneousFwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate = eLocal,
bool  Shuff = true,
bool  UnShuff = true 
)
inline

Definition at line 282 of file ExpListHomogeneous1D.h.

References v_HomogeneousFwdTrans().

Referenced by v_DealiasedProd(), v_FwdTrans(), v_FwdTrans_IterPerExp(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_HelmSolve(), Nektar::MultiRegions::DisContField3DHomogeneous1D::v_HelmSolve(), and v_PhysDeriv().

287  {
288  v_HomogeneousFwdTrans(inarray,outarray,coeffstate,Shuff,UnShuff);
289  }
virtual void v_HomogeneousFwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
void Nektar::MultiRegions::ExpListHomogeneous1D::PhysDeriv ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  out_d0,
Array< OneD, NekDouble > &  out_d1,
Array< OneD, NekDouble > &  out_d2 
)

Definition at line 1155 of file ExpListHomogeneous1D.cpp.

References v_PhysDeriv().

1160  {
1161  v_PhysDeriv(inarray,out_d0,out_d1,out_d2);
1162  }
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)
void Nektar::MultiRegions::ExpListHomogeneous1D::PhysDeriv ( Direction  edir,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  out_d 
)

Definition at line 1164 of file ExpListHomogeneous1D.cpp.

References v_PhysDeriv().

1167  {
1168  v_PhysDeriv(edir,inarray,out_d);
1169  }
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)
void Nektar::MultiRegions::ExpListHomogeneous1D::v_AppendFieldData ( LibUtilities::FieldDefinitionsSharedPtr fielddef,
std::vector< NekDouble > &  fielddata 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 722 of file ExpListHomogeneous1D.cpp.

References Nektar::MultiRegions::ExpList::m_coeffs.

723  {
724  v_AppendFieldData(fielddef,fielddata,m_coeffs);
725  }
Array< OneD, NekDouble > m_coeffs
Concatenation of all local expansion coefficients.
Definition: ExpList.h:939
virtual void v_AppendFieldData(LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata)
void Nektar::MultiRegions::ExpListHomogeneous1D::v_AppendFieldData ( LibUtilities::FieldDefinitionsSharedPtr fielddef,
std::vector< NekDouble > &  fielddata,
Array< OneD, NekDouble > &  coeffs 
)
protectedvirtual

This routine appends the data from the expansion list into the output format where each element is given by looping over its Fourier modes where as data in the expandion is stored with all consecutive elements and then the Fourier modes

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 697 of file ExpListHomogeneous1D.cpp.

References Nektar::MultiRegions::ExpList::m_coeff_offset, and m_planes.

698  {
699  int i,n;
700  int ncoeffs_per_plane = m_planes[0]->GetNcoeffs();
701 
702  // Determine mapping from element ids to location in
703  // expansion list
704  map<int, int> ElmtID_to_ExpID;
705  for(i = 0; i < m_planes[0]->GetExpSize(); ++i)
706  {
707  ElmtID_to_ExpID[(*m_exp)[i]->GetGeom()->GetGlobalID()] = i;
708  }
709 
710  for(i = 0; i < fielddef->m_elementIDs.size(); ++i)
711  {
712  int eid = ElmtID_to_ExpID[fielddef->m_elementIDs[i]];
713  int datalen = (*m_exp)[eid]->GetNcoeffs();
714 
715  for(n = 0; n < m_planes.num_elements(); ++n)
716  {
717  fielddata.insert(fielddata.end(),&coeffs[m_coeff_offset[eid]+n*ncoeffs_per_plane],&coeffs[m_coeff_offset[eid]+n*ncoeffs_per_plane]+datalen);
718  }
719  }
720  }
Array< OneD, int > m_coeff_offset
Offset of elemental data into the array m_coeffs.
Definition: ExpList.h:988
Array< OneD, ExpListSharedPtr > m_planes
void Nektar::MultiRegions::ExpListHomogeneous1D::v_BwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate 
)
protectedvirtual

Backward transform

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 299 of file ExpListHomogeneous1D.cpp.

References HomogeneousBwdTrans(), m_planes, and Nektar::MultiRegions::ExpList::m_WaveSpace.

300  {
301  int cnt = 0, cnt1 = 0;
302  Array<OneD, NekDouble> tmparray;
303 
304  for(int n = 0; n < m_planes.num_elements(); ++n)
305  {
306  m_planes[n]->BwdTrans(inarray+cnt, tmparray = outarray + cnt1,
307  coeffstate);
308  cnt += m_planes[n]->GetNcoeffs();
309  cnt1 += m_planes[n]->GetTotPoints();
310  }
311  if(!m_WaveSpace)
312  {
313  HomogeneousBwdTrans(outarray,outarray);
314  }
315  }
void HomogeneousBwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
Array< OneD, ExpListSharedPtr > m_planes
void Nektar::MultiRegions::ExpListHomogeneous1D::v_BwdTrans_IterPerExp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Backward transform element by element

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 320 of file ExpListHomogeneous1D.cpp.

References HomogeneousBwdTrans(), m_planes, and Nektar::MultiRegions::ExpList::m_WaveSpace.

321  {
322  int cnt = 0, cnt1 = 0;
323  Array<OneD, NekDouble> tmparray;
324 
325  for(int n = 0; n < m_planes.num_elements(); ++n)
326  {
327  m_planes[n]->BwdTrans_IterPerExp(inarray+cnt, tmparray = outarray + cnt1);
328 
329  cnt += m_planes[n]->GetNcoeffs();
330  cnt1 += m_planes[n]->GetTotPoints();
331  }
332  if(!m_WaveSpace)
333  {
334  HomogeneousBwdTrans(outarray,outarray);
335  }
336  }
void HomogeneousBwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
Array< OneD, ExpListSharedPtr > m_planes
void Nektar::MultiRegions::ExpListHomogeneous1D::v_DealiasedProd ( const Array< OneD, NekDouble > &  inarray1,
const Array< OneD, NekDouble > &  inarray2,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate = eLocal 
)
protectedvirtual

Dealiasing routine

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 170 of file ExpListHomogeneous1D.cpp.

References Nektar::LibUtilities::eXYtoZ, Nektar::LibUtilities::eZtoXY, HomogeneousBwdTrans(), HomogeneousFwdTrans(), Nektar::MultiRegions::ExpList::m_comm, m_FFT_deal, m_homogeneousBasis, m_padsize, m_planes, Nektar::MultiRegions::ExpList::m_session, m_StripZcomm, m_transposition, Vmath::Vcopy(), and Vmath::Vmul().

Referenced by DealiasedProd().

174  {
175  // inarray1 = first term of the product in full physical space
176  // inarray2 = second term of the product in full physical space
177  // dealiased product stored in outarray
178 
179  int num_dofs = inarray1.num_elements();
180 
181  int N = m_homogeneousBasis->GetNumPoints();
182 
183  Array<OneD, NekDouble> V1(num_dofs);
184  Array<OneD, NekDouble> V2(num_dofs);
185  Array<OneD, NekDouble> V1V2(num_dofs);
186 
187  HomogeneousFwdTrans(inarray1,V1,coeffstate);
188  HomogeneousFwdTrans(inarray2,V2,coeffstate);
189 
190  int num_points_per_plane = num_dofs/m_planes.num_elements();
191  int num_proc;
192  if(!m_session->DefinesSolverInfo("HomoStrip"))
193  {
194  num_proc = m_comm->GetColumnComm()->GetSize();
195  }
196  else
197  {
198  num_proc = m_StripZcomm->GetSize();
199  }
200  int num_dfts_per_proc = num_points_per_plane / num_proc
201  + (num_points_per_plane % num_proc > 0);
202 
203  Array<OneD, NekDouble> ShufV1(num_dfts_per_proc*N,0.0);
204  Array<OneD, NekDouble> ShufV2(num_dfts_per_proc*N,0.0);
205  Array<OneD, NekDouble> ShufV1V2(num_dfts_per_proc*N,0.0);
206 
207  Array<OneD, NekDouble> ShufV1_PAD_coef(m_padsize,0.0);
208  Array<OneD, NekDouble> ShufV2_PAD_coef(m_padsize,0.0);
209  Array<OneD, NekDouble> ShufV1_PAD_phys(m_padsize,0.0);
210  Array<OneD, NekDouble> ShufV2_PAD_phys(m_padsize,0.0);
211 
212  Array<OneD, NekDouble> ShufV1V2_PAD_coef(m_padsize,0.0);
213  Array<OneD, NekDouble> ShufV1V2_PAD_phys(m_padsize,0.0);
214 
215  m_transposition->Transpose(V1, ShufV1, false, LibUtilities::eXYtoZ);
216  m_transposition->Transpose(V2, ShufV2, false, LibUtilities::eXYtoZ);
217 
218  // Looping on the pencils
219  for(int i = 0 ; i < num_dfts_per_proc ; i++)
220  {
221  // Copying the i-th pencil pf lenght N into a bigger
222  // pencil of lenght 2N We are in Fourier space
223  Vmath::Vcopy(N, &(ShufV1[i*N]), 1, &(ShufV1_PAD_coef[0]), 1);
224  Vmath::Vcopy(N, &(ShufV2[i*N]), 1, &(ShufV2_PAD_coef[0]), 1);
225 
226  // Moving to physical space using the padded system
227  m_FFT_deal->FFTBwdTrans(ShufV1_PAD_coef, ShufV1_PAD_phys);
228  m_FFT_deal->FFTBwdTrans(ShufV2_PAD_coef, ShufV2_PAD_phys);
229 
230  // Perfroming the vectors multiplication in physical space on
231  // the padded system
232  Vmath::Vmul(m_padsize, ShufV1_PAD_phys, 1,
233  ShufV2_PAD_phys, 1,
234  ShufV1V2_PAD_phys, 1);
235 
236  // Moving back the result (V1*V2)_phys in Fourier space, padded
237  // system
238  m_FFT_deal->FFTFwdTrans(ShufV1V2_PAD_phys, ShufV1V2_PAD_coef);
239 
240  // Copying the first half of the padded pencil in the full
241  // vector (Fourier space)
242  Vmath::Vcopy(N, &(ShufV1V2_PAD_coef[0]), 1,
243  &(ShufV1V2[i*N]), 1);
244  }
245 
246  m_transposition->Transpose(ShufV1V2, V1V2, false,
248 
249  // Moving the results in physical space for the output
250  HomogeneousBwdTrans(V1V2, outarray, coeffstate);
251  }
void HomogeneousBwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
LibUtilities::TranspositionSharedPtr m_transposition
LibUtilities::NektarFFTSharedPtr m_FFT_deal
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
void HomogeneousFwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:907
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Vmul(int n, const T *x, const int incx, const T *y, const int incy, T *z, const int incz)
Multiply vector z = x*y.
Definition: Vmath.cpp:169
void Nektar::MultiRegions::ExpListHomogeneous1D::v_ExtractCoeffsToCoeffs ( const boost::shared_ptr< ExpList > &  fromExpList,
const Array< OneD, const NekDouble > &  fromCoeffs,
Array< OneD, NekDouble > &  toCoeffs 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 853 of file ExpListHomogeneous1D.cpp.

References m_planes.

855  {
856  int i;
857  int fromNcoeffs_per_plane = fromExpList->GetPlane(0)->GetNcoeffs();
858  int toNcoeffs_per_plane = m_planes[0]->GetNcoeffs();
859  Array<OneD, NekDouble> tocoeffs_tmp, fromcoeffs_tmp;
860 
861  for(i = 0; i < m_planes.num_elements(); ++i)
862  {
863  m_planes[i]->ExtractCoeffsToCoeffs(fromExpList->GetPlane(i),fromcoeffs_tmp = fromCoeffs + fromNcoeffs_per_plane*i, tocoeffs_tmp = toCoeffs + toNcoeffs_per_plane*i);
864  }
865  }
Array< OneD, ExpListSharedPtr > m_planes
void Nektar::MultiRegions::ExpListHomogeneous1D::v_ExtractDataToCoeffs ( LibUtilities::FieldDefinitionsSharedPtr fielddef,
std::vector< NekDouble > &  fielddata,
std::string &  field,
Array< OneD, NekDouble > &  coeffs 
)
protectedvirtual

Extract data from raw field data into expansion list.

Parameters
fielddefField definitions.
fielddataData for associated field.
fieldField variable name.
coeffsResulting coefficient array.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 728 of file ExpListHomogeneous1D.cpp.

References Nektar::LibUtilities::eFourierSingleMode, Nektar::LibUtilities::GetNumberOfCoefficients(), Nektar::iterator, Nektar::MultiRegions::ExpList::m_coeff_offset, Nektar::MultiRegions::ExpList::m_exp, m_homogeneousBasis, m_planes, m_transposition, and Vmath::Vcopy().

733  {
734  int i,n;
735  int offset = 0;
736  int nzmodes = 1;
737  int datalen = fielddata.size()/fielddef->m_fields.size();
738  std::vector<unsigned int> fieldDefHomoZids;
739 
740 
741  // Find data location according to field definition
742  for(i = 0; i < fielddef->m_fields.size(); ++i)
743  {
744  if(fielddef->m_fields[i] == field)
745  {
746  break;
747  }
748  offset += datalen;
749  }
750 
751  if(i == fielddef->m_fields.size())
752  {
753  cout << "Field "<< field<< "not found in data file. " << endl;
754  }
755  else
756  {
757 
758  int modes_offset = 0;
759  int planes_offset = 0;
760  Array<OneD, NekDouble> coeff_tmp;
762  int IDoffset = 0;
763 
764  // introduce a 2 plane offset for single mode case so can
765  // be post-processed or used in MultiMode expansion.
767  {
768  IDoffset = 2;
769  }
770 
771  // Build map of plane IDs lying on this processor.
772  std::map<int,int> homoZids;
773  for (i = 0; i < m_planes.num_elements(); ++i)
774  {
775  homoZids[m_transposition->GetPlaneID(i)+IDoffset] = i;
776  }
777 
778  if(fielddef->m_numHomogeneousDir)
779  {
780  nzmodes = fielddef->m_homogeneousZIDs.size();
781  fieldDefHomoZids = fielddef->m_homogeneousZIDs;
782  }
783  else // input file is 2D and so set nzmodes to 1
784  {
785  nzmodes = 1;
786  fieldDefHomoZids.push_back(0);
787  }
788 
789  // Determine mapping from element ids to location in expansion list.
790  map<int, int> ElmtID_to_ExpID;
791  for(i = 0; i < m_planes[0]->GetExpSize(); ++i)
792  {
793  ElmtID_to_ExpID[(*m_exp)[i]->GetGeom()->GetGlobalID()] = i;
794  }
795 
796 
797  // calculate number of modes in the current partition
798  int ncoeffs_per_plane = m_planes[0]->GetNcoeffs();
799 
800  for(i = 0; i < fielddef->m_elementIDs.size(); ++i)
801  {
802  if(fielddef->m_uniOrder == true) // reset modes_offset to zero
803  {
804  modes_offset = 0;
805  }
806 
807  int datalen = LibUtilities::GetNumberOfCoefficients(fielddef->m_shapeType,
808  fielddef->m_numModes,
809  modes_offset);
810 
811  it = ElmtID_to_ExpID.find(fielddef->m_elementIDs[i]);
812 
813  // ensure element is on this partition for parallel case.
814  if(it == ElmtID_to_ExpID.end())
815  {
816  // increase offset for correct FieldData access
817  offset += datalen*nzmodes;
818  modes_offset += (*m_exp)[0]->GetNumBases() +
819  fielddef->m_numHomogeneousDir;
820  continue;
821  }
822 
823  int eid = it->second;
824 
825 
826  for(n = 0; n < nzmodes; ++n, offset += datalen)
827  {
828 
829  it = homoZids.find(fieldDefHomoZids[n]);
830 
831  // Check to make sure this mode number lies in this field.
832  if (it == homoZids.end())
833  {
834  continue;
835  }
836 
837  planes_offset = it->second;
838  if(datalen == (*m_exp)[eid]->GetNcoeffs())
839  {
840  Vmath::Vcopy(datalen,&fielddata[offset],1,&coeffs[m_coeff_offset[eid]+planes_offset*ncoeffs_per_plane],1);
841  }
842  else // unpack data to new order
843  {
844  (*m_exp)[eid]->ExtractDataToCoeffs(&fielddata[offset], fielddef->m_numModes,modes_offset,&coeffs[m_coeff_offset[eid] + planes_offset*ncoeffs_per_plane]);
845  }
846  }
847  modes_offset += (*m_exp)[0]->GetNumBases() + fielddef->m_numHomogeneousDir;
848  }
849  }
850  }
LibUtilities::TranspositionSharedPtr m_transposition
int GetNumberOfCoefficients(ShapeType shape, std::vector< unsigned int > &modes, int offset)
Definition: ShapeType.hpp:312
Array< OneD, int > m_coeff_offset
Offset of elemental data into the array m_coeffs.
Definition: ExpList.h:988
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
boost::shared_ptr< LocalRegions::ExpansionVector > m_exp
The list of local expansions.
Definition: ExpList.h:977
Array< OneD, ExpListSharedPtr > m_planes
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
Fourier ModifiedExpansion with just the first mode .
Definition: BasisType.h:58
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Nektar::MultiRegions::ExpListHomogeneous1D::v_FwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate 
)
protectedvirtual

Forward transform

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 256 of file ExpListHomogeneous1D.cpp.

References HomogeneousFwdTrans(), m_planes, and Nektar::MultiRegions::ExpList::m_WaveSpace.

257  {
258  int cnt = 0, cnt1 = 0;
259  Array<OneD, NekDouble> tmparray;
260 
261  for(int n = 0; n < m_planes.num_elements(); ++n)
262  {
263  m_planes[n]->FwdTrans(inarray+cnt, tmparray = outarray + cnt1,
264  coeffstate);
265  cnt += m_planes[n]->GetTotPoints();
266 
267  cnt1 += m_planes[n]->GetNcoeffs(); // need to skip ncoeffs
268  }
269  if(!m_WaveSpace)
270  {
271  HomogeneousFwdTrans(outarray,outarray,coeffstate);
272  }
273  }
Array< OneD, ExpListSharedPtr > m_planes
void HomogeneousFwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
void Nektar::MultiRegions::ExpListHomogeneous1D::v_FwdTrans_IterPerExp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Forward transform element by element

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 278 of file ExpListHomogeneous1D.cpp.

References HomogeneousFwdTrans(), m_planes, and Nektar::MultiRegions::ExpList::m_WaveSpace.

279  {
280  int cnt = 0, cnt1 = 0;
281  Array<OneD, NekDouble> tmparray;
282 
283  //spectral element FwdTrans plane by plane
284  for(int n = 0; n < m_planes.num_elements(); ++n)
285  {
286  m_planes[n]->FwdTrans_IterPerExp(inarray+cnt, tmparray = outarray + cnt1);
287  cnt += m_planes[n]->GetTotPoints();
288  cnt1 += m_planes[n]->GetNcoeffs();
289  }
290  if(!m_WaveSpace)
291  {
292  HomogeneousFwdTrans(outarray,outarray);
293  }
294  }
Array< OneD, ExpListSharedPtr > m_planes
void HomogeneousFwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
std::vector< LibUtilities::FieldDefinitionsSharedPtr > Nektar::MultiRegions::ExpListHomogeneous1D::v_GetFieldDefinitions ( void  )
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 616 of file ExpListHomogeneous1D.cpp.

References Nektar::LibUtilities::eFourierSingleMode, m_homogeneousBasis, m_lhom, m_planes, Nektar::MultiRegions::ExpList::m_session, and m_transposition.

617  {
618  std::vector<LibUtilities::FieldDefinitionsSharedPtr> returnval;
619 
620  // Set up Homogeneous length details.
622 
623  std::vector<NekDouble> HomoLen;
624  HomoLen.push_back(m_lhom);
625 
626  std::vector<unsigned int> StripsIDs;
627 
628  bool strips;
629  m_session->MatchSolverInfo("HomoStrip","True",strips,false);
630  if (strips)
631  {
632  StripsIDs.push_back(m_transposition->GetStripID());
633  }
634 
635  std::vector<unsigned int> PlanesIDs;
636  int IDoffset = 0;
637 
638  // introduce a 2 plane offset for single mode case so can
639  // be post-processed or used in MultiMode expansion.
641  {
642  IDoffset = 2;
643  }
644 
645  for(int i = 0; i < m_planes.num_elements(); i++)
646  {
647  PlanesIDs.push_back(m_transposition->GetPlaneID(i)+IDoffset);
648  }
649 
650  m_planes[0]->GeneralGetFieldDefinitions(returnval, 1, HomoBasis,
651  HomoLen, strips, StripsIDs, PlanesIDs);
652  return returnval;
653  }
LibUtilities::TranspositionSharedPtr m_transposition
NekDouble m_lhom
Width of homogeneous direction.
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
Fourier ModifiedExpansion with just the first mode .
Definition: BasisType.h:58
void Nektar::MultiRegions::ExpListHomogeneous1D::v_GetFieldDefinitions ( std::vector< LibUtilities::FieldDefinitionsSharedPtr > &  fielddef)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 655 of file ExpListHomogeneous1D.cpp.

References Nektar::LibUtilities::eFourierSingleMode, m_homogeneousBasis, m_lhom, m_planes, Nektar::MultiRegions::ExpList::m_session, and m_transposition.

656  {
657  // Set up Homogeneous length details.
659 
660  std::vector<NekDouble> HomoLen;
661  HomoLen.push_back(m_lhom);
662 
663  std::vector<unsigned int> StripsIDs;
664 
665  bool strips;
666  m_session->MatchSolverInfo("HomoStrip","True",strips,false);
667  if (strips)
668  {
669  StripsIDs.push_back(m_transposition->GetStripID());
670  }
671 
672  std::vector<unsigned int> PlanesIDs;
673  int IDoffset = 0;
674 
676  {
677  IDoffset = 2;
678  }
679 
680  for(int i = 0; i < m_planes.num_elements(); i++)
681  {
682  PlanesIDs.push_back(m_transposition->GetPlaneID(i)+IDoffset);
683  }
684 
685  // enforce NumHomoDir == 1 by direct call
686  m_planes[0]->GeneralGetFieldDefinitions(fielddef, 1, HomoBasis,
687  HomoLen, strips, StripsIDs, PlanesIDs);
688  }
LibUtilities::TranspositionSharedPtr m_transposition
NekDouble m_lhom
Width of homogeneous direction.
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
Fourier ModifiedExpansion with just the first mode .
Definition: BasisType.h:58
virtual LibUtilities::BasisSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::v_GetHomogeneousBasis ( void  )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 188 of file ExpListHomogeneous1D.h.

References GetHomogeneousBasis().

189  {
190  return GetHomogeneousBasis();
191  }
LibUtilities::BasisSharedPtr GetHomogeneousBasis(void)
NekDouble Nektar::MultiRegions::ExpListHomogeneous1D::v_GetHomoLen ( void  )
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1176 of file ExpListHomogeneous1D.cpp.

References m_lhom.

1177  {
1178  return m_lhom;
1179  }
NekDouble m_lhom
Width of homogeneous direction.
virtual int Nektar::MultiRegions::ExpListHomogeneous1D::v_GetNumElmts ( void  )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 183 of file ExpListHomogeneous1D.h.

184  {
185  return m_planes[0]->GetExpSize();
186  }
Array< OneD, ExpListSharedPtr > m_planes
virtual ExpListSharedPtr& Nektar::MultiRegions::ExpListHomogeneous1D::v_GetPlane ( int  n)
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 265 of file ExpListHomogeneous1D.h.

References GetPlane().

266  {
267  return GetPlane(n);
268  }
LibUtilities::TranspositionSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::v_GetTransposition ( void  )
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1171 of file ExpListHomogeneous1D.cpp.

References m_transposition.

1172  {
1173  return m_transposition;
1174  }
LibUtilities::TranspositionSharedPtr m_transposition
Array< OneD, const unsigned int > Nektar::MultiRegions::ExpListHomogeneous1D::v_GetZIDs ( void  )
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1181 of file ExpListHomogeneous1D.cpp.

References m_transposition.

1182  {
1183  return m_transposition->GetPlanesIDs();
1184  }
LibUtilities::TranspositionSharedPtr m_transposition
void Nektar::MultiRegions::ExpListHomogeneous1D::v_HomogeneousBwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate = eLocal,
bool  Shuff = true,
bool  UnShuff = true 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 157 of file ExpListHomogeneous1D.cpp.

References Homogeneous1DTrans().

Referenced by HomogeneousBwdTrans().

162  {
163  // Backwards trans
164  Homogeneous1DTrans(inarray,outarray,false,coeffstate,Shuff,UnShuff);
165  }
void Homogeneous1DTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool IsForwards, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
void Nektar::MultiRegions::ExpListHomogeneous1D::v_HomogeneousFwdTrans ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate = eLocal,
bool  Shuff = true,
bool  UnShuff = true 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 147 of file ExpListHomogeneous1D.cpp.

References Homogeneous1DTrans().

Referenced by HomogeneousFwdTrans().

152  {
153  // Forwards trans
154  Homogeneous1DTrans(inarray,outarray,true,coeffstate,Shuff,UnShuff);
155  }
void Homogeneous1DTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool IsForwards, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
void Nektar::MultiRegions::ExpListHomogeneous1D::v_IProductWRTBase ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
CoeffState  coeffstate 
)
protectedvirtual

Inner product

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 341 of file ExpListHomogeneous1D.cpp.

References m_planes.

342  {
343  int cnt = 0, cnt1 = 0;
344  Array<OneD, NekDouble> tmparray;
345 
346  for(int n = 0; n < m_planes.num_elements(); ++n)
347  {
348  m_planes[n]->IProductWRTBase(inarray+cnt, tmparray = outarray + cnt1,coeffstate);
349 
350  cnt1 += m_planes[n]->GetNcoeffs();
351  cnt += m_planes[n]->GetTotPoints();
352  }
353  }
Array< OneD, ExpListSharedPtr > m_planes
void Nektar::MultiRegions::ExpListHomogeneous1D::v_IProductWRTBase_IterPerExp ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Inner product element by element

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 358 of file ExpListHomogeneous1D.cpp.

References m_planes.

359  {
360  int cnt = 0, cnt1 = 0;
361  Array<OneD, NekDouble> tmparray;
362 
363  for(int n = 0; n < m_planes.num_elements(); ++n)
364  {
365  m_planes[n]->IProductWRTBase_IterPerExp(inarray+cnt, tmparray = outarray + cnt1);
366 
367  cnt1 += m_planes[n]->GetNcoeffs();
368  cnt += m_planes[n]->GetTotPoints();
369  }
370  }
Array< OneD, ExpListSharedPtr > m_planes
void Nektar::MultiRegions::ExpListHomogeneous1D::v_PhysDeriv ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  out_d0,
Array< OneD, NekDouble > &  out_d1,
Array< OneD, NekDouble > &  out_d2 
)
protectedvirtual

Given a function $f(\boldsymbol{x})$ evaluated at the quadrature points, this function calculates the derivatives $\frac{d}{dx_1}$, $\frac{d}{dx_2}$ and $\frac{d}{dx_3}$ of the function $f(\boldsymbol{x})$ at the same quadrature points. The local distribution of the quadrature points allows an elemental evaluation of the derivative. This is done by a call to the function StdRegions::StdExpansion::PhysDeriv.

Parameters
inarrayAn array of size $Q_{\mathrm{tot}}$ containing the values of the function $f(\boldsymbol{x})$ at the quadrature points $\boldsymbol{x}_i$.
out_d0The discrete evaluation of the derivative $\frac{d}{dx_1}$ will be stored in this array of size $Q_{\mathrm{tot}}$.
out_d1The discrete evaluation of the derivative $\frac{d}{dx_2}$ will be stored in this array of size $Q_{\mathrm{tot}}$. Note that if no memory is allocated for out_d1, the derivative $\frac{d}{dx_2}$ will not be calculated.
out_d2The discrete evaluation of the derivative $\frac{d}{dx_3}$ will be stored in this array of size $Q_{\mathrm{tot}}$. Note that if no memory is allocated for out_d2, the derivative $\frac{d}{dx_3}$ will not be calculated.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 933 of file ExpListHomogeneous1D.cpp.

References ASSERTL0, Nektar::LibUtilities::eFourier, Nektar::LibUtilities::eFourierHalfModeIm, Nektar::LibUtilities::eFourierHalfModeRe, Nektar::LibUtilities::eFourierSingleMode, Nektar::LibUtilities::eXYtoZ, Nektar::LibUtilities::eZtoXY, HomogeneousBwdTrans(), HomogeneousFwdTrans(), Nektar::MultiRegions::ExpList::m_comm, m_homogeneousBasis, m_lhom, m_planes, Nektar::MultiRegions::ExpList::m_session, m_StripZcomm, m_transposition, Nektar::MultiRegions::ExpList::m_WaveSpace, Nektar::NullNekDouble1DArray, sign, and Vmath::Smul().

Referenced by PhysDeriv().

937  {
938  int nT_pts = inarray.num_elements(); //number of total points = n. of Fourier points * n. of points per plane (nT_pts)
939  int nP_pts = nT_pts/m_planes.num_elements(); //number of points per plane = n of Fourier transform required (nP_pts)
940 
941  Array<OneD, NekDouble> temparray(nT_pts);
942  Array<OneD, NekDouble> outarray(nT_pts);
946 
947  for(int i = 0; i < m_planes.num_elements(); i++)
948  {
949  m_planes[i]->PhysDeriv(inarray + i*nP_pts ,tmp2 = out_d0 + i*nP_pts , tmp3 = out_d1 + i*nP_pts );
950  }
951 
952  if(out_d2 != NullNekDouble1DArray)
953  {
956  {
957  if(m_WaveSpace)
958  {
959  temparray = inarray;
960  }
961  else
962  {
963  HomogeneousFwdTrans(inarray,temparray);
964  }
965 
966  NekDouble sign = -1.0;
967  NekDouble beta;
968 
969  //Half Mode
971  {
972  beta = sign*2*M_PI*(m_transposition->GetK(0))/m_lhom;
973 
974  Vmath::Smul(nP_pts,beta,temparray,1,outarray,1);
975  }
976  else if(m_homogeneousBasis->GetBasisType() == LibUtilities::eFourierHalfModeIm)
977  {
978  beta = -sign*2*M_PI*(m_transposition->GetK(0))/m_lhom;
979 
980  Vmath::Smul(nP_pts,beta,temparray,1,outarray,1);
981  }
982 
983  //Fully complex
984  else
985  {
986  for(int i = 0; i < m_planes.num_elements(); i++)
987  {
988  beta = -sign*2*M_PI*(m_transposition->GetK(i))/m_lhom;
989 
990  Vmath::Smul(nP_pts,beta,tmp1 = temparray + i*nP_pts,1,tmp2 = outarray + (i-int(sign))*nP_pts,1);
991 
992  sign = -1.0*sign;
993  }
994  }
995 
996  if(m_WaveSpace)
997  {
998  out_d2 = outarray;
999  }
1000  else
1001  {
1002  HomogeneousBwdTrans(outarray,out_d2);
1003  }
1004  }
1005  else
1006  {
1007  if(!m_session->DefinesSolverInfo("HomoStrip"))
1008  {
1009  ASSERTL0(m_comm->GetColumnComm()->GetSize() == 1,
1010  "Parallelisation in the homogeneous direction "
1011  "implemented just for Fourier basis");
1012  }
1013  else
1014  {
1015  ASSERTL0(m_StripZcomm->GetSize() == 1,
1016  "Parallelisation in the homogeneous direction "
1017  "implemented just for Fourier basis");
1018  }
1019 
1020  if(m_WaveSpace)
1021  {
1022  ASSERTL0(false, "Semi-phyisical time-stepping not "
1023  "implemented yet for non-Fourier "
1024  "basis");
1025  }
1026  else
1027  {
1028  StdRegions::StdSegExp StdSeg(m_homogeneousBasis->GetBasisKey());
1029 
1030  m_transposition->Transpose(inarray,temparray,false,LibUtilities::eXYtoZ);
1031 
1032  for(int i = 0; i < nP_pts; i++)
1033  {
1034  StdSeg.PhysDeriv(temparray + i*m_planes.num_elements(), tmp2 = outarray + i*m_planes.num_elements());
1035  }
1036 
1037  m_transposition->Transpose(outarray,out_d2,false,LibUtilities::eZtoXY);
1038 
1039  Vmath::Smul(nT_pts,2.0/m_lhom,out_d2,1,out_d2,1);
1040  }
1041  }
1042  }
1043  }
void HomogeneousBwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
static Array< OneD, NekDouble > NullNekDouble1DArray
#define sign(a, b)
return the sign(b)*a
Definition: Polylib.cpp:22
LibUtilities::TranspositionSharedPtr m_transposition
NekDouble m_lhom
Width of homogeneous direction.
Fourier Expansion .
Definition: BasisType.h:52
Class representing a segment element in reference space.
Definition: StdSegExp.h:54
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Fourier Modified expansions with just the real part of the first mode .
Definition: BasisType.h:59
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:199
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
double NekDouble
void HomogeneousFwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
Fourier Modified expansions with just the imaginary part of the first mode .
Definition: BasisType.h:60
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:907
Fourier ModifiedExpansion with just the first mode .
Definition: BasisType.h:58
void Nektar::MultiRegions::ExpListHomogeneous1D::v_PhysDeriv ( Direction  edir,
const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  out_d 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1045 of file ExpListHomogeneous1D.cpp.

References ASSERTL0, Nektar::LibUtilities::eFourier, Nektar::LibUtilities::eFourierHalfModeIm, Nektar::LibUtilities::eFourierHalfModeRe, Nektar::LibUtilities::eFourierSingleMode, Nektar::LibUtilities::eXYtoZ, Nektar::LibUtilities::eZtoXY, HomogeneousBwdTrans(), HomogeneousFwdTrans(), Nektar::MultiRegions::ExpList::m_comm, m_homogeneousBasis, m_lhom, m_planes, Nektar::MultiRegions::ExpList::m_session, m_StripZcomm, m_transposition, Nektar::MultiRegions::ExpList::m_WaveSpace, sign, and Vmath::Smul().

1048  {
1049  int nT_pts = inarray.num_elements(); //number of total points = n. of Fourier points * n. of points per plane (nT_pts)
1050  int nP_pts = nT_pts/m_planes.num_elements(); //number of points per plane = n of Fourier transform required (nP_pts)
1051 
1052  int dir= (int)edir;
1053 
1054  Array<OneD, NekDouble> temparray(nT_pts);
1055  Array<OneD, NekDouble> outarray(nT_pts);
1058 
1059  if (dir < 2)
1060  {
1061  for(int i=0; i < m_planes.num_elements(); i++)
1062  {
1063  m_planes[i]->PhysDeriv(edir, inarray + i*nP_pts ,tmp2 = out_d + i*nP_pts);
1064  }
1065  }
1066  else
1067  {
1070  {
1071  if(m_WaveSpace)
1072  {
1073  temparray = inarray;
1074  }
1075  else
1076  {
1077  HomogeneousFwdTrans(inarray,temparray);
1078  }
1079 
1080  NekDouble sign = -1.0;
1081  NekDouble beta;
1082 
1083  //HalfMode
1085  {
1086  beta = 2*sign*M_PI*(m_transposition->GetK(0))/m_lhom;
1087 
1088  Vmath::Smul(nP_pts,beta,temparray,1,outarray,1);
1089  }
1090  else if(m_homogeneousBasis->GetBasisType() == LibUtilities::eFourierHalfModeIm)
1091  {
1092  beta = -2*sign*M_PI*(m_transposition->GetK(0))/m_lhom;
1093 
1094  Vmath::Smul(nP_pts,beta,temparray,1,outarray,1);
1095  }
1096  //Fully complex
1097  else
1098  {
1099  for(int i = 0; i < m_planes.num_elements(); i++)
1100  {
1101  beta = -sign*2*M_PI*(m_transposition->GetK(i))/m_lhom;
1102 
1103  Vmath::Smul(nP_pts,beta,tmp1 = temparray + i*nP_pts,1,tmp2 = outarray + (i-int(sign))*nP_pts,1);
1104 
1105  sign = -1.0*sign;
1106  }
1107  }
1108  if(m_WaveSpace)
1109  {
1110  out_d = outarray;
1111  }
1112  else
1113  {
1114  HomogeneousBwdTrans(outarray,out_d);
1115  }
1116  }
1117  else
1118  {
1119  if(!m_session->DefinesSolverInfo("HomoStrip"))
1120  {
1121  ASSERTL0(m_comm->GetColumnComm()->GetSize() == 1,
1122  "Parallelisation in the homogeneous direction "
1123  "implemented just for Fourier basis");
1124  }
1125  else
1126  {
1127  ASSERTL0(m_StripZcomm->GetSize() == 1,
1128  "Parallelisation in the homogeneous direction "
1129  "implemented just for Fourier basis");
1130  }
1131 
1132  if(m_WaveSpace)
1133  {
1134  ASSERTL0(false,"Semi-phyisical time-stepping not implemented yet for non-Fourier basis");
1135  }
1136  else
1137  {
1138  StdRegions::StdSegExp StdSeg(m_homogeneousBasis->GetBasisKey());
1139 
1140  m_transposition->Transpose(inarray,temparray,false,LibUtilities::eXYtoZ);
1141 
1142  for(int i = 0; i < nP_pts; i++)
1143  {
1144  StdSeg.PhysDeriv(temparray + i*m_planes.num_elements(), tmp2 = outarray + i*m_planes.num_elements());
1145  }
1146 
1147  m_transposition->Transpose(outarray,out_d,false,LibUtilities::eZtoXY);
1148 
1149  Vmath::Smul(nT_pts,2.0/m_lhom,out_d,1,out_d,1);
1150  }
1151  }
1152  }
1153  }
void HomogeneousBwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:161
#define sign(a, b)
return the sign(b)*a
Definition: Polylib.cpp:22
LibUtilities::TranspositionSharedPtr m_transposition
NekDouble m_lhom
Width of homogeneous direction.
Fourier Expansion .
Definition: BasisType.h:52
Class representing a segment element in reference space.
Definition: StdSegExp.h:54
LibUtilities::BasisSharedPtr m_homogeneousBasis
Definition of the total number of degrees of freedom and quadrature points. Sets up the storage for m...
Fourier Modified expansions with just the real part of the first mode .
Definition: BasisType.h:59
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:199
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:910
Array< OneD, ExpListSharedPtr > m_planes
double NekDouble
void HomogeneousFwdTrans(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
Fourier Modified expansions with just the imaginary part of the first mode .
Definition: BasisType.h:60
LibUtilities::CommSharedPtr m_comm
Communicator.
Definition: ExpList.h:907
Fourier ModifiedExpansion with just the first mode .
Definition: BasisType.h:58
void Nektar::MultiRegions::ExpListHomogeneous1D::v_PhysGalerkinProjection1DScaled ( const NekDouble  scale,
const Array< OneD, NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 916 of file ExpListHomogeneous1D.cpp.

References ASSERTL1, and m_planes.

917  {
918  int cnt,cnt1;
919  Array<OneD, NekDouble> tmparray;
920  cnt = m_planes[0]->Get1DScaledTotPoints(scale);
921  cnt1 = m_planes[0]->GetTotPoints();
922 
923  ASSERTL1(m_planes.num_elements()*cnt <= inarray.num_elements(),"size of outarray does not match internal estimage");
924 
925 
926  for(int i = 0; i < m_planes.num_elements(); i++)
927  {
928  m_planes[i]->PhysGalerkinProjection1DScaled(scale,inarray+i*cnt,
929  tmparray = outarray+i*cnt1);
930  }
931 
932  }
Array< OneD, ExpListSharedPtr > m_planes
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
void Nektar::MultiRegions::ExpListHomogeneous1D::v_PhysInterp1DScaled ( const NekDouble  scale,
const Array< OneD, NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 897 of file ExpListHomogeneous1D.cpp.

References ASSERTL1, and m_planes.

898  {
899  int cnt,cnt1;
900  Array<OneD, NekDouble> tmparray;
901  cnt = m_planes[0]->GetTotPoints();
902  cnt1 = m_planes[0]->Get1DScaledTotPoints(scale);
903 
904  ASSERTL1(m_planes.num_elements()*cnt1 <= outarray.num_elements(),"size of outarray does not match internal estimage");
905 
906 
907  for(int i = 0; i < m_planes.num_elements(); i++)
908  {
909 
910  m_planes[i]->PhysInterp1DScaled(scale,inarray+i*cnt,
911  tmparray = outarray+i*cnt1);
912  }
913  }
Array< OneD, ExpListSharedPtr > m_planes
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:191
virtual void Nektar::MultiRegions::ExpListHomogeneous1D::v_SetHomo1DSpecVanVisc ( Array< OneD, NekDouble visc)
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 178 of file ExpListHomogeneous1D.h.

References m_specVanVisc.

179  {
180  m_specVanVisc = visc;
181  }
Array< OneD, NekDouble > m_specVanVisc
Spectral vanishing Viscosity coefficient for stabilisation.
void Nektar::MultiRegions::ExpListHomogeneous1D::v_WriteVtkPieceData ( std::ostream &  outfile,
int  expansion,
std::string  var 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 867 of file ExpListHomogeneous1D.cpp.

References Nektar::NekConstants::kNekZeroTol, Nektar::MultiRegions::ExpList::m_phys, Nektar::MultiRegions::ExpList::m_phys_offset, and m_planes.

869  {
870  // If there is only one plane (e.g. HalfMode), we write a 2D plane.
871  if (m_planes.num_elements() == 1)
872  {
873  m_planes[0]->WriteVtkPieceData(outfile, expansion, var);
874  return;
875  }
876 
877  int i;
878  int nq = (*m_exp)[expansion]->GetTotPoints();
879  int npoints_per_plane = m_planes[0]->GetTotPoints();
880 
881  // printing the fields of that zone
882  outfile << " <DataArray type=\"Float32\" Name=\""
883  << var << "\">" << endl;
884  outfile << " ";
885  for (int n = 0; n < m_planes.num_elements(); ++n)
886  {
887  const Array<OneD, NekDouble> phys = m_phys + m_phys_offset[expansion] + n*npoints_per_plane;
888  for(i = 0; i < nq; ++i)
889  {
890  outfile << (fabs(phys[i]) < NekConstants::kNekZeroTol ? 0 : phys[i]) << " ";
891  }
892  }
893  outfile << endl;
894  outfile << " </DataArray>" << endl;
895  }
Array< OneD, NekDouble > m_phys
The global expansion evaluated at the quadrature points.
Definition: ExpList.h:956
static const NekDouble kNekZeroTol
Array< OneD, int > m_phys_offset
Offset of elemental data into the array m_phys.
Definition: ExpList.h:991
Array< OneD, ExpListSharedPtr > m_planes

Member Data Documentation

bool Nektar::MultiRegions::ExpListHomogeneous1D::m_dealiasing
private

Definition at line 275 of file ExpListHomogeneous1D.h.

Referenced by ExpListHomogeneous1D().

LibUtilities::NektarFFTSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::m_FFT
protected

Definition at line 145 of file ExpListHomogeneous1D.h.

Referenced by ExpListHomogeneous1D(), and Homogeneous1DTrans().

LibUtilities::NektarFFTSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::m_FFT_deal
protected

Definition at line 147 of file ExpListHomogeneous1D.h.

Referenced by ExpListHomogeneous1D(), and v_DealiasedProd().

Homo1DBlockMatrixMapShPtr Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneous1DBlockMat
protected

Definition at line 157 of file ExpListHomogeneous1D.h.

Referenced by GetHomogeneous1DBlockMatrix().

LibUtilities::BasisSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::m_homogeneousBasis
protected
NekDouble Nektar::MultiRegions::ExpListHomogeneous1D::m_lhom
protected
int Nektar::MultiRegions::ExpListHomogeneous1D::m_padsize
private

Definition at line 276 of file ExpListHomogeneous1D.h.

Referenced by ExpListHomogeneous1D(), and v_DealiasedProd().

Array<OneD, ExpListSharedPtr> Nektar::MultiRegions::ExpListHomogeneous1D::m_planes
protected

Definition at line 158 of file ExpListHomogeneous1D.h.

Referenced by Nektar::MultiRegions::ContField3DHomogeneous1D::ContField3DHomogeneous1D(), Nektar::MultiRegions::DisContField3DHomogeneous1D::DisContField3DHomogeneous1D(), Nektar::MultiRegions::DisContField3DHomogeneous1D::EvaluateBoundaryConditions(), Nektar::MultiRegions::ExpList2DHomogeneous1D::ExpList2DHomogeneous1D(), Nektar::MultiRegions::ExpList3DHomogeneous1D::ExpList3DHomogeneous1D(), ExpListHomogeneous1D(), Nektar::MultiRegions::ExpList3DHomogeneous1D::GenExpList3DHomogeneous1D(), GenHomogeneous1DBlockMatrix(), Nektar::MultiRegions::DisContField3DHomogeneous1D::GetBCValues(), Nektar::MultiRegions::DisContField3DHomogeneous1D::GetBoundaryToElmtMap(), Nektar::MultiRegions::ExpList2DHomogeneous1D::GetCoords(), Nektar::MultiRegions::ExpList3DHomogeneous1D::GetCoords(), GetPlane(), Homogeneous1DTrans(), Nektar::MultiRegions::DisContField3DHomogeneous1D::NormVectorIProductWRTBase(), Nektar::MultiRegions::ExpList2DHomogeneous1D::SetCoeffPhys(), Nektar::MultiRegions::ExpList3DHomogeneous1D::SetCoeffPhys(), Nektar::MultiRegions::DisContField3DHomogeneous1D::SetupBoundaryConditions(), Nektar::MultiRegions::DisContField3DHomogeneous1D::SetUpDG(), v_AppendFieldData(), v_BwdTrans(), v_BwdTrans_IterPerExp(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_ClearGlobalLinSysManager(), v_DealiasedProd(), v_ExtractCoeffsToCoeffs(), v_ExtractDataToCoeffs(), Nektar::MultiRegions::DisContField3DHomogeneous1D::v_ExtractTracePhys(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_FillBndCondFromField(), v_FwdTrans(), v_FwdTrans_IterPerExp(), Nektar::MultiRegions::ExpList2DHomogeneous1D::v_GetCoords(), Nektar::MultiRegions::ExpList3DHomogeneous1D::v_GetCoords(), v_GetFieldDefinitions(), Nektar::MultiRegions::ExpList2DHomogeneous1D::v_GetNormals(), Nektar::MultiRegions::ExpList3DHomogeneous1D::v_GetPeriodicEntities(), Nektar::MultiRegions::DisContField3DHomogeneous1D::v_GetTraceMap(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_GlobalToLocal(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_HelmSolve(), Nektar::MultiRegions::DisContField3DHomogeneous1D::v_HelmSolve(), Nektar::MultiRegions::ExpList3DHomogeneous1D::v_HomogeneousEnergy(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_ImposeDirichletConditions(), v_IProductWRTBase(), v_IProductWRTBase_IterPerExp(), Nektar::MultiRegions::ExpList3DHomogeneous1D::v_L2(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_LocalToGlobal(), v_PhysDeriv(), v_PhysGalerkinProjection1DScaled(), v_PhysInterp1DScaled(), Nektar::MultiRegions::ContField3DHomogeneous1D::v_SmoothField(), Nektar::MultiRegions::ExpList3DHomogeneous1D::v_WriteTecplotConnectivity(), Nektar::MultiRegions::ExpList2DHomogeneous1D::v_WriteTecplotZone(), v_WriteVtkPieceData(), Nektar::MultiRegions::ExpList2DHomogeneous1D::v_WriteVtkPieceHeader(), and Nektar::MultiRegions::ExpList3DHomogeneous1D::v_WriteVtkPieceHeader().

Array<OneD, NekDouble> Nektar::MultiRegions::ExpListHomogeneous1D::m_specVanVisc
private

Spectral vanishing Viscosity coefficient for stabilisation.

Definition at line 279 of file ExpListHomogeneous1D.h.

Referenced by GetSpecVanVisc(), and v_SetHomo1DSpecVanVisc().

LibUtilities::CommSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::m_StripZcomm
Array<OneD,NekDouble> Nektar::MultiRegions::ExpListHomogeneous1D::m_tmpIN
protected

Definition at line 149 of file ExpListHomogeneous1D.h.

Referenced by Homogeneous1DTrans().

Array<OneD,NekDouble> Nektar::MultiRegions::ExpListHomogeneous1D::m_tmpOUT
protected

Definition at line 150 of file ExpListHomogeneous1D.h.

Referenced by Homogeneous1DTrans().

LibUtilities::TranspositionSharedPtr Nektar::MultiRegions::ExpListHomogeneous1D::m_transposition
bool Nektar::MultiRegions::ExpListHomogeneous1D::m_useFFT
protected