A global linear system. More...

#include <GlobalLinSys.h>

Inheritance diagram for Nektar::MultiRegions::GlobalLinSys:

Collaboration diagram for Nektar::MultiRegions::GlobalLinSys:

Public Member Functions
	GlobalLinSys (const GlobalLinSysKey &pKey, const boost::weak_ptr< ExpList > &pExpList, const boost::shared_ptr< AssemblyMap > &pLocToGloMap)
	Constructor for full direct matrix solve. More...

virtual	~GlobalLinSys ()

const GlobalLinSysKey &	GetKey (void) const
	Returns the key associated with the system. More...

const boost::weak_ptr< ExpList > &	GetLocMat (void) const

void	InitObject ()

void	Initialise (const boost::shared_ptr< AssemblyMap > &pLocToGloMap)

void	Solve (const Array< OneD, const NekDouble > &in, Array< OneD, NekDouble > &out, const AssemblyMapSharedPtr &locToGloMap, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)
	Solve the linear system for given input and output vectors using a specified local to global map. More...

boost::shared_ptr< GlobalLinSys >	GetSharedThisPtr ()
	Returns a shared pointer to the current object. More...

int	GetNumBlocks ()

DNekScalMatSharedPtr	GetBlock (unsigned int n)

DNekScalBlkMatSharedPtr	GetStaticCondBlock (unsigned int n)

void	DropStaticCondBlock (unsigned int n)

void	SolveLinearSystem (const int pNumRows, const Array< OneD, const NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const AssemblyMapSharedPtr &locToGloMap, const int pNumDir=0)
	Solve the linear system for given input and output vectors. More...

Protected Member Functions
virtual int	v_GetNumBlocks ()
	Get the number of blocks in this system. More...

virtual DNekScalMatSharedPtr	v_GetBlock (unsigned int n)
	Retrieves the block matrix from n-th expansion using the matrix key provided by the m_linSysKey. More...

virtual DNekScalBlkMatSharedPtr	v_GetStaticCondBlock (unsigned int n)
	Retrieves a the static condensation block matrices from n-th expansion using the matrix key provided by the m_linSysKey. More...

virtual void	v_DropStaticCondBlock (unsigned int n)
	Releases the static condensation block matrices from NekManager of n-th expansion using the matrix key provided by the m_linSysKey. More...

PreconditionerSharedPtr	CreatePrecon (AssemblyMapSharedPtr asmMap)
	Create a preconditioner object from the parameters defined in the supplied assembly map. More...

Protected Attributes
const GlobalLinSysKey	m_linSysKey
	Key associated with this linear system. More...

const boost::weak_ptr< ExpList >	m_expList
	Local Matrix System. More...

const std::map< int, RobinBCInfoSharedPtr >	m_robinBCInfo
	Robin boundary info. More...

bool	m_verbose

Private Member Functions
virtual void	v_Solve (const Array< OneD, const NekDouble > &in, Array< OneD, NekDouble > &out, const AssemblyMapSharedPtr &locToGloMap, const Array< OneD, const NekDouble > &dirForcing=NullNekDouble1DArray)=0
	Solve a linear system based on mapping. More...

virtual void	v_SolveLinearSystem (const int pNumRows, const Array< OneD, const NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const AssemblyMapSharedPtr &locToGloMap, const int pNumDir)=0
	Solve a basic matrix system. More...

virtual void	v_InitObject ()

virtual void	v_Initialise (const boost::shared_ptr< AssemblyMap > &pLocToGloMap)

Static Private Attributes
static std::string	lookupIds []

static std::string	def

Detailed Description

A global linear system.

Consider the linear system $\boldsymbol{M\hat{u}}_g=\boldsymbol{\hat{f}}$ . Distinguishing between the boundary and interior components of $\boldsymbol{\hat{u}}_g$ and $\boldsymbol{\hat{f}}$ using $\boldsymbol{\hat{u}}_b$ , $\boldsymbol{\hat{u}}_i$ and $\boldsymbol{\hat{f}}_b$ , $\boldsymbol{\hat{f}}_i$ respectively, this system can be split into its constituent parts as

$\left[\begin{array}{cc} \boldsymbol{M}_b&\boldsymbol{M}_{c1}\\ \boldsymbol{M}_{c2}&\boldsymbol{M}_i\\ \end{array}\right] \left[\begin{array}{c} \boldsymbol{\hat{u}_b}\\ \boldsymbol{\hat{u}_i}\\ \end{array}\right]= \left[\begin{array}{c} \boldsymbol{\hat{f}_b}\\ \boldsymbol{\hat{f}_i}\\ \end{array}\right]$

where $\boldsymbol{M}_b$ represents the components of $\boldsymbol{M}$ resulting from boundary-boundary mode interactions, $\boldsymbol{M}_{c1}$ and $\boldsymbol{M}_{c2}$ represent the components resulting from coupling between the boundary-interior modes, and $\boldsymbol{M}_i$ represents the components of $\boldsymbol{M}$ resulting from interior-interior mode interactions.

The solution of the linear system can now be determined in two steps:

$\begin{eqnarray*} \mathrm{step 1:}&\quad&(\boldsymbol{M}_b-\boldsymbol{M}_{c1} \boldsymbol{M}_i^{-1}\boldsymbol{M}_{c2}) \boldsymbol{\hat{u}_b} = \boldsymbol{\hat{f}}_b - \boldsymbol{M}_{c1}\boldsymbol{M}_i^{-1} \boldsymbol{\hat{f}}_i,\nonumber \\ \mathrm{step 2:}&\quad&\boldsymbol{\hat{u}_i}=\boldsymbol{M}_i^{-1} \left( \boldsymbol{\hat{f}}_i - \boldsymbol{M}_{c2}\boldsymbol{\hat{u}_b} \right). \nonumber \\ \end{eqnarray*}$

As the inverse of $\boldsymbol{M}_i^{-1}$ is

$\boldsymbol{M}_i^{-1} = \left [\underline{\boldsymbol{M}^e_i} \right ]^{-1} = \underline{[\boldsymbol{M}^e_i]}^{-1}$

and the following operations can be evaluated as,

$\begin{eqnarray*} \boldsymbol{M}_{c1}\boldsymbol{M}_i^{-1}\boldsymbol{\hat{f}}_i & =& \boldsymbol{\mathcal{A}}_b^T \underline{\boldsymbol{M}^e_{c1}} \underline{[\boldsymbol{M}^e_i]}^{-1} \boldsymbol{\hat{f}}_i \\ \boldsymbol{M}_{c2} \boldsymbol{\hat{u}_b} &=& \underline{\boldsymbol{M}^e_{c2}} \boldsymbol{\mathcal{A}}_b \boldsymbol{\hat{u}_b}.\end{eqnarray*}$

where $\boldsymbol{\mathcal{A}}_b$ is the permutation matrix which scatters from global to local degrees of freedom, only the following four matrices should be constructed:

$\underline{[\boldsymbol{M}^e_i]}^{-1}$
$\underline{\boldsymbol{M}^e_{c1}} \underline{[\boldsymbol{M}^e_i]}^{-1}$
$\underline{\boldsymbol{M}^e_{c2}}$
The Schur complement: $\boldsymbol{M}_{\mathrm{Schur}}= \quad\boldsymbol{M}_b-\boldsymbol{M}_{c1}\boldsymbol{M}_i^{-1} \boldsymbol{M}_{c2}$

The first three matrices are just a concatenation of the corresponding local matrices and they can be created as such. They also allow for an elemental evaluation of the operations concerned.

The global Schur complement however should be assembled from the concatenation of the local elemental Schur complements, that is,

$\boldsymbol{M}_{\mathrm{Schur}}=\boldsymbol{M}_b - \boldsymbol{M}_{c1} \boldsymbol{M}_i^{-1} \boldsymbol{M}_{c2} = \boldsymbol{\mathcal{A}}_b^T \left [\underline{\boldsymbol{M}^e_b - \boldsymbol{M}^e_{c1} [\boldsymbol{M}^e_i]^{-1} (\boldsymbol{M}^e_{c2})} \right ] \boldsymbol{\mathcal{A}}_b$

and it is the only matrix operation that need to be evaluated on a global level when using static condensation. However, due to the size and sparsity of the matrix $\boldsymbol{\mathcal{A}}_b$ , it is more efficient to assemble the global Schur matrix using the mapping array bmap $[e][i]$ contained in the input argument locToGloMap. The global Schur complement is then constructed as:

$\boldsymbol{M}_{\mathrm{Schur}}\left[\mathrm{\a bmap}[e][i]\right] \left[\mathrm{\a bmap}[e][j]\right]=\mathrm{\a bsign}[e][i]\cdot \mathrm{\a bsign}[e][j] \cdot\boldsymbol{M}^e_{\mathrm{Schur}}[i][j]$

All four matrices are stored in the GlobalLinSys returned by this function.

Definition at line 74 of file GlobalLinSys.h.

Constructor & Destructor Documentation

Nektar::MultiRegions::GlobalLinSys::GlobalLinSys	(	const GlobalLinSysKey &	pKey,
		const boost::weak_ptr< ExpList > &	pExpList,
		const boost::shared_ptr< AssemblyMap > &	pLocToGloMap
	)

Constructor for full direct matrix solve.

Given a block matrix, construct a global matrix system according to a local to global mapping. #m_linSys is constructed by AssembleFullMatrix().

Parameters

pkey	Associated linear system key.
locToGloMap	Local to global mapping.

Definition at line 189 of file GlobalLinSys.cpp.

                                                  :
             m_linSysKey(pKey),
             m_expList(pExpList),
             m_robinBCInfo(m_expList.lock()->GetRobinBCInfo()),
             m_verbose(m_expList.lock()->GetSession()->
                       DefinesCmdLineArgument("verbose"))
         {
         }

virtual Nektar::MultiRegions::GlobalLinSys::~GlobalLinSys ( )

inlinevirtual

Definition at line 84 of file GlobalLinSys.h.

84 {}

Member Function Documentation

PreconditionerSharedPtr Nektar::MultiRegions::GlobalLinSys::CreatePrecon ( AssemblyMapSharedPtr asmMap )

protected

Create a preconditioner object from the parameters defined in the supplied assembly map.

Parameters

asmMap Assembly map used to construct the global system.

Definition at line 219 of file GlobalLinSys.cpp.

References Nektar::LibUtilities::NekFactory< tKey, tBase, >::CreateInstance(), Nektar::MultiRegions::GetPreconFactory(), GetSharedThisPtr(), and Nektar::MultiRegions::PreconditionerTypeMap.

         {
             PreconditionerType pType = asmMap->GetPreconType();
             std::string PreconType = MultiRegions::PreconditionerTypeMap[pType];
             return GetPreconFactory().CreateInstance(
                 PreconType, GetSharedThisPtr(), asmMap);
         }

void Nektar::MultiRegions::GlobalLinSys::DropStaticCondBlock ( unsigned int n )

inline

Definition at line 232 of file GlobalLinSys.h.

References v_DropStaticCondBlock().

         {
             return v_DropStaticCondBlock(n);
         }

DNekScalMatSharedPtr Nektar::MultiRegions::GlobalLinSys::GetBlock ( unsigned int n )

inline

Definition at line 222 of file GlobalLinSys.h.

References v_GetBlock().

Referenced by Nektar::MultiRegions::GlobalLinSysDirectFull::AssembleFullMatrix(), Nektar::MultiRegions::GlobalLinSysXxtFull::AssembleMatrixArrays(), and Nektar::MultiRegions::GlobalLinSysPETScFull::GlobalLinSysPETScFull().

         {
             return v_GetBlock(n);
         }

const GlobalLinSysKey & Nektar::MultiRegions::GlobalLinSys::GetKey ( void ) const

inline

Returns the key associated with the system.

Definition at line 171 of file GlobalLinSys.h.

References m_linSysKey.

         {
             return m_linSysKey;
         }

const boost::weak_ptr< ExpList > & Nektar::MultiRegions::GlobalLinSys::GetLocMat ( void ) const

inline

Definition at line 179 of file GlobalLinSys.h.

References m_expList.

         {
             return m_expList;
         }

int Nektar::MultiRegions::GlobalLinSys::GetNumBlocks ( )

inline

Definition at line 237 of file GlobalLinSys.h.

References v_GetNumBlocks().

         {
             return v_GetNumBlocks();
         }

boost::shared_ptr<GlobalLinSys> Nektar::MultiRegions::GlobalLinSys::GetSharedThisPtr ( )

inline

Returns a shared pointer to the current object.

Definition at line 107 of file GlobalLinSys.h.

Referenced by CreatePrecon().

             {
                 return shared_from_this();
             }

DNekScalBlkMatSharedPtr Nektar::MultiRegions::GlobalLinSys::GetStaticCondBlock ( unsigned int n )

inline

Definition at line 227 of file GlobalLinSys.h.

References v_GetStaticCondBlock().

         {
             return v_GetStaticCondBlock(n);
         }

void Nektar::MultiRegions::GlobalLinSys::Initialise ( const boost::shared_ptr< AssemblyMap > & pLocToGloMap )

inline

Definition at line 216 of file GlobalLinSys.h.

References v_Initialise().

Referenced by Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_InitObject(), Nektar::MultiRegions::GlobalLinSysStaticCond::v_InitObject(), and Nektar::MultiRegions::GlobalLinSysIterativeStaticCond::v_InitObject().

         {
             v_Initialise(pLocToGloMap);
         }

void Nektar::MultiRegions::GlobalLinSys::InitObject ( )

inline

Definition at line 211 of file GlobalLinSys.h.

References v_InitObject().

         {
             v_InitObject();
         }

void Nektar::MultiRegions::GlobalLinSys::Solve	(	const Array< OneD, const NekDouble > &	in,
		Array< OneD, NekDouble > &	out,
		const AssemblyMapSharedPtr &	locToGloMap,
		const Array< OneD, const NekDouble > &	dirForcing = `NullNekDouble1DArray`
	)

inline

Solve the linear system for given input and output vectors using a specified local to global map.

Definition at line 188 of file GlobalLinSys.h.

References v_Solve().

         {
             v_Solve(in,out,locToGloMap,dirForcing);
         }

void Nektar::MultiRegions::GlobalLinSys::SolveLinearSystem	(	const int	pNumRows,
		const Array< OneD, const NekDouble > &	pInput,
		Array< OneD, NekDouble > &	pOutput,
		const AssemblyMapSharedPtr &	locToGloMap,
		const int	pNumDir = `0`
	)

inline

Solve the linear system for given input and output vectors.

Definition at line 201 of file GlobalLinSys.h.

References v_SolveLinearSystem().

Referenced by Nektar::MultiRegions::GlobalLinSysDirectFull::v_Solve(), Nektar::MultiRegions::GlobalLinSysXxtFull::v_Solve(), Nektar::MultiRegions::GlobalLinSysIterativeFull::v_Solve(), Nektar::MultiRegions::GlobalLinSysPETScFull::v_Solve(), and Nektar::MultiRegions::GlobalLinSysStaticCond::v_Solve().

         {
             v_SolveLinearSystem(pNumRows, pInput, pOutput, locToGloMap, pNumDir);
         }

void Nektar::MultiRegions::GlobalLinSys::v_DropStaticCondBlock ( unsigned int n )

protectedvirtual

Releases the static condensation block matrices from NekManager of n-th expansion using the matrix key provided by the m_linSysKey.

Parameters

n	Number of the expansion

Definition at line 396 of file GlobalLinSys.cpp.

References Nektar::MultiRegions::GlobalMatrixKey::GetConstFactors(), Nektar::MultiRegions::GlobalMatrixKey::GetMatrixType(), Nektar::MultiRegions::GlobalMatrixKey::GetNVarCoeffs(), Nektar::MultiRegions::GlobalMatrixKey::GetVarCoeffs(), m_expList, and m_linSysKey.

Referenced by DropStaticCondBlock(), and Nektar::MultiRegions::GlobalLinSysIterativeStaticCond::PrepareLocalSchurComplement().

         {
             boost::shared_ptr<MultiRegions::ExpList> expList = m_expList.lock();
 
             StdRegions::StdExpansionSharedPtr vExp = expList->GetExp(n);
 
             // need to be initialised with zero size for non variable
             // coefficient case
             StdRegions::VarCoeffMap vVarCoeffMap;
 
             // retrieve variable coefficients
             if(m_linSysKey.GetNVarCoeffs() > 0)
             {
                 StdRegions::VarCoeffMap::const_iterator x;
                 int cnt = expList->GetPhys_Offset(n);
                 for (x  = m_linSysKey.GetVarCoeffs().begin(); 
                      x != m_linSysKey.GetVarCoeffs().end  (); ++x)
                 {
                     vVarCoeffMap[x->first] = x->second + cnt;
                 }
             }
 
             LocalRegions::MatrixKey matkey(m_linSysKey.GetMatrixType(),
                                            vExp->DetShapeType(),
                                            *vExp,
                                            m_linSysKey.GetConstFactors(),
                                            vVarCoeffMap);
 
             vExp->DropLocStaticCondMatrix(matkey);
         }

DNekScalMatSharedPtr Nektar::MultiRegions::GlobalLinSys::v_GetBlock ( unsigned int n )

protectedvirtual

Retrieves the block matrix from n-th expansion using the matrix key provided by the m_linSysKey.

Parameters

n	Number of the expansion.

Returns: Block matrix for the specified expansion.

Definition at line 246 of file GlobalLinSys.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::GlobalMatrixKey::GetConstFactors(), Nektar::MultiRegions::GlobalMatrixKey::GetMatrixType(), Nektar::MultiRegions::GlobalMatrixKey::GetNVarCoeffs(), Nektar::MultiRegions::GlobalMatrixKey::GetVarCoeffs(), m_expList, m_linSysKey, and m_robinBCInfo.

Referenced by GetBlock(), and Nektar::MultiRegions::GlobalLinSysStaticCond::SetupTopLevel().

         {
             boost::shared_ptr<MultiRegions::ExpList> expList = m_expList.lock();
             int cnt = 0;
             DNekScalMatSharedPtr loc_mat;
 
             LocalRegions::ExpansionSharedPtr vExp = 
                 boost::dynamic_pointer_cast<LocalRegions::Expansion>(
                     expList->GetExp(n));
 
             // need to be initialised with zero size for non variable
             // coefficient case
             StdRegions::VarCoeffMap vVarCoeffMap;
 
             // retrieve variable coefficients
             if(m_linSysKey.GetNVarCoeffs() > 0)
             {
                 StdRegions::VarCoeffMap::const_iterator x;
                 cnt = expList->GetPhys_Offset(n);
                 
                 for (x = m_linSysKey.GetVarCoeffs().begin(); 
                      x != m_linSysKey.GetVarCoeffs().end(); ++x)
                 {
                     vVarCoeffMap[x->first] = x->second + cnt;
                 }
             }
 
             LocalRegions::MatrixKey matkey(m_linSysKey.GetMatrixType(),
                                            vExp->DetShapeType(),
                                            *vExp, m_linSysKey.GetConstFactors(),
                                            vVarCoeffMap);
             loc_mat = vExp->GetLocMatrix(matkey);
 
             // apply robin boundary conditions to the matrix.
             if(m_robinBCInfo.count(n) != 0) // add robin mass matrix
             {
                 RobinBCInfoSharedPtr rBC;
 
                 // declare local matrix from scaled matrix.
                 int rows = loc_mat->GetRows();
                 int cols = loc_mat->GetColumns();
                 const NekDouble *dat = loc_mat->GetRawPtr();
                 DNekMatSharedPtr new_mat = MemoryManager<DNekMat>::
                     AllocateSharedPtr(rows,cols,dat);
                 Blas::Dscal(rows*cols,loc_mat->Scale(),new_mat->GetRawPtr(),1);
 
                 // add local matrix contribution
                 for(rBC = m_robinBCInfo.find(n)->second;rBC; rBC = rBC->next)
                 {
                     vExp->AddRobinMassMatrix(
                         rBC->m_robinID, rBC->m_robinPrimitiveCoeffs, new_mat);
                 }
 
                 // redeclare loc_mat to point to new_mat plus the scalar.
                 loc_mat = MemoryManager<DNekScalMat>::AllocateSharedPtr(
                     1.0, new_mat);
             }
 
             // finally return the matrix.
             return loc_mat;
         }

int Nektar::MultiRegions::GlobalLinSys::v_GetNumBlocks ( )

protectedvirtual

Get the number of blocks in this system.

At the top level this corresponds to the number of elements in the expansion list.

Reimplemented in Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 234 of file GlobalLinSys.cpp.

References m_expList.

Referenced by GetNumBlocks().

         {
             return m_expList.lock()->GetExpSize();
         }

DNekScalBlkMatSharedPtr Nektar::MultiRegions::GlobalLinSys::v_GetStaticCondBlock ( unsigned int n )

protectedvirtual

Retrieves a the static condensation block matrices from n-th expansion using the matrix key provided by the m_linSysKey.

Parameters

n	Number of the expansion

Returns: 2x2 Block matrix holding the static condensation matrices for the n-th expansion.

Reimplemented in Nektar::MultiRegions::GlobalLinSysIterativeStaticCond, and Nektar::MultiRegions::GlobalLinSysPETScStaticCond.

Definition at line 316 of file GlobalLinSys.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::GlobalMatrixKey::GetConstFactors(), Nektar::MultiRegions::GlobalMatrixKey::GetMatrixType(), Nektar::MultiRegions::GlobalMatrixKey::GetNVarCoeffs(), Nektar::MultiRegions::GlobalMatrixKey::GetVarCoeffs(), m_expList, m_linSysKey, and m_robinBCInfo.

Referenced by GetStaticCondBlock(), and Nektar::MultiRegions::GlobalLinSysStaticCond::SetupTopLevel().

         {
             boost::shared_ptr<MultiRegions::ExpList> expList = m_expList.lock();
             int cnt = 0;
             DNekScalBlkMatSharedPtr loc_mat;
             DNekScalMatSharedPtr    tmp_mat;
 
             StdRegions::StdExpansionSharedPtr vExp = expList->GetExp(n);
 
             // need to be initialised with zero size for non variable
             // coefficient case
             StdRegions::VarCoeffMap vVarCoeffMap;
 
             // retrieve variable coefficients
             if(m_linSysKey.GetNVarCoeffs() > 0)
             {
                 StdRegions::VarCoeffMap::const_iterator x;
                 cnt = expList->GetPhys_Offset(n);
                 for (x  = m_linSysKey.GetVarCoeffs().begin(); 
                      x != m_linSysKey.GetVarCoeffs().end  (); ++x)
                 {
                     vVarCoeffMap[x->first] = x->second + cnt;
                 }
             }
 
             LocalRegions::MatrixKey matkey(m_linSysKey.GetMatrixType(),
                                            vExp->DetShapeType(),
                                            *vExp,
                                            m_linSysKey.GetConstFactors(),
                                            vVarCoeffMap);
 
             loc_mat = vExp->GetLocStaticCondMatrix(matkey);
 
             if(m_robinBCInfo.count(n) != 0) // add robin mass matrix
             {
                 RobinBCInfoSharedPtr rBC;
 
                 tmp_mat = loc_mat->GetBlock(0,0);
 
                 // declare local matrix from scaled matrix.
                 int rows = tmp_mat->GetRows();
                 int cols = tmp_mat->GetColumns();
                 const NekDouble *dat = tmp_mat->GetRawPtr();
                 DNekMatSharedPtr new_mat = MemoryManager<DNekMat>::
                     AllocateSharedPtr(rows, cols, dat);
                 Blas::Dscal(rows*cols,tmp_mat->Scale(),new_mat->GetRawPtr(),1);
 
                 // add local matrix contribution
                 for(rBC = m_robinBCInfo.find(n)->second;rBC; rBC = rBC->next)
                 {
                     vExp->AddRobinMassMatrix(
                         rBC->m_robinID, rBC->m_robinPrimitiveCoeffs, new_mat);
                 }
 
                 // redeclare loc_mat to point to new_mat plus the scalar.
                 tmp_mat = MemoryManager<DNekScalMat>::AllocateSharedPtr(
                     1.0, new_mat);
                 DNekScalBlkMatSharedPtr new_loc_mat;
                 unsigned int exp_size[] = {tmp_mat->GetRows(), loc_mat->GetBlock(1,1)->GetRows()};
                 unsigned int nblks = 2;
                 new_loc_mat = MemoryManager<DNekScalBlkMat>::AllocateSharedPtr(nblks, nblks, exp_size, exp_size);
 
 
                 new_loc_mat->SetBlock(0,0,tmp_mat);
                 new_loc_mat->SetBlock(0,1,loc_mat->GetBlock(0,1));
                 new_loc_mat->SetBlock(1,0,loc_mat->GetBlock(1,0));
                 new_loc_mat->SetBlock(1,1,loc_mat->GetBlock(1,1));
                 loc_mat = new_loc_mat;
             }
 
             return loc_mat;
         }

void Nektar::MultiRegions::GlobalLinSys::v_Initialise ( const boost::shared_ptr< AssemblyMap > & pLocToGloMap )

privatevirtual

Reimplemented in Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 432 of file GlobalLinSys.cpp.

References ErrorUtil::efatal, and NEKERROR.

Referenced by Initialise().

         {
             NEKERROR(ErrorUtil::efatal, "Method does not exist" );
     }

void Nektar::MultiRegions::GlobalLinSys::v_InitObject ( )

privatevirtual

Reimplemented in Nektar::MultiRegions::GlobalLinSysIterativeStaticCond, Nektar::MultiRegions::GlobalLinSysStaticCond, and Nektar::MultiRegions::GlobalLinSysPETScStaticCond.

Definition at line 427 of file GlobalLinSys.cpp.

References ErrorUtil::efatal, and NEKERROR.

Referenced by InitObject().

         {
             NEKERROR(ErrorUtil::efatal, "Method does not exist" );
     }

virtual void Nektar::MultiRegions::GlobalLinSys::v_Solve	(	const Array< OneD, const NekDouble > &	in,
		Array< OneD, NekDouble > &	out,
		const AssemblyMapSharedPtr &	locToGloMap,
		const Array< OneD, const NekDouble > &	dirForcing = `NullNekDouble1DArray`
	)

privatepure virtual

Solve a linear system based on mapping.

Implemented in Nektar::MultiRegions::GlobalLinSysStaticCond, Nektar::MultiRegions::GlobalLinSysPETScFull, Nektar::MultiRegions::GlobalLinSysIterativeFull, Nektar::MultiRegions::GlobalLinSysXxtFull, and Nektar::MultiRegions::GlobalLinSysDirectFull.

Referenced by Solve().

virtual void Nektar::MultiRegions::GlobalLinSys::v_SolveLinearSystem	(	const int	pNumRows,
		const Array< OneD, const NekDouble > &	pInput,
		Array< OneD, NekDouble > &	pOutput,
		const AssemblyMapSharedPtr &	locToGloMap,
		const int	pNumDir
	)

privatepure virtual

Solve a basic matrix system.

Implemented in Nektar::MultiRegions::GlobalLinSysIterative, Nektar::MultiRegions::GlobalLinSysPETSc, Nektar::MultiRegions::GlobalLinSysDirect, and Nektar::MultiRegions::GlobalLinSysXxt.

Referenced by SolveLinearSystem().

Member Data Documentation

std::string Nektar::MultiRegions::GlobalLinSys::def

staticprivate

Initial value:

= LibUtilities::SessionReader::
            RegisterDefaultSolverInfo("GlobalSysSoln",
                                      "DirectMultiLevelStaticCond")

Definition at line 164 of file GlobalLinSys.h.

std::string Nektar::MultiRegions::GlobalLinSys::lookupIds

staticprivate

Definition at line 163 of file GlobalLinSys.h.

const boost::weak_ptr<ExpList> Nektar::MultiRegions::GlobalLinSys::m_expList

protected

Local Matrix System.

Definition at line 129 of file GlobalLinSys.h.

const GlobalLinSysKey Nektar::MultiRegions::GlobalLinSys::m_linSysKey

protected

Key associated with this linear system.

Definition at line 127 of file GlobalLinSys.h.

const std::map<int, RobinBCInfoSharedPtr> Nektar::MultiRegions::GlobalLinSys::m_robinBCInfo

protected

Robin boundary info.

Definition at line 131 of file GlobalLinSys.h.

Referenced by v_GetBlock(), and v_GetStaticCondBlock().

bool Nektar::MultiRegions::GlobalLinSys::m_verbose

protected

Definition at line 133 of file GlobalLinSys.h.

Referenced by Nektar::MultiRegions::GlobalLinSysXxtFull::AssembleMatrixArrays(), Nektar::MultiRegions::GlobalLinSysIterative::DoConjugateGradient(), and Nektar::MultiRegions::GlobalLinSysXxtStaticCond::v_AssembleSchurComplement().

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation