Modified version of MultiRegions::AssemblyMapCG that allows for coupled fields [u,v,w] instead of individual scalar fields u, v and w. More...

#include <CoupledAssemblyMap.h>

Inheritance diagram for Nektar::CoupledAssemblyMap:

[legend]

Collaboration diagram for Nektar::CoupledAssemblyMap:

[legend]

Public Member Functions
	CoupledAssemblyMap (const LibUtilities::SessionReaderSharedPtr &pSession, const SpatialDomains::MeshGraphSharedPtr &graph, const MultiRegions::AssemblyMapCGSharedPtr &cgMap, const Array< OneD, const BoundaryCondShPtr > &boundaryConditions, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
	Take an existing assembly map and create a coupled version suitable for use in the linear elasticity solver. More...

Public Member Functions inherited from Nektar::MultiRegions::AssemblyMapCG
	AssemblyMapCG (const LibUtilities::SessionReaderSharedPtr &pSession, const std::string variable="DefaultVar")
	Default constructor. More...

	AssemblyMapCG (const LibUtilities::SessionReaderSharedPtr &pSession, const int numLocalCoeffs, const ExpList &locExp, const BndCondExp &bndCondExp=NullExpListSharedPtrArray, const BndCond &bndConditions=SpatialDomains::NullBoundaryConditionShPtrArray, const bool checkIfSingular=false, const std::string variable="defaultVar", const PeriodicMap &periodicVerts=NullPeriodicMap, const PeriodicMap &periodicEdges=NullPeriodicMap, const PeriodicMap &periodicFaces=NullPeriodicMap)
	General constructor for expansions of all dimensions without boundary conditions. More...

virtual	~AssemblyMapCG ()
	Destructor. More...

std::map< int, std::vector < ExtraDirDof > > &	GetExtraDirDofs ()

Public Member Functions inherited from Nektar::MultiRegions::AssemblyMap
	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 ()

boost::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

Private Types
typedef SpatialDomains::BoundaryConditionShPtr	BoundaryCondShPtr

Additional Inherited Members
Protected Member Functions inherited from Nektar::MultiRegions::AssemblyMapCG
int	CreateGraph (const ExpList &locExp, const BndCondExp &bndCondExp, const Array< OneD, const BndCond > &bndConditions, const bool checkIfSystemSingular, const PeriodicMap &periodicVerts, const PeriodicMap &periodicEdges, const PeriodicMap &periodicFaces, DofGraph &graph, BottomUpSubStructuredGraphSharedPtr &bottomUpGraph, std::set< int > &extraDirVerts, std::set< int > &extraDirEdges, int &firstNonDirGraphVertId, int &nExtraDirichlet, int mdswitch=1)

void	SetUpUniversalC0ContMap (const ExpList &locExp, const PeriodicMap &perVerts=NullPeriodicMap, const PeriodicMap &perEdges=NullPeriodicMap, const PeriodicMap &perFaces=NullPeriodicMap)

void	CalculateFullSystemBandWidth ()
	Calculate the bandwith of the full matrix system. More...

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 AssemblyMapSharedPtr	v_LinearSpaceMap (const ExpList &locexp, GlobalSysSolnType solnType)
	Construct an AssemblyMapCG object which corresponds to the linear space of the current object. More...

Protected Member Functions inherited from Nektar::MultiRegions::AssemblyMap
void	CalculateBndSystemBandWidth ()
	Calculates the bandwidth of the boundary system. More...

void	GlobalToLocalBndWithoutSign (const Array< OneD, const NekDouble > &global, Array< OneD, NekDouble > &loc)

Protected Attributes inherited from Nektar::MultiRegions::AssemblyMapCG
Array< OneD, int >	m_localToGlobalMap
	Integer map of local coeffs to global space. More...

Array< OneD, NekDouble >	m_localToGlobalSign
	Integer sign of local coeffs to global space. More...

int	m_fullSystemBandWidth
	Bandwith of the full matrix system (no static condensation). More...

Array< OneD, int >	m_globalToUniversalMap
	Integer map of process coeffs to universal space. More...

Array< OneD, int >	m_globalToUniversalMapUnique
	Integer map of unique process coeffs to universal space (signed) More...

int	m_numNonDirVertexModes
	Number of non Dirichlet vertex modes. More...

int	m_numNonDirEdgeModes
	Number of non Dirichlet edge modes. More...

int	m_numNonDirFaceModes
	Number of non Dirichlet face modes. More...

int	m_numDirEdges
	Number of Dirichlet edges. More...

int	m_numDirFaces
	Number of Dirichlet faces. More...

int	m_numNonDirEdges
	Number of Dirichlet edges. More...

int	m_numNonDirFaces
	Number of Dirichlet faces. More...

int	m_numLocalBndCondCoeffs
	Number of local boundary condition coefficients. More...

Array< OneD, int >	m_extraDirEdges
	Extra dirichlet edges in parallel. More...

int	m_numLocDirBndCondDofs
	Number of local boundary condition degrees of freedom. More...

int	m_maxStaticCondLevel
	Maximum static condensation level. More...

std::map< int, std::vector < ExtraDirDof > >	m_extraDirDofs
	Map indicating degrees of freedom which are Dirichlet but whose value is stored on another processor. More...

Protected Attributes inherited from Nektar::MultiRegions::AssemblyMap
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...

Detailed Description

Modified version of MultiRegions::AssemblyMapCG that allows for coupled fields [u,v,w] instead of individual scalar fields u, v and w.

Definition at line 50 of file CoupledAssemblyMap.h.

Member Typedef Documentation

typedef SpatialDomains::BoundaryConditionShPtr Nektar::CoupledAssemblyMap::BoundaryCondShPtr

private

Definition at line 52 of file CoupledAssemblyMap.h.

Constructor & Destructor Documentation

Nektar::CoupledAssemblyMap::CoupledAssemblyMap	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const SpatialDomains::MeshGraphSharedPtr &	graph,
		const MultiRegions::AssemblyMapCGSharedPtr &	cgMap,
		const Array< OneD, const BoundaryCondShPtr > &	boundaryConditions,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields
	)

Take an existing assembly map and create a coupled version suitable for use in the linear elasticity solver.

The linear elasticity solver requires a slight reordering of local and global coefficients to support problems of the form

[ A B C ] [ u ] = [ f_u ] [ D E F ] [ v ] [ f_v ] [ G H I ] [ w ] [ f_w ]

Definition at line 64 of file CoupledAssemblyMap.cpp.

                                                              :
     AssemblyMapCG(pSession)
 {
     int nVel = fields[0]->GetCoordim(0);
 
     // Multi-level static condensation doesn't work yet.
     ASSERTL0(m_solnType != MultiRegions::eDirectMultiLevelStaticCond    &&
              m_solnType != MultiRegions::eIterativeMultiLevelStaticCond &&
              m_solnType != MultiRegions::eXxtMultiLevelStaticCond,
              "Multi-level static condensation not supported.");
 
     // Copy various coefficient counts, and multiply by the dimension of the
     // problem to obtain our new values.
     m_numLocalDirBndCoeffs      = cgMap->GetNumLocalDirBndCoeffs()  * nVel;
     m_numLocalBndCoeffs         = cgMap->GetNumLocalBndCoeffs()     * nVel;
     m_numLocalCoeffs            = cgMap->GetNumLocalCoeffs()        * nVel;
     m_numGlobalBndCoeffs        = cgMap->GetNumGlobalBndCoeffs()    * nVel;
     m_numGlobalDirBndCoeffs     = cgMap->GetNumGlobalDirBndCoeffs() * nVel;
     m_numGlobalCoeffs           = cgMap->GetNumGlobalCoeffs()       * nVel;
     m_signChange                = cgMap->GetSignChange();
     m_systemSingular            = cgMap->GetSingularSystem();
 
     // Copy static condensation information. TODO: boundary and interior patches
     // need to be re-ordered in order to allow for multi-level static
     // condensation support.
     m_staticCondLevel           = cgMap->GetStaticCondLevel();
     m_lowestStaticCondLevel     = cgMap->GetLowestStaticCondLevel();
     m_numPatches                = cgMap->GetNumPatches();
     m_numLocalBndCoeffsPerPatch = cgMap->GetNumLocalBndCoeffsPerPatch();
     m_numLocalIntCoeffsPerPatch = cgMap->GetNumLocalIntCoeffsPerPatch();
 
     // Set up local to global and boundary condition maps.
     const int nLocBndCondDofs = cgMap->
         GetBndCondCoeffsToGlobalCoeffsMap().num_elements() * nVel;
 
     // Allocate storage for local to global maps.
     m_localToGlobalMap               =
         Array<OneD, int>(m_numLocalCoeffs,-1);
     m_localToGlobalBndMap            =
         Array<OneD, int>(m_numLocalBndCoeffs,-1);
     m_bndCondCoeffsToGlobalCoeffsMap =
         Array<OneD, int>(nLocBndCondDofs,-1);
 
     // Only require a sign map if we are using modal polynomials in the
     // expansion and the order is >= 3.
     if (m_signChange)
     {
         m_localToGlobalSign               =
             Array<OneD, NekDouble>(m_numLocalCoeffs,1.0);
         m_localToGlobalBndSign            =
             Array<OneD, NekDouble>(m_numLocalBndCoeffs,1.0);
         m_bndCondCoeffsToGlobalCoeffsSign =
             Array<OneD, NekDouble>(nLocBndCondDofs,1.0);
     }
     else
     {
         m_localToGlobalSign               = NullNekDouble1DArray;
         m_localToGlobalBndSign            = NullNekDouble1DArray;
         m_bndCondCoeffsToGlobalCoeffsSign = NullNekDouble1DArray;
     }
 
     const LocalRegions::ExpansionVector &locExpVector
         = *(fields[0]->GetExp());
 
     map<int, int> newGlobalIds;
     int i, j, n, cnt1, cnt2;
 
     // Order local boundary degrees of freedom. These are basically fine; we
     // reorder storage so that we loop over each element and then each component
     // of velocity, by applying a mapping l2g -> nVel*l2g + n, for 0 <= n <
     // nVel. Note that Dirichlet ordering is preserved under this
     // transformation.
     cnt1 = cnt2 = 0;
     for (i = 0; i < locExpVector.size(); ++i)
     {
         const int nBndCoeffs = locExpVector[i]->NumBndryCoeffs();
 
         for (n = 0; n < nVel; ++n)
         {
             for (j = 0; j < nBndCoeffs; ++j, ++cnt1)
             {
                 const int l2g = cgMap->GetLocalToGlobalBndMap()[cnt2+j];
                 m_localToGlobalBndMap[cnt1] = nVel * l2g + n;
 
                 if (m_signChange)
                 {
                     m_localToGlobalBndSign[cnt1] =
                         cgMap->GetLocalToGlobalBndSign()[cnt2+j];
                 }
 
                 if (n == 0)
                 {
                     const int l2gnew = m_localToGlobalBndMap[cnt1];
                     if (newGlobalIds.count(l2g))
                     {
                         ASSERTL1(newGlobalIds[l2g] == l2gnew,
                                  "Consistency error");
                     }
                     newGlobalIds[l2g] = l2gnew;
                 }
             }
         }
 
         cnt2 += nBndCoeffs;
     }
 
     // Grab map of extra Dirichlet degrees of freedom for parallel runs
     // (particularly in 3D).
     m_extraDirDofs = cgMap->GetExtraDirDofs();
 
     // Counter for remaining interior degrees of freedom.
     int globalId = m_numGlobalBndCoeffs;
 
     // Interior degrees of freedom are a bit more tricky -- global linear system
     // solve relies on them being in the same order as the BinvD, C and invD
     // matrices.
     cnt1 = cnt2 = 0;
     for (i = 0; i < locExpVector.size(); ++i)
     {
         const int nCoeffs    = locExpVector[i]->GetNcoeffs();
         const int nBndCoeffs = locExpVector[i]->NumBndryCoeffs();
 
         for (n = 0; n < nVel; ++n)
         {
             for (j = 0; j < nBndCoeffs; ++j, ++cnt1, ++cnt2)
             {
                 const int l2g = m_localToGlobalBndMap[cnt2];
                 m_localToGlobalMap[cnt1] = l2g;
 
                 if (m_signChange)
                 {
                     m_localToGlobalSign[cnt1] = m_localToGlobalBndSign[cnt2];
                 }
             }
         }
 
         for (n = 0; n < nVel; ++n)
         {
             for (j = 0; j < nCoeffs - nBndCoeffs; ++j, ++cnt1)
             {
                 m_localToGlobalMap[cnt1] = globalId++;
             }
         }
     }
 
     for (i = 0; i < m_localToGlobalMap.num_elements(); ++i)
     {
         ASSERTL1(m_localToGlobalMap[i] != -1, "Consistency error");
     }
 
     ASSERTL1(globalId == m_numGlobalCoeffs, "Consistency error");
 
     cnt1 = 0;
     for (n = 0; n < nVel; ++n)
     {
         for (i = 0; i < nLocBndCondDofs/nVel; ++i, ++cnt1)
         {
             const int l2g = cgMap->GetBndCondCoeffsToGlobalCoeffsMap()[i];
             int newId = newGlobalIds[l2g];
             m_bndCondCoeffsToGlobalCoeffsMap[cnt1] = newId + n;
 
             if (m_signChange)
             {
                 m_bndCondCoeffsToGlobalCoeffsSign[cnt1] =
                     cgMap->GetBndCondCoeffsToGlobalCoeffsSign(i);
             }
         }
     }
 
     // Finally, set up global to universal maps.
     m_globalToUniversalMap          = Array<OneD, int>(m_numGlobalCoeffs);
     m_globalToUniversalMapUnique    = Array<OneD, int>(m_numGlobalCoeffs);
     m_globalToUniversalBndMap       = Array<OneD, int>(m_numGlobalBndCoeffs);
     m_globalToUniversalBndMapUnique = Array<OneD, int>(m_numGlobalBndCoeffs);
 
     for (i = 0; i < cgMap->GetNumGlobalBndCoeffs(); ++i)
     {
         for (n = 0; n < nVel; ++n)
         {
             m_globalToUniversalBndMap[i*nVel + n] =
                 cgMap->GetGlobalToUniversalBndMap()[i]*nVel + n;
             m_globalToUniversalMap[i*nVel + n] =
                 cgMap->GetGlobalToUniversalBndMap()[i]*nVel + n;
         }
     }
 
     Array<OneD, long> tmp(m_numGlobalCoeffs);
     Vmath::Zero(m_numGlobalCoeffs, tmp, 1);
     Array<OneD, long> tmp2(m_numGlobalBndCoeffs, tmp);
     for (unsigned int i = 0; i < m_numGlobalBndCoeffs; ++i)
     {
         tmp[i] = m_globalToUniversalBndMap[i];
     }
 
     LibUtilities::CommSharedPtr vCommRow = m_comm->GetRowComm();
     m_gsh    = Gs::Init(tmp,  vCommRow);
     m_bndGsh = Gs::Init(tmp2, vCommRow);
     Gs::Unique(tmp, vCommRow);
     for (unsigned int i = 0; i < m_numGlobalCoeffs; ++i)
     {
         m_globalToUniversalMapUnique[i] = (tmp[i] >= 0 ? 1 : 0);
     }
     for (unsigned int i = 0; i < m_numGlobalBndCoeffs; ++i)
     {
         m_globalToUniversalBndMapUnique[i] = (tmp2[i] >= 0 ? 1 : 0);
     }
 
     m_hash = boost::hash_range(
         m_localToGlobalMap.begin(), m_localToGlobalMap.end());
 }

Public Member Functions

Private Types

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation