Nektar++
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Classes | Typedefs | Enumerations | Functions | Variables
Nektar::MultiRegions Namespace Reference

Classes

struct  _PeriodicEntity
 
class  AssemblyMap
 Base class for constructing local to global mapping of degrees of freedom. More...
 
class  AssemblyMapCG
 Constructs mappings for the C0 scalar continuous Galerkin formulation. More...
 
class  AssemblyMapDG
 
class  BottomUpSubStructuredGraph
 
struct  cmpop
 
class  ContField1D
 Abstraction of a global continuous one-dimensional spectral/hp element expansion which approximates the solution of a set of partial differential equations. More...
 
class  ContField2D
 This class is the abstraction of a global continuous two- dimensional spectral/hp element expansion which approximates the solution of a set of partial differential equations. More...
 
class  ContField3D
 
class  ContField3DHomogeneous1D
 
class  ContField3DHomogeneous2D
 
class  DisContField1D
 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...
 
class  DisContField2D
 
class  DisContField3D
 
class  DisContField3DHomogeneous1D
 
class  DisContField3DHomogeneous2D
 
class  ExpList
 Base class for all multi-elemental spectral/hp expansions. More...
 
class  ExpList0D
 This class is the abstraction of a collection of zero-dimensional expansions which is merely a collection of points/values. More...
 
class  ExpList1D
 This class is the abstraction of a one-dimensional multi-elemental expansions which is merely a collection of local expansions. More...
 
class  ExpList1DHomogeneous2D
 Abstraction of a one-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpList2D
 Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpList2DHomogeneous1D
 Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpList3D
 Abstraction of a three-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpList3DHomogeneous1D
 Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpList3DHomogeneous2D
 Abstraction of a one-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpListHomogeneous1D
 Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  ExpListHomogeneous2D
 Abstraction of a two-dimensional multi-elemental expansion which is merely a collection of local expansions. More...
 
class  GlobalLinSys
 A global linear system. More...
 
class  GlobalLinSysDirect
 A global linear system. More...
 
class  GlobalLinSysDirectFull
 A global linear system. More...
 
class  GlobalLinSysDirectStaticCond
 A global linear system. More...
 
class  GlobalLinSysIterative
 A global linear system. More...
 
class  GlobalLinSysIterativeFull
 A global linear system. More...
 
class  GlobalLinSysIterativeStaticCond
 A global linear system. More...
 
class  GlobalLinSysKey
 Describe a linear system. More...
 
class  GlobalLinSysPETSc
 A PETSc global linear system. More...
 
class  GlobalLinSysPETScFull
 A global linear system. More...
 
class  GlobalLinSysPETScStaticCond
 A global linear system. More...
 
class  GlobalLinSysStaticCond
 A global linear system. More...
 
class  GlobalLinSysXxt
 
class  GlobalLinSysXxtFull
 A global linear system. More...
 
class  GlobalLinSysXxtStaticCond
 A global linear system. More...
 
class  GlobalMatrix
 Represents a matrix of all degrees of freedom. More...
 
class  GlobalMatrixKey
 Describes a matrix with ordering defined by a local to global map. More...
 
class  LocTraceToTraceMap
 A helper class to deal with trace operations in the discontinuous Galerkin code. More...
 
class  MultiLevelBisectedGraph
 
class  PatchMap
 
class  Preconditioner
 
class  PreconditionerBlock
 
class  PreconditionerDiagonal
 
class  PreconditionerLinear
 
class  PreconditionerLinearWithBlock
 
class  PreconditionerLinearWithDiag
 
class  PreconditionerLinearWithLowEnergy
 
class  PreconditionerLowEnergy
 
class  PreconditionerNull
 
struct  RobinBCInfo
 
class  SubGraph
 

Typedefs

typedef boost::shared_ptr
< AssemblyMap
AssemblyMapSharedPtr
 
typedef boost::shared_ptr
< AssemblyMapCG
AssemblyMapCGSharedPtr
 
typedef boost::tuple< int, int,
NekDouble
ExtraDirDof
 
typedef std::vector< std::map
< int, int > > 
DofGraph
 
typedef boost::shared_ptr
< AssemblyMapDG
AssemblyMapDGSharedPtr
 
typedef boost::shared_ptr
< ExpList
ExpListSharedPtr
 Shared pointer to an ExpList object. More...
 
typedef boost::tuple
< LibUtilities::PointsKey,
LibUtilities::PointsKey,
LibUtilities::PointsKey,
LibUtilities::PointsKey
TraceInterpPoints
 Map holding points distributions required for interpolation of local traces onto global trace in two and three dimensions. More...
 
typedef boost::shared_ptr
< LocTraceToTraceMap
LocTraceToTraceMapSharedPtr
 
typedef boost::shared_ptr
< ContField1D
ContField1DSharedPtr
 
typedef boost::shared_ptr
< ContField2D
ContField2DSharedPtr
 
typedef boost::shared_ptr
< ContField3D
ContField3DSharedPtr
 
typedef boost::shared_ptr
< ContField3DHomogeneous1D
ContField3DHomogeneous1DSharedPtr
 
typedef boost::shared_ptr
< ContField3DHomogeneous2D
ContField3DHomogeneous2DSharedPtr
 
typedef boost::shared_ptr
< DisContField1D
DisContField1DSharedPtr
 
typedef boost::shared_ptr
< DisContField2D
DisContField2DSharedPtr
 
typedef boost::shared_ptr
< DisContField3D
DisContField3DSharedPtr
 
typedef boost::shared_ptr
< DisContField3DHomogeneous1D
DisContField3DHomogeneous1DSharedPtr
 
typedef boost::shared_ptr
< DisContField3DHomogeneous2D
DisContField3DHomogeneous2DSharedPtr
 
typedef std::map
< GlobalMatrixKey,
DNekScalBlkMatSharedPtr
BlockMatrixMap
 A map between global matrix keys and their associated block matrices. More...
 
typedef boost::shared_ptr
< BlockMatrixMap
BlockMatrixMapShPtr
 A shared pointer to a BlockMatrixMap. More...
 
typedef boost::shared_ptr
< ExpList0D
ExpList0DSharedPtr
 Shared pointer to an ExpList0D object. More...
 
typedef std::vector
< ExpList0DSharedPtr
ExpList0DVector
 Vector of pointers to ExpList0D objects. More...
 
typedef std::vector
< ExpList0DSharedPtr >
::iterator 
ExpList0DVectorIter
 Iterator for the vector of ExpList0D pointers. More...
 
typedef boost::shared_ptr
< ExpList1D
ExpList1DSharedPtr
 Shared pointer to an ExpList1D object. More...
 
typedef std::vector
< ExpList1DSharedPtr
ExpList1DVector
 Vector of pointers to ExpList1D objects. More...
 
typedef std::vector
< ExpList1DSharedPtr >
::iterator 
ExpList1DVectorIter
 Iterator for the vector of ExpList1D pointers. More...
 
typedef boost::shared_ptr
< ExpList1DHomogeneous2D
ExpList1DHomogeneous2DSharedPtr
 Shared pointer to an ExpList1DHomogeneous2D object. More...
 
typedef std::vector
< ExpList1DHomogeneous2DSharedPtr
ExpList1DHomogeneous2DVector
 Vector of pointers to ExpList1DHomogeneous2D objects. More...
 
typedef std::vector
< ExpList1DHomogeneous2DSharedPtr >
::iterator 
ExpList1DHomogeneous2DVectorIter
 Iterator for the vector of ExpList1DHomogeneous2D pointers. More...
 
typedef boost::shared_ptr
< ExpList2D
ExpList2DSharedPtr
 Shared pointer to an ExpList2D object. More...
 
typedef std::vector
< ExpList2DSharedPtr
ExpList2DVector
 Vector of pointers to ExpList2D objects. More...
 
typedef std::vector
< ExpList2DSharedPtr >
::iterator 
ExpList2DVectorIter
 Iterator for the vector of ExpList2D pointers. More...
 
typedef boost::shared_ptr
< ExpList2DHomogeneous1D
ExpList2DHomogeneous1DSharedPtr
 Shared pointer to an ExpList2DHomogeneous1D object. More...
 
typedef std::vector
< ExpList2DHomogeneous1DSharedPtr
ExpList2DHomogeneous1DVector
 Vector of pointers to ExpList2DHomogeneous1D objects. More...
 
typedef std::vector
< ExpList2DHomogeneous1DSharedPtr >
::iterator 
ExpList2DHomogeneous1DVectorIter
 Iterator for the vector of ExpList2DHomogeneous1D pointers. More...
 
typedef boost::shared_ptr
< ExpList3D
ExpList3DSharedPtr
 Shared pointer to an ExpList3D object. More...
 
typedef std::vector
< ExpList3DSharedPtr
ExpList3DVector
 Vector of pointers to ExpList3D objects. More...
 
typedef std::vector
< ExpList3DSharedPtr >
::iterator 
ExpList3DVectorIter
 Iterator over an ExpList3DVector. More...
 
typedef boost::shared_ptr
< ExpList3DHomogeneous1D
ExpList3DHomogeneous1DSharedPtr
 Shared pointer to an ExpList3DHomogeneous1D object. More...
 
typedef std::vector
< ExpList3DHomogeneous1DSharedPtr
ExpList3DHomogeneous1DVector
 Vector of pointers to ExpList3DHomogeneous1D objects. More...
 
typedef std::vector
< ExpList3DHomogeneous1DSharedPtr >
::iterator 
ExpList3DHomogeneous1DVectorIter
 Iterator for the vector of ExpList3DHomogeneous1D pointers. More...
 
typedef boost::shared_ptr
< ExpList3DHomogeneous2D
ExpList3DHomogeneous2DSharedPtr
 Shared pointer to an ExpList3DHomogeneous2D object. More...
 
typedef std::vector
< ExpList3DHomogeneous2DSharedPtr
ExpList3DHomogeneous2DVector
 Vector of pointers to ExpList3DHomogeneous2D objects. More...
 
typedef std::vector
< ExpList3DHomogeneous2DSharedPtr >
::iterator 
ExpList3DHomogeneous2DVectorIter
 Iterator for the vector of ExpList3DHomogeneous2D pointers. More...
 
typedef std::map
< Homogeneous1DMatType,
DNekBlkMatSharedPtr
Homo1DBlockMatrixMap
 A map between homo matrix keys and their associated block matrices. More...
 
typedef boost::shared_ptr
< Homo1DBlockMatrixMap
Homo1DBlockMatrixMapShPtr
 A shared pointer to a BlockMatrixMap. More...
 
typedef boost::shared_ptr
< ExpListHomogeneous1D
ExpListHomogeneous1DSharedPtr
 Shared pointer to an ExpList3DHomogeneous1D object. More...
 
typedef std::vector
< ExpListHomogeneous1DSharedPtr
ExpListHomogeneous1DVector
 Vector of pointers to ExpList3DHomogeneous1D objects. More...
 
typedef std::vector
< ExpListHomogeneous1DSharedPtr >
::iterator 
ExpListHomogeneous1DVectorIter
 Iterator for the vector of ExpList3DHomogeneous1D pointers. More...
 
typedef std::map
< Homogeneous2DMatType,
DNekBlkMatSharedPtr
Homo2DBlockMatrixMap
 A map between homo matrix keys and their associated block matrices. More...
 
typedef boost::shared_ptr
< Homo2DBlockMatrixMap
Homo2DBlockMatrixMapShPtr
 A shared pointer to a BlockMatrixMap. More...
 
typedef boost::shared_ptr
< ExpListHomogeneous2D
ExpListHomogeneous2DSharedPtr
 Shared pointer to an ExpList3DHomogeneous2D object. More...
 
typedef std::vector
< ExpListHomogeneous2DSharedPtr
ExpListHomogeneous2DVector
 Vector of pointers to ExpList3DHomogeneous2D objects. More...
 
typedef std::vector
< ExpListHomogeneous2DSharedPtr >
::iterator 
ExpListHomogeneous2DVectorIter
 Iterator for the vector of ExpList3DHomogeneous2D pointers. More...
 
typedef boost::shared_ptr
< GlobalLinSys
GlobalLinSysSharedPtr
 Pointer to a GlobalLinSys object. More...
 
typedef std::map
< GlobalLinSysKey,
GlobalLinSysSharedPtr
GlobalLinSysMap
 Mapping between GlobalLinSys objects and their associated keys. More...
 
typedef boost::shared_ptr
< GlobalLinSysMap
GlobalLinSysMapShPtr
 Pointer to a GlobalLinSys/key map. More...
 
typedef boost::shared_ptr
< Preconditioner
PreconditionerSharedPtr
 
typedef
LibUtilities::NekFactory
< std::string, GlobalLinSys,
const GlobalLinSysKey &, const
boost::weak_ptr< ExpList >
&, const boost::shared_ptr
< AssemblyMap > & > 
GlobalLinSysFactory
 Datatype of the NekFactory used to instantiate classes derived from the EquationSystem class. More...
 
typedef boost::shared_ptr
< GlobalLinSysDirectStaticCond
GlobalLinSysDirectStaticCondSharedPtr
 
typedef boost::shared_ptr
< GlobalLinSysIterativeStaticCond
GlobalLinSysIterativeStaticCondSharedPtr
 
typedef boost::shared_ptr
< GlobalLinSysPETScStaticCond
GlobalLinSysPETScStaticCondSharedPtr
 
typedef boost::shared_ptr
< GlobalLinSysStaticCond
GlobalLinSysStaticCondSharedPtr
 
typedef boost::shared_ptr
< GlobalLinSysXxtStaticCond
GlobalLinSysXxtStaticCondSharedPtr
 
typedef boost::shared_ptr
< GlobalMatrix
GlobalMatrixSharedPtr
 Shared pointer to a GlobalMatrix object. More...
 
typedef std::map
< GlobalMatrixKey,
GlobalMatrixSharedPtr
GlobalMatrixMap
 Mapping from global matrix keys to global matrices. More...
 
typedef boost::shared_ptr
< GlobalMatrixMap
GlobalMatrixMapShPtr
 Shared pointer to a global matrix map. More...
 
typedef boost::shared_ptr
< GlobalMatrixKey
GlobalMatrixKeySharedPtr
 A pointer to a GlobalMatrixKey. More...
 
typedef std::vector
< SpatialDomains::BoundaryConditionType
BndTypesVector
 
typedef std::vector
< SpatialDomains::BoundaryConditionType >
::iterator 
BndTypesVectorIter
 
typedef boost::shared_ptr
< RobinBCInfo
RobinBCInfoSharedPtr
 
typedef struct
Nektar::MultiRegions::_PeriodicEntity 
PeriodicEntity
 
typedef std::map< int,
std::vector< PeriodicEntity > > 
PeriodicMap
 
typedef
LibUtilities::NekFactory
< std::string, Preconditioner,
const boost::shared_ptr
< GlobalLinSys > &, const
boost::shared_ptr< AssemblyMap > & > 
PreconFactory
 
typedef boost::shared_ptr
< PreconditionerBlock
PreconditionerBlockSharedPtr
 
typedef boost::shared_ptr
< PreconditionerDiagonal
PreconditionerDiagonalSharedPtr
 
typedef boost::shared_ptr
< PreconditionerNull
PreconditionerNullSharedPtr
 
typedef boost::shared_ptr
< PreconditionerLinear
PreconditionerLinearSharedPtr
 
typedef boost::shared_ptr
< PreconditionerLinearWithBlock
PreconditionerLinearWithBlockSharedPtr
 
typedef boost::shared_ptr
< PreconditionerLinearWithDiag
PreconditionerLinearWithDiagSharedPtr
 
typedef boost::shared_ptr
< PreconditionerLinearWithLowEnergy
PreconditionerLinearWithLowEnergySharedPtr
 
typedef boost::shared_ptr
< PreconditionerLowEnergy
PreconditionerLowEnergySharedPtr
 
typedef boost::shared_ptr
< BottomUpSubStructuredGraph
BottomUpSubStructuredGraphSharedPtr
 
typedef boost::shared_ptr
< SubGraph
SubGraphSharedPtr
 
typedef boost::shared_ptr
< MultiLevelBisectedGraph
MultiLevelBisectedGraphSharedPtr
 
typedef boost::shared_ptr
< PatchMap
PatchMapSharedPtr
 

Enumerations

enum  InterpLocTraceToTrace {
  eNoInterp, eInterpDir0, eInterpEndPtDir0, eInterpDir1,
  eInterpEndPtDir1, eInterpBothDirs, eInterpEndPtDir0InterpDir1
}
 
enum  Direction {
  eX, eY, eZ, eS,
  eN
}
 
enum  ExpansionType {
  e0D, e1D, e2D, e3DH1D,
  e3DH2D, e3D, eNoType
}
 
enum  Homogeneous1DMatType { eForwardsCoeffSpace1D, eBackwardsCoeffSpace1D, eForwardsPhysSpace1D, eBackwardsPhysSpace1D }
 
enum  Homogeneous2DMatType {
  eForwardsCoeffSpaceY1D, eForwardsCoeffSpaceZ1D, eBackwardsCoeffSpaceY1D, eBackwardsCoeffSpaceZ1D,
  eForwardsPhysSpaceY1D, eForwardsPhysSpaceZ1D, eBackwardsPhysSpaceY1D, eBackwardsPhysSpaceZ1D
}
 
enum  LocalMatrixStorageStrategy { eNoStrategy, eContiguous, eNonContiguous, eSparse }
 
enum  PETScMatMult { ePETScMatMultSparse, ePETScMatMultShell }
 Enumerator. More...
 
enum  CoeffState { eNotSet, eLocal, eHybrid, eGlobal }
 
enum  AdjacentTraceOrientation { eAdjacentEdgeIsForwards, eAdjacentEdgeIsBackwards }
 
enum  AdjacentFaceOrientation {
  eAdjacentFaceDir1FwdDir1_Dir2FwdDir2, eAdjacentFaceDir1FwdDir1_Dir2BwdDir2, eAdjacentFaceDir1BwdDir1_Dir2FwdDir2, eAdjacentFaceDir1BwdDir1_Dir2BwdDir2,
  eAdjacentFaceDir1FwdDir2_Dir2FwdDir1, eAdjacentFaceDir1FwdDir2_Dir2BwdDir1, eAdjacentFaceDir1BwdDir2_Dir2FwdDir1, eAdjacentFaceDir1BwdDir2_Dir2BwdDir1
}
 
enum  GlobalSysSolnType {
  eNoSolnType, eDirectFullMatrix, eDirectStaticCond, eDirectMultiLevelStaticCond,
  eIterativeFull, eIterativeStaticCond, eIterativeMultiLevelStaticCond, eXxtFullMatrix,
  eXxtStaticCond, eXxtMultiLevelStaticCond, ePETScFullMatrix, ePETScStaticCond,
  ePETScMultiLevelStaticCond, eSIZE_GlobalSysSolnType
}
 
enum  ProjectionType { eGalerkin, eDiscontinuous, eMixed_CG_Discontinuous }
 Type of Galerkin projection. More...
 
enum  PreconditionerType {
  eNull, eDiagonal, eLinearWithDiagonal, eLinear,
  eLowEnergy, eLinearWithLowEnergy, eBlock, eLinearWithBlock
}
 
enum  MatrixStorageType { eSmvBSR }
 
enum  LinearPreconSolver { eLinearPreconXxt, eLinearPreconPETSc }
 

Functions

int RoundNekDoubleToInt (NekDouble x)
 Rounds a double precision number to an integer. More...
 
void RoundNekDoubleToInt (const Array< OneD, const NekDouble > inarray, Array< OneD, int > outarray)
 Rounds an array of double precision numbers to integers. More...
 
pair< int,
StdRegions::Orientation
DeterminePeriodicEdgeOrientId (int meshEdgeId, StdRegions::Orientation edgeOrient, const vector< PeriodicEntity > &periodicEdges)
 Determine orientation of an edge to its periodic equivalents, as well as the ID of the representative edge. More...
 
StdRegions::Orientation DeterminePeriodicFaceOrient (StdRegions::Orientation faceOrient, StdRegions::Orientation perFaceOrient)
 Determine relative orientation between two faces. More...
 
StdRegions::StdExpansionSharedPtr GetStdExp (StdRegions::StdExpansionSharedPtr exp)
 
void AlignFace (const StdRegions::Orientation orient, const int nquad1, const int nquad2, const Array< OneD, const NekDouble > &in, Array< OneD, NekDouble > &out)
 Helper function to re-align face to a given orientation. More...
 
GlobalLinSysFactoryGetGlobalLinSysFactory ()
 
bool operator< (const GlobalLinSysKey &lhs, const GlobalLinSysKey &rhs)
 
std::ostream & operator<< (std::ostream &os, const GlobalLinSysKey &rhs)
 Writes information about the object to a given stream. More...
 
bool operator< (const GlobalMatrixKey &lhs, const GlobalMatrixKey &rhs)
 
std::ostream & operator<< (std::ostream &os, const GlobalMatrixKey &rhs)
 Writes statistics about the matrix key to an output stream. More...
 
PreconFactoryGetPreconFactory ()
 
bool SubGraphWithoutVerts (const SubGraphSharedPtr g)
 
void CuthillMckeeReordering (const BoostGraph &graph, Array< OneD, int > &perm, Array< OneD, int > &iperm)
 
void MultiLevelBisectionReordering (const BoostGraph &graph, Array< OneD, int > &perm, Array< OneD, int > &iperm, BottomUpSubStructuredGraphSharedPtr &substructgraph, std::set< int > partVerts, int mdswitch)
 
void NoReordering (const BoostGraph &graph, Array< OneD, int > &perm, Array< OneD, int > &iperm)
 

Variables

static AssemblyMapSharedPtr NullAssemblyMapSharedPtr
 
static std::map< int, int > NullIntIntMap
 
static const std::vector
< std::map< int, int > > 
NullVecIntIntMap
 
MultiRegions::Direction const DirCartesianMap []
 
static ExpList NullExpList
 An empty ExpList object. More...
 
static ExpListSharedPtr NullExpListSharedPtr
 
static const Array< OneD,
ExpListSharedPtr
NullExpListSharedPtrArray
 
static const Array< OneD,
ExpList0DSharedPtr
NullExpList0DSharedPtrArray
 Empty ExpList0D object. More...
 
static const Array< OneD,
ExpList1DSharedPtr
NullExpList1DSharedPtrArray
 Empty ExpList1D object. More...
 
static const Array< OneD,
ExpList2DSharedPtr
NullExpList2DSharedPtrArray
 Empty ExpList2D object. More...
 
const char *const LocalMatrixStorageStrategyMap []
 
const char *const GlobalSysSolnTypeMap []
 
const char *const PreconditionerTypeMap []
 
const char *const MatrixStorageTypeMap []
 
static PeriodicMap NullPeriodicMap
 
static PreconditionerSharedPtr NullPreconditionerSharedPtr
 

Typedef Documentation

Definition at line 52 of file AssemblyMapCG.h.

Definition at line 49 of file AssemblyMapDG.h.

Definition at line 53 of file AssemblyMap.h.

A map between global matrix keys and their associated block matrices.

Definition at line 95 of file ExpList.h.

A shared pointer to a BlockMatrixMap.

Definition at line 97 of file ExpList.h.

Definition at line 159 of file MultiRegions.hpp.

Definition at line 160 of file MultiRegions.hpp.

Definition at line 60 of file SubStructuredGraph.h.

Definition at line 237 of file ContField1D.h.

Definition at line 293 of file ContField2D.h.

Definition at line 103 of file ContField3DHomogeneous1D.h.

Definition at line 98 of file ContField3DHomogeneous2D.h.

Definition at line 191 of file ContField3D.h.

Definition at line 275 of file DisContField1D.h.

Definition at line 298 of file DisContField2D.h.

Definition at line 297 of file DisContField3DHomogeneous1D.h.

Definition at line 154 of file DisContField3DHomogeneous2D.h.

Definition at line 266 of file DisContField3D.h.

typedef std::vector<std::map<int, int> > Nektar::MultiRegions::DofGraph

Definition at line 56 of file AssemblyMapCG.h.

Shared pointer to an ExpList0D object.

Definition at line 54 of file ExpList0D.h.

Vector of pointers to ExpList0D objects.

Definition at line 59 of file ExpList0D.h.

Iterator for the vector of ExpList0D pointers.

Definition at line 61 of file ExpList0D.h.

Shared pointer to an ExpList1DHomogeneous2D object.

Definition at line 50 of file ExpList1DHomogeneous2D.h.

Vector of pointers to ExpList1DHomogeneous2D objects.

Definition at line 55 of file ExpList1DHomogeneous2D.h.

Iterator for the vector of ExpList1DHomogeneous2D pointers.

Definition at line 57 of file ExpList1DHomogeneous2D.h.

Shared pointer to an ExpList1D object.

Definition at line 50 of file ExpList1D.h.

Vector of pointers to ExpList1D objects.

Definition at line 55 of file ExpList1D.h.

Iterator for the vector of ExpList1D pointers.

Definition at line 57 of file ExpList1D.h.

Shared pointer to an ExpList2DHomogeneous1D object.

Definition at line 50 of file ExpList2DHomogeneous1D.h.

Vector of pointers to ExpList2DHomogeneous1D objects.

Definition at line 57 of file ExpList2DHomogeneous1D.h.

Iterator for the vector of ExpList2DHomogeneous1D pointers.

Definition at line 60 of file ExpList2DHomogeneous1D.h.

Shared pointer to an ExpList2D object.

Definition at line 49 of file ExpList2D.h.

Vector of pointers to ExpList2D objects.

Definition at line 54 of file ExpList2D.h.

Iterator for the vector of ExpList2D pointers.

Definition at line 56 of file ExpList2D.h.

Shared pointer to an ExpList3DHomogeneous1D object.

Definition at line 50 of file ExpList3DHomogeneous1D.h.

Vector of pointers to ExpList3DHomogeneous1D objects.

Definition at line 55 of file ExpList3DHomogeneous1D.h.

Iterator for the vector of ExpList3DHomogeneous1D pointers.

Definition at line 57 of file ExpList3DHomogeneous1D.h.

Shared pointer to an ExpList3DHomogeneous2D object.

Definition at line 50 of file ExpList3DHomogeneous2D.h.

Vector of pointers to ExpList3DHomogeneous2D objects.

Definition at line 55 of file ExpList3DHomogeneous2D.h.

Iterator for the vector of ExpList3DHomogeneous2D pointers.

Definition at line 57 of file ExpList3DHomogeneous2D.h.

Shared pointer to an ExpList3D object.

Definition at line 114 of file ExpList3D.h.

Vector of pointers to ExpList3D objects.

Definition at line 116 of file ExpList3D.h.

Iterator over an ExpList3DVector.

Definition at line 118 of file ExpList3D.h.

Shared pointer to an ExpList3DHomogeneous1D object.

Definition at line 66 of file ExpListHomogeneous1D.h.

Vector of pointers to ExpList3DHomogeneous1D objects.

Definition at line 71 of file ExpListHomogeneous1D.h.

Iterator for the vector of ExpList3DHomogeneous1D pointers.

Definition at line 73 of file ExpListHomogeneous1D.h.

Shared pointer to an ExpList3DHomogeneous2D object.

Definition at line 70 of file ExpListHomogeneous2D.h.

Vector of pointers to ExpList3DHomogeneous2D objects.

Definition at line 75 of file ExpListHomogeneous2D.h.

Iterator for the vector of ExpList3DHomogeneous2D pointers.

Definition at line 77 of file ExpListHomogeneous2D.h.

typedef boost::shared_ptr< ExpList > Nektar::MultiRegions::ExpListSharedPtr

Shared pointer to an ExpList object.

Definition at line 56 of file LocTraceToTraceMap.h.

typedef boost::tuple<int, int, NekDouble> Nektar::MultiRegions::ExtraDirDof

Definition at line 54 of file AssemblyMapCG.h.

Definition at line 48 of file GlobalLinSysDirectStaticCond.h.

typedef LibUtilities::NekFactory< std::string, GlobalLinSys, const GlobalLinSysKey&, const boost::weak_ptr<ExpList>&, const boost::shared_ptr<AssemblyMap>& > Nektar::MultiRegions::GlobalLinSysFactory

Datatype of the NekFactory used to instantiate classes derived from the EquationSystem class.

Definition at line 69 of file GlobalLinSys.h.

Definition at line 50 of file GlobalLinSysIterativeStaticCond.h.

Mapping between GlobalLinSys objects and their associated keys.

Definition at line 57 of file GlobalLinSys.h.

Pointer to a GlobalLinSys/key map.

Definition at line 59 of file GlobalLinSys.h.

Definition at line 49 of file GlobalLinSysPETScStaticCond.h.

Pointer to a GlobalLinSys object.

Definition at line 52 of file GlobalLinSys.h.

Definition at line 49 of file GlobalLinSysStaticCond.h.

Definition at line 52 of file GlobalLinSysXxtStaticCond.h.

A pointer to a GlobalMatrixKey.

Definition at line 113 of file GlobalMatrixKey.h.

Mapping from global matrix keys to global matrices.

Definition at line 91 of file GlobalMatrix.h.

Shared pointer to a global matrix map.

Definition at line 93 of file GlobalMatrix.h.

Shared pointer to a GlobalMatrix object.

Definition at line 89 of file GlobalMatrix.h.

A map between homo matrix keys and their associated block matrices.

Definition at line 61 of file ExpListHomogeneous1D.h.

A shared pointer to a BlockMatrixMap.

Definition at line 63 of file ExpListHomogeneous1D.h.

A map between homo matrix keys and their associated block matrices.

Definition at line 65 of file ExpListHomogeneous2D.h.

A shared pointer to a BlockMatrixMap.

Definition at line 67 of file ExpListHomogeneous2D.h.

Definition at line 273 of file LocTraceToTraceMap.h.

Definition at line 67 of file SubStructuredGraph.h.

Definition at line 69 of file SubStructuredGraph.h.

typedef std::map<int, std::vector<PeriodicEntity> > Nektar::MultiRegions::PeriodicMap

Definition at line 201 of file MultiRegions.hpp.

Definition at line 48 of file PreconditionerBlock.h.

Definition at line 47 of file PreconditionerDiagonal.h.

Definition at line 55 of file PreconditionerLinear.h.

Definition at line 46 of file PreconditionerLinearWithBlock.h.

Definition at line 48 of file PreconditionerLinearWithDiag.h.

Definition at line 46 of file PreconditionerLinearWithLowEnergy.h.

Definition at line 50 of file PreconditionerLowEnergy.h.

Definition at line 99 of file PreconditionerDiagonal.h.

Definition at line 62 of file GlobalLinSys.h.

typedef LibUtilities::NekFactory< std::string, Preconditioner, const boost::shared_ptr<GlobalLinSys>&, const boost::shared_ptr<AssemblyMap>& > Nektar::MultiRegions::PreconFactory

Definition at line 63 of file Preconditioner.h.

Definition at line 181 of file MultiRegions.hpp.

Definition at line 65 of file SubStructuredGraph.h.

Map holding points distributions required for interpolation of local traces onto global trace in two and three dimensions.

The tuple contains four PointsKey entries. All four are used in 3D, but for 2D only two are used. The first two denote the local elemental points distributions, the latter the target trace distributions. In the 2D case, unused trace directions are set to LibUtilities::eNoPointsType.

Definition at line 82 of file LocTraceToTraceMap.h.

Enumeration Type Documentation

Enumerator
eAdjacentEdgeIsForwards 
eAdjacentEdgeIsBackwards 

Definition at line 57 of file MultiRegions.hpp.

Enumerator
eNotSet 

No state set.

eLocal 

Local coefficients.

eHybrid 

Hybrid coefficients.

eGlobal 

Global coefficients.

Definition at line 47 of file MultiRegions.hpp.

48  {
49  eNotSet, ///< No state set
50  eLocal, ///< Local coefficients
51  eHybrid, ///< Hybrid coefficients
52  eGlobal, ///< Global coefficients
53  };
Local coefficients.
Hybrid coefficients.
Global coefficients.
Enumerator
eX 
eY 
eZ 
eS 
eN 

Definition at line 66 of file ExpList.h.

Enumerator
e0D 
e1D 
e2D 
e3DH1D 
e3DH2D 
e3D 
eNoType 

Definition at line 75 of file ExpList.h.

Enumerator
eNoSolnType 

No Solution type specified.

eDirectFullMatrix 
eDirectStaticCond 
eDirectMultiLevelStaticCond 
eIterativeFull 
eIterativeStaticCond 
eIterativeMultiLevelStaticCond 
eXxtFullMatrix 
eXxtStaticCond 
eXxtMultiLevelStaticCond 
ePETScFullMatrix 
ePETScStaticCond 
ePETScMultiLevelStaticCond 
eSIZE_GlobalSysSolnType 

Definition at line 77 of file MultiRegions.hpp.

Enumerator
eLinearPreconXxt 
eLinearPreconPETSc 

Definition at line 49 of file PreconditionerLinear.h.

Enumerator
eSmvBSR 

Definition at line 148 of file MultiRegions.hpp.

149  {
150  eSmvBSR
151  };

Enumerator.

Enumerator
ePETScMatMultSparse 
ePETScMatMultShell 

Definition at line 52 of file GlobalLinSysPETSc.h.

Enumerator
eNull 

No Solution type specified.

eDiagonal 
eLinearWithDiagonal 
eLinear 
eLowEnergy 
eLinearWithLowEnergy 
eBlock 
eLinearWithBlock 

Definition at line 121 of file MultiRegions.hpp.

Type of Galerkin projection.

Enumerator
eGalerkin 
eDiscontinuous 
eMixed_CG_Discontinuous 

Definition at line 114 of file MultiRegions.hpp.

Function Documentation

void Nektar::MultiRegions::AlignFace ( const StdRegions::Orientation  orient,
const int  nquad1,
const int  nquad2,
const Array< OneD, const NekDouble > &  in,
Array< OneD, NekDouble > &  out 
)

Helper function to re-align face to a given orientation.

Definition at line 967 of file ExpList2D.cpp.

References Nektar::StdRegions::eDir1BwdDir1_Dir2BwdDir2, Nektar::StdRegions::eDir1BwdDir1_Dir2FwdDir2, Nektar::StdRegions::eDir1BwdDir2_Dir2BwdDir1, Nektar::StdRegions::eDir1BwdDir2_Dir2FwdDir1, Nektar::StdRegions::eDir1FwdDir1_Dir2BwdDir2, Nektar::StdRegions::eDir1FwdDir2_Dir2BwdDir1, and Nektar::StdRegions::eDir1FwdDir2_Dir2FwdDir1.

Referenced by Nektar::MultiRegions::ExpList2D::v_GetNormals().

972  {
973  // Copy transpose.
974  if (orient == StdRegions::eDir1FwdDir2_Dir2FwdDir1 ||
978  {
979  for (int i = 0; i < nquad2; ++i)
980  {
981  for (int j = 0; j < nquad1; ++j)
982  {
983  out[i*nquad1 + j] = in[j*nquad2 + i];
984  }
985  }
986  }
987  else
988  {
989  for (int i = 0; i < nquad2; ++i)
990  {
991  for (int j = 0; j < nquad1; ++j)
992  {
993  out[i*nquad1 + j] = in[i*nquad1 + j];
994  }
995  }
996  }
997 
998  if (orient == StdRegions::eDir1BwdDir1_Dir2FwdDir2 ||
1002  {
1003  // Reverse x direction
1004  for (int i = 0; i < nquad2; ++i)
1005  {
1006  for (int j = 0; j < nquad1/2; ++j)
1007  {
1008  swap(out[i*nquad1 + j],
1009  out[i*nquad1 + nquad1-j-1]);
1010  }
1011  }
1012  }
1013 
1014  if (orient == StdRegions::eDir1FwdDir1_Dir2BwdDir2 ||
1018  {
1019  // Reverse y direction
1020  for (int j = 0; j < nquad1; ++j)
1021  {
1022  for (int i = 0; i < nquad2/2; ++i)
1023  {
1024  swap(out[i*nquad1 + j],
1025  out[(nquad2-i-1)*nquad1 + j]);
1026  }
1027  }
1028  }
1029  }
void Nektar::MultiRegions::CuthillMckeeReordering ( const BoostGraph &  graph,
Array< OneD, int > &  perm,
Array< OneD, int > &  iperm 
)

Definition at line 839 of file SubStructuredGraph.cpp.

References ASSERTL1.

Referenced by Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), and Nektar::MultiRegions::AssemblyMapCG::CreateGraph().

842  {
843  int nGraphVerts = boost::num_vertices(graph);
844 
845  ASSERTL1(perm. num_elements() >= nGraphVerts &&
846  iperm.num_elements() >= nGraphVerts,
847  "Non-matching dimensions");
848 
849  // Call boost::cuthill_mckee_ordering to reorder the graph-vertices
850  // using the reverse Cuthill-Mckee algorithm
851  std::vector<BoostVertex> reorderedVerts(nGraphVerts);
852  boost::cuthill_mckee_ordering(graph, reorderedVerts.rbegin());
853 
854  //copy the reordering to the Arrays perm and iperm
855  for(int i = 0; i < nGraphVerts; i++)
856  {
857  perm[i] = reorderedVerts[i];
858  iperm[ reorderedVerts[i] ] = i;
859  }
860  }
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:218
std::pair< int, StdRegions::Orientation > Nektar::MultiRegions::DeterminePeriodicEdgeOrientId ( int  meshEdgeId,
StdRegions::Orientation  edgeOrient,
const vector< PeriodicEntity > &  periodicEdges 
)

Determine orientation of an edge to its periodic equivalents, as well as the ID of the representative edge.

Since an edge may be periodic with more than one other edge (e.g. a periodic cube has sets of four periodic edges in each coordinate direction), we have to define a 'representative' edge. In this assembly map we define it to be the one with the minimum ID. This routine is set up to calculate the orientation of a given edge with ID meshEdgeId with respect to the edge ID.

Parameters
meshEdgeIdID of a periodic edge.
edgeOrientEdge orientation of meshEdgeId with respect to its parent element.
periodicEdgesThe map of all periodic edges.
Returns
Pair containing the ID of the periodic edge and the orientation of meshEdgeID with respect to this edge.

Definition at line 1904 of file AssemblyMapCG.cpp.

References Nektar::StdRegions::eBackwards, and Nektar::StdRegions::eForwards.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerCG(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), and Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap().

1908  {
1909  int minId = periodicEdges[0].id;
1910  int minIdK = 0;
1911  int k;
1912 
1913  for (k = 1; k < periodicEdges.size(); ++k)
1914  {
1915  if (periodicEdges[k].id < minId)
1916  {
1917  minId = min(minId, periodicEdges[k].id);
1918  minIdK = k;
1919  }
1920  }
1921 
1922  minId = min(minId, meshEdgeId);
1923 
1924  if (meshEdgeId != minId)
1925  {
1926  if (periodicEdges[minIdK].orient == StdRegions::eBackwards)
1927  {
1928  // Swap edge orientation
1929  edgeOrient = (edgeOrient == StdRegions::eForwards) ?
1931  }
1932  }
1933 
1934  return make_pair(minId, edgeOrient);
1935  }
StdRegions::Orientation Nektar::MultiRegions::DeterminePeriodicFaceOrient ( StdRegions::Orientation  faceOrient,
StdRegions::Orientation  perFaceOrient 
)

Determine relative orientation between two faces.

Given the orientation of a local element to its local face, defined as faceOrient, and perFaceOrient which states the alignment of one periodic face to the other global face, this routine determines the orientation that takes this local element face to the global/unique face.

Parameters
faceOrientOrientation of the face with respect to its parent element.
perFaceOrientOrientation of the representative/global face.
Returns
Orientation between the two faces.

Definition at line 1952 of file AssemblyMapCG.cpp.

References Nektar::StdRegions::eDir1FwdDir1_Dir2FwdDir2.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerCG(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), and Nektar::MultiRegions::PreconditionerLowEnergy::v_BuildPreconditioner().

1955  {
1956  StdRegions::Orientation returnval = faceOrient;
1957 
1958  if(perFaceOrient != StdRegions::eDir1FwdDir1_Dir2FwdDir2)
1959  {
1960  int tmp1 = (int)faceOrient - 5;
1961  int tmp2 = (int)perFaceOrient - 5;
1962 
1963  int flipDir1Map [8] = {2,3,0,1,6,7,4,5};
1964  int flipDir2Map [8] = {1,0,3,2,5,4,7,6};
1965  int transposeMap[8] = {4,5,6,7,0,2,1,3};
1966 
1967  // Transpose orientation
1968  if (tmp2 > 3)
1969  {
1970  tmp1 = transposeMap[tmp1];
1971  }
1972 
1973  // Reverse orientation in direction 1.
1974  if (tmp2 == 2 || tmp2 == 3 || tmp2 == 6 || tmp2 == 7)
1975  {
1976  tmp1 = flipDir1Map[tmp1];
1977  }
1978 
1979  // Reverse orientation in direction 2
1980  if (tmp2 % 2 == 1)
1981  {
1982  tmp1 = flipDir2Map[tmp1];
1983  }
1984 
1985  returnval = (StdRegions::Orientation)(tmp1+5);
1986  }
1987  return returnval;
1988  }
GlobalLinSysFactory & Nektar::MultiRegions::GetGlobalLinSysFactory ( )

Definition at line 204 of file GlobalLinSys.cpp.

Referenced by Nektar::MultiRegions::ExpList::GenGlobalBndLinSys(), and Nektar::MultiRegions::ExpList::GenGlobalLinSys().

205  {
206  typedef Loki::SingletonHolder<GlobalLinSysFactory,
207  Loki::CreateUsingNew,
208  Loki::NoDestroy,
209  Loki::SingleThreaded> Type;
210  return Type::Instance();
211  }
LibUtilities::NekFactory< std::string, GlobalLinSys, const GlobalLinSysKey &, const boost::weak_ptr< ExpList > &, const boost::shared_ptr< AssemblyMap > & > GlobalLinSysFactory
Datatype of the NekFactory used to instantiate classes derived from the EquationSystem class...
Definition: GlobalLinSys.h:69
PreconFactory & Nektar::MultiRegions::GetPreconFactory ( )

Definition at line 88 of file Preconditioner.cpp.

Referenced by Nektar::MultiRegions::GlobalLinSys::CreatePrecon(), Nektar::MultiRegions::PreconditionerLinearWithBlock::v_InitObject(), Nektar::MultiRegions::PreconditionerLinearWithLowEnergy::v_InitObject(), and Nektar::MultiRegions::PreconditionerLinearWithDiag::v_InitObject().

89  {
90  typedef Loki::SingletonHolder<PreconFactory,
91  Loki::CreateUsingNew,
92  Loki::NoDestroy,
93  Loki::SingleThreaded> Type;
94  return Type::Instance();
95  }
LibUtilities::NekFactory< std::string, Preconditioner, const boost::shared_ptr< GlobalLinSys > &, const boost::shared_ptr< AssemblyMap > & > PreconFactory
StdRegions::StdExpansionSharedPtr Nektar::MultiRegions::GetStdExp ( StdRegions::StdExpansionSharedPtr  exp)

Definition at line 2886 of file ExpList.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, Nektar::LibUtilities::eHexahedron, Nektar::LibUtilities::ePrism, Nektar::LibUtilities::ePyramid, Nektar::LibUtilities::eQuadrilateral, Nektar::LibUtilities::eSegment, Nektar::LibUtilities::eTetrahedron, and Nektar::LibUtilities::eTriangle.

2887  {
2888 
2890 
2891  switch(exp->DetShapeType())
2892  {
2894  stdExp = MemoryManager<StdRegions::StdSegExp>
2895  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey());
2896  break;
2898  {
2900  if((nexp = exp->as<StdRegions::StdNodalTriExp>()))
2901  {
2902  stdExp = MemoryManager<StdRegions::StdNodalTriExp>
2903  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2904  exp->GetBasis(1)->GetBasisKey(),
2905  nexp->GetNodalPointsKey().GetPointsType());
2906  }
2907  else
2908  {
2909  stdExp = MemoryManager<StdRegions::StdTriExp>
2910  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2911  exp->GetBasis(1)->GetBasisKey());
2912  }
2913  }
2914  break;
2916  stdExp = MemoryManager<StdRegions::StdQuadExp>
2917  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2918  exp->GetBasis(1)->GetBasisKey());
2919  break;
2921  stdExp = MemoryManager<StdRegions::StdTetExp>
2922  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2923  exp->GetBasis(1)->GetBasisKey(),
2924  exp->GetBasis(2)->GetBasisKey());
2925  break;
2927  stdExp = MemoryManager<StdRegions::StdPyrExp>
2928  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2929  exp->GetBasis(1)->GetBasisKey(),
2930  exp->GetBasis(2)->GetBasisKey());
2931  break;
2932  case LibUtilities::ePrism:
2933  stdExp = MemoryManager<StdRegions::StdPrismExp>
2934  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2935  exp->GetBasis(1)->GetBasisKey(),
2936  exp->GetBasis(2)->GetBasisKey());
2937  break;
2939  stdExp = MemoryManager<StdRegions::StdHexExp>
2940  ::AllocateSharedPtr(exp->GetBasis(0)->GetBasisKey(),
2941  exp->GetBasis(1)->GetBasisKey(),
2942  exp->GetBasis(2)->GetBasisKey());
2943  break;
2944  default:
2945  ASSERTL0(false,"Shape type not setup");
2946  break;
2947  }
2948 
2949  return stdExp;
2950  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:188
boost::shared_ptr< StdNodalTriExp > StdNodalTriExpSharedPtr
boost::shared_ptr< StdExpansion > StdExpansionSharedPtr
void Nektar::MultiRegions::MultiLevelBisectionReordering ( const BoostGraph &  graph,
Array< OneD, int > &  perm,
Array< OneD, int > &  iperm,
BottomUpSubStructuredGraphSharedPtr &  substructgraph,
std::set< int >  partVerts,
int  mdswitch 
)

Definition at line 862 of file SubStructuredGraph.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Metis::as_onmetis(), ASSERTL0, ASSERTL1, ErrorUtil::efatal, Nektar::iterator, NEKERROR, and Vmath::Vcopy().

Referenced by Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), and Nektar::MultiRegions::AssemblyMapCG::CreateGraph().

869  {
870  int nGraphVerts = boost::num_vertices(graph);
871  int nGraphEdges = boost::num_edges (graph);
872 
873  ASSERTL1(perm. num_elements() >= nGraphVerts &&
874  iperm.num_elements() >= nGraphVerts,
875  "Non-matching dimensions");
876 
877  // We will now use METIS to reorder the graph. For the purpose of
878  // multi-level static condensation, we will use a METIS routine that
879  // partitions the graph recursively using a multi-level bisection
880  // algorithm. The name of this routine is METIS_NodeND and it was
881  // originally designed to reorder the DOFs in a matrix in order to
882  // minimise the fill-in when applying a factorisation technique
883  // (such as Cholesky). However, this reordering of DOFs also seems
884  // to be perfectly suited in the context of multilevel
885  // substructering. Therefore, we will use this metis routine instead
886  // of the more well-known graph-partitioning routines. However, as
887  // the standard metis implementation of METIS_NodeND only gives the
888  // resulting re-ordering as an output, we we will use an modified
889  // version of this routine that also returns information about the
890  // structure of the multi-level bisected partitioning.
891 
892  // This modified implementation has been written by W. GAO and
893  // collaborators and it additionally returns the separator tree
894  // compared to the standard implementation.
895 
896  // The name of this modified routine AS_METIS_NodeND (where AS
897  // stands for automated substructering) More information can be
898  // found in the paper:
899  //
900  // W. Gao, S. Li Xiaoye, C. Yang and Z. Bai
901  // 'An implementation and evaluation of the AMLS method
902  // for sparse eigenvalue problems'
903  // - ACM Trans. Math. Softw. 34, 4, Article 20 (July 2008)
904 
905  if(nGraphEdges)
906  {
907  // Step 1: Convert boost graph to a graph in adjncy-list format
908  // as required by METIS
909  int acnt = 0;
910  int vcnt = 0;
911  int i, cnt;
912  int nPartition = partVerts.size();
913  int nNonPartition = nGraphVerts - partVerts.size();
914  BoostVertexIterator vertit, vertit_end;
915  BoostAdjacencyIterator adjvertit, adjvertit_end;
916  Array<OneD, int> xadj(nNonPartition+1,0);
917  Array<OneD, int> adjncy(2*nGraphEdges);
918  Array<OneD, int> initial_perm(nGraphVerts);
919  Array<OneD, int> iinitial_perm(nGraphVerts);
920  Array<OneD, int> perm_tmp (nNonPartition);
921  Array<OneD, int> iperm_tmp(nNonPartition);
922 
924 
925  // First reorder vertices so that partition nodes are at the
926  // end. This allows METIS to partition the interior nodes.
927  for (i = cnt = 0; i < nGraphVerts; ++i)
928  {
929  if (partVerts.count(i) == 0)
930  {
931  initial_perm [i] = cnt;
932  iinitial_perm[cnt++] = i;
933  }
934  }
935 
936  for (i = 0; i < nGraphVerts; ++i)
937  {
938  if (partVerts.count(i) > 0)
939  {
940  initial_perm [i] = cnt;
941  iinitial_perm[cnt++] = i;
942  }
943  }
944 
945  // Apply this reordering to the graph.
946  boost::property_map<BoostGraph, boost::vertex_index_t>::type
947  index = get(boost::vertex_index, graph);
948 
949  for (boost::tie(vertit, vertit_end) = boost::vertices(graph);
950  vertit != vertit_end; ++vertit)
951  {
952  if (partVerts.count(index[*vertit]) > 0)
953  {
954  continue;
955  }
956 
957  for (boost::tie(adjvertit, adjvertit_end) =
958  boost::adjacent_vertices(*vertit,graph);
959  adjvertit != adjvertit_end;
960  ++adjvertit)
961  {
962  if (partVerts.count(index[*adjvertit]) > 0)
963  {
964  continue;
965  }
966  adjncy[acnt++] = initial_perm[*adjvertit];
967  }
968  xadj[++vcnt] = acnt;
969  }
970 
971  // Step 2: use metis to reorder the dofs. We do not know on
972  // forehand the size of the separator tree that METIS will
973  // return, so we just assume a really big value and try with
974  // that
975  int sizeSeparatorTree = nGraphVerts*10;
976  Array<OneD,int> septreeTmp(sizeSeparatorTree,-1);
977 
978  // The separator tree returned by metis has the following
979  // structure: It is a one dimensional array and information per
980  // level is contained per 5 elements:
981  //
982  // m_septree[i*5 + 0]: the level of recursion (top-level = 1)
983  // m_septree[i*5 + 1]: is this substructure a left or right
984  // branch? 1 = left branch, 2 = right branch
985  // m_septree[i*5 + 2]: the number of 'interior' DOFs in left
986  // branch
987  // m_septree[i*5 + 3]: the number of 'interior' DOFs in right
988  // branch
989  // m_septree[i*5 + 4]: the number of 'boundary' DOFs
990 
991  // Now try to call Call METIS.
992  try
993  {
995  nNonPartition,xadj,adjncy,perm_tmp,iperm_tmp,
996  septreeTmp, mdswitch);
997  }
998  catch(...)
999  {
1001  "Error in calling metis (the size of the separator"
1002  " tree might not be sufficient)");
1003  }
1004 
1005  // Change permutations from METIS to account for initial offset.
1006  for (i = 0; i < nGraphVerts; ++i)
1007  {
1008  if (partVerts.count(i) == 0)
1009  {
1010  iperm[i] = iperm_tmp[initial_perm[i]];
1011  perm[iperm[i]] = i;
1012  }
1013  }
1014 
1015  for (i = nNonPartition, it = partVerts.begin();
1016  it != partVerts.end(); ++it, ++i)
1017  {
1018  perm [i] = *it;
1019  iperm[*it] = i;
1020  }
1021 
1022  for (i = 0; i < nGraphVerts; ++i)
1023  {
1024  ASSERTL0(perm[iperm[i]] == i,
1025  "Perm error " + boost::lexical_cast<std::string>(i));
1026  }
1027 
1028  // Post-process the separator tree
1029  int trueSizeSepTree = 0;
1030  for (i = 0; septreeTmp[i] != -1; i++)
1031  {
1032  trueSizeSepTree++;
1033  }
1034  Array<OneD,int> septree(trueSizeSepTree);
1035  Vmath::Vcopy(trueSizeSepTree,septreeTmp,1,septree,1);
1036 
1037  // Based upon the separator tree, where are going to set up an
1038  // object of the class BottomUpSubStructuredGraph. The
1039  // constructor will read the separatortree and will interprete
1040  // the information from a bottom-up point of view.
1041  substructgraph = MemoryManager<BottomUpSubStructuredGraph>::
1042  AllocateSharedPtr(septree, nPartition);
1043 
1044  // Important, we cannot simply use the ordering given by metis
1045  // as it does not order the different blocks as we would like
1046  // it. Therefore, we use following command to re-order them
1047  // again in the context of the bottom-up substructering. As a
1048  // result, we will now obtain an ordering where the interior
1049  // degrees of freedom of the first (=bottom) level will be
1050  // ordered last (block by block ofcoarse). The interior degrees
1051  // of freedom of the second level will be ordered second to
1052  // last, etc ... As a result, the boundary degrees of freedom of
1053  // the last level (i.e. the dofs that will have to solved
1054  // non-recursively) will be ordered first (after the Dirichlet
1055  // Dofs that is). (this way, we actually follow the same idea
1056  // and convention in the standard (non-multi-level) static
1057  // condensation approach)
1058  substructgraph->UpdateBottomUpReordering(perm,iperm);
1059  }
1060  else
1061  {
1062  // This is the very special case of a graph without any
1063  // connectivity i.e. a collection of vertices without any edges
1064  for(int i = 0; i < nGraphVerts; i++)
1065  {
1066  perm[i] = i;
1067  iperm[i] = i;
1068  }
1069  substructgraph = MemoryManager<BottomUpSubStructuredGraph>::
1070  AllocateSharedPtr(nGraphVerts);
1071  }
1072  }
#define ASSERTL0(condition, msg)
Definition: ErrorUtil.hpp:188
#define NEKERROR(type, msg)
Assert Level 0 – Fundamental assert which is used whether in FULLDEBUG, DEBUG or OPT compilation mod...
Definition: ErrorUtil.hpp:185
static void as_onmetis(int nVerts, Nektar::Array< Nektar::OneD, int > xadj, Nektar::Array< Nektar::OneD, int > adjncy, Nektar::Array< Nektar::OneD, int > perm, Nektar::Array< Nektar::OneD, int > iperm, Nektar::Array< Nektar::OneD, int > map, int mdswitch=1)
Definition: Metis.hpp:53
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:218
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
Definition: Vmath.cpp:1047
void Nektar::MultiRegions::NoReordering ( const BoostGraph &  graph,
Array< OneD, int > &  perm,
Array< OneD, int > &  iperm 
)

Definition at line 1074 of file SubStructuredGraph.cpp.

References ASSERTL1.

Referenced by Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), and Nektar::MultiRegions::AssemblyMapCG::CreateGraph().

1077  {
1078  int nGraphVerts = boost::num_vertices(graph);
1079 
1080  ASSERTL1(perm. num_elements() >= nGraphVerts &&
1081  iperm.num_elements() >= nGraphVerts,
1082  "Non-matching dimensions");
1083 
1084  for (int i = 0; i < nGraphVerts; i++)
1085  {
1086  perm [i] = i;
1087  iperm[i] = i;
1088  }
1089  }
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:218
bool Nektar::MultiRegions::operator< ( const GlobalMatrixKey &  lhs,
const GlobalMatrixKey &  rhs 
)

Definition at line 80 of file GlobalMatrixKey.cpp.

References Nektar::MultiRegions::GlobalMatrixKey::m_constFactors, Nektar::MultiRegions::GlobalMatrixKey::m_locToGloMap, Nektar::MultiRegions::GlobalMatrixKey::m_matrixType, Nektar::MultiRegions::GlobalMatrixKey::m_shapeType, and Nektar::MultiRegions::GlobalMatrixKey::m_varCoeffs.

81  {
82  if(lhs.m_matrixType < rhs.m_matrixType)
83  {
84  return true;
85  }
86 
87  if(lhs.m_matrixType > rhs.m_matrixType)
88  {
89  return false;
90  }
91 
92 
93  if(lhs.m_shapeType < rhs.m_shapeType)
94  {
95  return true;
96  }
97 
98 
99  if(lhs.m_shapeType > rhs.m_shapeType)
100  {
101  return false;
102  }
103 
104  if(lhs.m_constFactors.size() < rhs.m_constFactors.size())
105  {
106  return true;
107  }
108  else if(lhs.m_constFactors.size() > rhs.m_constFactors.size())
109  {
110  return false;
111  }
112  else
113  {
114  StdRegions::ConstFactorMap::const_iterator x, y;
115  for(x = lhs.m_constFactors.begin(), y = rhs.m_constFactors.begin();
116  x != lhs.m_constFactors.end(); ++x, ++y)
117  {
118  if (x->second < y->second)
119  {
120  return true;
121  }
122  if (x->second > y->second)
123  {
124  return false;
125  }
126  }
127  }
128 
129  if(lhs.m_varCoeffs.size() < rhs.m_varCoeffs.size())
130  {
131  return true;
132  }
133  else if(lhs.m_varCoeffs.size() > rhs.m_varCoeffs.size())
134  {
135  return false;
136  }
137 // else
138 // {
139 // StdRegions::VarCoeffMap::const_iterator x, y;
140 // for (x = lhs.m_varCoeffs.begin(), y = rhs.m_varCoeffs.begin();
141 // x != lhs.m_varCoeffs.end(); ++x, ++y)
142 // {
143 // if (x->second.get() < y->second.get())
144 // {
145 // return true;
146 // }
147 // if (x->second.get() > y->second.get())
148 // {
149 // return false;
150 // }
151 // }
152 // }
153 
154  if(!rhs.m_locToGloMap.get())
155  {
156  return false;
157  }
158  else if(!lhs.m_locToGloMap.get() && rhs.m_locToGloMap.get() )
159  {
160  return true;
161  }
162  if(lhs.m_locToGloMap->GetHash() < rhs.m_locToGloMap->GetHash())
163  {
164  return true;
165  }
166 
167  return false;
168  }
StandardMatrixTag & lhs
bool Nektar::MultiRegions::operator< ( const GlobalLinSysKey &  lhs,
const GlobalLinSysKey &  rhs 
)

Compares two GlobalLinSysKeys by comparing their solution types and matrix keys.

Parameters
lhsFirst operand.
rhsSecond operand.
Returns
true if the first operand is considered less than the second operand.

Definition at line 96 of file GlobalLinSysKey.cpp.

References Nektar::MultiRegions::GlobalLinSysKey::m_solnType.

97  {
98  if(lhs.m_solnType < rhs.m_solnType)
99  {
100  return true;
101  }
102 
103  if(lhs.m_solnType > rhs.m_solnType)
104  {
105  return false;
106  }
107 
108  return (*dynamic_cast<const GlobalMatrixKey*>(&lhs)
109  < *dynamic_cast<const GlobalMatrixKey*>(&rhs));
110  }
StandardMatrixTag & lhs
std::ostream & Nektar::MultiRegions::operator<< ( std::ostream &  os,
const GlobalLinSysKey &  rhs 
)

Writes information about the object to a given stream.

Writes the vital statistics of a global linear system to a stream.

Parameters
osOutput stream.
rhsGlobalLinSys object to use.
Returns
Reference to the output stream os.

Definition at line 119 of file GlobalLinSysKey.cpp.

References Nektar::StdRegions::ConstFactorTypeMap, Nektar::MultiRegions::GlobalMatrixKey::GetConstFactors(), Nektar::MultiRegions::GlobalLinSysKey::GetGlobalSysSolnType(), Nektar::MultiRegions::GlobalMatrixKey::GetMatrixType(), Nektar::MultiRegions::GlobalMatrixKey::GetNConstFactors(), Nektar::MultiRegions::GlobalMatrixKey::GetNVarCoeffs(), Nektar::MultiRegions::GlobalMatrixKey::GetShapeType(), GlobalSysSolnTypeMap, Nektar::StdRegions::MatrixTypeMap, and Nektar::LibUtilities::ShapeTypeMap.

120  {
121  os << "MatrixType: " << StdRegions::MatrixTypeMap[rhs.GetMatrixType()] << ", ShapeType: "
122  << LibUtilities::ShapeTypeMap[rhs.GetShapeType()]
123  << std::endl;
124  os << "Solution Type: "
125  << GlobalSysSolnTypeMap[rhs.GetGlobalSysSolnType()] << endl;
126  os << "Number of constants: " << rhs.GetNConstFactors() << endl;
127  StdRegions::ConstFactorMap::const_iterator x;
128  for (x = rhs.GetConstFactors().begin(); x != rhs.GetConstFactors().end(); ++x)
129  {
130  os << " Constant " << StdRegions::ConstFactorTypeMap[x->first]
131  << ": " << x->second << endl;
132  }
133  os << "Number of variable coefficients: "
134  << rhs.GetNVarCoeffs() << endl;
135 
136  return os;
137  }
const char *const ConstFactorTypeMap[]
Definition: StdRegions.hpp:241
const char *const ShapeTypeMap[]
Definition: ShapeType.hpp:66
const char *const GlobalSysSolnTypeMap[]
const char *const MatrixTypeMap[]
Definition: StdRegions.hpp:152
std::ostream & Nektar::MultiRegions::operator<< ( std::ostream &  os,
const GlobalMatrixKey &  rhs 
)

Writes statistics about the matrix key to an output stream.

Definition at line 170 of file GlobalMatrixKey.cpp.

References Nektar::StdRegions::ConstFactorTypeMap, Nektar::MultiRegions::GlobalMatrixKey::GetConstFactors(), Nektar::MultiRegions::GlobalMatrixKey::GetMatrixType(), Nektar::MultiRegions::GlobalMatrixKey::GetNConstFactors(), and Nektar::MultiRegions::GlobalMatrixKey::GetNVarCoeffs().

171  {
172  os << "MatrixType: " << rhs.GetMatrixType() << endl;
173  os << "Number of constants: " << rhs.GetNConstFactors() << endl;
174  StdRegions::ConstFactorMap::const_iterator x;
175  for(x = rhs.GetConstFactors().begin(); x != rhs.GetConstFactors().end(); ++x)
176  {
177  os << " Constant " << StdRegions::ConstFactorTypeMap[x->first]
178  << ": " << x->second << endl;
179  }
180  os << "Number of variable coefficients: "
181  << rhs.GetNVarCoeffs() << endl;
182 
183  return os;
184  }
const char *const ConstFactorTypeMap[]
Definition: StdRegions.hpp:241
int Nektar::MultiRegions::RoundNekDoubleToInt ( NekDouble  x)

Rounds a double precision number to an integer.

Definition at line 57 of file AssemblyMap.cpp.

Referenced by Nektar::MultiRegions::AssemblyMap::AssemblyMap().

58  {
59  return int(x > 0.0 ? x + 0.5 : x - 0.5);
60  }
void Nektar::MultiRegions::RoundNekDoubleToInt ( const Array< OneD, const NekDouble >  inarray,
Array< OneD, int >  outarray 
)

Rounds an array of double precision numbers to integers.

Definition at line 63 of file AssemblyMap.cpp.

References ASSERTL1.

64  {
65  int size = inarray.num_elements();
66  ASSERTL1(outarray.num_elements()>=size,"Array sizes not compatible");
67 
68  NekDouble x;
69  for(int i = 0; i < size; i++)
70  {
71  x = inarray[i];
72  outarray[i] = int(x > 0.0 ? x + 0.5 : x - 0.5);
73  }
74  }
double NekDouble
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...
Definition: ErrorUtil.hpp:218
bool Nektar::MultiRegions::SubGraphWithoutVerts ( const SubGraphSharedPtr  g)

Definition at line 817 of file SubStructuredGraph.cpp.

Referenced by Nektar::MultiRegions::BottomUpSubStructuredGraph::ExpandGraphWithVertexWeights().

818  {
819  return g->GetNverts() == 0;
820  }

Variable Documentation

MultiRegions::Direction const Nektar::MultiRegions::DirCartesianMap[]
Initial value:

Definition at line 86 of file ExpList.h.

Referenced by Nektar::Extrapolate::CalcOutflowBCs(), Nektar::GlobalMapping::MappingGeneral::CalculateChristoffel(), Nektar::GlobalMapping::MappingXYofXY::CalculateChristoffel(), Nektar::SolverUtils::FilterAeroForces::CalculateForcesMapping(), Nektar::GlobalMapping::MappingGeneral::CalculateMetricTerms(), Nektar::Extrapolate::CurlCurl(), Nektar::NonlinearSWE::DoOdeRhs(), Nektar::LinearSWE::DoOdeRhs(), main(), Nektar::VCSMapping::MappingAccelerationCorrection(), Nektar::Utilities::ProcessGrad::Process(), Nektar::CoupledLinearNS::SetUpCoupledMatrix(), Nektar::AlternateSkewAdvection::v_Advect(), Nektar::SkewSymmetricAdvection::v_Advect(), Nektar::NavierStokesAdvection::v_Advect(), Nektar::MappingExtrapolate::v_CalcNeumannPressureBCs(), Nektar::MappingExtrapolate::v_CorrectPressureBCs(), Nektar::GlobalMapping::Mapping::v_CurlCurlField(), Nektar::GlobalMapping::Mapping::v_Divergence(), Nektar::VCSMapping::v_DoInitialise(), Nektar::GlobalMapping::Mapping::v_DotGradJacobian(), Nektar::ShallowWaterSystem::v_InitObject(), Nektar::MultiRegions::ExpList::v_PhysDeriv(), Nektar::VCSMapping::v_SetUpPressureForcing(), Nektar::VelocityCorrectionScheme::v_SetUpPressureForcing(), Nektar::VCSMapping::v_SolvePressure(), Nektar::GlobalMapping::MappingXYofZ::v_UpdateGeomInfo(), Nektar::GlobalMapping::MappingXYofXY::v_UpdateGeomInfo(), Nektar::GlobalMapping::MappingXofXZ::v_UpdateGeomInfo(), Nektar::GlobalMapping::MappingXofZ::v_UpdateGeomInfo(), Nektar::GlobalMapping::Mapping::v_VelocityLaplacian(), and Nektar::SolverUtils::EquationSystem::WeakAdvectionDivergenceForm().

const char* const Nektar::MultiRegions::GlobalSysSolnTypeMap[]
Initial value:
=
{
"No Solution Type",
"DirectFull",
"DirectStaticCond",
"DirectMultiLevelStaticCond",
"IterativeFull",
"IterativeStaticCond",
"IterativeMultiLevelStaticCond",
"XxtFull",
"XxtStaticCond",
"XxtMultiLevelStaticCond",
"PETScFull",
"PETScStaticCond",
"PETScMultiLevelStaticCond"
}

Definition at line 96 of file MultiRegions.hpp.

Referenced by Nektar::MultiRegions::AssemblyMapCG::CreateGraph(), Nektar::MultiRegions::ExpList::GenGlobalBndLinSys(), Nektar::MultiRegions::ExpList::GenGlobalLinSys(), and operator<<().

const char* const Nektar::MultiRegions::LocalMatrixStorageStrategyMap[]
Initial value:
=
{
"Contiguous",
"Non-contiguous",
"Sparse"
}

Definition at line 63 of file GlobalLinSysIterativeStaticCond.h.

const char* const Nektar::MultiRegions::MatrixStorageTypeMap[]
Initial value:
=
{
"SmvBSR"
}

Definition at line 153 of file MultiRegions.hpp.

Referenced by Nektar::MultiRegions::GlobalMatrix::GlobalMatrix().

AssemblyMapSharedPtr Nektar::MultiRegions::NullAssemblyMapSharedPtr
static
ExpList Nektar::MultiRegions::NullExpList
static

An empty ExpList object.

Definition at line 1397 of file ExpList.h.

const Array<OneD, ExpList0DSharedPtr> Nektar::MultiRegions::NullExpList0DSharedPtrArray
static

Empty ExpList0D object.

Definition at line 107 of file ExpList0D.h.

const Array<OneD, ExpList1DSharedPtr> Nektar::MultiRegions::NullExpList1DSharedPtrArray
static

Empty ExpList1D object.

Definition at line 197 of file ExpList1D.h.

const Array<OneD, ExpList2DSharedPtr> Nektar::MultiRegions::NullExpList2DSharedPtrArray
static

Empty ExpList2D object.

Definition at line 161 of file ExpList2D.h.

ExpListSharedPtr Nektar::MultiRegions::NullExpListSharedPtr
static
const Array<OneD, ExpListSharedPtr> Nektar::MultiRegions::NullExpListSharedPtrArray
static

Definition at line 2179 of file ExpList.h.

std::map<int,int> Nektar::MultiRegions::NullIntIntMap
static

Definition at line 49 of file AssemblyMapCG.h.

PeriodicMap Nektar::MultiRegions::NullPeriodicMap
static

Definition at line 202 of file MultiRegions.hpp.

PreconditionerSharedPtr Nektar::MultiRegions::NullPreconditionerSharedPtr
static

Definition at line 59 of file Preconditioner.h.

Referenced by Nektar::LinearElasticSystem::v_DoSolve().

const std::vector<std::map<int,int> > Nektar::MultiRegions::NullVecIntIntMap
static

Definition at line 50 of file AssemblyMapCG.h.

const char* const Nektar::MultiRegions::PreconditionerTypeMap[]
Initial value:
=
{
"Null",
"Diagonal",
"FullLinearSpaceWithDiagonal",
"FullLinearSpace",
"LowEnergyBlock",
"FullLinearSpaceWithLowEnergyBlock",
"Block",
"FullLinearSpaceWithBlock"
}

Definition at line 133 of file MultiRegions.hpp.

Referenced by Nektar::MultiRegions::GlobalLinSys::CreatePrecon().