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

This class is the abstraction of a global discontinuous two- dimensional spectral/hp element expansion which approximates the solution of a set of partial differential equations. More...

#include <DisContField1D.h>

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

Public Member Functions

 DisContField1D ()
 Default constructor. More...
 
 DisContField1D (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &graph1D, const std::string &variable, const bool SetUpJustDG=true)
 Constructs a 1D discontinuous field based on a mesh and boundary conditions. More...
 
 DisContField1D (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::CompositeMap &domain, const SpatialDomains::BoundaryConditions &Allbcs, const std::string &variable, bool SetToOneSpaceDimensions=false)
 Constructor for a DisContField1D from a List of subdomains New Constructor for arterial network. More...
 
 DisContField1D (const DisContField1D &In)
 Constructs a 1D discontinuous field based on an existing field. More...
 
 DisContField1D (const ExpList1D &In)
 Constructs a 1D discontinuous field based on an existing field. (needed in order to use ContField( const ExpList1D &In) constructor. More...
 
virtual ~DisContField1D ()
 Destructor. More...
 
GlobalLinSysSharedPtr GetGlobalBndLinSys (const GlobalLinSysKey &mkey)
 For a given key, returns the associated global linear system. More...
 
std::vector< bool > & GetNegatedFluxNormal (void)
 
- Public Member Functions inherited from Nektar::MultiRegions::ExpList1D
 ExpList1D ()
 The default constructor. More...
 
 ExpList1D (const ExpList1D &In, const bool DeclareCoeffPhysArrays=true)
 The copy constructor. More...
 
 ExpList1D (const ExpList1D &In, const std::vector< unsigned int > &eIDs, const bool DeclareCoeffPhysArrays=true)
 Constructor copying only elements defined in eIds. More...
 
 ExpList1D (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::BasisKey &Ba, const SpatialDomains::MeshGraphSharedPtr &graph1D)
 Construct an ExpList1D from a given graph. More...
 
 ExpList1D (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &graph1D, const bool DeclareCoeffPhysArrays=true)
 This constructor sets up a list of local expansions based on an input graph1D. More...
 
 ExpList1D (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::CompositeMap &domain, const bool DeclareCoeffPhysArrays=true, const std::string var="DefaultVar", bool SetToOneSpaceDimension=false)
 This constructor sets up a list of local expansions based on an input compositeMap. More...
 
 ExpList1D (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::CompositeMap &domain, const SpatialDomains::MeshGraphSharedPtr &graph2D, const bool DeclareCoeffPhysArrays=true, const std::string variable="DefaultVar")
 Specialised constructor for Neumann boundary conditions in DisContField2D and ContField2D. More...
 
 ExpList1D (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::CompositeMap &domain, const SpatialDomains::MeshGraphSharedPtr &graph1D, int i, const bool DeclareCoeffPhysArrays=true)
 
 ExpList1D (const LibUtilities::SessionReaderSharedPtr &pSession, const Array< OneD, const ExpListSharedPtr > &bndConstraint, const Array< OneD, const SpatialDomains::BoundaryConditionShPtr > &bndCond, const LocalRegions::ExpansionVector &locexp, const SpatialDomains::MeshGraphSharedPtr &graph2D, const PeriodicMap &periodicEdges, const bool DeclareCoeffPhysArrays=true, const std::string variable="DefaultVar")
 Specialised constructor for trace expansions. More...
 
virtual ~ExpList1D ()
 Destructor. More...
 
void PostProcess (LibUtilities::KernelSharedPtr kernel, Array< OneD, NekDouble > &inarray, Array< OneD, NekDouble > &outarray, NekDouble h, int elmId=0)
 Performs the post-processing on a specified element. More...
 
void PeriodicEval (Array< OneD, NekDouble > &inarray1, Array< OneD, NekDouble > &inarray2, NekDouble h, int nmodes, Array< OneD, NekDouble > &outarray)
 Evaluates the global spectral/hp expansion at some arbitray set of points. More...
 
void ParNormalSign (Array< OneD, NekDouble > &normsign)
 Set up the normals on each expansion. More...
 
- Public Member Functions inherited from Nektar::MultiRegions::ExpList
 ExpList ()
 The default constructor. More...
 
 ExpList (const LibUtilities::SessionReaderSharedPtr &pSession)
 The default constructor. More...
 
 ExpList (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &pGraph)
 The default constructor. More...
 
 ExpList (const ExpList &in, const std::vector< unsigned int > &eIDs, const bool DeclareCoeffPhysArrays=true)
 Constructor copying only elements defined in eIds. More...
 
 ExpList (const ExpList &in, const bool DeclareCoeffPhysArrays=true)
 The copy constructor. More...
 
virtual ~ExpList ()
 The default destructor. More...
 
int GetNcoeffs (void) const
 Returns the total number of local degrees of freedom $N_{\mathrm{eof}}=\sum_{e=1}^{{N_{\mathrm{el}}}}N^{e}_m$. More...
 
int GetNcoeffs (const int eid) const
 Returns the total number of local degrees of freedom for element eid. More...
 
ExpansionType GetExpType (void)
 Returns the type of the expansion. More...
 
void SetExpType (ExpansionType Type)
 Returns the type of the expansion. More...
 
int EvalBasisNumModesMax (void) const
 Evaulates the maximum number of modes in the elemental basis order over all elements. More...
 
const Array< OneD, int > EvalBasisNumModesMaxPerExp (void) const
 Returns the vector of the number of modes in the elemental basis order over all elements. More...
 
int GetTotPoints (void) const
 Returns the total number of quadrature points m_npoints $=Q_{\mathrm{tot}}$. More...
 
int GetTotPoints (const int eid) const
 Returns the total number of quadrature points for eid's element $=Q_{\mathrm{tot}}$. More...
 
int GetNpoints (void) const
 Returns the total number of quadrature points m_npoints $=Q_{\mathrm{tot}}$. More...
 
int Get1DScaledTotPoints (const NekDouble scale) const
 Returns the total number of qudature points scaled by the factor scale on each 1D direction. More...
 
void SetWaveSpace (const bool wavespace)
 Sets the wave space to the one of the possible configuration true or false. More...
 
void SetModifiedBasis (const bool modbasis)
 Set Modified Basis for the stability analysis. More...
 
void SetPhys (int i, NekDouble val)
 Set the i th value of m_phys to value val. More...
 
bool GetWaveSpace (void) const
 This function returns the third direction expansion condition, which can be in wave space (coefficient) or not It is stored in the variable m_WaveSpace. More...
 
void SetPhys (const Array< OneD, const NekDouble > &inarray)
 Fills the array m_phys. More...
 
void SetPhysArray (Array< OneD, NekDouble > &inarray)
 Sets the array m_phys. More...
 
void SetPhysState (const bool physState)
 This function manually sets whether the array of physical values $\boldsymbol{u}_l$ (implemented as m_phys) is filled or not. More...
 
bool GetPhysState (void) const
 This function indicates whether the array of physical values $\boldsymbol{u}_l$ (implemented as m_phys) is filled or not. More...
 
NekDouble PhysIntegral (void)
 This function integrates a function $f(\boldsymbol{x})$ over the domain consisting of all the elements of the expansion. More...
 
NekDouble PhysIntegral (const Array< OneD, const NekDouble > &inarray)
 This function integrates a function $f(\boldsymbol{x})$ over the domain consisting of all the elements of the expansion. More...
 
void IProductWRTBase_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function calculates the inner product of a function $f(\boldsymbol{x})$ with respect to all {local} expansion modes $\phi_n^e(\boldsymbol{x})$. More...
 
void IProductWRTBase (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void IProductWRTDerivBase (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function calculates the inner product of a function $f(\boldsymbol{x})$ with respect to the derivative (in direction. More...
 
void IProductWRTDerivBase (const Array< OneD, const Array< OneD, NekDouble > > &inarray, Array< OneD, NekDouble > &outarray)
 This function calculates the inner product of a function $f(\boldsymbol{x})$ with respect to the derivative (in direction. More...
 
void FwdTrans_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function elementally evaluates the forward transformation of a function $u(\boldsymbol{x})$ onto the global spectral/hp expansion. More...
 
void FwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void MultiplyByElmtInvMass (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function elementally mulplies the coefficient space of Sin my the elemental inverse of the mass matrix. More...
 
void MultiplyByInvMassMatrix (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void SmoothField (Array< OneD, NekDouble > &field)
 Smooth a field across elements. More...
 
void HelmSolve (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const FlagList &flags, const StdRegions::ConstFactorMap &factors, const StdRegions::VarCoeffMap &varcoeff=StdRegions::NullVarCoeffMap, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray, const bool PhysSpaceForcing=true)
 Solve helmholtz problem. More...
 
void LinearAdvectionDiffusionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 Solve Advection Diffusion Reaction. More...
 
void LinearAdvectionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 Solve Advection Diffusion Reaction. More...
 
void FwdTrans_BndConstrained (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void BwdTrans_IterPerExp (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This function elementally evaluates the backward transformation of the global spectral/hp element expansion. More...
 
void BwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void GetCoords (Array< OneD, NekDouble > &coord_0, Array< OneD, NekDouble > &coord_1=NullNekDouble1DArray, Array< OneD, NekDouble > &coord_2=NullNekDouble1DArray)
 This function calculates the coordinates of all the elemental quadrature points $\boldsymbol{x}_i$. More...
 
void HomogeneousFwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
void HomogeneousBwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
void DealiasedProd (const Array< OneD, NekDouble > &inarray1, const Array< OneD, NekDouble > &inarray2, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
void DealiasedDotProd (const Array< OneD, Array< OneD, NekDouble > > &inarray1, const Array< OneD, Array< OneD, NekDouble > > &inarray2, Array< OneD, Array< OneD, NekDouble > > &outarray, CoeffState coeffstate=eLocal)
 
void GetBCValues (Array< OneD, NekDouble > &BndVals, const Array< OneD, NekDouble > &TotField, int BndID)
 
void NormVectorIProductWRTBase (Array< OneD, const NekDouble > &V1, Array< OneD, const NekDouble > &V2, Array< OneD, NekDouble > &outarray, int BndID)
 
void NormVectorIProductWRTBase (Array< OneD, Array< OneD, NekDouble > > &V, Array< OneD, NekDouble > &outarray)
 
void ApplyGeomInfo ()
 Apply geometry information to each expansion. More...
 
void Reset ()
 Reset geometry information and reset matrices. More...
 
void WriteTecplotHeader (std::ostream &outfile, std::string var="")
 
void WriteTecplotZone (std::ostream &outfile, int expansion=-1)
 
void WriteTecplotField (std::ostream &outfile, int expansion=-1)
 
void WriteTecplotConnectivity (std::ostream &outfile, int expansion=-1)
 
void WriteVtkHeader (std::ostream &outfile)
 
void WriteVtkFooter (std::ostream &outfile)
 
void WriteVtkPieceHeader (std::ostream &outfile, int expansion, int istrip=0)
 
void WriteVtkPieceFooter (std::ostream &outfile, int expansion)
 
void WriteVtkPieceData (std::ostream &outfile, int expansion, std::string var="v")
 
int GetCoordim (int eid)
 This function returns the dimension of the coordinates of the element eid. More...
 
void SetCoeff (int i, NekDouble val)
 Set the i th coefficiient in m_coeffs to value val. More...
 
void SetCoeffs (int i, NekDouble val)
 Set the i th coefficiient in m_coeffs to value val. More...
 
void SetCoeffsArray (Array< OneD, NekDouble > &inarray)
 Set the m_coeffs array to inarray. More...
 
const Array< OneD, const
NekDouble > & 
GetCoeffs () const
 This function returns (a reference to) the array $\boldsymbol{\hat{u}}_l$ (implemented as m_coeffs) containing all local expansion coefficients. More...
 
void ImposeDirichletConditions (Array< OneD, NekDouble > &outarray)
 Impose Dirichlet Boundary Conditions onto Array. More...
 
void FillBndCondFromField (void)
 Fill Bnd Condition expansion from the values stored in expansion. More...
 
void FillBndCondFromField (const int nreg)
 Fill Bnd Condition expansion in nreg from the values stored in expansion. More...
 
void LocalToGlobal (bool useComm=true)
 Gathers the global coefficients $\boldsymbol{\hat{u}}_g$ from the local coefficients $\boldsymbol{\hat{u}}_l$. More...
 
void LocalToGlobal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool useComm=true)
 
void GlobalToLocal (void)
 Scatters from the global coefficients $\boldsymbol{\hat{u}}_g$ to the local coefficients $\boldsymbol{\hat{u}}_l$. More...
 
void GlobalToLocal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
NekDouble GetCoeff (int i)
 Get the i th value (coefficient) of m_coeffs. More...
 
NekDouble GetCoeffs (int i)
 Get the i th value (coefficient) of m_coeffs. More...
 
const Array< OneD, const
NekDouble > & 
GetPhys () const
 This function returns (a reference to) the array $\boldsymbol{u}_l$ (implemented as m_phys) containing the function $u^{\delta}(\boldsymbol{x})$ evaluated at the quadrature points. More...
 
NekDouble Linf (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 This function calculates the $L_\infty$ error of the global spectral/hp element approximation. More...
 
NekDouble L2 (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 This function calculates the $L_2$ error with respect to soln of the global spectral/hp element approximation. More...
 
NekDouble H1 (const Array< OneD, const NekDouble > &inarray, const Array< OneD, const NekDouble > &soln=NullNekDouble1DArray)
 Calculates the $H^1$ error of the global spectral/hp element approximation. More...
 
NekDouble Integral (const Array< OneD, const NekDouble > &inarray)
 
Array< OneD, const 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)
 
void CurlCurl (Array< OneD, Array< OneD, NekDouble > > &Vel, Array< OneD, Array< OneD, NekDouble > > &Q)
 
const Array< OneD, const
boost::shared_ptr< ExpList > > & 
GetBndCondExpansions ()
 
boost::shared_ptr< ExpList > & UpdateBndCondExpansion (int i)
 
void Upwind (const Array< OneD, const Array< OneD, NekDouble > > &Vec, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &Upwind)
 
void Upwind (const Array< OneD, const NekDouble > &Vn, const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &Upwind)
 
boost::shared_ptr< ExpList > & GetTrace ()
 
boost::shared_ptr
< AssemblyMapDG > & 
GetTraceMap (void)
 
const Array< OneD, const int > & GetTraceBndMap (void)
 
void GetNormals (Array< OneD, Array< OneD, NekDouble > > &normals)
 
void AddTraceIntegral (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, Array< OneD, NekDouble > &outarray)
 
void AddTraceIntegral (const Array< OneD, const NekDouble > &Fn, Array< OneD, NekDouble > &outarray)
 
void AddFwdBwdTraceIntegral (const Array< OneD, const NekDouble > &Fwd, const Array< OneD, const NekDouble > &Bwd, Array< OneD, NekDouble > &outarray)
 
void GetFwdBwdTracePhys (Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
 
void GetFwdBwdTracePhys (const Array< OneD, const NekDouble > &field, Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
 
const std::vector< bool > & GetLeftAdjacentFaces (void) const
 
void ExtractTracePhys (Array< OneD, NekDouble > &outarray)
 
void ExtractTracePhys (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
const Array< OneD, const
SpatialDomains::BoundaryConditionShPtr > & 
GetBndConditions ()
 
Array< OneD,
SpatialDomains::BoundaryConditionShPtr > & 
UpdateBndConditions ()
 
void EvaluateBoundaryConditions (const NekDouble time=0.0, const std::string varName="", const NekDouble=NekConstants::kNekUnsetDouble, const NekDouble=NekConstants::kNekUnsetDouble)
 
void GeneralMatrixOp (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 This function calculates the result of the multiplication of a matrix of type specified by mkey with a vector given by inarray. More...
 
void GeneralMatrixOp_IterPerExp (const GlobalMatrixKey &gkey, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
void SetUpPhysNormals ()
 
void GetBoundaryToElmtMap (Array< OneD, int > &ElmtID, Array< OneD, int > &EdgeID)
 
void GetBndElmtExpansion (int i, boost::shared_ptr< ExpList > &result, const bool DeclareCoeffPhysArrays=true)
 
void ExtractElmtToBndPhys (int i, Array< OneD, NekDouble > &elmt, Array< OneD, NekDouble > &boundary)
 
void ExtractPhysToBndElmt (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bndElmt)
 
void ExtractPhysToBnd (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bnd)
 
void GetBoundaryNormals (int i, Array< OneD, Array< OneD, NekDouble > > &normals)
 
void GeneralGetFieldDefinitions (std::vector< LibUtilities::FieldDefinitionsSharedPtr > &fielddef, int NumHomoDir=0, Array< OneD, LibUtilities::BasisSharedPtr > &HomoBasis=LibUtilities::NullBasisSharedPtr1DArray, std::vector< NekDouble > &HomoLen=LibUtilities::NullNekDoubleVector, bool homoStrips=false, std::vector< unsigned int > &HomoSIDs=LibUtilities::NullUnsignedIntVector, std::vector< unsigned int > &HomoZIDs=LibUtilities::NullUnsignedIntVector, std::vector< unsigned int > &HomoYIDs=LibUtilities::NullUnsignedIntVector)
 
const
NekOptimize::GlobalOptParamSharedPtr
GetGlobalOptParam (void)
 
std::map< int,
RobinBCInfoSharedPtr
GetRobinBCInfo ()
 
void GetPeriodicEntities (PeriodicMap &periodicVerts, PeriodicMap &periodicEdges, PeriodicMap &periodicFaces=NullPeriodicMap)
 
std::vector
< LibUtilities::FieldDefinitionsSharedPtr
GetFieldDefinitions ()
 
void GetFieldDefinitions (std::vector< LibUtilities::FieldDefinitionsSharedPtr > &fielddef)
 
void AppendFieldData (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata)
 Append the element data listed in elements fielddef->m_ElementIDs onto fielddata. More...
 
void AppendFieldData (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, Array< OneD, NekDouble > &coeffs)
 Append the data in coeffs listed in elements fielddef->m_ElementIDs onto fielddata. More...
 
void ExtractElmtDataToCoeffs (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, std::string &field, Array< OneD, NekDouble > &coeffs)
 Extract the data in fielddata into the coeffs using the basic ExpList Elemental expansions rather than planes in homogeneous case. More...
 
void ExtractCoeffsToCoeffs (const boost::shared_ptr< ExpList > &fromExpList, const Array< OneD, const NekDouble > &fromCoeffs, Array< OneD, NekDouble > &toCoeffs)
 Extract the data from fromField using fromExpList the coeffs using the basic ExpList Elemental expansions rather than planes in homogeneous case. More...
 
void ExtractDataToCoeffs (LibUtilities::FieldDefinitionsSharedPtr &fielddef, std::vector< NekDouble > &fielddata, std::string &field, Array< OneD, NekDouble > &coeffs)
 Extract the data in fielddata into the coeffs. More...
 
boost::shared_ptr< ExpListGetSharedThisPtr ()
 Returns a shared pointer to the current object. More...
 
boost::shared_ptr
< LibUtilities::SessionReader
GetSession () const
 Returns the session object. More...
 
boost::shared_ptr
< LibUtilities::Comm
GetComm ()
 Returns the comm object. More...
 
SpatialDomains::MeshGraphSharedPtr GetGraph ()
 
LibUtilities::BasisSharedPtr GetHomogeneousBasis (void)
 
boost::shared_ptr< ExpList > & GetPlane (int n)
 
void CreateCollections (Collections::ImplementationType ImpType=Collections::eNoImpType)
 Construct collections of elements containing a single element type and polynomial order from the list of expansions. More...
 
void ClearGlobalLinSysManager (void)
 

Protected Member Functions

void GenerateBoundaryConditionExpansion (const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::BoundaryConditions &bcs, const std::string variable)
 Discretises the boundary conditions. More...
 
void FindPeriodicVertices (const SpatialDomains::BoundaryConditions &bcs, const std::string variable)
 Generate a associative map of periodic vertices in a mesh. More...
 
virtual ExpListSharedPtrv_GetTrace ()
 
virtual AssemblyMapDGSharedPtrv_GetTraceMap (void)
 
virtual void v_AddTraceIntegral (const Array< OneD, const NekDouble > &Fn, 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)
 This method extracts the "forward" and "backward" trace data from the array field and puts the data into output vectors Fwd and Bwd. More...
 
virtual void v_ExtractTracePhys (Array< OneD, NekDouble > &outarray)
 
virtual void v_ExtractTracePhys (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 This method extracts the trace (verts in 1D) from the field inarray and puts the values in outarray. More...
 
void SetBoundaryConditionExpansion (const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::BoundaryConditions &bcs, const std::string variable, Array< OneD, MultiRegions::ExpListSharedPtr > &bndCondExpansions, Array< OneD, SpatialDomains::BoundaryConditionShPtr > &bndConditions)
 Populates the list of boundary condition expansions. More...
 
void SetMultiDomainBoundaryConditionExpansion (const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::BoundaryConditions &bcs, const std::string variable, Array< OneD, MultiRegions::ExpListSharedPtr > &bndCondExpansions, Array< OneD, SpatialDomains::BoundaryConditionShPtr > &bndConditions, int subdomain)
 Populates the list of boundary condition expansions in multidomain case. More...
 
void GenerateFieldBnd1D (SpatialDomains::BoundaryConditions &bcs, const std::string variable)
 
virtual std::map< int,
RobinBCInfoSharedPtr
v_GetRobinBCInfo ()
 
virtual const Array< OneD,
const
MultiRegions::ExpListSharedPtr > & 
v_GetBndCondExpansions ()
 
virtual const Array< OneD,
const
SpatialDomains::BoundaryConditionShPtr > & 
v_GetBndConditions ()
 
virtual
MultiRegions::ExpListSharedPtr
v_UpdateBndCondExpansion (int i)
 
virtual Array< OneD,
SpatialDomains::BoundaryConditionShPtr > & 
v_UpdateBndConditions ()
 
virtual void v_GetBoundaryToElmtMap (Array< OneD, int > &ElmtID, Array< OneD, int > &VertID)
 
virtual void v_GetBndElmtExpansion (int i, boost::shared_ptr< ExpList > &result, const bool DeclareCoeffPhysArrays)
 
virtual void v_Reset ()
 Reset this field, so that geometry information can be updated. More...
 
virtual void v_EvaluateBoundaryConditions (const NekDouble time=0.0, const std::string varName="", const NekDouble x2_in=NekConstants::kNekUnsetDouble, const NekDouble x3_in=NekConstants::kNekUnsetDouble)
 Evaluate all boundary conditions at a given time.. More...
 
virtual void v_HelmSolve (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const FlagList &flags, const StdRegions::ConstFactorMap &factors, const StdRegions::VarCoeffMap &varcoeff, const Array< OneD, const NekDouble > &dirForcing, const bool PhysSpaceForcing)
 Solve the Helmholtz equation. More...
 
- Protected Member Functions inherited from Nektar::MultiRegions::ExpList1D
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)
 Upwind the Fwd and Bwd states based on the velocity field given by Vec. More...
 
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)
 Upwind the Fwd and Bwd states based on the one- dimensional normal velocity field given by Vn. More...
 
void v_GetNormals (Array< OneD, Array< OneD, NekDouble > > &normals)
 Populate normals with the normals of all expansions. More...
 
- 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 int v_GetNumElmts (void)
 
virtual const Array< OneD,
const int > & 
v_GetTraceBndMap ()
 
virtual void v_AddTraceIntegral (const Array< OneD, const NekDouble > &Fx, const Array< OneD, const NekDouble > &Fy, 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 const std::vector< bool > & v_GetLeftAdjacentFaces (void) const
 
virtual void v_MultiplyByInvMassMatrix (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate)
 
virtual void v_LinearAdvectionDiffusionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 
virtual void v_LinearAdvectionReactionSolve (const Array< OneD, Array< OneD, NekDouble > > &velocity, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const NekDouble lambda, CoeffState coeffstate=eLocal, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
 
virtual void v_ImposeDirichletConditions (Array< OneD, NekDouble > &outarray)
 
virtual void v_FillBndCondFromField ()
 
virtual void v_FillBndCondFromField (const int nreg)
 
virtual void v_LocalToGlobal (bool UseComm)
 
virtual void v_LocalToGlobal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, bool UseComm)
 
virtual void v_GlobalToLocal (void)
 
virtual void v_GlobalToLocal (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray)
 
virtual void v_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_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_SmoothField (Array< OneD, NekDouble > &field)
 
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 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 Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d0, Array< OneD, NekDouble > &out_d1, Array< OneD, NekDouble > &out_d2)
 
virtual void v_PhysDeriv (const int dir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
virtual void v_PhysDeriv (Direction edir, const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &out_d)
 
virtual void v_CurlCurl (Array< OneD, Array< OneD, NekDouble > > &Vel, Array< OneD, Array< OneD, NekDouble > > &Q)
 
virtual void v_HomogeneousFwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
virtual void v_HomogeneousBwdTrans (const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal, bool Shuff=true, bool UnShuff=true)
 
virtual void v_DealiasedProd (const Array< OneD, NekDouble > &inarray1, const Array< OneD, NekDouble > &inarray2, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
 
virtual void v_DealiasedDotProd (const Array< OneD, Array< OneD, NekDouble > > &inarray1, const Array< OneD, Array< OneD, NekDouble > > &inarray2, Array< OneD, Array< OneD, NekDouble > > &outarray, CoeffState coeffstate=eLocal)
 
virtual void v_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_ExtractElmtToBndPhys (int i, Array< OneD, NekDouble > &elmt, Array< OneD, NekDouble > &boundary)
 
virtual void v_ExtractPhysToBndElmt (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bndElmt)
 
virtual void v_ExtractPhysToBnd (int i, const Array< OneD, const NekDouble > &phys, Array< OneD, NekDouble > &bnd)
 
virtual void v_GetBoundaryNormals (int i, Array< OneD, Array< OneD, NekDouble > > &normals)
 
virtual 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_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_WriteVtkPieceData (std::ostream &outfile, int expansion, std::string var)
 
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
LibUtilities::TranspositionSharedPtr 
v_GetTransposition (void)
 
virtual NekDouble v_GetHomoLen (void)
 
virtual Array< OneD, const
unsigned int > 
v_GetZIDs (void)
 
virtual Array< OneD, const
unsigned int > 
v_GetYIDs (void)
 
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_ClearGlobalLinSysManager (void)
 
void ExtractFileBCs (const std::string &fileName, LibUtilities::CommSharedPtr comm, const std::string &varName, const boost::shared_ptr< ExpList > locExpList)
 

Protected Attributes

int m_numDirBndCondExpansions
 The number of boundary segments on which Dirichlet boundary conditions are imposed. More...
 
Array< OneD,
MultiRegions::ExpListSharedPtr
m_bndCondExpansions
 Discretised boundary conditions. More...
 
Array< OneD,
SpatialDomains::BoundaryConditionShPtr
m_bndConditions
 An array which contains the information about the boundary condition on the different boundary regions. More...
 
GlobalLinSysMapShPtr m_globalBndMat
 Global boundary matrix. More...
 
ExpListSharedPtr m_trace
 Trace space storage for points between elements. More...
 
AssemblyMapDGSharedPtr m_traceMap
 Local to global DG mapping for trace space. More...
 
std::set< int > m_boundaryVerts
 A set storing the global IDs of any boundary edges. More...
 
PeriodicMap m_periodicVerts
 A map which identifies groups of periodic vertices. More...
 
std::vector< int > m_periodicFwdCopy
 A vector indicating degress of freedom which need to be copied from forwards to backwards space in case of a periodic boundary condition. More...
 
std::vector< int > m_periodicBwdCopy
 
std::vector< bool > m_leftAdjacentVerts
 
- 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
 
boost::unordered_map< int, int > m_elmtToExpId
 Mapping from geometry ID of element to index inside m_exp. More...
 

Private Member Functions

void SetUpDG (const std::string &variable)
 
bool IsLeftAdjacentVertex (const int n, const int e)
 
SpatialDomains::BoundaryConditionsSharedPtr GetDomainBCs (const SpatialDomains::CompositeMap &domain, const SpatialDomains::BoundaryConditions &Allbcs, const std::string &variable)
 

Private Attributes

std::vector< bool > m_negatedFluxNormal
 

Additional Inherited Members

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

Detailed Description

This class is the abstraction of a global discontinuous two- dimensional spectral/hp element expansion which approximates the solution of a set of partial differential equations.

This class augments the list of local expansions inherited from ExpList1D with boundary conditions. Inter-element boundaries are handled using an discontinuous Galerkin scheme.

Definition at line 55 of file DisContField1D.h.

Constructor & Destructor Documentation

Nektar::MultiRegions::DisContField1D::DisContField1D ( )

Default constructor.

Constructs an empty expansion list with no boundary conditions.

Definition at line 58 of file DisContField1D.cpp.

58  :
59  ExpList1D(),
62  {
63  }
ExpList1D()
The default constructor.
Definition: ExpList1D.cpp:77
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
Nektar::MultiRegions::DisContField1D::DisContField1D ( const LibUtilities::SessionReaderSharedPtr pSession,
const SpatialDomains::MeshGraphSharedPtr graph1D,
const std::string &  variable,
const bool  SetUpJustDG = true 
)

Constructs a 1D discontinuous field based on a mesh and boundary conditions.

An expansion list for the boundary expansions is generated first for the field. These are subsequently evaluated for time zero. The trace map is then constructed.

Parameters
graph1DA mesh, containing information about the domain and the spectral/hp element expansions.
bcsInformation about the enforced boundary conditions.
variableThe session variable associated with the boundary conditions to enforce.
solnTypeType of global system to use.

Definition at line 78 of file DisContField1D.cpp.

References Nektar::MultiRegions::ExpList::ApplyGeomInfo(), Nektar::MultiRegions::ExpList::EvaluateBoundaryConditions(), FindPeriodicVertices(), GenerateBoundaryConditionExpansion(), Nektar::MultiRegions::ExpList::GetBoundaryToElmtMap(), m_bndCondExpansions, Nektar::MultiRegions::ExpList::m_session, SetUpDG(), and Nektar::MultiRegions::ExpList::SetUpPhysNormals().

83  : ExpList1D(pSession,graph1D),
86  {
87  if (variable.compare("DefaultVar") != 0)
88  {
89  SpatialDomains::BoundaryConditions bcs(m_session, graph1D);
90 
91  GenerateBoundaryConditionExpansion(graph1D,bcs,variable);
92  EvaluateBoundaryConditions(0.0, variable);
93  ApplyGeomInfo();
94  FindPeriodicVertices(bcs,variable);
95  }
96 
97  if(SetUpJustDG)
98  {
99  SetUpDG(variable);
100  }
101  else
102  {
103  int i;
104  Array<OneD, int> ElmtID, VertexID;
105  GetBoundaryToElmtMap(ElmtID, VertexID);
106 
107  for(i = 0; i < m_bndCondExpansions.num_elements(); ++i)
108  {
110  locExpList = m_bndCondExpansions[i];
111 
113  (*m_exp)[ElmtID[i]]->
114  as<LocalRegions::Expansion1D>();
116  locExpList->GetExp(0)->
117  as<LocalRegions::Expansion0D>();
118 
119  exp0d->SetAdjacentElementExp(VertexID[i], exp1d);
120  }
121 
123  }
124 
125  }
ExpList1D()
The default constructor.
Definition: ExpList1D.cpp:77
void EvaluateBoundaryConditions(const NekDouble time=0.0, const std::string varName="", const NekDouble=NekConstants::kNekUnsetDouble, const NekDouble=NekConstants::kNekUnsetDouble)
Definition: ExpList.h:2249
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
void FindPeriodicVertices(const SpatialDomains::BoundaryConditions &bcs, const std::string variable)
Generate a associative map of periodic vertices in a mesh.
void SetUpDG(const std::string &variable)
void GetBoundaryToElmtMap(Array< OneD, int > &ElmtID, Array< OneD, int > &EdgeID)
Definition: ExpList.h:2302
void ApplyGeomInfo()
Apply geometry information to each expansion.
Definition: ExpList.cpp:1500
boost::shared_ptr< Expansion0D > Expansion0DSharedPtr
Definition: Expansion0D.h:49
boost::shared_ptr< ExpList > ExpListSharedPtr
Shared pointer to an ExpList object.
void GenerateBoundaryConditionExpansion(const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::BoundaryConditions &bcs, const std::string variable)
Discretises the boundary conditions.
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:969
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
boost::shared_ptr< Expansion1D > Expansion1DSharedPtr
Definition: Expansion1D.h:53
Nektar::MultiRegions::DisContField1D::DisContField1D ( const LibUtilities::SessionReaderSharedPtr pSession,
const SpatialDomains::MeshGraphSharedPtr graph1D,
const SpatialDomains::CompositeMap domain,
const SpatialDomains::BoundaryConditions Allbcs,
const std::string &  variable,
bool  SetToOneSpaceDimension = false 
)

Constructor for a DisContField1D from a List of subdomains New Constructor for arterial network.

Constructor for use in multidomain computations where a domain list can be passed instead of graph1D

Parameters
domainSubdomain specified in the inputfile from which the DisContField1D is set up

Definition at line 451 of file DisContField1D.cpp.

References Nektar::MultiRegions::ExpList::ApplyGeomInfo(), Nektar::MultiRegions::ExpList::EvaluateBoundaryConditions(), FindPeriodicVertices(), GenerateBoundaryConditionExpansion(), GetDomainBCs(), and SetUpDG().

457  :
458  ExpList1D(pSession,graph1D,domain, true,variable,SetToOneSpaceDimension),
461  {
462  if (variable.compare("DefaultVar") != 0)
463  {
464  SpatialDomains::BoundaryConditionsSharedPtr DomBCs = GetDomainBCs(domain,Allbcs,variable);
465 
466  GenerateBoundaryConditionExpansion(graph1D,*DomBCs,variable);
467  EvaluateBoundaryConditions(0.0, variable);
468  ApplyGeomInfo();
469  FindPeriodicVertices(*DomBCs,variable);
470  }
471 
472  SetUpDG(variable);
473  }
ExpList1D()
The default constructor.
Definition: ExpList1D.cpp:77
void EvaluateBoundaryConditions(const NekDouble time=0.0, const std::string varName="", const NekDouble=NekConstants::kNekUnsetDouble, const NekDouble=NekConstants::kNekUnsetDouble)
Definition: ExpList.h:2249
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
void FindPeriodicVertices(const SpatialDomains::BoundaryConditions &bcs, const std::string variable)
Generate a associative map of periodic vertices in a mesh.
SpatialDomains::BoundaryConditionsSharedPtr GetDomainBCs(const SpatialDomains::CompositeMap &domain, const SpatialDomains::BoundaryConditions &Allbcs, const std::string &variable)
void SetUpDG(const std::string &variable)
void ApplyGeomInfo()
Apply geometry information to each expansion.
Definition: ExpList.cpp:1500
void GenerateBoundaryConditionExpansion(const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::BoundaryConditions &bcs, const std::string variable)
Discretises the boundary conditions.
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
boost::shared_ptr< BoundaryConditions > BoundaryConditionsSharedPtr
Definition: Conditions.h:284
Nektar::MultiRegions::DisContField1D::DisContField1D ( const DisContField1D In)

Constructs a 1D discontinuous field based on an existing field.

Constructs a field as a copy of an existing field.

Parameters
InExisting DisContField1D object to copy.

Definition at line 479 of file DisContField1D.cpp.

479  :
480  ExpList1D(In),
481  m_bndCondExpansions(In.m_bndCondExpansions),
482  m_bndConditions(In.m_bndConditions),
483  m_globalBndMat(In.m_globalBndMat),
484  m_trace(In.m_trace),
485  m_traceMap(In.m_traceMap),
486  m_boundaryVerts(In.m_boundaryVerts),
487  m_periodicVerts(In.m_periodicVerts),
488  m_periodicFwdCopy(In.m_periodicFwdCopy),
489  m_periodicBwdCopy(In.m_periodicBwdCopy),
490  m_leftAdjacentVerts(In.m_leftAdjacentVerts)
491  {
492  }
ExpList1D()
The default constructor.
Definition: ExpList1D.cpp:77
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
GlobalLinSysMapShPtr m_globalBndMat
Global boundary matrix.
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
std::set< int > m_boundaryVerts
A set storing the global IDs of any boundary edges.
ExpListSharedPtr m_trace
Trace space storage for points between elements.
std::vector< int > m_periodicFwdCopy
A vector indicating degress of freedom which need to be copied from forwards to backwards space in ca...
PeriodicMap m_periodicVerts
A map which identifies groups of periodic vertices.
Nektar::MultiRegions::DisContField1D::DisContField1D ( const ExpList1D In)

Constructs a 1D discontinuous field based on an existing field. (needed in order to use ContField( const ExpList1D &In) constructor.

Constructs a field as a copy of an existing explist1D field.

Parameters
InExisting ExpList1D object to copy.

Definition at line 499 of file DisContField1D.cpp.

499  :
500  ExpList1D(In)
501  {
502  }
ExpList1D()
The default constructor.
Definition: ExpList1D.cpp:77
Nektar::MultiRegions::DisContField1D::~DisContField1D ( )
virtual

Destructor.

Definition at line 507 of file DisContField1D.cpp.

508  {
509  }

Member Function Documentation

void Nektar::MultiRegions::DisContField1D::FindPeriodicVertices ( const SpatialDomains::BoundaryConditions bcs,
const std::string  variable 
)
protected

Generate a associative map of periodic vertices in a mesh.

Parameters
graph1DA mesh containing information about the domain and the spectral/hp element expansion.
bcsInformation about the boundary conditions.
variableSpecifies the field.
periodicVerticesMap into which the list of periodic vertices is placed.

Definition at line 571 of file DisContField1D.cpp.

References ASSERTL0, Nektar::StdRegions::eNoOrientation, Nektar::SpatialDomains::ePeriodic, Nektar::MultiRegions::ExpList::GetBoundaryCondition(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryConditions(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryRegions(), Nektar::iterator, Nektar::MultiRegions::ExpList::m_graph, m_periodicVerts, Nektar::MultiRegions::ExpList::m_session, CellMLToNektar.cellml_metadata::p, Nektar::LibUtilities::ReduceSum, and Vmath::Vsum().

Referenced by DisContField1D().

574  {
576  = bcs.GetBoundaryRegions();
578  = bcs.GetBoundaryConditions();
579 
581  = boost::dynamic_pointer_cast<
582  SpatialDomains::MeshGraph1D>(m_graph);
583  SpatialDomains::BoundaryRegionCollection::const_iterator it;
584 
586  m_session->GetComm()->GetRowComm();
587 
588  int i, region1ID, region2ID;
589 
591 
592  map<int,int> BregionToVertMap;
593 
594  // Construct list of all periodic Region and their global vertex on
595  // this process.
596  for (it = bregions.begin(); it != bregions.end(); ++it)
597  {
598  locBCond = GetBoundaryCondition(bconditions, it->first, variable);
599 
600  if (locBCond->GetBoundaryConditionType()
602  {
603  continue;
604  }
605 
606  int id = (*(it->second->begin()->second))[0]->GetGlobalID();
607 
608  BregionToVertMap[it->first] = id;
609  }
610 
612  set<int> islocal;
613 
614  int n = vComm->GetSize();
615  int p = vComm->GetRank();
616 
617  Array<OneD, int> nregions(n, 0);
618  nregions[p] = BregionToVertMap.size();
619  vComm->AllReduce(nregions, LibUtilities::ReduceSum);
620 
621  int totRegions = Vmath::Vsum(n, nregions, 1);
622 
623  Array<OneD, int> regOffset(n, 0);
624 
625  for (i = 1; i < n; ++i)
626  {
627  regOffset[i] = regOffset[i-1] + nregions[i-1];
628  }
629 
630  Array<OneD, int> bregmap(totRegions, 0);
631  Array<OneD, int> bregid (totRegions, 0);
632  for(i = regOffset[p], iit = BregionToVertMap.begin();
633  iit != BregionToVertMap.end(); ++iit, ++i)
634  {
635  bregid [i] = iit->first;
636  bregmap[i] = iit->second;
637  islocal.insert(iit->first);
638  }
639 
640  vComm->AllReduce(bregmap, LibUtilities::ReduceSum);
641  vComm->AllReduce(bregid, LibUtilities::ReduceSum);
642 
643  for (int i = 0; i < totRegions; ++i)
644  {
645  BregionToVertMap[bregid[i]] = bregmap[i];
646  }
647 
648  // Construct list of all periodic pairs local to this process.
649  for (it = bregions.begin(); it != bregions.end(); ++it)
650  {
651  locBCond = GetBoundaryCondition(bconditions, it->first, variable);
652 
653  if (locBCond->GetBoundaryConditionType()
655  {
656  continue;
657  }
658 
659  // Identify periodic boundary region IDs.
660  region1ID = it->first;
661  region2ID = boost::static_pointer_cast<
662  SpatialDomains::PeriodicBoundaryCondition>(
663  locBCond)->m_connectedBoundaryRegion;
664 
665  ASSERTL0(BregionToVertMap.count(region1ID) != 0,
666  "Cannot determine vertex of region1ID");
667 
668  ASSERTL0(BregionToVertMap.count(region2ID) != 0,
669  "Cannot determine vertex of region2ID");
670 
671  PeriodicEntity ent(BregionToVertMap[region2ID],
673  islocal.count(region2ID) != 0);
674 
675  m_periodicVerts[BregionToVertMap[region1ID]].push_back(ent);
676  }
677  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
boost::shared_ptr< Comm > CommSharedPtr
Pointer to a Communicator object.
Definition: Comm.h:55
std::map< int, BoundaryRegionShPtr > BoundaryRegionCollection
Definition: Conditions.h:217
SpatialDomains::MeshGraphSharedPtr m_graph
Mesh associated with this expansion list.
Definition: ExpList.h:972
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:969
std::map< int, BoundaryConditionMapShPtr > BoundaryConditionCollection
Definition: Conditions.h:226
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
boost::shared_ptr< MeshGraph1D > MeshGraph1DSharedPtr
Definition: MeshGraph1D.h:86
static SpatialDomains::BoundaryConditionShPtr GetBoundaryCondition(const SpatialDomains::BoundaryConditionCollection &collection, unsigned int index, const std::string &variable)
Definition: ExpList.cpp:3037
boost::shared_ptr< BoundaryConditionBase > BoundaryConditionShPtr
Definition: Conditions.h:219
PeriodicMap m_periodicVerts
A map which identifies groups of periodic vertices.
T Vsum(int n, const T *x, const int incx)
Subtract return sum(x)
Definition: Vmath.cpp:737
struct Nektar::MultiRegions::_PeriodicEntity PeriodicEntity
void Nektar::MultiRegions::DisContField1D::GenerateBoundaryConditionExpansion ( const SpatialDomains::MeshGraphSharedPtr graph1D,
const SpatialDomains::BoundaryConditions bcs,
const std::string  variable 
)
protected

Discretises the boundary conditions.

Generate the boundary condition expansion list

Parameters
graph1DA mesh, containing information about the domain and the spectral/hp element expansions.
bcsInformation about the enforced boundary conditions.
variableThe session variable associated with the boundary conditions to enforce.

Definition at line 521 of file DisContField1D.cpp.

References Nektar::SpatialDomains::ePeriodic, Nektar::MultiRegions::ExpList::GetBoundaryCondition(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryConditions(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryRegions(), Nektar::iterator, m_bndCondExpansions, m_bndConditions, and SetBoundaryConditionExpansion().

Referenced by DisContField1D().

525  {
526  int cnt = 0;
527 
529  = bcs.GetBoundaryRegions();
531  = bcs.GetBoundaryConditions();
532  SpatialDomains::BoundaryRegionCollection::const_iterator it;
533 
534  // count the number of non-periodic boundary points
535  for (it = bregions.begin(); it != bregions.end(); ++it)
536  {
537  const SpatialDomains::BoundaryConditionShPtr boundaryCondition =
538  GetBoundaryCondition(bconditions, it->first, variable);
539  if (boundaryCondition->GetBoundaryConditionType() !=
541  {
543  for (bregionIt = it->second->begin();
544  bregionIt != it->second->end(); bregionIt++)
545  {
546  cnt += bregionIt->second->size();
547  }
548  }
549  }
550 
552  = Array<OneD,MultiRegions::ExpListSharedPtr>(cnt);
553 
555  = Array<OneD,SpatialDomains::BoundaryConditionShPtr>(cnt);
556 
557  SetBoundaryConditionExpansion(graph1D,bcs,variable,
559  m_bndConditions);
560  }
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
void SetBoundaryConditionExpansion(const SpatialDomains::MeshGraphSharedPtr &graph1D, const SpatialDomains::BoundaryConditions &bcs, const std::string variable, Array< OneD, MultiRegions::ExpListSharedPtr > &bndCondExpansions, Array< OneD, SpatialDomains::BoundaryConditionShPtr > &bndConditions)
Populates the list of boundary condition expansions.
std::map< int, BoundaryRegionShPtr > BoundaryRegionCollection
Definition: Conditions.h:217
std::map< int, BoundaryConditionMapShPtr > BoundaryConditionCollection
Definition: Conditions.h:226
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
static SpatialDomains::BoundaryConditionShPtr GetBoundaryCondition(const SpatialDomains::BoundaryConditionCollection &collection, unsigned int index, const std::string &variable)
Definition: ExpList.cpp:3037
boost::shared_ptr< BoundaryConditionBase > BoundaryConditionShPtr
Definition: Conditions.h:219
void Nektar::MultiRegions::DisContField1D::GenerateFieldBnd1D ( SpatialDomains::BoundaryConditions bcs,
const std::string  variable 
)
protected
SpatialDomains::BoundaryConditionsSharedPtr Nektar::MultiRegions::DisContField1D::GetDomainBCs ( const SpatialDomains::CompositeMap domain,
const SpatialDomains::BoundaryConditions Allbcs,
const std::string &  variable 
)
private

Definition at line 341 of file DisContField1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL1, Nektar::SpatialDomains::BoundaryConditions::GetBoundaryConditions(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryRegions(), Nektar::iterator, and Nektar::MultiRegions::ExpList::m_session.

Referenced by DisContField1D().

345  {
347 
349 
351  map<int,int> GeometryToRegionsMap;
352 
353  SpatialDomains::BoundaryRegionCollection::const_iterator it;
354 
356  = Allbcs.GetBoundaryRegions();
358  = Allbcs.GetBoundaryConditions();
359 
360  // Set up a map of all boundary regions
361  for(it = bregions.begin(); it != bregions.end(); ++it)
362  {
364  for (bregionIt = it->second->begin();
365  bregionIt != it->second->end(); bregionIt++)
366  {
367  // can assume that all regions only contain one point in 1D
368  // Really do not need loop above
369  int id = (*(bregionIt->second))[0]->GetGlobalID();
370  GeometryToRegionsMap[id] = it->first;
371  }
372  }
373 
375  map<int,SpatialDomains::GeometrySharedPtr> EndOfDomain;
376 
377  // Now find out which points in domain have only one vertex
378  for(domIt = domain.begin(); domIt != domain.end(); ++domIt)
379  {
380  SpatialDomains::Composite geomvector = domIt->second;
381  for(int i = 0; i < geomvector->size(); ++i)
382  {
383  for(int j = 0; j < 2; ++j)
384  {
385  int vid = (*geomvector)[i]->GetVid(j);
386  if(EndOfDomain.count(vid) == 0)
387  {
388  EndOfDomain[vid] = (*geomvector)[i]->GetVertex(j);
389  }
390  else
391  {
392  EndOfDomain.erase(vid);
393  }
394  }
395  }
396  }
397  ASSERTL1(EndOfDomain.size() == 2,"Did not find two ends of domain");
398 
400  int numNewBc = 1;
401  for(regIt = EndOfDomain.begin(); regIt != EndOfDomain.end(); ++regIt)
402  {
403  if(GeometryToRegionsMap.count(regIt->first) != 0) // Set up boundary condition up
404  {
405  map<int,int>::iterator iter = GeometryToRegionsMap.find(regIt->first);
406  ASSERTL1(iter != GeometryToRegionsMap.end(),"Failied to find GeometryToRegionMap");
407  int regionId = iter->second;
408  SpatialDomains::BoundaryRegionCollection::const_iterator bregionsIter = bregions.find(regionId);
409  ASSERTL1(bregionsIter != bregions.end(),"Failed to find boundary region");
410  SpatialDomains::BoundaryRegionShPtr breg = bregionsIter->second;
411  returnval->AddBoundaryRegions (regionId,breg);
412 
413  SpatialDomains::BoundaryConditionCollection::const_iterator bconditionsIter = bconditions.find(regionId);
414  ASSERTL1(bconditionsIter != bconditions.end(),"Failed to find boundary collection");
415  SpatialDomains::BoundaryConditionMapShPtr bcond = bconditionsIter->second;
416  returnval->AddBoundaryConditions(regionId,bcond);
417  }
418  else // Set up an undefined region.
419  {
420  SpatialDomains::BoundaryRegionShPtr breg(MemoryManager<SpatialDomains::BoundaryRegion>::AllocateSharedPtr());
421 
422  // Set up Composite (GemetryVector) to contain vertex and put into bRegion
423  SpatialDomains::Composite gvec(MemoryManager<SpatialDomains::GeometryVector>::AllocateSharedPtr());
424  gvec->push_back(regIt->second);
425  (*breg)[regIt->first] = gvec;
426 
427  returnval->AddBoundaryRegions(bregions.size()+numNewBc,breg);
428 
430 
431  // Set up just boundary condition for this variable.
432  SpatialDomains::BoundaryConditionShPtr notDefinedCondition(MemoryManager<SpatialDomains::NotDefinedBoundaryCondition>::AllocateSharedPtr(m_session, "0"));
433  (*bCondition)[variable] = notDefinedCondition;
434 
435  returnval->AddBoundaryConditions(bregions.size()+numNewBc,bCondition);
436  ++numNewBc;
437 
438  }
439  }
440 
441  return returnval;
442  }
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
std::map< int, Composite >::const_iterator CompositeMapConstIter
Definition: MeshGraph.h:117
std::map< int, BoundaryRegionShPtr > BoundaryRegionCollection
Definition: Conditions.h:217
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:969
std::map< int, BoundaryConditionMapShPtr > BoundaryConditionCollection
Definition: Conditions.h:226
boost::shared_ptr< GeometryVector > Composite
Definition: MeshGraph.h:114
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
boost::shared_ptr< BoundaryConditionMap > BoundaryConditionMapShPtr
Definition: Conditions.h:225
boost::shared_ptr< BoundaryConditions > BoundaryConditionsSharedPtr
Definition: Conditions.h:284
boost::shared_ptr< BoundaryConditionBase > BoundaryConditionShPtr
Definition: Conditions.h:219
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:228
boost::shared_ptr< BoundaryRegion > BoundaryRegionShPtr
Definition: Conditions.h:215
GlobalLinSysSharedPtr Nektar::MultiRegions::DisContField1D::GetGlobalBndLinSys ( const GlobalLinSysKey mkey)

For a given key, returns the associated global linear system.

Definition at line 752 of file DisContField1D.cpp.

References ASSERTL0, ASSERTL1, Nektar::MultiRegions::eDirectFullMatrix, Nektar::StdRegions::eHybridDGHelmBndLam, Nektar::MultiRegions::ExpList::GenGlobalBndLinSys(), Nektar::MultiRegions::GlobalLinSysKey::GetGlobalSysSolnType(), Nektar::MultiRegions::GlobalMatrixKey::GetMatrixType(), Nektar::iterator, m_globalBndMat, and m_traceMap.

Referenced by v_HelmSolve().

754  {
755  ASSERTL0(mkey.GetMatrixType() == StdRegions::eHybridDGHelmBndLam,
756  "Routine currently only tested for HybridDGHelmholtz");
757 
758  ASSERTL1(mkey.GetGlobalSysSolnType() != eDirectFullMatrix,
759  "Full matrix global systems are not supported for HDG "
760  "expansions");
761 
762  ASSERTL1(mkey.GetGlobalSysSolnType()
763  ==m_traceMap->GetGlobalSysSolnType(),
764  "The local to global map is not set up for the requested "
765  "solution type");
766 
767  GlobalLinSysSharedPtr glo_matrix;
768  GlobalLinSysMap::iterator matrixIter = m_globalBndMat->find(mkey);
769 
770  if (matrixIter == m_globalBndMat->end())
771  {
772  glo_matrix = GenGlobalBndLinSys(mkey,m_traceMap);
773  (*m_globalBndMat)[mkey] = glo_matrix;
774  }
775  else
776  {
777  glo_matrix = matrixIter->second;
778  }
779 
780  return glo_matrix;
781  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
boost::shared_ptr< GlobalLinSys > GenGlobalBndLinSys(const GlobalLinSysKey &mkey, const AssemblyMapSharedPtr &locToGloMap)
Generate a GlobalLinSys from information provided by the key "mkey" and the mapping provided in LocTo...
Definition: ExpList.cpp:1270
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
GlobalLinSysMapShPtr m_globalBndMat
Global boundary matrix.
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
boost::shared_ptr< GlobalLinSys > GlobalLinSysSharedPtr
Pointer to a GlobalLinSys object.
Definition: GlobalLinSys.h:52
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:228
vector< bool > & Nektar::MultiRegions::DisContField1D::GetNegatedFluxNormal ( void  )

Definition at line 784 of file DisContField1D.cpp.

References Nektar::MultiRegions::ExpList::GetExpSize(), m_negatedFluxNormal, and m_traceMap.

Referenced by v_AddTraceIntegral().

785  {
786  if(m_negatedFluxNormal.size() == 0)
787  {
788  Array<OneD, Array<OneD, LocalRegions::ExpansionSharedPtr> >
789  &elmtToTrace = m_traceMap->GetElmtToTrace();
790 
791  m_negatedFluxNormal.resize(2*GetExpSize());
792 
793  for(int i = 0; i < GetExpSize(); ++i)
794  {
795 
796  for(int v = 0; v < 2; ++v)
797  {
798 
800  elmtToTrace[i][v]->as<LocalRegions::Expansion0D>();
801 
802  if(vertExp->GetLeftAdjacentElementExp()->GetGeom()->GetGlobalID() != (*m_exp)[i]->GetGeom()->GetGlobalID())
803  {
804  m_negatedFluxNormal[2*i+v] = true;
805  }
806  else
807  {
808  m_negatedFluxNormal[2*i+v] = false;
809  }
810  }
811  }
812 
813  }
814 
815  return m_negatedFluxNormal;
816  }
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2054
boost::shared_ptr< Expansion0D > Expansion0DSharedPtr
Definition: Expansion0D.h:49
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
bool Nektar::MultiRegions::DisContField1D::IsLeftAdjacentVertex ( const int  n,
const int  e 
)
private

Definition at line 286 of file DisContField1D.cpp.

References ASSERTL2, Nektar::iterator, m_boundaryVerts, Nektar::MultiRegions::ExpList::m_exp, m_periodicVerts, m_trace, and m_traceMap.

Referenced by SetUpDG().

287  {
290  m_traceMap->GetElmtToTrace()[n][e]->as<LocalRegions::Expansion0D>();
291 
292 
293  bool fwd = true;
294  if (traceEl->GetLeftAdjacentElementVertex () == -1 ||
295  traceEl->GetRightAdjacentElementVertex() == -1)
296  {
297  // Boundary edge (1 connected element). Do nothing in
298  // serial.
299  it = m_boundaryVerts.find(traceEl->GetElmtId());
300 
301  // If the edge does not have a boundary condition set on
302  // it, then assume it is a partition edge or periodic.
303  if (it == m_boundaryVerts.end())
304  {
305  int traceGeomId = traceEl->GetGeom0D()->GetGlobalID();
307  traceGeomId);
308 
309  if (pIt != m_periodicVerts.end() && !pIt->second[0].isLocal)
310  {
311  fwd = traceGeomId == min(traceGeomId,pIt->second[0].id);
312  }
313  else
314  {
315  int offset = m_trace->GetPhys_Offset(traceEl->GetElmtId());
316  fwd = m_traceMap->
317  GetTraceToUniversalMapUnique(offset) >= 0;
318  }
319  }
320  }
321  else if (traceEl->GetLeftAdjacentElementVertex () != -1 &&
322  traceEl->GetRightAdjacentElementVertex() != -1)
323  {
324  // Non-boundary edge (2 connected elements).
325  fwd = dynamic_cast<Nektar::StdRegions::StdExpansion*>
326  (traceEl->GetLeftAdjacentElementExp().get()) ==
327  (*m_exp)[n].get();
328  }
329  else
330  {
331  ASSERTL2(false, "Unconnected trace element!");
332  }
333 
334  return fwd;
335  }
boost::shared_ptr< Expansion0D > Expansion0DSharedPtr
Definition: Expansion0D.h:49
The base class for all shapes.
Definition: StdExpansion.h:69
boost::shared_ptr< LocalRegions::ExpansionVector > m_exp
The list of local expansions.
Definition: ExpList.h:1036
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
#define ASSERTL2(condition, msg)
Assert Level 2 – Debugging which is used FULLDEBUG compilation mode. This level assert is designed t...
Definition: ErrorUtil.hpp:250
std::set< int > m_boundaryVerts
A set storing the global IDs of any boundary edges.
ExpListSharedPtr m_trace
Trace space storage for points between elements.
PeriodicMap m_periodicVerts
A map which identifies groups of periodic vertices.
void Nektar::MultiRegions::DisContField1D::SetBoundaryConditionExpansion ( const SpatialDomains::MeshGraphSharedPtr graph1D,
const SpatialDomains::BoundaryConditions bcs,
const std::string  variable,
Array< OneD, MultiRegions::ExpListSharedPtr > &  bndCondExpansions,
Array< OneD, SpatialDomains::BoundaryConditionShPtr > &  bndConditions 
)
protected

Populates the list of boundary condition expansions.

Parameters
graph1DA mesh containing information about the domain and the Spectral/hp element expansion.
bcsInformation about the boundary conditions.
variableSpecifies the field.
bndCondExpansionsArray of ExpList1D objects each containing a 1D spectral/hp element expansion on a single boundary region.
bncConditionsArray of BoundaryCondition objects which contain information about the boundary conditions on the different boundary regions.

Definition at line 692 of file DisContField1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::SpatialDomains::ePeriodic, Nektar::MultiRegions::ExpList::GetBoundaryCondition(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryConditions(), Nektar::SpatialDomains::BoundaryConditions::GetBoundaryRegions(), and Nektar::iterator.

Referenced by GenerateBoundaryConditionExpansion().

700  {
701  int k;
702  int cnt = 0;
703 
705  = bcs.GetBoundaryRegions();
707  = bcs.GetBoundaryConditions();
708  SpatialDomains::BoundaryRegionCollection::const_iterator it;
709 
713 
714  cnt = 0;
715  for (it = bregions.begin(); it != bregions.end(); ++it)
716  {
717  locBCond = GetBoundaryCondition(bconditions, it->first, variable);
718 
719  if (locBCond->GetBoundaryConditionType() !=
721  {
723  for (bregionIt = it->second->begin();
724  bregionIt != it->second->end(); bregionIt++)
725  {
726  for (k = 0; k < bregionIt->second->size(); k++)
727  {
728  if((vert = boost::dynamic_pointer_cast
729  <SpatialDomains::PointGeom>(
730  (*bregionIt->second)[k])))
731  {
732  locPointExp
735  bndCondExpansions[cnt] = locPointExp;
736  bndConditions[cnt++] = locBCond;
737  }
738  else
739  {
740  ASSERTL0(false,
741  "dynamic cast to a vertex failed");
742  }
743  }
744  }
745  }
746  }
747  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
std::map< int, BoundaryRegionShPtr > BoundaryRegionCollection
Definition: Conditions.h:217
std::map< int, BoundaryConditionMapShPtr > BoundaryConditionCollection
Definition: Conditions.h:226
boost::shared_ptr< ExpList0D > ExpList0DSharedPtr
Shared pointer to an ExpList0D object.
Definition: ExpList0D.h:54
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
static SpatialDomains::BoundaryConditionShPtr GetBoundaryCondition(const SpatialDomains::BoundaryConditionCollection &collection, unsigned int index, const std::string &variable)
Definition: ExpList.cpp:3037
boost::shared_ptr< BoundaryConditionBase > BoundaryConditionShPtr
Definition: Conditions.h:219
boost::shared_ptr< PointGeom > PointGeomSharedPtr
Definition: Geometry.h:60
void Nektar::MultiRegions::DisContField1D::SetMultiDomainBoundaryConditionExpansion ( const SpatialDomains::MeshGraphSharedPtr graph1D,
const SpatialDomains::BoundaryConditions bcs,
const std::string  variable,
Array< OneD, MultiRegions::ExpListSharedPtr > &  bndCondExpansions,
Array< OneD, SpatialDomains::BoundaryConditionShPtr > &  bndConditions,
int  subdomain 
)
protected

Populates the list of boundary condition expansions in multidomain case.

void Nektar::MultiRegions::DisContField1D::SetUpDG ( const std::string &  variable)
private

Definition at line 127 of file DisContField1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, ASSERTL1, Nektar::SpatialDomains::ePeriodic, Nektar::MultiRegions::_PeriodicEntity::id, IsLeftAdjacentVertex(), Nektar::MultiRegions::_PeriodicEntity::isLocal, Nektar::iterator, m_bndCondExpansions, m_bndConditions, m_boundaryVerts, Nektar::MultiRegions::ExpList::m_exp, m_globalBndMat, Nektar::MultiRegions::ExpList::m_graph, m_leftAdjacentVerts, m_periodicBwdCopy, m_periodicFwdCopy, m_periodicVerts, Nektar::MultiRegions::ExpList::m_session, m_trace, m_traceMap, Nektar::MultiRegions::NullExpListSharedPtr, Nektar::LocalRegions::Expansion0D::SetAdjacentElementExp(), and Nektar::MultiRegions::ExpList::SetUpPhysNormals().

Referenced by DisContField1D().

128  {
129  // Check for multiple calls
131  {
132  return;
133  }
134 
136  boost::dynamic_pointer_cast<SpatialDomains::MeshGraph1D>(
137  m_graph);
138 
140 
145  *m_exp,graph1D,
147 
148  m_trace = boost::dynamic_pointer_cast<ExpList>(trace);
149 
151  AllocateSharedPtr(m_session, graph1D, trace, *this,
153  m_periodicVerts, variable);
154 
155  if (m_session->DefinesCmdLineArgument("verbose"))
156  {
157  m_traceMap->PrintStats(std::cout, variable);
158  }
159 
160  // Scatter trace points to 1D elements. For each element, we find
161  // the trace point associated to each vertex. The element then
162  // retains a pointer to the trace space points, to ensure
163  // uniqueness of normals when retrieving from two adjoining
164  // elements which do not lie in a plane.
165 
166  int ElmtPointGeom = 0;
167  int TracePointGeom = 0;
170  for (int i = 0; i < m_exp->size(); ++i)
171  {
172  exp1d = (*m_exp)[i]->as<LocalRegions::Expansion1D>();
173  for (int j = 0; j < exp1d->GetNverts(); ++j)
174  {
175  ElmtPointGeom = (exp1d->GetGeom1D())->GetVid(j);
176 
177  for (int k = 0; k < m_trace->GetExpSize(); ++k)
178  {
179  TracePointGeom = m_trace->GetExp(k)->GetGeom()->GetVid(0);
180 
181  if (TracePointGeom == ElmtPointGeom)
182  {
183  exp0d = m_trace->GetExp(k)->as<LocalRegions::Expansion0D>();
184  exp0d->SetAdjacentElementExp(j,exp1d);
185  break;
186  }
187  }
188  }
189  }
190 
192 
193  int cnt, n, e;
194 
195  // Identify boundary verts
196  for(cnt = 0, n = 0; n < m_bndCondExpansions.num_elements(); ++n)
197  {
198  if (m_bndConditions[n]->GetBoundaryConditionType() !=
200  {
201  for(e = 0; e < m_bndCondExpansions[n]->GetExpSize(); ++e)
202  {
203  m_boundaryVerts.insert(
204  m_traceMap->GetBndCondTraceToGlobalTraceMap(cnt+e));
205  }
206  }
207  else
208  {
209  ASSERTL0(false,"Periodic verts need setting up");
210  }
211  cnt += m_bndCondExpansions[n]->GetExpSize();
212  }
213 
214  // Set up left-adjacent edge list.
215  m_leftAdjacentVerts.resize(2*((*m_exp).size()));
216 
217  // count size of trace
218  for (cnt = n = 0; n < m_exp->size(); ++n)
219  {
220  for (int v = 0; v < (*m_exp)[n]->GetNverts(); ++v, ++cnt)
221  {
223  }
224  }
225 
226 
227  boost::unordered_map<int,pair<int,int> > perVertToExpMap;
228  boost::unordered_map<int,pair<int,int> >::iterator it2;
229  for (n = 0; n < m_exp->size(); ++n)
230  {
231  for (int v = 0; v < (*m_exp)[n]->GetNverts(); ++v)
232  {
234  (*m_exp)[n]->GetGeom()->GetVid(v));
235 
236  if (it != m_periodicVerts.end())
237  {
238  perVertToExpMap[it->first] = make_pair(n,v);
239  }
240  }
241  }
242 
243 
244  // Set up mapping to copy Fwd of periodic bcs to Bwd of other edge.
245  for (n = 0; n < m_exp->size(); ++n)
246  {
247  for (int v = 0; v < (*m_exp)[n]->GetNverts(); ++v)
248  {
249  int vertGeomId = (*m_exp)[n]->GetGeom()->GetVid(v);
250 
251  // Check to see if this face is periodic.
252  PeriodicMap::iterator it = m_periodicVerts.find(vertGeomId);
253 
254  if (it != m_periodicVerts.end())
255  {
256  const PeriodicEntity &ent = it->second[0];
257  it2 = perVertToExpMap.find(ent.id);
258 
259  if (it2 == perVertToExpMap.end())
260  {
261  if (m_session->GetComm()->GetRowComm()->GetSize() > 1 &&
262  !ent.isLocal)
263  {
264  continue;
265  }
266  else
267  {
268  ASSERTL1(false, "Periodic vert not found!");
269  }
270  }
271 
272  int offset = m_trace->GetPhys_Offset((m_traceMap->GetElmtToTrace())
273  [n][v]->GetElmtId());
274  m_periodicFwdCopy.push_back(offset);
275 
276  int offset2 = m_trace->GetPhys_Offset((m_traceMap->GetElmtToTrace())
277  [it2->second.first]
278  [it2->second.second]->GetElmtId());
279  m_periodicBwdCopy.push_back(offset2);
280  }
281  }
282  }
283  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
ExpList()
The default constructor.
Definition: ExpList.cpp:95
static ExpListSharedPtr NullExpListSharedPtr
Definition: ExpList.h:1485
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
boost::shared_ptr< Expansion0D > Expansion0DSharedPtr
Definition: Expansion0D.h:49
boost::shared_ptr< LocalRegions::ExpansionVector > m_exp
The list of local expansions.
Definition: ExpList.h:1036
SpatialDomains::MeshGraphSharedPtr m_graph
Mesh associated with this expansion list.
Definition: ExpList.h:972
LibUtilities::SessionReaderSharedPtr m_session
Session.
Definition: ExpList.h:969
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
bool IsLeftAdjacentVertex(const int n, const int e)
GlobalLinSysMapShPtr m_globalBndMat
Global boundary matrix.
boost::shared_ptr< ExpList0D > ExpList0DSharedPtr
Shared pointer to an ExpList0D object.
Definition: ExpList0D.h:54
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
boost::shared_ptr< MeshGraph1D > MeshGraph1DSharedPtr
Definition: MeshGraph1D.h:86
std::set< int > m_boundaryVerts
A set storing the global IDs of any boundary edges.
ExpListSharedPtr m_trace
Trace space storage for points between elements.
std::vector< int > m_periodicFwdCopy
A vector indicating degress of freedom which need to be copied from forwards to backwards space in ca...
PeriodicMap m_periodicVerts
A map which identifies groups of periodic vertices.
struct Nektar::MultiRegions::_PeriodicEntity PeriodicEntity
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:228
boost::shared_ptr< Expansion1D > Expansion1DSharedPtr
Definition: Expansion1D.h:53
void Nektar::MultiRegions::DisContField1D::v_AddTraceIntegral ( const Array< OneD, const NekDouble > &  Fn,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 998 of file DisContField1D.cpp.

References Nektar::LibUtilities::BasisManager(), Nektar::LibUtilities::eGauss_Lagrange, Nektar::LibUtilities::eGaussGaussLegendre, Nektar::MultiRegions::ExpList::GetCoeff_Offset(), Nektar::MultiRegions::ExpList::GetExpSize(), GetNegatedFluxNormal(), Nektar::MultiRegions::ExpList::GetTrace(), Nektar::MultiRegions::ExpList::m_exp, m_traceMap, and CellMLToNektar.cellml_metadata::p.

1001  {
1002  int n,offset, t_offset;
1003 
1004  Array<OneD, Array<OneD, LocalRegions::ExpansionSharedPtr> >
1005  &elmtToTrace = m_traceMap->GetElmtToTrace();
1006 
1007  vector<bool> negatedFluxNormal = GetNegatedFluxNormal();
1008 
1009  for (n = 0; n < GetExpSize(); ++n)
1010  {
1011  // Number of coefficients on each element
1012  int e_ncoeffs = (*m_exp)[n]->GetNcoeffs();
1013 
1014  offset = GetCoeff_Offset(n);
1015 
1016  // Implementation for every points except Gauss points
1017  if ((*m_exp)[n]->GetBasis(0)->GetBasisType() !=
1019  {
1020  t_offset = GetTrace()->GetCoeff_Offset(elmtToTrace[n][0]->GetElmtId());
1021  if(negatedFluxNormal[2*n])
1022  {
1023  outarray[offset] -= Fn[t_offset];
1024  }
1025  else
1026  {
1027  outarray[offset] += Fn[t_offset];
1028  }
1029 
1030  t_offset = GetTrace()->GetCoeff_Offset(elmtToTrace[n][1]->GetElmtId());
1031 
1032  if(negatedFluxNormal[2*n+1])
1033  {
1034  outarray[offset+(*m_exp)[n]->GetVertexMap(1)] -= Fn[t_offset];
1035  }
1036  else
1037  {
1038  outarray[offset+(*m_exp)[n]->GetVertexMap(1)] += Fn[t_offset];
1039  }
1040 
1041  }
1042  else
1043  {
1044 #if 0
1045  DNekMatSharedPtr m_Ixm;
1047  const LibUtilities::PointsKey
1048  BS_p(e_ncoeffs,LibUtilities::eGaussGaussLegendre);
1049  const LibUtilities::BasisKey
1050  BS_k(LibUtilities::eGauss_Lagrange,e_ncoeffs,BS_p);
1051 
1052  BASE = LibUtilities::BasisManager()[BS_k];
1053 
1054  Array<OneD, NekDouble> coords(3, 0.0);
1055 
1056  int j;
1057 
1058  for(p = 0; p < 2; ++p)
1059  {
1060  NekDouble vertnorm = 0.0;
1061  for (int i=0; i<((*m_exp)[n]->
1062  GetVertexNormal(p)).num_elements(); i++)
1063  {
1064  vertnorm += ((*m_exp)[n]->GetVertexNormal(p))[i][0];
1065  coords[0] = vertnorm ;
1066  }
1067 
1068  t_offset = GetTrace()->GetPhys_Offset(n+p);
1069 
1070  if (vertnorm >= 0.0)
1071  {
1072  m_Ixm = BASE->GetI(coords);
1073 
1074 
1075  for (j = 0; j < e_ncoeffs; j++)
1076  {
1077  outarray[offset + j] +=
1078  (m_Ixm->GetPtr())[j] * Fn[t_offset];
1079  }
1080  }
1081 
1082  if (vertnorm < 0.0)
1083  {
1084  m_Ixm = BASE->GetI(coords);
1085 
1086  for (j = 0; j < e_ncoeffs; j++)
1087  {
1088  outarray[offset + j] -=
1089  (m_Ixm->GetPtr())[j] * Fn[t_offset];
1090  }
1091  }
1092  }
1093 #else
1094  int j;
1095  static DNekMatSharedPtr m_Ixm, m_Ixp;
1096  static int sav_ncoeffs = 0;
1097  if(!m_Ixm.get() || e_ncoeffs != sav_ncoeffs)
1098  {
1100  const LibUtilities::PointsKey
1101  BS_p(e_ncoeffs,LibUtilities::eGaussGaussLegendre);
1102  const LibUtilities::BasisKey
1103  BS_k(LibUtilities::eGauss_Lagrange,e_ncoeffs,BS_p);
1104 
1105  BASE = LibUtilities::BasisManager()[BS_k];
1106 
1107  Array<OneD, NekDouble> coords(1, 0.0);
1108 
1109  coords[0] = -1.0;
1110  m_Ixm = BASE->GetI(coords);
1111 
1112  coords[0] = 1.0;
1113  m_Ixp = BASE->GetI(coords);
1114 
1115  sav_ncoeffs = e_ncoeffs;
1116  }
1117 
1118  t_offset = GetTrace()->GetCoeff_Offset(elmtToTrace[n][0]->GetElmtId());
1119  if(negatedFluxNormal[2*n])
1120  {
1121  for (j = 0; j < e_ncoeffs; j++)
1122  {
1123  outarray[offset + j] -=
1124  (m_Ixm->GetPtr())[j] * Fn[t_offset];
1125  }
1126  }
1127  else
1128  {
1129  for (j = 0; j < e_ncoeffs; j++)
1130  {
1131  outarray[offset + j] +=
1132  (m_Ixm->GetPtr())[j] * Fn[t_offset];
1133  }
1134  }
1135 
1136  t_offset = GetTrace()->GetCoeff_Offset(elmtToTrace[n][1]->GetElmtId());
1137  if (negatedFluxNormal[2*n+1])
1138  {
1139  for (j = 0; j < e_ncoeffs; j++)
1140  {
1141  outarray[offset + j] -=
1142  (m_Ixp->GetPtr())[j] * Fn[t_offset];
1143  }
1144  }
1145  else
1146  {
1147  for (j = 0; j < e_ncoeffs; j++)
1148  {
1149  outarray[offset + j] +=
1150  (m_Ixp->GetPtr())[j] * Fn[t_offset];
1151  }
1152  }
1153 #endif
1154  }
1155  }
1156  }
int GetCoeff_Offset(int n) const
Get the start offset position for a global list of m_coeffs correspoinding to element n...
Definition: ExpList.h:2084
Lagrange Polynomials using the Gauss points .
Definition: BasisType.h:54
boost::shared_ptr< DNekMat > DNekMatSharedPtr
Definition: NekTypeDefs.hpp:70
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2054
1D Gauss-Gauss-Legendre quadrature points
Definition: PointsType.h:49
BasisManagerT & BasisManager(void)
boost::shared_ptr< ExpList > & GetTrace()
Definition: ExpList.h:2159
boost::shared_ptr< LocalRegions::ExpansionVector > m_exp
The list of local expansions.
Definition: ExpList.h:1036
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
double NekDouble
std::vector< bool > & GetNegatedFluxNormal(void)
boost::shared_ptr< Basis > BasisSharedPtr
void Nektar::MultiRegions::DisContField1D::v_EvaluateBoundaryConditions ( const NekDouble  time = 0.0,
const std::string  varName = "",
const NekDouble  x2_in = NekConstants::kNekUnsetDouble,
const NekDouble  x3_in = NekConstants::kNekUnsetDouble 
)
protectedvirtual

Evaluate all boundary conditions at a given time..

Based on the expression $g(x,t)$ for the boundary conditions, this function evaluates the boundary conditions for all boundaries at time-level t.

Parameters
timeThe time at which the boundary conditions should be evaluated.
bndCondExpansionsList of boundary expansions.
bndConditionsInformation about the boundary conditions.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1295 of file DisContField1D.cpp.

References ASSERTL0, Nektar::SpatialDomains::eDirichlet, Nektar::SpatialDomains::eNeumann, Nektar::SpatialDomains::eNotDefined, Nektar::SpatialDomains::eRobin, Nektar::MultiRegions::ExpList::GetCoeff(), Nektar::NekConstants::kNekUnsetDouble, m_bndCondExpansions, and m_bndConditions.

1300  {
1301  int i;
1302 
1303  Array<OneD, NekDouble> x0(1);
1304  Array<OneD, NekDouble> x1(1);
1305  Array<OneD, NekDouble> x2(1);
1306 
1307  for (i = 0; i < m_bndCondExpansions.num_elements(); ++i)
1308  {
1309  if (time == 0.0 || m_bndConditions[i]->IsTimeDependent())
1310  {
1311  m_bndCondExpansions[i]->GetCoords(x0, x1, x2);
1312 
1313  if (x2_in != NekConstants::kNekUnsetDouble && x3_in !=
1315  {
1316  x1[0] = x2_in;
1317  x2[0] = x3_in;
1318  }
1319 
1320  if (m_bndConditions[i]->GetBoundaryConditionType() ==
1322  {
1323  m_bndCondExpansions[i]->SetCoeff(0,
1324  (boost::static_pointer_cast<SpatialDomains
1325  ::DirichletBoundaryCondition>(m_bndConditions[i])
1326  ->m_dirichletCondition).Evaluate(x0[0],x1[0],x2[0],time));
1327  m_bndCondExpansions[i]->SetPhys(0,m_bndCondExpansions[i]->GetCoeff(0));
1328  }
1329  else if (m_bndConditions[i]->GetBoundaryConditionType()
1331  {
1332  m_bndCondExpansions[i]->SetCoeff(0,
1333  (boost::static_pointer_cast<SpatialDomains
1334  ::NeumannBoundaryCondition>(m_bndConditions[i])
1335  ->m_neumannCondition).Evaluate(x0[0],x1[0],x2[0],time));
1336  }
1337  else if (m_bndConditions[i]->GetBoundaryConditionType()
1339  {
1340  m_bndCondExpansions[i]->SetCoeff(0,
1341  (boost::static_pointer_cast<SpatialDomains
1342  ::RobinBoundaryCondition>(m_bndConditions[i])
1343  ->m_robinFunction).Evaluate(x0[0],x1[0],x2[0],time));
1344 
1345  }
1346  else if (m_bndConditions[i]->GetBoundaryConditionType()
1348  {
1349  }
1350  else
1351  {
1352  ASSERTL0(false, "This type of BC not implemented yet");
1353  }
1354  }
1355  }
1356  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
static const NekDouble kNekUnsetDouble
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
NekDouble GetCoeff(int i)
Get the i th value (coefficient) of m_coeffs.
Definition: ExpList.h:2023
void Nektar::MultiRegions::DisContField1D::v_ExtractTracePhys ( Array< OneD, NekDouble > &  outarray)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 953 of file DisContField1D.cpp.

References ASSERTL1, Nektar::MultiRegions::ExpList::m_phys, and Nektar::MultiRegions::ExpList::m_physState.

955  {
956  ASSERTL1(m_physState == true,"local physical space is not true ");
957  v_ExtractTracePhys(m_phys, outarray);
958  }
Array< OneD, NekDouble > m_phys
The global expansion evaluated at the quadrature points.
Definition: ExpList.h:1015
virtual void v_ExtractTracePhys(Array< OneD, NekDouble > &outarray)
bool m_physState
The state of the array m_phys.
Definition: ExpList.h:1024
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:228
void Nektar::MultiRegions::DisContField1D::v_ExtractTracePhys ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray 
)
protectedvirtual

This method extracts the trace (verts in 1D) from the field inarray and puts the values in outarray.

It assumes the field is C0 continuous so that it can overwrite the edge data when visited by the two adjacent elements.

Parameters
inarrayAn array containing the 1D data from which we wish to extract the edge data.
outarrayThe resulting edge information.

This will not work for non-boundary expansions

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 973 of file DisContField1D.cpp.

References ASSERTL1, Nektar::MultiRegions::ExpList::GetExpSize(), Nektar::MultiRegions::ExpList::GetPhys_Offset(), m_trace, m_traceMap, and CellMLToNektar.cellml_metadata::p.

976  {
977  // Loop over elemente and collect forward expansion
978  int nexp = GetExpSize();
979  int n,p,offset,phys_offset;
980 
981  ASSERTL1(outarray.num_elements() >= m_trace->GetExpSize(),
982  "input array is of insufficient length");
983 
984  for (n = 0; n < nexp; ++n)
985  {
986  phys_offset = GetPhys_Offset(n);
987 
988  for (p = 0; p < (*m_exp)[n]->GetNverts(); ++p)
989  {
990  offset = m_trace->GetPhys_Offset(
991  (m_traceMap->GetElmtToTrace())[n][p]->GetElmtId());
992  (*m_exp)[n]->GetVertexPhysVals(p, inarray + phys_offset,
993  outarray[offset]);
994  }
995  }
996  }
int GetPhys_Offset(int n) const
Get the start offset position for a global list of m_phys correspoinding to element n...
Definition: ExpList.h:2092
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2054
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
ExpListSharedPtr m_trace
Trace space storage for points between elements.
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:228
virtual const Array<OneD,const MultiRegions::ExpListSharedPtr>& Nektar::MultiRegions::DisContField1D::v_GetBndCondExpansions ( void  )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 218 of file DisContField1D.h.

References m_bndCondExpansions.

219  {
220  return m_bndCondExpansions;
221  }
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
virtual const Array<OneD,const SpatialDomains::BoundaryConditionShPtr>& Nektar::MultiRegions::DisContField1D::v_GetBndConditions ( void  )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Reimplemented in Nektar::MultiRegions::ContField1D.

Definition at line 224 of file DisContField1D.h.

References m_bndConditions.

225  {
226  return m_bndConditions;
227  }
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
void Nektar::MultiRegions::DisContField1D::v_GetBndElmtExpansion ( int  i,
boost::shared_ptr< ExpList > &  result,
const bool  DeclareCoeffPhysArrays 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1413 of file DisContField1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::ExpList::GetBoundaryToElmtMap(), Nektar::MultiRegions::ExpList::GetCoeff_Offset(), Nektar::MultiRegions::ExpList::GetCoeffs(), Nektar::MultiRegions::ExpList::GetExp(), Nektar::MultiRegions::ExpList::GetPhys(), Nektar::MultiRegions::ExpList::GetPhys_Offset(), m_bndCondExpansions, and Vmath::Vcopy().

1416  {
1417  int n, cnt, nq;
1418  int offsetOld, offsetNew;
1419  std::vector<unsigned int> eIDs;
1420 
1421  Array<OneD, int> ElmtID,EdgeID;
1422  GetBoundaryToElmtMap(ElmtID,EdgeID);
1423 
1424  // Skip other boundary regions
1425  for (cnt = n = 0; n < i; ++n)
1426  {
1427  cnt += m_bndCondExpansions[n]->GetExpSize();
1428  }
1429 
1430  // Populate eIDs with information from BoundaryToElmtMap
1431  for (n = 0; n < m_bndCondExpansions[i]->GetExpSize(); ++n)
1432  {
1433  eIDs.push_back(ElmtID[cnt+n]);
1434  }
1435 
1436  // Create expansion list
1437  result =
1439  (*this, eIDs, DeclareCoeffPhysArrays);
1440 
1441  // Copy phys and coeffs to new explist
1442  if( DeclareCoeffPhysArrays)
1443  {
1444  Array<OneD, NekDouble> tmp1, tmp2;
1445  for (n = 0; n < result->GetExpSize(); ++n)
1446  {
1447  nq = GetExp(ElmtID[cnt+n])->GetTotPoints();
1448  offsetOld = GetPhys_Offset(ElmtID[cnt+n]);
1449  offsetNew = result->GetPhys_Offset(n);
1450  Vmath::Vcopy(nq, tmp1 = GetPhys()+ offsetOld, 1,
1451  tmp2 = result->UpdatePhys()+ offsetNew, 1);
1452 
1453  nq = GetExp(ElmtID[cnt+n])->GetNcoeffs();
1454  offsetOld = GetCoeff_Offset(ElmtID[cnt+n]);
1455  offsetNew = result->GetCoeff_Offset(n);
1456  Vmath::Vcopy(nq, tmp1 = GetCoeffs()+ offsetOld, 1,
1457  tmp2 = result->UpdateCoeffs()+ offsetNew, 1);
1458  }
1459  }
1460  }
const Array< OneD, const NekDouble > & GetCoeffs() const
This function returns (a reference to) the array (implemented as m_coeffs) containing all local expa...
Definition: ExpList.h:1946
int GetCoeff_Offset(int n) const
Get the start offset position for a global list of m_coeffs correspoinding to element n...
Definition: ExpList.h:2084
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
int GetPhys_Offset(int n) const
Get the start offset position for a global list of m_phys correspoinding to element n...
Definition: ExpList.h:2092
const boost::shared_ptr< LocalRegions::ExpansionVector > GetExp() const
This function returns the vector of elements in the expansion.
Definition: ExpList.h:2075
void GetBoundaryToElmtMap(Array< OneD, int > &ElmtID, Array< OneD, int > &EdgeID)
Definition: ExpList.h:2302
const Array< OneD, const NekDouble > & GetPhys() const
This function returns (a reference to) the array (implemented as m_phys) containing the function ev...
Definition: ExpList.h:2045
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1061
void Nektar::MultiRegions::DisContField1D::v_GetBoundaryToElmtMap ( Array< OneD, int > &  ElmtID,
Array< OneD, int > &  VertID 
)
protectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1360 of file DisContField1D.cpp.

References ASSERTL1, Nektar::MultiRegions::ExpList::GetExpSize(), m_bndCondExpansions, and m_bndConditions.

1362  {
1363  map<int, int> VertGID;
1364  int i,n,id;
1365  int bid,cnt,Vid;
1366  int nbcs = m_bndConditions.num_elements();
1367 
1368  // make sure arrays are of sufficient length
1369  if (ElmtID.num_elements() != nbcs)
1370  {
1371  ElmtID = Array<OneD, int>(nbcs,-1);
1372  }
1373  else
1374  {
1375  fill(ElmtID.get(), ElmtID.get()+nbcs, -1);
1376  }
1377 
1378  if (VertID.num_elements() != nbcs)
1379  {
1380  VertID = Array<OneD, int>(nbcs);
1381  }
1382 
1383  // setup map of all global ids along boundary
1384  for (cnt = n = 0; n < m_bndCondExpansions.num_elements(); ++n)
1385  {
1386  Vid = m_bndCondExpansions[n]->GetExp(0)->GetGeom()->GetVertex(0)->GetVid();
1387  VertGID[Vid] = cnt++;
1388  }
1389 
1390  // Loop over elements and find verts that match;
1392  for(cnt = n = 0; n < GetExpSize(); ++n)
1393  {
1394  exp = (*m_exp)[n];
1395  for(i = 0; i < exp->GetNverts(); ++i)
1396  {
1397  id = exp->GetGeom()->GetVid(i);
1398 
1399  if (VertGID.count(id) > 0)
1400  {
1401  bid = VertGID.find(id)->second;
1402  ASSERTL1(ElmtID[bid] == -1,"Edge already set");
1403  ElmtID[bid] = n;
1404  VertID[bid] = i;
1405  cnt ++;
1406  }
1407  }
1408  }
1409 
1410  ASSERTL1(cnt == nbcs,"Failed to visit all boundary condtiions");
1411  }
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2054
boost::shared_ptr< Expansion > ExpansionSharedPtr
Definition: Expansion.h:68
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:228
void Nektar::MultiRegions::DisContField1D::v_GetFwdBwdTracePhys ( Array< OneD, NekDouble > &  Fwd,
Array< OneD, NekDouble > &  Bwd 
)
protectedvirtual

Generate the forward or backward state for each trace point.

Parameters
FwdForward state.
BwdBackward state.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 823 of file DisContField1D.cpp.

References Nektar::MultiRegions::ExpList::m_phys.

825  {
826  v_GetFwdBwdTracePhys(m_phys,Fwd,Bwd);
827  }
Array< OneD, NekDouble > m_phys
The global expansion evaluated at the quadrature points.
Definition: ExpList.h:1015
virtual void v_GetFwdBwdTracePhys(Array< OneD, NekDouble > &Fwd, Array< OneD, NekDouble > &Bwd)
void Nektar::MultiRegions::DisContField1D::v_GetFwdBwdTracePhys ( const Array< OneD, const NekDouble > &  field,
Array< OneD, NekDouble > &  Fwd,
Array< OneD, NekDouble > &  Bwd 
)
protectedvirtual

This method extracts the "forward" and "backward" trace data from the array field and puts the data into output vectors Fwd and Bwd.

We first define the convention which defines "forwards" and "backwards". First an association is made between the vertex of each element and its corresponding vertex in the trace space using the mapping m_traceMap. The element can either be left-adjacent or right-adjacent to this trace edge (see Expansion0D::GetLeftAdjacentElementExp). Boundary edges are never left-adjacent since elemental left-adjacency is populated first.

If the element is left-adjacent we extract the vertex trace data from field into the forward trace space Fwd; otherwise, we place it in the backwards trace space Bwd. In this way, we form a unique set of trace normals since these are always extracted from left-adjacent elements.

Parameters
fieldis a NekDouble array which contains the 1D data from which we wish to extract the backward and forward orientated trace/edge arrays.
FwdThe resulting forwards space.
BwdThe resulting backwards space.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 856 of file DisContField1D.cpp.

References ASSERTL0, Nektar::SpatialDomains::eDirichlet, Nektar::SpatialDomains::eNeumann, Nektar::SpatialDomains::eNotDefined, Nektar::SpatialDomains::ePeriodic, Nektar::SpatialDomains::eRobin, Nektar::MultiRegions::ExpList::GetExpSize(), Nektar::MultiRegions::ExpList::GetPhys(), Nektar::MultiRegions::ExpList::GetPhys_Offset(), m_bndCondExpansions, m_bndConditions, m_leftAdjacentVerts, m_periodicBwdCopy, m_periodicFwdCopy, m_trace, m_traceMap, and Vmath::Zero().

860  {
861  // Counter variables
862  int n, v;
863 
864  // Number of elements
865  int nElements = GetExpSize();
866 
867  // Initial index of each element
868  int phys_offset;
869 
870  Array<OneD, Array<OneD, LocalRegions::ExpansionSharedPtr> >
871  &elmtToTrace = m_traceMap->GetElmtToTrace();
872 
873  // Set forward and backard state to zero
874  Vmath::Zero(Fwd.num_elements(), Fwd, 1);
875  Vmath::Zero(Bwd.num_elements(), Bwd, 1);
876 
877  int cnt;
878 
879  // Loop on the elements
880  for (cnt = n = 0; n < nElements; ++n)
881  {
882  // Set the offset of each element
883  phys_offset = GetPhys_Offset(n);
884 
885  for(v = 0; v < 2; ++v, ++cnt)
886  {
887  int offset = m_trace->GetPhys_Offset(elmtToTrace[n][v]->GetElmtId());
888 
889  if (m_leftAdjacentVerts[cnt])
890  {
891  (*m_exp)[n]->GetVertexPhysVals(v, field + phys_offset,
892  Fwd[offset]);
893  }
894  else
895  {
896  (*m_exp)[n]->GetVertexPhysVals(v, field + phys_offset,
897  Bwd[offset]);
898  }
899  }
900  }
901 
902  // Fill boundary conditions into missing elements.
903  int id = 0;
904 
905  for(cnt = n = 0; n < m_bndCondExpansions.num_elements(); ++n)
906  {
907  if (m_bndConditions[n]->GetBoundaryConditionType() ==
909  {
910  id = m_trace->GetPhys_Offset(m_traceMap->GetBndCondTraceToGlobalTraceMap(cnt));
911  Bwd[id] = m_bndCondExpansions[n]->GetPhys()[0]; //this is not getting the correct value?
912  cnt++;
913  }
914  else if (m_bndConditions[n]->GetBoundaryConditionType() ==
916  m_bndConditions[n]->GetBoundaryConditionType() ==
918  {
919  ASSERTL0((m_bndCondExpansions[n]->GetPhys())[0]==0.0,
920  "Method not set up for non-zero Neumann "
921  "boundary condition");
922  id = m_trace->GetPhys_Offset(m_traceMap->GetBndCondTraceToGlobalTraceMap(cnt));
923  Bwd[id] = Fwd[id];
924 
925  cnt++;
926  }
927  else if (m_bndConditions[n]->GetBoundaryConditionType() ==
929  {
930  // Do nothing
931  }
932  else if (m_bndConditions[n]->GetBoundaryConditionType() !=
934  {
935  ASSERTL0(false,
936  "Method not set up for this boundary condition.");
937  }
938  }
939 
940  // Copy any periodic boundary conditions.
941  for (n = 0; n < m_periodicFwdCopy.size(); ++n)
942  {
943  Bwd[m_periodicBwdCopy[n]] = Fwd[m_periodicFwdCopy[n]];
944  }
945 
946  // Do parallel exchange for forwards/backwards spaces.
947  m_traceMap->UniversalTraceAssemble(Fwd);
948  m_traceMap->UniversalTraceAssemble(Bwd);
949 
950  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:198
int GetPhys_Offset(int n) const
Get the start offset position for a global list of m_phys correspoinding to element n...
Definition: ExpList.h:2092
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2054
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
const Array< OneD, const NekDouble > & GetPhys() const
This function returns (a reference to) the array (implemented as m_phys) containing the function ev...
Definition: ExpList.h:2045
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
ExpListSharedPtr m_trace
Trace space storage for points between elements.
std::vector< int > m_periodicFwdCopy
A vector indicating degress of freedom which need to be copied from forwards to backwards space in ca...
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:373
map< int, RobinBCInfoSharedPtr > Nektar::MultiRegions::DisContField1D::v_GetRobinBCInfo ( void  )
protectedvirtual

Search through the edge expansions and identify which ones have Robin/Mixed type boundary conditions. If find a Robin boundary then store the edge id of the boundary condition and the array of points of the physical space boundary condition which are hold the boundary condition primitive variable coefficient at the quatrature points

Returns
std map containing the robin boundary condition info using a key of the element id

There is a next member to allow for more than one Robin boundary condition per element

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1490 of file DisContField1D.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::SpatialDomains::eRobin, Nektar::MultiRegions::ExpList::GetBoundaryToElmtMap(), m_bndCondExpansions, and m_bndConditions.

1491  {
1492  int i;
1493  map<int, RobinBCInfoSharedPtr> returnval;
1494  Array<OneD, int> ElmtID,VertID;
1495  GetBoundaryToElmtMap(ElmtID,VertID);
1496 
1497  for (i = 0; i < m_bndCondExpansions.num_elements(); ++i)
1498  {
1499  if (m_bndConditions[i]->GetBoundaryConditionType() ==
1501  {
1502  int elmtid;
1503 
1504  Array<OneD, NekDouble> x0(1);
1505  Array<OneD, NekDouble> x1(1);
1506  Array<OneD, NekDouble> x2(1);
1507  Array<OneD, NekDouble> coeffphys(1);
1508 
1509  m_bndCondExpansions[i]->GetCoords(x0, x1, x2);
1510 
1511  coeffphys[0] = (boost::static_pointer_cast<SpatialDomains
1512  ::RobinBoundaryCondition>(m_bndConditions[i])
1513  ->m_robinPrimitiveCoeff).Evaluate(x0[0],x1[0],x2[0],0.0);
1514 
1515  RobinBCInfoSharedPtr rInfo =
1517  AllocateSharedPtr(VertID[i],coeffphys);
1518 
1519  elmtid = ElmtID[i];
1520  // make link list if necessary (not likely in
1521  // 1D but needed in 2D & 3D)
1522  if(returnval.count(elmtid) != 0)
1523  {
1524  rInfo->next = returnval.find(elmtid)->second;
1525  }
1526  returnval[elmtid] = rInfo;
1527  }
1528  }
1529 
1530  return returnval;
1531  }
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
boost::shared_ptr< RobinBCInfo > RobinBCInfoSharedPtr
void GetBoundaryToElmtMap(Array< OneD, int > &ElmtID, Array< OneD, int > &EdgeID)
Definition: ExpList.h:2302
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
virtual ExpListSharedPtr& Nektar::MultiRegions::DisContField1D::v_GetTrace ( )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 164 of file DisContField1D.h.

References m_trace.

165  {
166  return m_trace;
167  }
ExpListSharedPtr m_trace
Trace space storage for points between elements.
virtual AssemblyMapDGSharedPtr& Nektar::MultiRegions::DisContField1D::v_GetTraceMap ( void  )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 169 of file DisContField1D.h.

References m_traceMap.

170  {
171  return m_traceMap;
172  }
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
void Nektar::MultiRegions::DisContField1D::v_HelmSolve ( const Array< OneD, const NekDouble > &  inarray,
Array< OneD, NekDouble > &  outarray,
const FlagList flags,
const StdRegions::ConstFactorMap factors,
const StdRegions::VarCoeffMap varcoeff,
const Array< OneD, const NekDouble > &  dirForcing,
const bool  PhysSpaceForcing 
)
protectedvirtual

Solve the Helmholtz equation.

Reimplemented from Nektar::MultiRegions::ExpList.

Reimplemented in Nektar::MultiRegions::ContField1D.

Definition at line 1159 of file DisContField1D.cpp.

References Nektar::SpatialDomains::eDirichlet, Nektar::StdRegions::eHybridDGHelmBndLam, Nektar::StdRegions::eHybridDGLamToU, Nektar::StdRegions::eInvHybridDGHelmholtz, Nektar::eWrapper, Nektar::MultiRegions::ExpList::GetBlockMatrix(), Nektar::MultiRegions::ExpList::GetExpSize(), GetGlobalBndLinSys(), Nektar::MultiRegions::ExpList::IProductWRTBase(), m_bndCondExpansions, m_bndConditions, Nektar::MultiRegions::ExpList::m_coeff_offset, Nektar::MultiRegions::ExpList::m_ncoeffs, m_traceMap, Vmath::Neg(), Nektar::MultiRegions::NullAssemblyMapSharedPtr, Vmath::Smul(), Nektar::Transpose(), and Vmath::Zero().

1167  {
1168  int i,n,cnt,nbndry;
1169  int nexp = GetExpSize();
1170  Array<OneD,NekDouble> f(m_ncoeffs);
1171  DNekVec F(m_ncoeffs,f,eWrapper);
1172  Array<OneD,NekDouble> e_f, e_l;
1173 
1174  //----------------------------------
1175  // Setup RHS Inner product if required
1176  //----------------------------------
1177  if(PhysSpaceForcing)
1178  {
1179  IProductWRTBase(inarray,f);
1180  Vmath::Neg(m_ncoeffs,f,1);
1181  }
1182  else
1183  {
1184  Vmath::Smul(m_ncoeffs,-1.0,inarray,1,f,1);
1185  }
1186 
1187  //----------------------------------
1188  // Solve continuous Boundary System
1189  //----------------------------------
1190  int GloBndDofs = m_traceMap->GetNumGlobalBndCoeffs();
1191  int NumDirBCs = m_traceMap->GetNumLocalDirBndCoeffs();
1192  int e_ncoeffs,id;
1193 
1194  GlobalMatrixKey HDGLamToUKey(
1197  factors,
1198  varcoeff);
1199 
1200  const DNekScalBlkMatSharedPtr &HDGLamToU =
1201  GetBlockMatrix(HDGLamToUKey);
1202 
1203  // Retrieve global trace space storage, \Lambda, from trace expansion
1204  Array<OneD,NekDouble> BndSol = Array<OneD,NekDouble>
1205  (m_traceMap->GetNumLocalBndCoeffs());
1206 
1207 
1208  Array<OneD,NekDouble> BndRhs(GloBndDofs,0.0);
1209  // Zero trace space
1210  Vmath::Zero(GloBndDofs,BndSol,1);
1211 
1212  int LocBndCoeffs = m_traceMap->GetNumLocalBndCoeffs();
1213  Array<OneD, NekDouble> loc_lambda(LocBndCoeffs);
1214  DNekVec LocLambda(LocBndCoeffs,loc_lambda,eWrapper);
1215 
1216  //----------------------------------
1217  // Evaluate Trace Forcing
1218  //----------------------------------
1219  // Determing <u_lam,f> terms using HDGLamToU matrix
1220  for (cnt = n = 0; n < nexp; ++n)
1221  {
1222  nbndry = (*m_exp)[n]->NumDGBndryCoeffs();
1223 
1224  e_ncoeffs = (*m_exp)[n]->GetNcoeffs();
1225  e_f = f+m_coeff_offset[n];
1226  e_l = loc_lambda + cnt;
1227 
1228  // use outarray as tmp space
1229  DNekVec Floc (nbndry, e_l, eWrapper);
1230  DNekVec ElmtFce (e_ncoeffs, e_f, eWrapper);
1231  Floc = Transpose(*(HDGLamToU->GetBlock(n,n)))*ElmtFce;
1232 
1233  cnt += nbndry;
1234  }
1235 
1236  // Assemble into global operator
1237  m_traceMap->AssembleBnd(loc_lambda,BndRhs);
1238 
1239  cnt = 0;
1240  // Copy Dirichlet boundary conditions into trace space
1241  for (i = 0; i < m_bndCondExpansions.num_elements(); ++i)
1242  {
1243  if (m_bndConditions[i]->GetBoundaryConditionType() ==
1245  {
1246  id = m_traceMap->GetBndCondCoeffsToGlobalCoeffsMap(i);
1247  BndSol[id] = m_bndCondExpansions[i]->GetCoeff(0);
1248  }
1249  else
1250  {
1251  id = m_traceMap->GetBndCondCoeffsToGlobalCoeffsMap(i);
1252  BndRhs[id] += m_bndCondExpansions[i]->GetCoeff(0);
1253  }
1254  }
1255 
1256  //----------------------------------
1257  // Solve trace problem
1258  //----------------------------------
1259  if (GloBndDofs - NumDirBCs > 0)
1260  {
1261  GlobalLinSysKey key(StdRegions::eHybridDGHelmBndLam,
1262  m_traceMap,factors);
1264  LinSys->Solve(BndRhs,BndSol,m_traceMap);
1265  }
1266 
1267  //----------------------------------
1268  // Internal element solves
1269  //----------------------------------
1270  GlobalMatrixKey invHDGhelmkey(StdRegions::eInvHybridDGHelmholtz,
1272  factors);
1273 
1274  const DNekScalBlkMatSharedPtr& InvHDGHelm =
1275  GetBlockMatrix(invHDGhelmkey);
1276  DNekVec out(m_ncoeffs,outarray,eWrapper);
1277  Vmath::Zero(m_ncoeffs,outarray,1);
1278 
1279  // get local trace solution from BndSol
1280  m_traceMap->GlobalToLocalBnd(BndSol,loc_lambda);
1281 
1282  // out = u_f + u_lam = (*InvHDGHelm)*f + (LamtoU)*Lam
1283  out = (*InvHDGHelm)*F + (*HDGLamToU)*LocLambda;
1284  }
const DNekScalBlkMatSharedPtr & GetBlockMatrix(const GlobalMatrixKey &gkey)
Definition: ExpList.cpp:890
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...
int GetExpSize(void)
This function returns the number of elements in the expansion.
Definition: ExpList.h:2054
GlobalLinSysSharedPtr GetGlobalBndLinSys(const GlobalLinSysKey &mkey)
For a given key, returns the associated global linear system.
Array< OneD, int > m_coeff_offset
Offset of elemental data into the array m_coeffs.
Definition: ExpList.h:1047
void Smul(int n, const T alpha, const T *x, const int incx, T *y, const int incy)
Scalar multiply y = alpha*y.
Definition: Vmath.cpp:213
int m_ncoeffs
The total number of local degrees of freedom. m_ncoeffs .
Definition: ExpList.h:976
boost::shared_ptr< DNekScalBlkMat > DNekScalBlkMatSharedPtr
Definition: NekTypeDefs.hpp:74
NekMatrix< InnerMatrixType, BlockMatrixTag > Transpose(NekMatrix< InnerMatrixType, BlockMatrixTag > &rhs)
AssemblyMapDGSharedPtr m_traceMap
Local to global DG mapping for trace space.
void Neg(int n, T *x, const int incx)
Negate x = -x.
Definition: Vmath.cpp:396
static AssemblyMapSharedPtr NullAssemblyMapSharedPtr
Definition: AssemblyMap.h:55
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
NekVector< NekDouble > DNekVec
Definition: NekTypeDefs.hpp:49
void IProductWRTBase(const Array< OneD, const NekDouble > &inarray, Array< OneD, NekDouble > &outarray, CoeffState coeffstate=eLocal)
Definition: ExpList.h:1649
boost::shared_ptr< GlobalLinSys > GlobalLinSysSharedPtr
Pointer to a GlobalLinSys object.
Definition: GlobalLinSys.h:52
void Zero(int n, T *x, const int incx)
Zero vector.
Definition: Vmath.cpp:373
void Nektar::MultiRegions::DisContField1D::v_Reset ( )
protectedvirtual

Reset this field, so that geometry information can be updated.

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 1465 of file DisContField1D.cpp.

References m_bndCondExpansions, and Nektar::MultiRegions::ExpList::v_Reset().

1466  {
1467  ExpList::v_Reset();
1468 
1469  // Reset boundary condition expansions.
1470  for (int n = 0; n < m_bndCondExpansions.num_elements(); ++n)
1471  {
1472  m_bndCondExpansions[n]->Reset();
1473  }
1474  }
virtual void v_Reset()
Reset geometry information, metrics, matrix managers and geometry information.
Definition: ExpList.cpp:1513
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
virtual MultiRegions::ExpListSharedPtr& Nektar::MultiRegions::DisContField1D::v_UpdateBndCondExpansion ( int  i)
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 230 of file DisContField1D.h.

231  {
232  return m_bndCondExpansions[i];
233  }
Array< OneD, MultiRegions::ExpListSharedPtr > m_bndCondExpansions
Discretised boundary conditions.
virtual Array<OneD, SpatialDomains::BoundaryConditionShPtr>& Nektar::MultiRegions::DisContField1D::v_UpdateBndConditions ( )
inlineprotectedvirtual

Reimplemented from Nektar::MultiRegions::ExpList.

Definition at line 236 of file DisContField1D.h.

References m_bndConditions.

237  {
238  return m_bndConditions;
239  }
Array< OneD, SpatialDomains::BoundaryConditionShPtr > m_bndConditions
An array which contains the information about the boundary condition on the different boundary region...

Member Data Documentation

Array<OneD,MultiRegions::ExpListSharedPtr> Nektar::MultiRegions::DisContField1D::m_bndCondExpansions
protected
Array<OneD,SpatialDomains::BoundaryConditionShPtr> Nektar::MultiRegions::DisContField1D::m_bndConditions
protected
std::set<int> Nektar::MultiRegions::DisContField1D::m_boundaryVerts
protected

A set storing the global IDs of any boundary edges.

Definition at line 128 of file DisContField1D.h.

Referenced by IsLeftAdjacentVertex(), and SetUpDG().

GlobalLinSysMapShPtr Nektar::MultiRegions::DisContField1D::m_globalBndMat
protected

Global boundary matrix.

Definition at line 117 of file DisContField1D.h.

Referenced by GetGlobalBndLinSys(), and SetUpDG().

std::vector<bool> Nektar::MultiRegions::DisContField1D::m_leftAdjacentVerts
protected

Definition at line 150 of file DisContField1D.h.

Referenced by SetUpDG(), and v_GetFwdBwdTracePhys().

std::vector<bool> Nektar::MultiRegions::DisContField1D::m_negatedFluxNormal
private

Definition at line 270 of file DisContField1D.h.

Referenced by GetNegatedFluxNormal().

int Nektar::MultiRegions::DisContField1D::m_numDirBndCondExpansions
protected

The number of boundary segments on which Dirichlet boundary conditions are imposed.

Definition at line 102 of file DisContField1D.h.

std::vector<int> Nektar::MultiRegions::DisContField1D::m_periodicBwdCopy
protected

Definition at line 143 of file DisContField1D.h.

Referenced by SetUpDG(), and v_GetFwdBwdTracePhys().

std::vector<int> Nektar::MultiRegions::DisContField1D::m_periodicFwdCopy
protected

A vector indicating degress of freedom which need to be copied from forwards to backwards space in case of a periodic boundary condition.

Definition at line 142 of file DisContField1D.h.

Referenced by SetUpDG(), and v_GetFwdBwdTracePhys().

PeriodicMap Nektar::MultiRegions::DisContField1D::m_periodicVerts
protected

A map which identifies groups of periodic vertices.

Definition at line 134 of file DisContField1D.h.

Referenced by Nektar::MultiRegions::ContField1D::ContField1D(), FindPeriodicVertices(), IsLeftAdjacentVertex(), and SetUpDG().

ExpListSharedPtr Nektar::MultiRegions::DisContField1D::m_trace
protected

Trace space storage for points between elements.

Definition at line 120 of file DisContField1D.h.

Referenced by IsLeftAdjacentVertex(), SetUpDG(), v_ExtractTracePhys(), v_GetFwdBwdTracePhys(), and v_GetTrace().

AssemblyMapDGSharedPtr Nektar::MultiRegions::DisContField1D::m_traceMap
protected