Nektar::MultiRegions::AssemblyMap Class Reference

Base class for constructing local to global mapping of degrees of freedom. More...

#include <AssemblyMap.h>

Inheritance diagram for Nektar::MultiRegions::AssemblyMap:

Public Member Functions
	AssemblyMap ()
	Default constructor. More...
	AssemblyMap (const LibUtilities::SessionReaderSharedPtr &pSession, const std::string variable="DefaultVar")
	Constructor with a communicator. More...
	AssemblyMap (AssemblyMap *oldLevelMap, const BottomUpSubStructuredGraphSharedPtr &multiLevelGraph)
	Constructor for next level in multi-level static condensation. More...
virtual	~AssemblyMap ()
	Destructor. More...
LibUtilities::CommSharedPtr	GetComm ()
	Retrieves the communicator. More...
size_t	GetHash () const
	Retrieves the hash of this map. More...
int	GetLocalToGlobalMap (const int i) const
int	GetGlobalToUniversalMap (const int i) const
int	GetGlobalToUniversalMapUnique (const int i) const
const Array< OneD, const int > &	GetLocalToGlobalMap ()
const Array< OneD, const int > &	GetGlobalToUniversalMap ()
const Array< OneD, const int > &	GetGlobalToUniversalMapUnique ()
NekDouble	GetLocalToGlobalSign (const int i) const
const Array< OneD, NekDouble > &	GetLocalToGlobalSign () const
void	LocalToGlobal (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global, bool useComm=true) const
void	LocalToGlobal (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global, bool useComm=true) const
void	GlobalToLocal (const Array< OneD, const NekDouble > &global, Array< OneD, NekDouble > &loc) const
void	GlobalToLocal (const NekVector< NekDouble > &global, NekVector< NekDouble > &loc) const
void	Assemble (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global) const
void	Assemble (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global) const
void	UniversalAssemble (Array< OneD, NekDouble > &pGlobal) const
void	UniversalAssemble (NekVector< NekDouble > &pGlobal) const
void	UniversalAssemble (Array< OneD, NekDouble > &pGlobal, int offset) const
int	GetLocalToGlobalBndMap (const int i) const
	Retrieve the global index of a given local boundary mode. More...
const Array< OneD, const int > &	GetLocalToGlobalBndMap ()
	Retrieve the global indices of the local boundary modes. More...
const Array< OneD, const int > &	GetGlobalToUniversalBndMap ()
const Array< OneD, const int > &	GetGlobalToUniversalBndMapUnique ()
bool	GetSignChange ()
	Returns true if using a modal expansion requiring a change of sign of some modes. More...
NekDouble	GetLocalToGlobalBndSign (const int i) const
	Retrieve the sign change of a given local boundary mode. More...
Array< OneD, const NekDouble >	GetLocalToGlobalBndSign () const
	Retrieve the sign change for all local boundary modes. More...
int	GetBndCondCoeffsToGlobalCoeffsMap (const int i)
	Retrieves the global index corresponding to a boundary expansion mode. More...
const Array< OneD, const int > &	GetBndCondCoeffsToGlobalCoeffsMap ()
	Retrieves the global indices corresponding to the boundary expansion modes. More...
NekDouble	GetBndCondCoeffsToGlobalCoeffsSign (const int i)
	Returns the modal sign associated with a given boundary expansion mode. More...
int	GetBndCondTraceToGlobalTraceMap (const int i)
	Returns the global index of the boundary trace giving the index on the boundary expansion. More...
const Array< OneD, const int > &	GetBndCondTraceToGlobalTraceMap ()
int	GetNumGlobalDirBndCoeffs () const
	Returns the number of global Dirichlet boundary coefficients. More...
int	GetNumLocalDirBndCoeffs () const
	Returns the number of local Dirichlet boundary coefficients. More...
int	GetNumGlobalBndCoeffs () const
	Returns the total number of global boundary coefficients. More...
int	GetNumLocalBndCoeffs () const
	Returns the total number of local boundary coefficients. More...
int	GetNumLocalCoeffs () const
	Returns the total number of local coefficients. More...
int	GetNumGlobalCoeffs () const
	Returns the total number of global coefficients. More...
bool	GetSingularSystem () const
	Retrieves if the system is singular (true) or not (false) More...
void	GlobalToLocalBnd (const NekVector< NekDouble > &global, NekVector< NekDouble > &loc, int offset) const
void	GlobalToLocalBnd (const NekVector< NekDouble > &global, NekVector< NekDouble > &loc) const
void	GlobalToLocalBnd (const Array< OneD, const NekDouble > &global, Array< OneD, NekDouble > &loc, int offset) const
void	GlobalToLocalBnd (const Array< OneD, const NekDouble > &global, Array< OneD, NekDouble > &loc) const
void	LocalBndToGlobal (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global, int offset) const
void	LocalBndToGlobal (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global) const
void	LocalBndToGlobal (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global, int offset) const
void	LocalBndToGlobal (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global) const
void	AssembleBnd (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global, int offset) const
void	AssembleBnd (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global) const
void	AssembleBnd (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global, int offset) const
void	AssembleBnd (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global) const
void	UniversalAssembleBnd (Array< OneD, NekDouble > &pGlobal) const
void	UniversalAssembleBnd (NekVector< NekDouble > &pGlobal) const
void	UniversalAssembleBnd (Array< OneD, NekDouble > &pGlobal, int offset) const
int	GetFullSystemBandWidth () const
int	GetNumNonDirVertexModes () const
int	GetNumNonDirEdgeModes () const
int	GetNumNonDirFaceModes () const
int	GetNumDirEdges () const
int	GetNumDirFaces () const
int	GetNumNonDirEdges () const
int	GetNumNonDirFaces () const
void	PrintStats (std::ostream &out, std::string variable, bool printHeader=true) const
const Array< OneD, const int > &	GetExtraDirEdges ()
std::shared_ptr< AssemblyMap >	LinearSpaceMap (const ExpList &locexp, GlobalSysSolnType solnType)
int	GetBndSystemBandWidth () const
	Returns the bandwidth of the boundary system. More...
int	GetStaticCondLevel () const
	Returns the level of static condensation for this map. More...
int	GetNumPatches () const
	Returns the number of patches in this static condensation level. More...
const Array< OneD, const unsigned int > &	GetNumLocalBndCoeffsPerPatch ()
	Returns the number of local boundary coefficients in each patch. More...
const Array< OneD, const unsigned int > &	GetNumLocalIntCoeffsPerPatch ()
	Returns the number of local interior coefficients in each patch. More...
const AssemblyMapSharedPtr	GetNextLevelLocalToGlobalMap () const
	Returns the local to global mapping for the next level in the multi-level static condensation. More...
void	SetNextLevelLocalToGlobalMap (AssemblyMapSharedPtr pNextLevelLocalToGlobalMap)
const PatchMapSharedPtr &	GetPatchMapFromPrevLevel (void) const
	Returns the patch map from the previous level of the multi-level static condensation. More...
bool	AtLastLevel () const
	Returns true if this is the last level in the multi-level static condensation. More...
GlobalSysSolnType	GetGlobalSysSolnType () const
	Returns the method of solving global systems. More...
PreconditionerType	GetPreconType () const
NekDouble	GetIterativeTolerance () const
int	GetMaxIterations () const
int	GetSuccessiveRHS () const
int	GetLowestStaticCondLevel () const

Protected Member Functions
void	CalculateBndSystemBandWidth ()
	Calculates the bandwidth of the boundary system. More...
void	GlobalToLocalBndWithoutSign (const Array< OneD, const NekDouble > &global, Array< OneD, NekDouble > &loc)

Protected Attributes
LibUtilities::SessionReaderSharedPtr	m_session
	Session object. More...
LibUtilities::CommSharedPtr	m_comm
	Communicator. More...
size_t	m_hash
	Hash for map. More...
int	m_numLocalBndCoeffs
	Number of local boundary coefficients. More...
int	m_numGlobalBndCoeffs
	Total number of global boundary coefficients. More...
int	m_numLocalDirBndCoeffs
	Number of Local Dirichlet Boundary Coefficients. More...
int	m_numGlobalDirBndCoeffs
	Number of Global Dirichlet Boundary Coefficients. More...
bool	m_systemSingular
	Flag indicating if the system is singular or not. More...
int	m_numLocalCoeffs
	Total number of local coefficients. More...
int	m_numGlobalCoeffs
	Total number of global coefficients. More...
bool	m_signChange
	Flag indicating if modes require sign reversal. More...
Array< OneD, int >	m_localToGlobalBndMap
	Integer map of local boundary coeffs to global space. More...
Array< OneD, NekDouble >	m_localToGlobalBndSign
	Integer sign of local boundary coeffs to global space. More...
Array< OneD, int >	m_bndCondCoeffsToGlobalCoeffsMap
	Integer map of bnd cond coeffs to global coefficients. More...
Array< OneD, NekDouble >	m_bndCondCoeffsToGlobalCoeffsSign
	Integer map of bnd cond coeffs to global coefficients. More...
Array< OneD, int >	m_bndCondTraceToGlobalTraceMap
	Integer map of bnd cond trace number to global trace number. More...
Array< OneD, int >	m_globalToUniversalBndMap
	Integer map of process coeffs to universal space. More...
Array< OneD, int >	m_globalToUniversalBndMapUnique
	Integer map of unique process coeffs to universal space (signed) More...
GlobalSysSolnType	m_solnType
	The solution type of the global system. More...
int	m_bndSystemBandWidth
	The bandwith of the global bnd system. More...
PreconditionerType	m_preconType
	Type type of preconditioner to use in iterative solver. More...
int	m_maxIterations
	Maximum iterations for iterative solver. More...
NekDouble	m_iterativeTolerance
	Tolerance for iterative solver. More...
int	m_successiveRHS
	sucessive RHS for iterative solver More...
Gs::gs_data *	m_gsh
Gs::gs_data *	m_bndGsh
int	m_staticCondLevel
	The level of recursion in the case of multi-level static condensation. More...
int	m_numPatches
	The number of patches (~elements) in the current level. More...
Array< OneD, unsigned int >	m_numLocalBndCoeffsPerPatch
	The number of bnd dofs per patch. More...
Array< OneD, unsigned int >	m_numLocalIntCoeffsPerPatch
	The number of int dofs per patch. More...
AssemblyMapSharedPtr	m_nextLevelLocalToGlobalMap
	Map from the patches of the previous level to the patches of the current level. More...
int	m_lowestStaticCondLevel
	Lowest static condensation level. More...

Private Member Functions
virtual int	v_GetLocalToGlobalMap (const int i) const
virtual int	v_GetGlobalToUniversalMap (const int i) const
virtual int	v_GetGlobalToUniversalMapUnique (const int i) const
virtual const Array< OneD, const int > &	v_GetLocalToGlobalMap ()
virtual const Array< OneD, const int > &	v_GetGlobalToUniversalMap ()
virtual const Array< OneD, const int > &	v_GetGlobalToUniversalMapUnique ()
virtual NekDouble	v_GetLocalToGlobalSign (const int i) const
virtual const Array< OneD, NekDouble > &	v_GetLocalToGlobalSign () const
virtual void	v_LocalToGlobal (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global, bool useComm) const
virtual void	v_LocalToGlobal (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global, bool useComm) const
virtual void	v_GlobalToLocal (const Array< OneD, const NekDouble > &global, Array< OneD, NekDouble > &loc) const
virtual void	v_GlobalToLocal (const NekVector< NekDouble > &global, NekVector< NekDouble > &loc) const
virtual void	v_Assemble (const Array< OneD, const NekDouble > &loc, Array< OneD, NekDouble > &global) const
virtual void	v_Assemble (const NekVector< NekDouble > &loc, NekVector< NekDouble > &global) const
virtual void	v_UniversalAssemble (Array< OneD, NekDouble > &pGlobal) const
virtual void	v_UniversalAssemble (NekVector< NekDouble > &pGlobal) const
virtual void	v_UniversalAssemble (Array< OneD, NekDouble > &pGlobal, int offset) const
virtual int	v_GetFullSystemBandWidth () const
virtual int	v_GetNumNonDirVertexModes () const
virtual int	v_GetNumNonDirEdgeModes () const
virtual int	v_GetNumNonDirFaceModes () const
virtual int	v_GetNumDirEdges () const
virtual int	v_GetNumDirFaces () const
virtual int	v_GetNumNonDirEdges () const
virtual int	v_GetNumNonDirFaces () const
virtual const Array< OneD, const int > &	v_GetExtraDirEdges ()
virtual std::shared_ptr< AssemblyMap >	v_LinearSpaceMap (const ExpList &locexp, GlobalSysSolnType solnType)
	Generate a linear space mapping from existing mapping. More...

Private Attributes
PatchMapSharedPtr	m_patchMapFromPrevLevel
	Mapping information for previous level in MultiLevel Solver. More...

Detailed Description

Base class for constructing local to global mapping of degrees of freedom.

This class acts as a base class for constructing mappings between local, global and boundary degrees of freedom. It holds the storage for the maps and provides the accessors needed to retrieve them.

There are two derived classes: AssemblyMapCG and AssemblyMapDG. These perform the actual construction of the maps within their specific contexts.

Definition at line 58 of file AssemblyMap.h.

Constructor & Destructor Documentation

AssemblyMap() [1/3]

Nektar::MultiRegions::AssemblyMap::AssemblyMap

(

)

Default constructor.

Initialises an empty mapping.

Definition at line 80 of file AssemblyMap.cpp.

                                :
            m_session(),
            m_comm(),
            m_hash(0),
            m_numLocalBndCoeffs(0),
            m_numGlobalBndCoeffs(0),
            m_numLocalDirBndCoeffs(0),
            m_numGlobalDirBndCoeffs(0),
            m_solnType(eNoSolnType),
            m_bndSystemBandWidth(0),
            m_successiveRHS(0),
            m_gsh(0),
            m_bndGsh(0)
        {
        }

AssemblyMap() [2/3]

Nektar::MultiRegions::AssemblyMap::AssemblyMap	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const std::string	variable = "DefaultVar"
	)

Constructor with a communicator.

Definition at line 96 of file AssemblyMap.cpp.

References Nektar::NekConstants::kNekIterativeTol, m_iterativeTolerance, m_maxIterations, m_preconType, m_solnType, and m_successiveRHS.

                                         :
            m_session(pSession),
            m_comm(pSession->GetComm()),
            m_hash(0),
            m_numLocalBndCoeffs(0),
            m_numGlobalBndCoeffs(0),
            m_numLocalDirBndCoeffs(0),
            m_numGlobalDirBndCoeffs(0),
            m_bndSystemBandWidth(0),
            m_successiveRHS(0),
            m_gsh(0),
            m_bndGsh(0)
        {
            // Default value from Solver Info
            m_solnType = pSession->GetSolverInfoAsEnum<GlobalSysSolnType>(
                                                            "GlobalSysSoln");
            m_preconType = pSession->GetSolverInfoAsEnum<PreconditionerType>(
                                                            "Preconditioner");

            // Override values with data from GlobalSysSolnInfo section 
            if(pSession->DefinesGlobalSysSolnInfo(variable, "GlobalSysSoln"))
            {
                std::string sysSoln = pSession->GetGlobalSysSolnInfo(variable,
                                                            "GlobalSysSoln");
                m_solnType = pSession->GetValueAsEnum<GlobalSysSolnType>(
                                                    "GlobalSysSoln", sysSoln);
            }

            if(pSession->DefinesGlobalSysSolnInfo(variable, "Preconditioner"))
            {
                std::string precon = pSession->GetGlobalSysSolnInfo(variable,
                                                            "Preconditioner");
                m_preconType = pSession->GetValueAsEnum<PreconditionerType>(
                                                    "Preconditioner", precon);
            }

            if(pSession->DefinesGlobalSysSolnInfo(variable,
                                                  "IterativeSolverTolerance"))
            {
                m_iterativeTolerance = boost::lexical_cast<NekDouble>(
                        pSession->GetGlobalSysSolnInfo(variable,
                                "IterativeSolverTolerance").c_str());
            }
            else
            {
                pSession->LoadParameter("IterativeSolverTolerance",
                                        m_iterativeTolerance,
                                        NekConstants::kNekIterativeTol);
            }


            if(pSession->DefinesGlobalSysSolnInfo(variable,
                                                  "MaxIterations"))
            {
                m_maxIterations = boost::lexical_cast<int>(
                        pSession->GetGlobalSysSolnInfo(variable,
                                "MaxIterations").c_str());
            }
            else
            {
                pSession->LoadParameter("MaxIterations",
                                        m_maxIterations,
                                        5000);
            }


            if(pSession->DefinesGlobalSysSolnInfo(variable,"SuccessiveRHS"))
            {
                m_successiveRHS = boost::lexical_cast<int>(
                        pSession->GetGlobalSysSolnInfo(variable,
                                "SuccessiveRHS").c_str());
            }
            else
            {
                pSession->LoadParameter("SuccessiveRHS",
                                        m_successiveRHS,0);
            }

        }

AssemblyMap() [3/3]

Nektar::MultiRegions::AssemblyMap::AssemblyMap	(	AssemblyMap *	oldLevelMap,
		const BottomUpSubStructuredGraphSharedPtr &	multiLevelGraph
	)

Constructor for next level in multi-level static condensation.

Create a new level of mapping using the information in multiLevelGraph and performing the following steps:

• STEP 1: setup a mask array to determine to which patch of the new level every patch of the current level belongs. To do so we make four arrays, #gloPatchMask, #globHomPatchMask, #locPatchMask_NekDouble and #locPatchMask. These arrays are then used to check which local dofs of the old level belong to which patch of the new level
• STEP 2: We calculate how many local bnd dofs of the old level belong to the boundaries of each patch at the new level. We need this information to set up the mapping between different levels.

Definition at line 182 of file AssemblyMap.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), ASSERTL0, ASSERTL1, CalculateBndSystemBandWidth(), Nektar::MultiRegions::eDirectMultiLevelStaticCond, Nektar::MultiRegions::eIterativeMultiLevelStaticCond, Nektar::MultiRegions::ePETScMultiLevelStaticCond, Nektar::MultiRegions::eXxtMultiLevelStaticCond, GetGlobalSysSolnType(), GetGlobalToUniversalBndMap(), GetGlobalToUniversalBndMapUnique(), GetLocalToGlobalBndMap(), GetLocalToGlobalBndSign(), GetNumGlobalBndCoeffs(), GetNumGlobalDirBndCoeffs(), GetNumLocalBndCoeffs(), GetNumLocalBndCoeffsPerPatch(), GetNumLocalDirBndCoeffs(), GetNumPatches(), GetSignChange(), GetStaticCondLevel(), GlobalToLocalBndWithoutSign(), m_globalToUniversalBndMap, m_globalToUniversalBndMapUnique, m_localToGlobalBndMap, m_localToGlobalBndSign, m_nextLevelLocalToGlobalMap, m_numGlobalBndCoeffs, m_numGlobalCoeffs, m_numGlobalDirBndCoeffs, m_numLocalBndCoeffs, m_numLocalBndCoeffsPerPatch, m_numLocalDirBndCoeffs, m_numLocalIntCoeffsPerPatch, m_numPatches, m_patchMapFromPrevLevel, m_signChange, m_solnType, m_staticCondLevel, Nektar::MultiRegions::RoundNekDoubleToInt(), and sign.

                                                                               :
            m_session(oldLevelMap->m_session),
            m_comm(oldLevelMap->GetComm()),
            m_hash(0),
            m_solnType(oldLevelMap->m_solnType),
            m_preconType(oldLevelMap->m_preconType),
            m_maxIterations(oldLevelMap->m_maxIterations),
            m_iterativeTolerance(oldLevelMap->m_iterativeTolerance),
            m_successiveRHS(oldLevelMap->m_successiveRHS),
            m_gsh(oldLevelMap->m_gsh),
            m_bndGsh(oldLevelMap->m_bndGsh),
            m_lowestStaticCondLevel(oldLevelMap->m_lowestStaticCondLevel)
        {
            int i;
            int j;
            int cnt;

            //--------------------------------------------------------------
            // -- Extract information from the input argument
            int numGlobalBndCoeffsOld    = oldLevelMap->GetNumGlobalBndCoeffs();
            int numGlobalDirBndCoeffsOld = oldLevelMap->GetNumGlobalDirBndCoeffs();
            int numLocalBndCoeffsOld     = oldLevelMap->GetNumLocalBndCoeffs();
            int numLocalDirBndCoeffsOld  = oldLevelMap->GetNumLocalDirBndCoeffs();
            bool signChangeOld           = oldLevelMap->GetSignChange();

            int staticCondLevelOld       = oldLevelMap->GetStaticCondLevel();
            int numPatchesOld            = oldLevelMap->GetNumPatches();
            GlobalSysSolnType solnTypeOld = oldLevelMap->GetGlobalSysSolnType();
            const Array<OneD, const unsigned int>& numLocalBndCoeffsPerPatchOld = oldLevelMap->GetNumLocalBndCoeffsPerPatch();
            //--------------------------------------------------------------

            //--------------------------------------------------------------
            int newLevel = staticCondLevelOld+1;
            /** - STEP 1: setup a mask array to determine to which patch
             *          of the new level every patch of the current
             *          level belongs.  To do so we make four arrays,
             *          #gloPatchMask, #globHomPatchMask,
             *          #locPatchMask_NekDouble and #locPatchMask.
             *          These arrays are then used to check which local
             *          dofs of the old level belong to which patch of
             *          the new level
             */
            Array<OneD, NekDouble> globPatchMask         (numGlobalBndCoeffsOld,-1.0);
            Array<OneD, NekDouble> globHomPatchMask      (globPatchMask+numGlobalDirBndCoeffsOld);
            Array<OneD, NekDouble> locPatchMask_NekDouble(numLocalBndCoeffsOld,-3.0);
            Array<OneD, int>       locPatchMask          (numLocalBndCoeffsOld);

            // Fill the array globPatchMask as follows:
            // - The first part (i.e. the glob bnd dofs) is filled with the
            //   value -1
            // - The second part (i.e. the glob interior dofs) is numbered
            //   according to the patch it belongs to (i.e. dofs in first block
            //   all are numbered 0, the second block numbered are 1, etc...)
            multiLevelGraph->MaskPatches(newLevel,globHomPatchMask);
            
            // Map from Global Dofs to Local Dofs
            // As a result, we know for each local dof whether
            // it is mapped to the boundary of the next level, or to which
            // patch. Based upon this, we can than later associate every patch
            // of the current level with a patch in the next level.
            oldLevelMap->GlobalToLocalBndWithoutSign(globPatchMask,locPatchMask_NekDouble);

            // Convert the result to an array of integers rather than doubles
            RoundNekDoubleToInt(locPatchMask_NekDouble,locPatchMask);

            /** - STEP 2: We calculate how many local bnd dofs of the
             *  old level belong to the boundaries of each patch at
             *  the new level. We need this information to set up the
             *  mapping between different levels.
             */

            // Retrieve the number of patches at the next level
            int numPatchesWithIntNew = multiLevelGraph->GetNpatchesWithInterior(newLevel);
            int numPatchesNew        = numPatchesWithIntNew;
            
            // Allocate memory to store the number of local dofs associated to
            // each of elemental boundaries of these patches
            std::map<int, int> numLocalBndCoeffsPerPatchNew;
            for(int i = 0; i < numPatchesNew; i++)
            {
                numLocalBndCoeffsPerPatchNew[i] = 0;
            }

            int minval;
            int maxval;
            int curPatch;
            for(i = cnt = 0; i < numPatchesOld; i++)
            {
                // For every patch at the current level, the mask array
                // locPatchMask should be filled with
                // - the same (positive) number for each entry (which will
                //   correspond to the patch at the next level it belongs to)
                // - the same (positive) number for each entry, except some
                //   entries that are -1 (the enties correspond to -1, will be
                //   mapped to the local boundary of the next level patch given
                //   by the positive number)
                // - -1 for all entries. In this case, we will make an
                //   additional patch only consisting of boundaries at the next
                //   level
                minval = *min_element(&locPatchMask[cnt],
                                      &locPatchMask[cnt]+numLocalBndCoeffsPerPatchOld[i]);
                maxval = *max_element(&locPatchMask[cnt],
                                      &locPatchMask[cnt]+numLocalBndCoeffsPerPatchOld[i]);
                ASSERTL0((minval==maxval)||(minval==-1),"These values should never be the same");
                
                if(maxval == -1)
                {
                    curPatch = numPatchesNew;
                    numLocalBndCoeffsPerPatchNew[curPatch] = 0;
                    numPatchesNew++;
                }
                else
                {
                    curPatch = maxval;
                }
                
                for(j = 0; j < numLocalBndCoeffsPerPatchOld[i]; j++ )
                {
                    ASSERTL0((locPatchMask[cnt]==maxval)||(locPatchMask[cnt]==minval),
                             "These values should never be the same");
                    if(locPatchMask[cnt] == -1)
                    {
                        ++numLocalBndCoeffsPerPatchNew[curPatch];
                    }
                    cnt++;
                }
            }

            // Count how many local dofs of the old level are mapped
            // to the local boundary dofs of the new level
            m_numLocalBndCoeffs  = 0;
            m_numPatches         = numLocalBndCoeffsPerPatchNew.size();
            m_numLocalBndCoeffsPerPatch = Array<OneD, unsigned int>(m_numPatches);
            m_numLocalIntCoeffsPerPatch = Array<OneD, unsigned int>(m_numPatches,0u);
            for(int i = 0; i < m_numPatches; i++)
            {
                m_numLocalBndCoeffsPerPatch[i] = (unsigned int) numLocalBndCoeffsPerPatchNew[i];
                m_numLocalBndCoeffs           += numLocalBndCoeffsPerPatchNew[i];
            }
            multiLevelGraph->GetNintDofsPerPatch(newLevel,m_numLocalIntCoeffsPerPatch);

            // Also initialise some more data members
            m_solnType              = solnTypeOld;
            ASSERTL1(m_solnType==eDirectMultiLevelStaticCond
                    ||m_solnType==eIterativeMultiLevelStaticCond
                    ||m_solnType==eXxtMultiLevelStaticCond
                    ||m_solnType==ePETScMultiLevelStaticCond,
                     "This method should only be called for in "
                     "case of multi-level static condensation.");
            m_staticCondLevel       = newLevel;
            m_signChange            = signChangeOld;
            m_numLocalDirBndCoeffs  = numLocalDirBndCoeffsOld;
            m_numGlobalDirBndCoeffs = numGlobalDirBndCoeffsOld;
            m_numGlobalBndCoeffs    = multiLevelGraph->GetInteriorOffset(newLevel) +
                m_numGlobalDirBndCoeffs;
            m_numGlobalCoeffs       = multiLevelGraph->GetNumGlobalDofs(newLevel) +
                m_numGlobalDirBndCoeffs;
            m_localToGlobalBndMap   = Array<OneD,int>(m_numLocalBndCoeffs);
            if(m_signChange)
            {
                m_localToGlobalBndSign    = Array<OneD,NekDouble>(m_numLocalBndCoeffs);
            }
            
            m_patchMapFromPrevLevel = MemoryManager<PatchMap>::AllocateSharedPtr(numLocalBndCoeffsOld);

            m_globalToUniversalBndMap = Array<OneD, int>(
                m_numGlobalBndCoeffs, oldLevelMap->GetGlobalToUniversalBndMap());
            m_globalToUniversalBndMapUnique = Array<OneD, int>(
                m_numGlobalBndCoeffs, oldLevelMap->GetGlobalToUniversalBndMapUnique());

            // Set up an offset array that denotes the offset of the local
            // boundary degrees of freedom of the next level
            Array<OneD, int> numLocalBndCoeffsPerPatchOffset(m_numPatches+1,0);
            for(int i = 1; i < m_numPatches+1; i++)
            {
                numLocalBndCoeffsPerPatchOffset[i] += numLocalBndCoeffsPerPatchOffset[i-1] + numLocalBndCoeffsPerPatchNew[i-1];
            }

            int additionalPatchCnt = numPatchesWithIntNew;
            int newid;
            int blockid;
            bool isBndDof;
            NekDouble sign;
            Array<OneD, int> bndDofPerPatchCnt(m_numPatches,0);
            for(i = cnt = 0; i < numPatchesOld; i++)
            {
                minval = *min_element(&locPatchMask[cnt],&locPatchMask[cnt]+numLocalBndCoeffsPerPatchOld[i]);
                maxval = *max_element(&locPatchMask[cnt],&locPatchMask[cnt]+numLocalBndCoeffsPerPatchOld[i]);
                ASSERTL0((minval==maxval)||(minval==-1),"These values should never be the same");

                if(maxval == -1)
                {
                    curPatch = additionalPatchCnt;
                    additionalPatchCnt++;
                }
                else
                {
                    curPatch = maxval;
                }

                for(j = 0; j < numLocalBndCoeffsPerPatchOld[i]; j++ )
                {
                    ASSERTL0((locPatchMask[cnt]==maxval)||(locPatchMask[cnt]==minval),
                             "These values should never be the same");

                    sign = oldLevelMap->GetLocalToGlobalBndSign(cnt);

                    if(locPatchMask[cnt] == -1)
                    {
                        newid = numLocalBndCoeffsPerPatchOffset[curPatch];

                        m_localToGlobalBndMap[newid] = oldLevelMap->GetLocalToGlobalBndMap(cnt);
                        if(m_signChange)
                        {
                            m_localToGlobalBndSign[ newid ] = sign;
                        }

                        blockid = bndDofPerPatchCnt[curPatch];
                        isBndDof = true;


                        numLocalBndCoeffsPerPatchOffset[curPatch]++;
                        bndDofPerPatchCnt[curPatch]++;
                    }
                    else
                    {
                        newid = oldLevelMap->GetLocalToGlobalBndMap(cnt) -
                            m_numGlobalBndCoeffs+m_numLocalBndCoeffs;

                        blockid = oldLevelMap->GetLocalToGlobalBndMap(cnt)-
                            m_numGlobalDirBndCoeffs - multiLevelGraph->GetInteriorOffset(newLevel,curPatch);
                        isBndDof = false;
                    }

                    sign = isBndDof?1.0:sign;

                    m_patchMapFromPrevLevel->SetPatchMap(cnt,curPatch,blockid,isBndDof,sign);

                    cnt++;
                }
            }
            CalculateBndSystemBandWidth();

            // Postprocess the computed information - Update the old
            // level with the mapping to new evel
            // oldLevelMap->SetLocalBndToNextLevelMap(oldLocalBndToNewLevelMap,oldLocalBndToNewLevelSign);
            // - Construct the next level mapping object
            if(m_staticCondLevel < (multiLevelGraph->GetNlevels()-1))
            {
                m_nextLevelLocalToGlobalMap = MemoryManager<AssemblyMap>::AllocateSharedPtr(this,multiLevelGraph);
            }
        }

~AssemblyMap()

Nektar::MultiRegions::AssemblyMap::~AssemblyMap ( void  )

virtual

Destructor.

Definition at line 438 of file AssemblyMap.cpp.

        {
        }

Member Function Documentation

Assemble() [1/2]

void Nektar::MultiRegions::AssemblyMap::Assemble	(	const Array< OneD, const NekDouble > &	loc,
		Array< OneD, NekDouble > &	global
	)		const

Definition at line 793 of file AssemblyMap.cpp.

References v_Assemble().

Referenced by Nektar::MultiRegions::AssemblyMapCG::v_Assemble().

        {
            v_Assemble(loc,global);
        }

Assemble() [2/2]

void Nektar::MultiRegions::AssemblyMap::Assemble	(	const NekVector< NekDouble > &	loc,
		NekVector< NekDouble > &	global
	)		const

Definition at line 800 of file AssemblyMap.cpp.

References v_Assemble().

        {
            v_Assemble(loc,global);
        }

AssembleBnd() [1/4]

void Nektar::MultiRegions::AssemblyMap::AssembleBnd	(	const NekVector< NekDouble > &	loc,
		NekVector< NekDouble > &	global,
		int	offset
	)		const

Definition at line 1113 of file AssemblyMap.cpp.

References Nektar::NekVector< DataType >::GetPtr().

Referenced by AssembleBnd(), Nektar::MultiRegions::AssemblyMapDG::v_Assemble(), and Nektar::MultiRegions::AssemblyMapDG::v_LocalToGlobal().

        {
            AssembleBnd(loc.GetPtr(), global.GetPtr(), offset);
        }

AssembleBnd() [2/4]

void Nektar::MultiRegions::AssemblyMap::AssembleBnd	(	const NekVector< NekDouble > &	loc,
		NekVector< NekDouble > &	global
	)		const

Definition at line 1121 of file AssemblyMap.cpp.

References AssembleBnd(), and Nektar::NekVector< DataType >::GetPtr().

        {
            AssembleBnd(loc.GetPtr(), global.GetPtr());
        }

AssembleBnd() [3/4]

void Nektar::MultiRegions::AssemblyMap::AssembleBnd	(	const Array< OneD, const NekDouble > &	loc,
		Array< OneD, NekDouble > &	global,
		int	offset
	)		const

Definition at line 1129 of file AssemblyMap.cpp.

References ASSERTL1, Vmath::Assmb(), m_localToGlobalBndMap, m_localToGlobalBndSign, m_numGlobalBndCoeffs, m_numLocalBndCoeffs, m_signChange, UniversalAssembleBnd(), and Vmath::Vcopy().

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local array is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs-offset,"Global array is not of correct dimension");
            Array<OneD,NekDouble> tmp(m_numGlobalBndCoeffs,0.0);

            if(m_signChange)
            {
                Vmath::Assmb(m_numLocalBndCoeffs, m_localToGlobalBndSign.get(),loc.get(), m_localToGlobalBndMap.get(), tmp.get());
            }
            else
            {
                Vmath::Assmb(m_numLocalBndCoeffs,loc.get(), m_localToGlobalBndMap.get(), tmp.get());
            }
            UniversalAssembleBnd(tmp);
            Vmath::Vcopy(m_numGlobalBndCoeffs-offset, tmp.get() + offset, 1, global.get(), 1);
        }

AssembleBnd() [4/4]

void Nektar::MultiRegions::AssemblyMap::AssembleBnd	(	const Array< OneD, const NekDouble > &	loc,
		Array< OneD, NekDouble > &	global
	)		const

Definition at line 1150 of file AssemblyMap.cpp.

References ASSERTL1, Vmath::Assmb(), m_localToGlobalBndMap, m_localToGlobalBndSign, m_numGlobalBndCoeffs, m_numLocalBndCoeffs, m_signChange, UniversalAssembleBnd(), and Vmath::Zero().

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local vector is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs,"Global vector is not of correct dimension");

            Vmath::Zero(m_numGlobalBndCoeffs, global.get(), 1);

            if(m_signChange)
            {
                Vmath::Assmb(m_numLocalBndCoeffs,m_localToGlobalBndSign.get(),
                             loc.get(), m_localToGlobalBndMap.get(), global.get());
            }
            else
            {
                Vmath::Assmb(m_numLocalBndCoeffs,loc.get(), m_localToGlobalBndMap.get(), global.get());
            }
            UniversalAssembleBnd(global);
        }

AtLastLevel()

bool Nektar::MultiRegions::AssemblyMap::AtLastLevel

(

)

const

Returns true if this is the last level in the multi-level static condensation.

Definition at line 1236 of file AssemblyMap.cpp.

References m_nextLevelLocalToGlobalMap.

        {
            return !( m_nextLevelLocalToGlobalMap.get() );
        }

CalculateBndSystemBandWidth()

void Nektar::MultiRegions::AssemblyMap::CalculateBndSystemBandWidth ( )

protected

Calculates the bandwidth of the boundary system.

The bandwidth calculated corresponds to what is referred to as half-bandwidth. If the elements of the matrix are designated as a_ij, it corresponds to the maximum value of |i-j| for non-zero a_ij. As a result, the value also corresponds to the number of sub- or super-diagonals.

The bandwith can be calculated elementally as it corresponds to the maximal elemental bandwith (i.e. the maximal difference in global DOF index for every element).

We here calculate the bandwith of the global boundary system (as used for static condensation).

Definition at line 457 of file AssemblyMap.cpp.

References m_bndSystemBandWidth, m_localToGlobalBndMap, m_numGlobalDirBndCoeffs, m_numLocalBndCoeffs, m_numLocalBndCoeffsPerPatch, and m_numPatches.

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

        {
            int i,j;
            int cnt = 0;
            int locSize;
            int maxId;
            int minId;
            int bwidth = -1;
            for(i = 0; i < m_numPatches; ++i)
            {
                locSize = m_numLocalBndCoeffsPerPatch[i];
                maxId = -1;
                minId = m_numLocalBndCoeffs+1;
                for(j = 0; j < locSize; j++)
                {
                    if(m_localToGlobalBndMap[cnt+j] >= m_numGlobalDirBndCoeffs)
                    {
                        if(m_localToGlobalBndMap[cnt+j] > maxId)
                        {
                            maxId = m_localToGlobalBndMap[cnt+j];
                        }

                        if(m_localToGlobalBndMap[cnt+j] < minId)
                        {
                            minId = m_localToGlobalBndMap[cnt+j];
                        }
                    }
                }
                bwidth = (bwidth>(maxId-minId))?bwidth:(maxId-minId);

                cnt+=locSize;
            }

            m_bndSystemBandWidth = bwidth;
        }

GetBndCondCoeffsToGlobalCoeffsMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::GetBndCondCoeffsToGlobalCoeffsMap

(

const int

i

)

Retrieves the global index corresponding to a boundary expansion mode.

Definition at line 922 of file AssemblyMap.cpp.

References m_bndCondCoeffsToGlobalCoeffsMap.

        {
            return m_bndCondCoeffsToGlobalCoeffsMap[i];
        }

GetBndCondCoeffsToGlobalCoeffsMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetBndCondCoeffsToGlobalCoeffsMap

(

)

Retrieves the global indices corresponding to the boundary expansion modes.

Definition at line 956 of file AssemblyMap.cpp.

References m_bndCondCoeffsToGlobalCoeffsMap.

Referenced by Nektar::CoupledAssemblyMap::CoupledAssemblyMap().

        {
            return m_bndCondCoeffsToGlobalCoeffsMap;
        }

GetBndCondCoeffsToGlobalCoeffsSign()

NekDouble Nektar::MultiRegions::AssemblyMap::GetBndCondCoeffsToGlobalCoeffsSign

(

const int

i

)

Returns the modal sign associated with a given boundary expansion mode.

Definition at line 943 of file AssemblyMap.cpp.

References m_bndCondCoeffsToGlobalCoeffsSign, and m_signChange.

        {
            if(m_signChange)
            {
                return m_bndCondCoeffsToGlobalCoeffsSign[i];
            }
            else
            {
                return 1.0;
            }
        }

GetBndCondTraceToGlobalTraceMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::GetBndCondTraceToGlobalTraceMap

(

const int

i

)

Returns the global index of the boundary trace giving the index on the boundary expansion.

Definition at line 929 of file AssemblyMap.cpp.

References ASSERTL1, GetBndCondTraceToGlobalTraceMap(), and m_bndCondTraceToGlobalTraceMap.

        {
            ASSERTL1(i < m_bndCondTraceToGlobalTraceMap.num_elements(),
                     "Index out of range.");
            return m_bndCondTraceToGlobalTraceMap[i];
        }

GetBndCondTraceToGlobalTraceMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetBndCondTraceToGlobalTraceMap

(

)

Definition at line 938 of file AssemblyMap.cpp.

References m_bndCondTraceToGlobalTraceMap.

Referenced by GetBndCondTraceToGlobalTraceMap().

        {
            return m_bndCondTraceToGlobalTraceMap;
        }

GetBndSystemBandWidth()

int Nektar::MultiRegions::AssemblyMap::GetBndSystemBandWidth

(

)

const

Returns the bandwidth of the boundary system.

Definition at line 1195 of file AssemblyMap.cpp.

References m_bndSystemBandWidth.

Referenced by Nektar::MultiRegions::AssemblyMapDG::v_GetFullSystemBandWidth().

        {
            return m_bndSystemBandWidth;
        }

GetComm()

LibUtilities::CommSharedPtr Nektar::MultiRegions::AssemblyMap::GetComm

(

)

Retrieves the communicator.

Definition at line 713 of file AssemblyMap.cpp.

References m_comm.

        {
            return m_comm;
        }

GetExtraDirEdges()

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetExtraDirEdges

(

)

Definition at line 866 of file AssemblyMap.cpp.

References v_GetExtraDirEdges().

        {
            return v_GetExtraDirEdges();
        }

GetFullSystemBandWidth()

int Nektar::MultiRegions::AssemblyMap::GetFullSystemBandWidth

(

)

const

Definition at line 826 of file AssemblyMap.cpp.

References v_GetFullSystemBandWidth().

        {
            return v_GetFullSystemBandWidth();
        }

GetGlobalSysSolnType()

GlobalSysSolnType Nektar::MultiRegions::AssemblyMap::GetGlobalSysSolnType

(

)

const

Returns the method of solving global systems.

Definition at line 1242 of file AssemblyMap.cpp.

References m_solnType.

Referenced by AssemblyMap().

        {
            return m_solnType;
        }

GetGlobalToUniversalBndMap()

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetGlobalToUniversalBndMap

(

)

Definition at line 899 of file AssemblyMap.cpp.

References m_globalToUniversalBndMap.

Referenced by AssemblyMap().

        {
            return m_globalToUniversalBndMap;
        }

GetGlobalToUniversalBndMapUnique()

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetGlobalToUniversalBndMapUnique

(

)

Definition at line 904 of file AssemblyMap.cpp.

References m_globalToUniversalBndMapUnique.

Referenced by AssemblyMap().

        {
            return m_globalToUniversalBndMapUnique;
        }

GetGlobalToUniversalMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::GetGlobalToUniversalMap

(

const int

i

)

const

Definition at line 728 of file AssemblyMap.cpp.

References v_GetGlobalToUniversalMap().

        {
            return v_GetGlobalToUniversalMap(i);
        }

GetGlobalToUniversalMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetGlobalToUniversalMap

(

)

Definition at line 743 of file AssemblyMap.cpp.

References v_GetGlobalToUniversalMap().

        {
            return v_GetGlobalToUniversalMap();
        }

GetGlobalToUniversalMapUnique() [1/2]

int Nektar::MultiRegions::AssemblyMap::GetGlobalToUniversalMapUnique

(

const int

i

)

const

Definition at line 733 of file AssemblyMap.cpp.

References v_GetGlobalToUniversalMapUnique().

        {
            return v_GetGlobalToUniversalMapUnique(i);
        }

GetGlobalToUniversalMapUnique() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetGlobalToUniversalMapUnique

(

)

Definition at line 748 of file AssemblyMap.cpp.

References v_GetGlobalToUniversalMapUnique().

        {
            return v_GetGlobalToUniversalMapUnique();
        }

GetHash()

size_t Nektar::MultiRegions::AssemblyMap::GetHash

(

)

const

Retrieves the hash of this map.

Definition at line 718 of file AssemblyMap.cpp.

References m_hash.

        {
            return m_hash;
        }

GetIterativeTolerance()

NekDouble Nektar::MultiRegions::AssemblyMap::GetIterativeTolerance

(

)

const

Definition at line 1252 of file AssemblyMap.cpp.

References m_iterativeTolerance.

        {
            return m_iterativeTolerance;
        }

GetLocalToGlobalBndMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalBndMap

(

const int

i

)

const

Retrieve the global index of a given local boundary mode.

Definition at line 876 of file AssemblyMap.cpp.

References m_localToGlobalBndMap.

Referenced by AssemblyMap().

        {
            return m_localToGlobalBndMap[i];
        }

GetLocalToGlobalBndMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalBndMap

(

void

)

Retrieve the global indices of the local boundary modes.

Definition at line 882 of file AssemblyMap.cpp.

References m_localToGlobalBndMap.

        {
            return m_localToGlobalBndMap;
        }

GetLocalToGlobalBndSign() [1/2]

NekDouble Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalBndSign

(

const int

i

)

const

Retrieve the sign change of a given local boundary mode.

Definition at line 909 of file AssemblyMap.cpp.

References m_localToGlobalBndSign, and m_signChange.

Referenced by AssemblyMap().

        {
            if(m_signChange)
            {
                return m_localToGlobalBndSign[i];
            }
            else
            {
                return 1.0;
            }
        }

GetLocalToGlobalBndSign() [2/2]

Array< OneD, const NekDouble > Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalBndSign

(

void

)

const

Retrieve the sign change for all local boundary modes.

Definition at line 894 of file AssemblyMap.cpp.

References m_localToGlobalBndSign.

Referenced by Nektar::MultiRegions::AssemblyMapDG::v_GetLocalToGlobalSign().

        {
            return m_localToGlobalBndSign;
        }

GetLocalToGlobalMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalMap

(

const int

i

)

const

Definition at line 723 of file AssemblyMap.cpp.

References v_GetLocalToGlobalMap().

        {
            return v_GetLocalToGlobalMap(i);
        }

GetLocalToGlobalMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalMap

(

)

Definition at line 738 of file AssemblyMap.cpp.

References v_GetLocalToGlobalMap().

        {
            return v_GetLocalToGlobalMap();
        }

GetLocalToGlobalSign() [1/2]

NekDouble Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalSign

(

const int

i

)

const

Definition at line 753 of file AssemblyMap.cpp.

References v_GetLocalToGlobalSign().

        {
            return v_GetLocalToGlobalSign(i);
        }

GetLocalToGlobalSign() [2/2]

const Array< OneD, NekDouble > & Nektar::MultiRegions::AssemblyMap::GetLocalToGlobalSign

(

)

const

Definition at line 758 of file AssemblyMap.cpp.

References v_GetLocalToGlobalSign().

        {
            return v_GetLocalToGlobalSign();
        }

GetLowestStaticCondLevel()

int Nektar::MultiRegions::AssemblyMap::GetLowestStaticCondLevel ( ) const

inline

Definition at line 297 of file AssemblyMap.h.

References m_lowestStaticCondLevel.

            {
                return m_lowestStaticCondLevel;
            }

GetMaxIterations()

int Nektar::MultiRegions::AssemblyMap::GetMaxIterations

(

)

const

Definition at line 1257 of file AssemblyMap.cpp.

References m_maxIterations.

        {
            return m_maxIterations;
        }

GetNextLevelLocalToGlobalMap()

const AssemblyMapSharedPtr Nektar::MultiRegions::AssemblyMap::GetNextLevelLocalToGlobalMap

(

)

const

Returns the local to global mapping for the next level in the multi-level static condensation.

Definition at line 1224 of file AssemblyMap.cpp.

References m_nextLevelLocalToGlobalMap.

        {
            return  m_nextLevelLocalToGlobalMap;
        }

GetNumDirEdges()

int Nektar::MultiRegions::AssemblyMap::GetNumDirEdges

(

)

const

Definition at line 846 of file AssemblyMap.cpp.

References v_GetNumDirEdges().

        {
            return v_GetNumDirEdges();
        }

GetNumDirFaces()

int Nektar::MultiRegions::AssemblyMap::GetNumDirFaces

(

)

const

Definition at line 851 of file AssemblyMap.cpp.

References v_GetNumDirFaces().

        {
            return v_GetNumDirFaces();
        }

GetNumGlobalBndCoeffs()

int Nektar::MultiRegions::AssemblyMap::GetNumGlobalBndCoeffs

(

)

const

Returns the total number of global boundary coefficients.

Definition at line 978 of file AssemblyMap.cpp.

References m_numGlobalBndCoeffs.

Referenced by AssemblyMap().

        {
            return m_numGlobalBndCoeffs;
        }

GetNumGlobalCoeffs()

int Nektar::MultiRegions::AssemblyMap::GetNumGlobalCoeffs

(

)

const

Returns the total number of global coefficients.

Definition at line 988 of file AssemblyMap.cpp.

References m_numGlobalCoeffs.

        {
            return m_numGlobalCoeffs;
        }

GetNumGlobalDirBndCoeffs()

int Nektar::MultiRegions::AssemblyMap::GetNumGlobalDirBndCoeffs

(

)

const

Returns the number of global Dirichlet boundary coefficients.

Definition at line 962 of file AssemblyMap.cpp.

References m_numGlobalDirBndCoeffs.

Referenced by AssemblyMap(), and Nektar::MultiRegions::PreconditionerBlock::BlockPreconditionerHDG().

        {
            return m_numGlobalDirBndCoeffs;
        }

GetNumLocalBndCoeffs()

int Nektar::MultiRegions::AssemblyMap::GetNumLocalBndCoeffs

(

)

const

Returns the total number of local boundary coefficients.

Definition at line 973 of file AssemblyMap.cpp.

References m_numLocalBndCoeffs.

Referenced by AssemblyMap().

        {
            return m_numLocalBndCoeffs;
        }

GetNumLocalBndCoeffsPerPatch()

const Array< OneD, const unsigned int > & Nektar::MultiRegions::AssemblyMap::GetNumLocalBndCoeffsPerPatch

(

)

Returns the number of local boundary coefficients in each patch.

Definition at line 1211 of file AssemblyMap.cpp.

References m_numLocalBndCoeffsPerPatch.

Referenced by AssemblyMap().

        {
            return m_numLocalBndCoeffsPerPatch;
        }

GetNumLocalCoeffs()

int Nektar::MultiRegions::AssemblyMap::GetNumLocalCoeffs

(

)

const

Returns the total number of local coefficients.

Definition at line 983 of file AssemblyMap.cpp.

References m_numLocalCoeffs.

        {
            return m_numLocalCoeffs;
        }

GetNumLocalDirBndCoeffs()

int Nektar::MultiRegions::AssemblyMap::GetNumLocalDirBndCoeffs

(

)

const

Returns the number of local Dirichlet boundary coefficients.

Definition at line 968 of file AssemblyMap.cpp.

References m_numLocalDirBndCoeffs.

Referenced by AssemblyMap().

        {
            return m_numLocalDirBndCoeffs;
        }

GetNumLocalIntCoeffsPerPatch()

const Array< OneD, const unsigned int > & Nektar::MultiRegions::AssemblyMap::GetNumLocalIntCoeffsPerPatch

(

)

Returns the number of local interior coefficients in each patch.

Definition at line 1218 of file AssemblyMap.cpp.

References m_numLocalIntCoeffsPerPatch.

        {
            return m_numLocalIntCoeffsPerPatch;
        }

GetNumNonDirEdgeModes()

int Nektar::MultiRegions::AssemblyMap::GetNumNonDirEdgeModes

(

)

const

Definition at line 836 of file AssemblyMap.cpp.

References v_GetNumNonDirEdgeModes().

        {
            return v_GetNumNonDirEdgeModes();
        }

GetNumNonDirEdges()

int Nektar::MultiRegions::AssemblyMap::GetNumNonDirEdges

(

)

const

Definition at line 856 of file AssemblyMap.cpp.

References v_GetNumNonDirEdges().

        {
            return v_GetNumNonDirEdges();
        }

GetNumNonDirFaceModes()

int Nektar::MultiRegions::AssemblyMap::GetNumNonDirFaceModes

(

)

const

Definition at line 841 of file AssemblyMap.cpp.

References v_GetNumNonDirFaceModes().

        {
            return v_GetNumNonDirFaceModes();
        }

GetNumNonDirFaces()

int Nektar::MultiRegions::AssemblyMap::GetNumNonDirFaces

(

)

const

Definition at line 861 of file AssemblyMap.cpp.

References v_GetNumNonDirFaces().

        {
            return v_GetNumNonDirFaces();
        }

GetNumNonDirVertexModes()

int Nektar::MultiRegions::AssemblyMap::GetNumNonDirVertexModes

(

)

const

Definition at line 831 of file AssemblyMap.cpp.

References v_GetNumNonDirVertexModes().

        {
            return v_GetNumNonDirVertexModes();
        }

GetNumPatches()

int Nektar::MultiRegions::AssemblyMap::GetNumPatches

(

)

const

Returns the number of patches in this static condensation level.

Definition at line 1205 of file AssemblyMap.cpp.

References m_numPatches.

Referenced by AssemblyMap().

        {
            return m_numPatches;
        }

GetPatchMapFromPrevLevel()

const PatchMapSharedPtr & Nektar::MultiRegions::AssemblyMap::GetPatchMapFromPrevLevel

(

void

)

const

Returns the patch map from the previous level of the multi-level static condensation.

Definition at line 1230 of file AssemblyMap.cpp.

References m_patchMapFromPrevLevel.

        {
            return m_patchMapFromPrevLevel;
        }

GetPreconType()

PreconditionerType Nektar::MultiRegions::AssemblyMap::GetPreconType

(

)

const

Definition at line 1247 of file AssemblyMap.cpp.

References m_preconType.

        {
            return m_preconType;
        }

GetSignChange()

bool Nektar::MultiRegions::AssemblyMap::GetSignChange

(

)

Returns true if using a modal expansion requiring a change of sign of some modes.

Definition at line 887 of file AssemblyMap.cpp.

References m_signChange.

Referenced by AssemblyMap().

        {
            return m_signChange;
        }

GetSingularSystem()

bool Nektar::MultiRegions::AssemblyMap::GetSingularSystem

(

)

const

Retrieves if the system is singular (true) or not (false)

Definition at line 993 of file AssemblyMap.cpp.

References m_systemSingular.

        {
            return m_systemSingular;
        }

GetStaticCondLevel()

int Nektar::MultiRegions::AssemblyMap::GetStaticCondLevel

(

)

const

Returns the level of static condensation for this map.

Definition at line 1200 of file AssemblyMap.cpp.

References m_staticCondLevel.

Referenced by AssemblyMap().

        {
            return m_staticCondLevel;
        }

GetSuccessiveRHS()

int Nektar::MultiRegions::AssemblyMap::GetSuccessiveRHS

(

)

const

Definition at line 1262 of file AssemblyMap.cpp.

References m_successiveRHS.

        {
            return m_successiveRHS;
        }

GlobalToLocal() [1/2]

void Nektar::MultiRegions::AssemblyMap::GlobalToLocal	(	const Array< OneD, const NekDouble > &	global,
		Array< OneD, NekDouble > &	loc
	)		const

Definition at line 779 of file AssemblyMap.cpp.

References v_GlobalToLocal().

Referenced by Nektar::MultiRegions::AssemblyMapCG::v_GlobalToLocal().

        {
            v_GlobalToLocal(global,loc);
        }

GlobalToLocal() [2/2]

void Nektar::MultiRegions::AssemblyMap::GlobalToLocal	(	const NekVector< NekDouble > &	global,
		NekVector< NekDouble > &	loc
	)		const

Definition at line 786 of file AssemblyMap.cpp.

References v_GlobalToLocal().

        {
            v_GlobalToLocal(global,loc);
        }

GlobalToLocalBnd() [1/4]

void Nektar::MultiRegions::AssemblyMap::GlobalToLocalBnd	(	const NekVector< NekDouble > &	global,
		NekVector< NekDouble > &	loc,
		int	offset
	)		const

Definition at line 998 of file AssemblyMap.cpp.

References Nektar::NekVector< DataType >::GetPtr().

Referenced by GlobalToLocalBnd(), and Nektar::MultiRegions::AssemblyMapDG::v_GlobalToLocal().

        {
            GlobalToLocalBnd(global.GetPtr(), loc.GetPtr(), offset);
        }

GlobalToLocalBnd() [2/4]

void Nektar::MultiRegions::AssemblyMap::GlobalToLocalBnd	(	const NekVector< NekDouble > &	global,
		NekVector< NekDouble > &	loc
	)		const

Definition at line 1007 of file AssemblyMap.cpp.

References Nektar::NekVector< DataType >::GetPtr(), and GlobalToLocalBnd().

        {
            GlobalToLocalBnd(global.GetPtr(), loc.GetPtr());
        }

GlobalToLocalBnd() [3/4]

void Nektar::MultiRegions::AssemblyMap::GlobalToLocalBnd	(	const Array< OneD, const NekDouble > &	global,
		Array< OneD, NekDouble > &	loc,
		int	offset
	)		const

Definition at line 1015 of file AssemblyMap.cpp.

References ASSERTL1, Vmath::Gathr(), m_localToGlobalBndMap, m_localToGlobalBndSign, m_numGlobalBndCoeffs, m_numLocalBndCoeffs, m_signChange, and Vmath::Vcopy().

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local vector is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs-offset,"Global vector is not of correct dimension");

            // offset input data by length "offset" for Dirichlet boundary conditions.
            Array<OneD,NekDouble> tmp(m_numGlobalBndCoeffs,0.0);
            Vmath::Vcopy(m_numGlobalBndCoeffs-offset, global.get(), 1, tmp.get() + offset, 1);

            if(m_signChange)
            {
                Vmath::Gathr(m_numLocalBndCoeffs, m_localToGlobalBndSign.get(), tmp.get(), m_localToGlobalBndMap.get(), loc.get());
            }
            else
            {
                Vmath::Gathr(m_numLocalBndCoeffs, tmp.get(), m_localToGlobalBndMap.get(), loc.get());
            }
        }

GlobalToLocalBnd() [4/4]

void Nektar::MultiRegions::AssemblyMap::GlobalToLocalBnd	(	const Array< OneD, const NekDouble > &	global,
		Array< OneD, NekDouble > &	loc
	)		const

Definition at line 1037 of file AssemblyMap.cpp.

References ASSERTL1, Vmath::Gathr(), m_localToGlobalBndMap, m_localToGlobalBndSign, m_numGlobalBndCoeffs, m_numLocalBndCoeffs, and m_signChange.

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local vector is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs,"Global vector is not of correct dimension");

            if(m_signChange)
            {
                Vmath::Gathr(m_numLocalBndCoeffs, m_localToGlobalBndSign.get(), global.get(), m_localToGlobalBndMap.get(), loc.get());
            }
            else
            {
                Vmath::Gathr(m_numLocalBndCoeffs, global.get(), m_localToGlobalBndMap.get(), loc.get());
            }
        }

GlobalToLocalBndWithoutSign()

void Nektar::MultiRegions::AssemblyMap::GlobalToLocalBndWithoutSign ( const Array< OneD, const NekDouble > &  global, Array< OneD, NekDouble > &  loc  )

protected

Definition at line 1267 of file AssemblyMap.cpp.

References ASSERTL1, Vmath::Gathr(), m_localToGlobalBndMap, m_numGlobalBndCoeffs, and m_numLocalBndCoeffs.

Referenced by AssemblyMap().

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local vector is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs,"Global vector is not of correct dimension");

            Vmath::Gathr(m_numLocalBndCoeffs, global.get(), m_localToGlobalBndMap.get(), loc.get());
        }

LinearSpaceMap()

std::shared_ptr< AssemblyMap > Nektar::MultiRegions::AssemblyMap::LinearSpaceMap	(	const ExpList &	locexp,
		GlobalSysSolnType	solnType
	)

Definition at line 871 of file AssemblyMap.cpp.

References v_LinearSpaceMap().

        {
            return v_LinearSpaceMap(locexp, solnType);
        }

LocalBndToGlobal() [1/4]

void Nektar::MultiRegions::AssemblyMap::LocalBndToGlobal	(	const NekVector< NekDouble > &	loc,
		NekVector< NekDouble > &	global,
		int	offset
	)		const

Definition at line 1054 of file AssemblyMap.cpp.

References Nektar::NekVector< DataType >::GetPtr().

Referenced by LocalBndToGlobal().

        {
            LocalBndToGlobal(loc.GetPtr(), global.GetPtr(), offset);
        }

LocalBndToGlobal() [2/4]

void Nektar::MultiRegions::AssemblyMap::LocalBndToGlobal	(	const NekVector< NekDouble > &	loc,
		NekVector< NekDouble > &	global
	)		const

Definition at line 1087 of file AssemblyMap.cpp.

References Nektar::NekVector< DataType >::GetPtr(), and LocalBndToGlobal().

        {
            LocalBndToGlobal(loc.GetPtr(), global.GetPtr());
        }

LocalBndToGlobal() [3/4]

void Nektar::MultiRegions::AssemblyMap::LocalBndToGlobal	(	const Array< OneD, const NekDouble > &	loc,
		Array< OneD, NekDouble > &	global,
		int	offset
	)		const

Definition at line 1063 of file AssemblyMap.cpp.

References ASSERTL1, m_localToGlobalBndMap, m_localToGlobalBndSign, m_numGlobalBndCoeffs, m_numLocalBndCoeffs, m_signChange, Vmath::Scatr(), UniversalAssembleBnd(), and Vmath::Vcopy().

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local vector is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs-offset,"Global vector is not of correct dimension");

            // offset input data by length "offset" for Dirichlet boundary conditions.
            Array<OneD,NekDouble> tmp(m_numGlobalBndCoeffs,0.0);

            if(m_signChange)
            {
                Vmath::Scatr(m_numLocalBndCoeffs, m_localToGlobalBndSign.get(), loc.get(), m_localToGlobalBndMap.get(), tmp.get());
            }
            else
            {
                Vmath::Scatr(m_numLocalBndCoeffs, loc.get(), m_localToGlobalBndMap.get(), tmp.get());
            }

            UniversalAssembleBnd(tmp);
            Vmath::Vcopy(m_numGlobalBndCoeffs-offset, tmp.get()+offset, 1, global.get(), 1);
        }

LocalBndToGlobal() [4/4]

void Nektar::MultiRegions::AssemblyMap::LocalBndToGlobal	(	const Array< OneD, const NekDouble > &	loc,
		Array< OneD, NekDouble > &	global
	)		const

Definition at line 1095 of file AssemblyMap.cpp.

References ASSERTL1, m_localToGlobalBndMap, m_localToGlobalBndSign, m_numGlobalBndCoeffs, m_numLocalBndCoeffs, m_signChange, and Vmath::Scatr().

        {
            ASSERTL1(loc.num_elements() >= m_numLocalBndCoeffs,"Local vector is not of correct dimension");
            ASSERTL1(global.num_elements() >= m_numGlobalBndCoeffs,"Global vector is not of correct dimension");

            if(m_signChange)
            {
                Vmath::Scatr(m_numLocalBndCoeffs, m_localToGlobalBndSign.get(), loc.get(), m_localToGlobalBndMap.get(), global.get());
            }
            else
            {
                Vmath::Scatr(m_numLocalBndCoeffs, loc.get(), m_localToGlobalBndMap.get(), global.get());
            }
        }

LocalToGlobal() [1/2]

void Nektar::MultiRegions::AssemblyMap::LocalToGlobal	(	const Array< OneD, const NekDouble > &	loc,
		Array< OneD, NekDouble > &	global,
		bool	useComm = true
	)		const

Definition at line 763 of file AssemblyMap.cpp.

References v_LocalToGlobal().

Referenced by Nektar::MultiRegions::AssemblyMapCG::v_LocalToGlobal().

        {
            v_LocalToGlobal(loc,global,useComm);
        }

LocalToGlobal() [2/2]

void Nektar::MultiRegions::AssemblyMap::LocalToGlobal	(	const NekVector< NekDouble > &	loc,
		NekVector< NekDouble > &	global,
		bool	useComm = true
	)		const

Definition at line 771 of file AssemblyMap.cpp.

References v_LocalToGlobal().

        {
            v_LocalToGlobal(loc,global,useComm);
        }

PrintStats()

void Nektar::MultiRegions::AssemblyMap::PrintStats	(	std::ostream &	out,
		std::string	variable,
		bool	printHeader = true
	)		const

Definition at line 1277 of file AssemblyMap.cpp.

References Vmath::Assmb(), CellMLToNektar.pycml::level, m_globalToUniversalBndMapUnique, m_localToGlobalBndMap, m_nextLevelLocalToGlobalMap, m_numGlobalBndCoeffs, m_numGlobalCoeffs, m_numGlobalDirBndCoeffs, m_numLocalBndCoeffs, m_numLocalCoeffs, m_numLocalDirBndCoeffs, m_session, Nektar::LibUtilities::ReduceMax, Nektar::LibUtilities::ReduceMin, Nektar::LibUtilities::ReduceSum, and UniversalAssembleBnd().

        {
            LibUtilities::CommSharedPtr vRowComm
                = m_session->GetComm()->GetRowComm();
            bool isRoot = vRowComm->GetRank() == 0;
            int n = vRowComm->GetSize();
            int i;

            // Determine number of global degrees of freedom.
            int globBndCnt = 0, globDirCnt = 0;

            for (i = 0; i < m_numGlobalBndCoeffs; ++i)
            {
                if (m_globalToUniversalBndMapUnique[i] > 0)
                {
                    globBndCnt++;

                    if (i < m_numGlobalDirBndCoeffs)
                    {
                        globDirCnt++;
                    }
                }
            }

            int globCnt = m_numGlobalCoeffs - m_numGlobalBndCoeffs + globBndCnt;

            // Calculate maximum valency
            Array<OneD, NekDouble> tmpLoc (m_numLocalBndCoeffs,  1.0);
            Array<OneD, NekDouble> tmpGlob(m_numGlobalBndCoeffs, 0.0);

            Vmath::Assmb(
                m_numLocalBndCoeffs, tmpLoc.get(),
                m_localToGlobalBndMap.get(), tmpGlob.get());
            UniversalAssembleBnd(tmpGlob);

            int totGlobDof     = globCnt;
            int totGlobBndDof  = globBndCnt;
            int totGlobDirDof  = globDirCnt;
            int totLocalDof    = m_numLocalCoeffs;
            int totLocalBndDof = m_numLocalBndCoeffs;
            int totLocalDirDof = m_numLocalDirBndCoeffs;

            int meanValence = 0;
            int maxValence = 0;
            int minValence = 10000000;
            for (int i = 0; i < m_numGlobalBndCoeffs; ++i)
            {
                if (!m_globalToUniversalBndMapUnique[i])
                {
                    continue;
                }

                if (tmpGlob[i] > maxValence)
                {
                    maxValence = tmpGlob[i];
                }
                if (tmpGlob[i] < minValence)
                {
                    minValence = tmpGlob[i];
                }
                meanValence += tmpGlob[i];
            }

            vRowComm->AllReduce(maxValence,     LibUtilities::ReduceMax);
            vRowComm->AllReduce(minValence,     LibUtilities::ReduceMin);
            vRowComm->AllReduce(meanValence,    LibUtilities::ReduceSum);
            vRowComm->AllReduce(totGlobDof,     LibUtilities::ReduceSum);
            vRowComm->AllReduce(totGlobBndDof,  LibUtilities::ReduceSum);
            vRowComm->AllReduce(totGlobDirDof,  LibUtilities::ReduceSum);
            vRowComm->AllReduce(totLocalDof,    LibUtilities::ReduceSum);
            vRowComm->AllReduce(totLocalBndDof, LibUtilities::ReduceSum);
            vRowComm->AllReduce(totLocalDirDof, LibUtilities::ReduceSum);

            meanValence /= totGlobBndDof;

            if (isRoot)
            {
                if (printHeader)
                {
                    out << "Assembly map statistics for field " << variable
                        << ":" << endl;
                }

                out << "  - Number of local/global dof             : "
                    << totLocalDof << " " << totGlobDof << endl;
                out << "  - Number of local/global boundary dof    : "
                    << totLocalBndDof << " " << totGlobBndDof << endl;
                out << "  - Number of local/global Dirichlet dof   : "
                    << totLocalDirDof << " " << totGlobDirDof << endl;
                out << "  - dof valency (min/max/mean)             : "
                    << minValence << " " << maxValence << " " << meanValence
                    << endl;

                if (n > 1)
                {
                    NekDouble mean = m_numLocalCoeffs, mean2 = mean * mean;
                    NekDouble minval = mean, maxval = mean;
                    Array<OneD, NekDouble> tmp(1);

                    for (i = 1; i < n; ++i)
                    {
                        vRowComm->Recv(i, tmp);
                        mean     += tmp[0];
                        mean2    += tmp[0]*tmp[0];

                        if (tmp[0] > maxval)
                        {
                            maxval = tmp[0];
                        }
                        if (tmp[0] < minval)
                        {
                            minval = tmp[0];
                        }
                    }

                    if (maxval > 0.1)
                    {
                        out << "  - Local dof dist. (min/max/mean/dev)     : "
                            << minval << " " << maxval << " " << (mean / n)
                            << " " << sqrt(mean2/n - mean*mean/n/n) << endl;
                    }

                    vRowComm->Block();

                    mean = minval = maxval = m_numLocalBndCoeffs;
                    mean2 = mean * mean;

                    for (i = 1; i < n; ++i)
                    {
                        vRowComm->Recv(i, tmp);
                        mean     += tmp[0];
                        mean2    += tmp[0]*tmp[0];

                        if (tmp[0] > maxval)
                        {
                            maxval = tmp[0];
                        }
                        if (tmp[0] < minval)
                        {
                            minval = tmp[0];
                        }
                    }

                    out << "  - Local bnd dof dist. (min/max/mean/dev) : "
                        << minval << " " << maxval << " " << (mean / n) << " "
                        << sqrt(mean2/n - mean*mean/n/n) << endl;
                }
            }
            else
            {
                Array<OneD, NekDouble> tmp(1);
                tmp[0] = m_numLocalCoeffs;
                vRowComm->Send(0, tmp);
                vRowComm->Block();
                tmp[0] = m_numLocalBndCoeffs;
                vRowComm->Send(0, tmp);
            }

            // Either we have no more levels in the static condensation, or we
            // are not multi-level.
            if (!m_nextLevelLocalToGlobalMap)
            {
                return;
            }

            int level = 2;
            AssemblyMapSharedPtr tmp = m_nextLevelLocalToGlobalMap;
            while (tmp->m_nextLevelLocalToGlobalMap)
            {
                tmp = tmp->m_nextLevelLocalToGlobalMap;
                ++level;
            }

            // Print out multi-level static condensation information.
            if (n > 1)
            {
                if (isRoot)
                {
                    NekDouble mean = level, mean2 = mean * mean;
                    int minval = level, maxval = level;

                    Array<OneD, NekDouble> tmpRecv(1);
                    for (i = 1; i < n; ++i)
                    {
                        vRowComm->Recv(i, tmpRecv);
                        mean  += tmpRecv[0];
                        mean2 += tmpRecv[0]*tmpRecv[0];

                        if (tmpRecv[0] > maxval)
                        {
                            maxval = (int)(tmpRecv[0] + 0.5);
                        }
                        if (tmpRecv[0] < minval)
                        {
                            minval = (int)(tmpRecv[0] + 0.5);
                        }
                    }

                    out << "  - M-level sc. dist. (min/max/mean/dev)   : "
                        << minval << " " << maxval << " " << (mean / n) << " "
                        << sqrt(mean2/n - mean*mean/n/n) << endl;
                }
                else
                {
                    Array<OneD, NekDouble> tmpSend(1);
                    tmpSend[0] = level;
                    vRowComm->Send(0, tmpSend);
                }
            }
            else
            {
                out << "  - Number of static cond. levels          : "
                    << level << endl;
            }

            if (isRoot)
            {
                out << "Stats at lowest static cond. level:" << endl;
            }
            tmp->PrintStats(out, variable, false);
        }

SetNextLevelLocalToGlobalMap()

void Nektar::MultiRegions::AssemblyMap::SetNextLevelLocalToGlobalMap

(

AssemblyMapSharedPtr

pNextLevelLocalToGlobalMap

)

UniversalAssemble() [1/3]

void Nektar::MultiRegions::AssemblyMap::UniversalAssemble

(

Array< OneD, NekDouble > &

pGlobal

)

const

Definition at line 807 of file AssemblyMap.cpp.

References v_UniversalAssemble().

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

        {
            v_UniversalAssemble(pGlobal);
        }

UniversalAssemble() [2/3]

void Nektar::MultiRegions::AssemblyMap::UniversalAssemble

(

NekVector< NekDouble > &

pGlobal

)

const

Definition at line 813 of file AssemblyMap.cpp.

References v_UniversalAssemble().

        {
            v_UniversalAssemble(pGlobal);
        }

UniversalAssemble() [3/3]

void Nektar::MultiRegions::AssemblyMap::UniversalAssemble	(	Array< OneD, NekDouble > &	pGlobal,
		int	offset
	)		const

Definition at line 819 of file AssemblyMap.cpp.

References v_UniversalAssemble().

        {
            v_UniversalAssemble(pGlobal, offset);
        }

UniversalAssembleBnd() [1/3]

void Nektar::MultiRegions::AssemblyMap::UniversalAssembleBnd

(

Array< OneD, NekDouble > &

pGlobal

)

const

Definition at line 1171 of file AssemblyMap.cpp.

References ASSERTL1, Gs::Gather(), Gs::gs_add, m_bndGsh, and m_numGlobalBndCoeffs.

Referenced by AssembleBnd(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), LocalBndToGlobal(), PrintStats(), and UniversalAssembleBnd().

        {
            ASSERTL1(pGlobal.num_elements() == m_numGlobalBndCoeffs,
                     "Wrong size.");
            Gs::Gather(pGlobal, Gs::gs_add, m_bndGsh);
        }

UniversalAssembleBnd() [2/3]

void Nektar::MultiRegions::AssemblyMap::UniversalAssembleBnd

(

NekVector< NekDouble > &

pGlobal

)

const

Definition at line 1179 of file AssemblyMap.cpp.

References Nektar::NekVector< DataType >::GetPtr(), and UniversalAssembleBnd().

        {
            UniversalAssembleBnd(pGlobal.GetPtr());
        }

UniversalAssembleBnd() [3/3]

void Nektar::MultiRegions::AssemblyMap::UniversalAssembleBnd	(	Array< OneD, NekDouble > &	pGlobal,
		int	offset
	)		const

Definition at line 1185 of file AssemblyMap.cpp.

References UniversalAssembleBnd(), and Vmath::Vcopy().

        {
            Array<OneD, NekDouble> tmp(offset);
            if (offset > 0)  Vmath::Vcopy(offset, pGlobal, 1, tmp, 1);
            UniversalAssembleBnd(pGlobal);
            if (offset > 0)  Vmath::Vcopy(offset, tmp, 1, pGlobal, 1);
        }

v_Assemble() [1/2]

void Nektar::MultiRegions::AssemblyMap::v_Assemble ( const Array< OneD, const NekDouble > &  loc, Array< OneD, NekDouble > &  global  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 596 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by Assemble().

        {
            boost::ignore_unused(loc, global);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
        }

v_Assemble() [2/2]

void Nektar::MultiRegions::AssemblyMap::v_Assemble ( const NekVector< NekDouble > &  loc, NekVector< NekDouble > &  global  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 605 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(loc, global);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
        }

v_GetExtraDirEdges()

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::v_GetExtraDirEdges ( )

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 695 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetExtraDirEdges().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            static Array<OneD, const int> result;
            return result;
        }

v_GetFullSystemBandWidth()

int Nektar::MultiRegions::AssemblyMap::v_GetFullSystemBandWidth ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 639 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetFullSystemBandWidth().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetGlobalToUniversalMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::v_GetGlobalToUniversalMap ( const int  i ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 502 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(i);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetGlobalToUniversalMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::v_GetGlobalToUniversalMap ( void  )

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 526 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetGlobalToUniversalMap().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            static Array<OneD, const int> result;
            return result;
        }

v_GetGlobalToUniversalMapUnique() [1/2]

int Nektar::MultiRegions::AssemblyMap::v_GetGlobalToUniversalMapUnique ( const int  i ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 510 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(i);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetGlobalToUniversalMapUnique() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::v_GetGlobalToUniversalMapUnique ( void  )

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 534 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetGlobalToUniversalMapUnique().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            static Array<OneD, const int> result;
            return result;
        }

v_GetLocalToGlobalMap() [1/2]

int Nektar::MultiRegions::AssemblyMap::v_GetLocalToGlobalMap ( const int  i ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 494 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(i);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetLocalToGlobalMap() [2/2]

const Array< OneD, const int > & Nektar::MultiRegions::AssemblyMap::v_GetLocalToGlobalMap ( void  )

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 518 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetLocalToGlobalMap().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            static Array<OneD,const int> result;
            return result;
        }

v_GetLocalToGlobalSign() [1/2]

NekDouble Nektar::MultiRegions::AssemblyMap::v_GetLocalToGlobalSign ( const int  i ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 542 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(i);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0.0;
        }

v_GetLocalToGlobalSign() [2/2]

const Array< OneD, NekDouble > & Nektar::MultiRegions::AssemblyMap::v_GetLocalToGlobalSign ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 550 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetLocalToGlobalSign().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            static Array<OneD, NekDouble> result;
            return result;
        }

v_GetNumDirEdges()

int Nektar::MultiRegions::AssemblyMap::v_GetNumDirEdges ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 667 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumDirEdges().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetNumDirFaces()

int Nektar::MultiRegions::AssemblyMap::v_GetNumDirFaces ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 674 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumDirFaces().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetNumNonDirEdgeModes()

int Nektar::MultiRegions::AssemblyMap::v_GetNumNonDirEdgeModes ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 653 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumNonDirEdgeModes().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetNumNonDirEdges()

int Nektar::MultiRegions::AssemblyMap::v_GetNumNonDirEdges ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 681 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumNonDirEdges().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetNumNonDirFaceModes()

int Nektar::MultiRegions::AssemblyMap::v_GetNumNonDirFaceModes ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 660 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumNonDirFaceModes().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetNumNonDirFaces()

int Nektar::MultiRegions::AssemblyMap::v_GetNumNonDirFaces ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 688 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumNonDirFaces().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GetNumNonDirVertexModes()

int Nektar::MultiRegions::AssemblyMap::v_GetNumNonDirVertexModes ( ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 646 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GetNumNonDirVertexModes().

        {
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            return 0;
        }

v_GlobalToLocal() [1/2]

void Nektar::MultiRegions::AssemblyMap::v_GlobalToLocal ( const Array< OneD, const NekDouble > &  global, Array< OneD, NekDouble > &  loc  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 578 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by GlobalToLocal().

        {
            boost::ignore_unused(loc, global);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
        }

v_GlobalToLocal() [2/2]

void Nektar::MultiRegions::AssemblyMap::v_GlobalToLocal ( const NekVector< NekDouble > &  global, NekVector< NekDouble > &  loc  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 587 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(loc, global);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
        }

v_LinearSpaceMap()

std::shared_ptr< AssemblyMap > Nektar::MultiRegions::AssemblyMap::v_LinearSpaceMap ( const ExpList &  locexp, GlobalSysSolnType  solnType  )

privatevirtual

Generate a linear space mapping from existing mapping.

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 703 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by LinearSpaceMap().

        {
            boost::ignore_unused(locexp, solnType);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
            static std::shared_ptr<AssemblyMap> result;
            return result;
        }

v_LocalToGlobal() [1/2]

void Nektar::MultiRegions::AssemblyMap::v_LocalToGlobal ( const Array< OneD, const NekDouble > &  loc, Array< OneD, NekDouble > &  global, bool  useComm  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 558 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

Referenced by LocalToGlobal().

        {
            boost::ignore_unused(loc, global, useComm);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
        }

v_LocalToGlobal() [2/2]

void Nektar::MultiRegions::AssemblyMap::v_LocalToGlobal ( const NekVector< NekDouble > &  loc, NekVector< NekDouble > &  global, bool  useComm  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 568 of file AssemblyMap.cpp.

References Nektar::ErrorUtil::efatal, and NEKERROR.

        {
            boost::ignore_unused(loc, global, useComm);
            NEKERROR(ErrorUtil::efatal,
                     "Not defined for this type of mapping.");
        }

v_UniversalAssemble() [1/3]

void Nektar::MultiRegions::AssemblyMap::v_UniversalAssemble ( Array< OneD, NekDouble > &  pGlobal ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 614 of file AssemblyMap.cpp.

Referenced by UniversalAssemble().

        {
            boost::ignore_unused(pGlobal);
            // Do nothing here since multi-level static condensation uses a
            // AssemblyMap and thus will call this routine in serial.
        }

v_UniversalAssemble() [2/3]

void Nektar::MultiRegions::AssemblyMap::v_UniversalAssemble ( NekVector< NekDouble > &  pGlobal ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG, and Nektar::MultiRegions::AssemblyMapDG.

Definition at line 622 of file AssemblyMap.cpp.

        {
            boost::ignore_unused(pGlobal);
            // Do nothing here since multi-level static condensation uses a
            // AssemblyMap and thus will call this routine in serial.
        }

v_UniversalAssemble() [3/3]

void Nektar::MultiRegions::AssemblyMap::v_UniversalAssemble ( Array< OneD, NekDouble > &  pGlobal, int  offset  ) const

privatevirtual

Reimplemented in Nektar::MultiRegions::AssemblyMapCG.

Definition at line 630 of file AssemblyMap.cpp.

        {
            boost::ignore_unused(pGlobal, offset);
            // Do nothing here since multi-level static condensation uses a
            // AssemblyMap and thus will call this routine in serial.
        }

Member Data Documentation

m_bndCondCoeffsToGlobalCoeffsMap

Array<OneD,int> Nektar::MultiRegions::AssemblyMap::m_bndCondCoeffsToGlobalCoeffsMap

protected

Integer map of bnd cond coeffs to global coefficients.

Definition at line 353 of file AssemblyMap.h.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), and GetBndCondCoeffsToGlobalCoeffsMap().

m_bndCondCoeffsToGlobalCoeffsSign

Array<OneD,NekDouble> Nektar::MultiRegions::AssemblyMap::m_bndCondCoeffsToGlobalCoeffsSign

protected

Integer map of bnd cond coeffs to global coefficients.

Definition at line 355 of file AssemblyMap.h.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), and GetBndCondCoeffsToGlobalCoeffsSign().

m_bndCondTraceToGlobalTraceMap

Array<OneD,int> Nektar::MultiRegions::AssemblyMap::m_bndCondTraceToGlobalTraceMap

protected

Integer map of bnd cond trace number to global trace number.

Definition at line 357 of file AssemblyMap.h.

Referenced by Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), and GetBndCondTraceToGlobalTraceMap().

m_bndGsh

Gs::gs_data* Nektar::MultiRegions::AssemblyMap::m_bndGsh

protected

Definition at line 381 of file AssemblyMap.h.

Referenced by Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), UniversalAssembleBnd(), and Nektar::MultiRegions::AssemblyMapCG::~AssemblyMapCG().

m_bndSystemBandWidth

int Nektar::MultiRegions::AssemblyMap::m_bndSystemBandWidth

protected

The bandwith of the global bnd system.

Definition at line 366 of file AssemblyMap.h.

Referenced by CalculateBndSystemBandWidth(), and GetBndSystemBandWidth().

m_comm

LibUtilities::CommSharedPtr Nektar::MultiRegions::AssemblyMap::m_comm

protected

Communicator.

Definition at line 307 of file AssemblyMap.h.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::CreateGraph(), GetComm(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalTraceMap(), and Nektar::MultiRegions::AssemblyMapCG::v_LinearSpaceMap().

m_globalToUniversalBndMap

Array<OneD,int> Nektar::MultiRegions::AssemblyMap::m_globalToUniversalBndMap

protected

Integer map of process coeffs to universal space.

Definition at line 359 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), GetGlobalToUniversalBndMap(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), and Nektar::MultiRegions::AssemblyMapDG::v_GetGlobalToUniversalMap().

m_globalToUniversalBndMapUnique

Array<OneD,int> Nektar::MultiRegions::AssemblyMap::m_globalToUniversalBndMapUnique

protected

Integer map of unique process coeffs to universal space (signed)

Definition at line 361 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), GetGlobalToUniversalBndMapUnique(), PrintStats(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), and Nektar::MultiRegions::AssemblyMapDG::v_GetGlobalToUniversalMapUnique().

m_gsh

Gs::gs_data* Nektar::MultiRegions::AssemblyMap::m_gsh

protected

Definition at line 380 of file AssemblyMap.h.

Referenced by Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), Nektar::MultiRegions::AssemblyMapCG::v_LocalToGlobal(), Nektar::MultiRegions::AssemblyMapDG::v_UniversalAssemble(), Nektar::MultiRegions::AssemblyMapCG::v_UniversalAssemble(), and Nektar::MultiRegions::AssemblyMapCG::~AssemblyMapCG().

m_hash

size_t Nektar::MultiRegions::AssemblyMap::m_hash

protected

Hash for map.

Definition at line 310 of file AssemblyMap.h.

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

m_iterativeTolerance

NekDouble Nektar::MultiRegions::AssemblyMap::m_iterativeTolerance

protected

Tolerance for iterative solver.

Definition at line 375 of file AssemblyMap.h.

Referenced by AssemblyMap(), and GetIterativeTolerance().

m_localToGlobalBndMap

Array<OneD,int> Nektar::MultiRegions::AssemblyMap::m_localToGlobalBndMap

protected

Integer map of local boundary coeffs to global space.

Definition at line 349 of file AssemblyMap.h.

Referenced by AssembleBnd(), AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), CalculateBndSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetLocalToGlobalBndMap(), GlobalToLocalBnd(), GlobalToLocalBndWithoutSign(), LocalBndToGlobal(), PrintStats(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), and Nektar::MultiRegions::AssemblyMapDG::v_GetLocalToGlobalMap().

m_localToGlobalBndSign

Array<OneD,NekDouble> Nektar::MultiRegions::AssemblyMap::m_localToGlobalBndSign

protected

Integer sign of local boundary coeffs to global space.

Definition at line 351 of file AssemblyMap.h.

Referenced by AssembleBnd(), AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetLocalToGlobalBndSign(), GlobalToLocalBnd(), and LocalBndToGlobal().

m_lowestStaticCondLevel

int Nektar::MultiRegions::AssemblyMap::m_lowestStaticCondLevel

protected

Lowest static condensation level.

Definition at line 398 of file AssemblyMap.h.

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

m_maxIterations

int Nektar::MultiRegions::AssemblyMap::m_maxIterations

protected

Maximum iterations for iterative solver.

Definition at line 372 of file AssemblyMap.h.

Referenced by AssemblyMap(), and GetMaxIterations().

m_nextLevelLocalToGlobalMap

AssemblyMapSharedPtr Nektar::MultiRegions::AssemblyMap::m_nextLevelLocalToGlobalMap

protected

Map from the patches of the previous level to the patches of the current level.

The local to global mapping of the next level of recursion

Definition at line 396 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), AtLastLevel(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetNextLevelLocalToGlobalMap(), and PrintStats().

m_numGlobalBndCoeffs

int Nektar::MultiRegions::AssemblyMap::m_numGlobalBndCoeffs

protected

Total number of global boundary coefficients.

Definition at line 315 of file AssemblyMap.h.

Referenced by AssembleBnd(), AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetNumGlobalBndCoeffs(), GlobalToLocalBnd(), GlobalToLocalBndWithoutSign(), LocalBndToGlobal(), PrintStats(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapDG::SetUpUniversalDGMap(), and UniversalAssembleBnd().

m_numGlobalCoeffs

int Nektar::MultiRegions::AssemblyMap::m_numGlobalCoeffs

protected

Total number of global coefficients.

This corresponds to the number of total number of coefficients

• For CG this corresponds to the total of bnd + int DOFs.
• For DG this corresponds to the number of bnd DOFs. This means that m_numGlobalCoeffs = m_numGlobalBndCoeffs This way, we can consider the trace-system solve as a statically condensed solve without interior DOFs. This allows us to use the same global system solver for both cases.

Definition at line 343 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetNumGlobalCoeffs(), PrintStats(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapCG::v_Assemble(), and Nektar::MultiRegions::AssemblyMapCG::v_LinearSpaceMap().

m_numGlobalDirBndCoeffs

int Nektar::MultiRegions::AssemblyMap::m_numGlobalDirBndCoeffs

protected

Number of Global Dirichlet Boundary Coefficients.

Definition at line 319 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), CalculateBndSystemBandWidth(), Nektar::MultiRegions::AssemblyMapCG::CalculateFullSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetNumGlobalDirBndCoeffs(), PrintStats(), and Nektar::MultiRegions::AssemblyMapCG::v_LinearSpaceMap().

m_numLocalBndCoeffs

int Nektar::MultiRegions::AssemblyMap::m_numLocalBndCoeffs

protected

Number of local boundary coefficients.

Definition at line 313 of file AssemblyMap.h.

Referenced by AssembleBnd(), AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), CalculateBndSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), Nektar::MultiRegions::AssemblyMapCG::CreateGraph(), GetNumLocalBndCoeffs(), GlobalToLocalBnd(), GlobalToLocalBndWithoutSign(), LocalBndToGlobal(), and PrintStats().

m_numLocalBndCoeffsPerPatch

Array<OneD, unsigned int> Nektar::MultiRegions::AssemblyMap::m_numLocalBndCoeffsPerPatch

protected

The number of bnd dofs per patch.

Definition at line 389 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), CalculateBndSystemBandWidth(), Nektar::MultiRegions::AssemblyMapCG::CalculateFullSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), and GetNumLocalBndCoeffsPerPatch().

m_numLocalCoeffs

int Nektar::MultiRegions::AssemblyMap::m_numLocalCoeffs

protected

Total number of local coefficients.

This corresponds to the number of total number of coefficients

• For CG this corresponds to the total of bnd + int DOFs
• For DG this corresponds to the number of bnd DOFs. This means that m_numLocalCoeffs = m_numLocalBndCoeffs This way, we can consider the trace-system solve as a statically condensed solve without interior DOFs. This allows us to use the same global system solver for both cases.

Definition at line 332 of file AssemblyMap.h.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::MultiRegions::AssemblyMapCG::CalculateFullSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetNumLocalCoeffs(), PrintStats(), Nektar::MultiRegions::AssemblyMapCG::v_Assemble(), Nektar::MultiRegions::AssemblyMapCG::v_GlobalToLocal(), and Nektar::MultiRegions::AssemblyMapCG::v_LocalToGlobal().

m_numLocalDirBndCoeffs

int Nektar::MultiRegions::AssemblyMap::m_numLocalDirBndCoeffs

protected

Number of Local Dirichlet Boundary Coefficients.

Definition at line 317 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), Nektar::MultiRegions::AssemblyMapCG::CreateGraph(), GetNumLocalDirBndCoeffs(), and PrintStats().

m_numLocalIntCoeffsPerPatch

Array<OneD, unsigned int> Nektar::MultiRegions::AssemblyMap::m_numLocalIntCoeffsPerPatch

protected

The number of int dofs per patch.

Definition at line 391 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::MultiRegions::AssemblyMapCG::CalculateFullSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), and GetNumLocalIntCoeffsPerPatch().

m_numPatches

int Nektar::MultiRegions::AssemblyMap::m_numPatches

protected

The number of patches (~elements) in the current level.

Definition at line 387 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), CalculateBndSystemBandWidth(), Nektar::MultiRegions::AssemblyMapCG::CalculateFullSystemBandWidth(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), and GetNumPatches().

m_patchMapFromPrevLevel

PatchMapSharedPtr Nektar::MultiRegions::AssemblyMap::m_patchMapFromPrevLevel

private

Mapping information for previous level in MultiLevel Solver.

Definition at line 409 of file AssemblyMap.h.

Referenced by AssemblyMap(), and GetPatchMapFromPrevLevel().

m_preconType

PreconditionerType Nektar::MultiRegions::AssemblyMap::m_preconType

protected

Type type of preconditioner to use in iterative solver.

Definition at line 369 of file AssemblyMap.h.

Referenced by AssemblyMap(), and GetPreconType().

m_session

LibUtilities::SessionReaderSharedPtr Nektar::MultiRegions::AssemblyMap::m_session

protected

Session object.

Definition at line 304 of file AssemblyMap.h.

Referenced by Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), Nektar::MultiRegions::AssemblyMapCG::CreateGraph(), PrintStats(), and Nektar::MultiRegions::AssemblyMapCG::v_LinearSpaceMap().

m_signChange

bool Nektar::MultiRegions::AssemblyMap::m_signChange

protected

Flag indicating if modes require sign reversal.

Definition at line 346 of file AssemblyMap.h.

Referenced by AssembleBnd(), AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), GetBndCondCoeffsToGlobalCoeffsSign(), GetLocalToGlobalBndSign(), GetSignChange(), GlobalToLocalBnd(), LocalBndToGlobal(), Nektar::MultiRegions::AssemblyMapCG::SetUpUniversalC0ContMap(), Nektar::MultiRegions::AssemblyMapCG::v_Assemble(), Nektar::MultiRegions::AssemblyMapCG::v_GetLocalToGlobalSign(), Nektar::MultiRegions::AssemblyMapCG::v_GlobalToLocal(), and Nektar::MultiRegions::AssemblyMapCG::v_LocalToGlobal().

m_solnType

GlobalSysSolnType Nektar::MultiRegions::AssemblyMap::m_solnType

protected

The solution type of the global system.

Definition at line 364 of file AssemblyMap.h.

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

m_staticCondLevel

int Nektar::MultiRegions::AssemblyMap::m_staticCondLevel

protected

The level of recursion in the case of multi-level static condensation.

Definition at line 385 of file AssemblyMap.h.

Referenced by AssemblyMap(), Nektar::MultiRegions::AssemblyMapCG::AssemblyMapCG(), Nektar::MultiRegions::AssemblyMapDG::AssemblyMapDG(), Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), and GetStaticCondLevel().

m_successiveRHS

int Nektar::MultiRegions::AssemblyMap::m_successiveRHS

protected

sucessive RHS for iterative solver

Definition at line 378 of file AssemblyMap.h.

Referenced by AssemblyMap(), and GetSuccessiveRHS().

m_systemSingular

bool Nektar::MultiRegions::AssemblyMap::m_systemSingular

protected

Flag indicating if the system is singular or not.

Definition at line 321 of file AssemblyMap.h.

Referenced by Nektar::CoupledAssemblyMap::CoupledAssemblyMap(), Nektar::CoupledLocalToGlobalC0ContMap::CoupledLocalToGlobalC0ContMap(), Nektar::MultiRegions::AssemblyMapCG::CreateGraph(), and GetSingularSystem().