A global linear system. More...

#include <GlobalLinSysPETScStaticCond.h>

Inheritance diagram for Nektar::MultiRegions::GlobalLinSysPETScStaticCond:

Collaboration diagram for Nektar::MultiRegions::GlobalLinSysPETScStaticCond:

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

	GlobalLinSysPETScStaticCond (const GlobalLinSysKey &mkey, const boost::weak_ptr< ExpList > &pExpList, const DNekScalBlkMatSharedPtr pSchurCompl, const DNekScalBlkMatSharedPtr pBinvD, const DNekScalBlkMatSharedPtr pC, const DNekScalBlkMatSharedPtr pInvD, const boost::shared_ptr< AssemblyMap > &locToGloMap, const PreconditionerSharedPtr pPrecon=PreconditionerSharedPtr())
	Constructor for full direct matrix solve. More...

virtual	~GlobalLinSysPETScStaticCond ()

Public Member Functions inherited from Nektar::MultiRegions::GlobalLinSysPETSc
	GlobalLinSysPETSc (const GlobalLinSysKey &pKey, const boost::weak_ptr< ExpList > &pExp, const boost::shared_ptr< AssemblyMap > &pLocToGloMap)
	Constructor for full direct matrix solve. More...

virtual	~GlobalLinSysPETSc ()
	Clean up PETSc objects. More...

virtual void	v_SolveLinearSystem (const int pNumRows, const Array< OneD, const NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const AssemblyMapSharedPtr &locToGloMap, const int pNumDir)
	Solve linear system using PETSc. More...

Public Member Functions inherited from 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. 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...

Public Member Functions inherited from Nektar::MultiRegions::GlobalLinSysStaticCond
	GlobalLinSysStaticCond (const GlobalLinSysKey &mkey, const boost::weak_ptr< ExpList > &pExpList, const boost::shared_ptr< AssemblyMap > &locToGloMap)
	Constructor for full direct matrix solve. More...

virtual	~GlobalLinSysStaticCond ()

Static Public Member Functions
static GlobalLinSysSharedPtr	create (const GlobalLinSysKey &pLinSysKey, const boost::weak_ptr< ExpList > &pExpList, const boost::shared_ptr< AssemblyMap > &pLocToGloMap)
	Creates an instance of this class. More...

Static Public Attributes
static std::string	className
	Name of class. More...

static std::string	className2

Protected Member Functions
virtual void	v_InitObject ()

virtual void	v_AssembleSchurComplement (boost::shared_ptr< AssemblyMap > locToGloMap)
	Assemble the Schur complement matrix. More...

virtual void	v_DoMatrixMultiply (const Array< OneD, const NekDouble > &input, Array< OneD, NekDouble > &output)
	Apply matrix-vector multiplication using local approach and the assembly map. 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 DNekScalBlkMatSharedPtr	v_PreSolve (int scLevel, NekVector< NekDouble > &F_GlobBnd)

virtual void	v_BasisTransform (Array< OneD, NekDouble > &pInOut, int offset)

virtual void	v_BasisInvTransform (Array< OneD, NekDouble > &pInOut)

virtual GlobalLinSysStaticCondSharedPtr	v_Recurse (const GlobalLinSysKey &mkey, const boost::weak_ptr< ExpList > &pExpList, const DNekScalBlkMatSharedPtr pSchurCompl, const DNekScalBlkMatSharedPtr pBinvD, const DNekScalBlkMatSharedPtr pC, const DNekScalBlkMatSharedPtr pInvD, const boost::shared_ptr< AssemblyMap > &locToGloMap)

Protected Member Functions inherited from Nektar::MultiRegions::GlobalLinSysPETSc
void	SetUpScatter ()
	Set up PETSc local (equivalent to Nektar++ global) and global (equivalent to universal) scatter maps. More...

void	SetUpMatVec (int nGlobal, int nDir)
	Construct PETSc matrix and vector handles. More...

void	SetUpSolver (NekDouble tolerance)
	Set up KSP solver object. More...

void	CalculateReordering (const Array< OneD, const int > &glo2uniMap, const Array< OneD, const int > &glo2unique, const AssemblyMapSharedPtr &pLocToGloMap)
	Calculate a reordering of universal IDs for PETSc. More...

Protected Member Functions inherited from Nektar::MultiRegions::GlobalLinSys
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 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 Member Functions inherited from Nektar::MultiRegions::GlobalLinSysStaticCond
virtual int	v_GetNumBlocks ()
	Get the number of blocks in this system. More...

virtual void	v_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...

virtual void	v_Initialise (const boost::shared_ptr< AssemblyMap > &locToGloMap)
	Initialise this object. More...

void	SetupTopLevel (const boost::shared_ptr< AssemblyMap > &locToGloMap)
	Set up the storage for the Schur complement or the top level of the multi-level Schur complement. More...

void	ConstructNextLevelCondensedSystem (const boost::shared_ptr< AssemblyMap > &locToGloMap)

Protected Attributes
DNekScalBlkMatSharedPtr	m_S1Blk

Protected Attributes inherited from Nektar::MultiRegions::GlobalLinSysPETSc
Mat	m_matrix
	PETSc matrix object. More...

Vec	m_x
	PETSc vector objects used for local storage. More...

Vec	m_b

Vec	m_locVec

KSP	m_ksp
	KSP object that represents solver system. More...

PC	m_pc
	PCShell for preconditioner. More...

PETScMatMult	m_matMult
	Enumerator to select matrix multiplication type. More...

std::vector< int >	m_reorderedMap
	Reordering that takes universal IDs to a unique row in the PETSc matrix. More...

VecScatter	m_ctx
	PETSc scatter context that takes us between Nektar++ global ordering and PETSc vector ordering. More...

int	m_nLocal
	Number of unique degrees of freedom on this process. More...

PreconditionerSharedPtr	m_precon

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

Protected Attributes inherited from Nektar::MultiRegions::GlobalLinSysStaticCond
GlobalLinSysStaticCondSharedPtr	m_recursiveSchurCompl
	Schur complement for Direct Static Condensation. More...

DNekScalBlkMatSharedPtr	m_schurCompl
	Block Schur complement matrix. More...

DNekScalBlkMatSharedPtr	m_BinvD
	Block $BD^{-1}$ matrix. More...

DNekScalBlkMatSharedPtr	m_C
	Block matrix. More...

DNekScalBlkMatSharedPtr	m_invD
	Block $D^{-1}$ matrix. More...

boost::shared_ptr< AssemblyMap >	m_locToGloMap
	Local to global map. More...

Array< OneD, NekDouble >	m_wsp
	Workspace array for matrix multiplication. More...

Detailed Description

A global linear system.

Definition at line 55 of file GlobalLinSysPETScStaticCond.h.

Constructor & Destructor Documentation

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

Constructor for full direct matrix solve.

For a matrix system of the form

$\left[ \begin{array}{cc} \boldsymbol{A} & \boldsymbol{B}\\ \boldsymbol{C} & \boldsymbol{D} \end{array} \right] \left[ \begin{array}{c} \boldsymbol{x_1}\\ \boldsymbol{x_2} \end{array}\right] = \left[ \begin{array}{c} \boldsymbol{y_1}\\ \boldsymbol{y_2} \end{array}\right],$

where $\boldsymbol{D}$ and $(\boldsymbol{A-BD^{-1}C})$ are invertible, store and assemble a static condensation system, according to a given local to global mapping. #m_linSys is constructed by AssembleSchurComplement().

Parameters

mKey	Associated matrix key.
pLocMatSys	LocalMatrixSystem
locToGloMap	Local to global mapping.

Definition at line 89 of file GlobalLinSysPETScStaticCond.cpp.

References ASSERTL1, Nektar::MultiRegions::ePETScMultiLevelStaticCond, Nektar::MultiRegions::ePETScStaticCond, and Nektar::MultiRegions::GlobalLinSysKey::GetGlobalSysSolnType().

             : GlobalLinSys          (pKey, pExpList, pLocToGloMap),
               GlobalLinSysPETSc     (pKey, pExpList, pLocToGloMap),
               GlobalLinSysStaticCond(pKey, pExpList, pLocToGloMap)
         {
             ASSERTL1((pKey.GetGlobalSysSolnType()==ePETScStaticCond)||
                      (pKey.GetGlobalSysSolnType()==ePETScMultiLevelStaticCond),
                      "This constructor is only valid when using static "
                      "condensation");
             ASSERTL1(pKey.GetGlobalSysSolnType()
                         == pLocToGloMap->GetGlobalSysSolnType(),
                      "The local to global map is not set up for the requested "
                      "solution type");
         }

Nektar::MultiRegions::GlobalLinSysPETScStaticCond::GlobalLinSysPETScStaticCond	(	const GlobalLinSysKey &	mkey,
		const boost::weak_ptr< ExpList > &	pExpList,
		const DNekScalBlkMatSharedPtr	pSchurCompl,
		const DNekScalBlkMatSharedPtr	pBinvD,
		const DNekScalBlkMatSharedPtr	pC,
		const DNekScalBlkMatSharedPtr	pInvD,
		const boost::shared_ptr< AssemblyMap > &	locToGloMap,
		const PreconditionerSharedPtr	pPrecon = `PreconditionerSharedPtr()`
	)

Constructor for full direct matrix solve.

Definition at line 110 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::MultiRegions::GlobalLinSysStaticCond::m_BinvD, Nektar::MultiRegions::GlobalLinSysStaticCond::m_C, Nektar::MultiRegions::GlobalLinSysStaticCond::m_invD, Nektar::MultiRegions::GlobalLinSysPETSc::m_precon, and Nektar::MultiRegions::GlobalLinSysStaticCond::m_schurCompl.

             : GlobalLinSys          (pKey, pExpList, pLocToGloMap),
               GlobalLinSysPETSc     (pKey, pExpList, pLocToGloMap),
               GlobalLinSysStaticCond(pKey, pExpList, pLocToGloMap)
         {
             m_schurCompl = pSchurCompl;
             m_BinvD      = pBinvD;
             m_C          = pC;
             m_invD       = pInvD;
             m_precon     = pPrecon;
         }

Nektar::MultiRegions::GlobalLinSysPETScStaticCond::~GlobalLinSysPETScStaticCond ( )

virtual

Definition at line 133 of file GlobalLinSysPETScStaticCond.cpp.

         {
 
         }

Member Function Documentation

static GlobalLinSysSharedPtr Nektar::MultiRegions::GlobalLinSysPETScStaticCond::create	(	const GlobalLinSysKey &	pLinSysKey,
		const boost::weak_ptr< ExpList > &	pExpList,
		const boost::shared_ptr< AssemblyMap > &	pLocToGloMap
	)

inlinestatic

Creates an instance of this class.

Definition at line 60 of file GlobalLinSysPETScStaticCond.h.

             {
                 GlobalLinSysSharedPtr p = MemoryManager<
                     GlobalLinSysPETScStaticCond>::AllocateSharedPtr(
                         pLinSysKey, pExpList, pLocToGloMap);
                 p->InitObject();
                 return p;
             }

void Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_AssembleSchurComplement ( boost::shared_ptr< AssemblyMap > locToGloMap )

protectedvirtual

Assemble the Schur complement matrix.

Assemble the schur complement matrix from the block matrices stored in #m_blkMatrices and the given local to global mapping information.

Parameters

locToGloMap Local to global mapping information.

Reimplemented from Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 187 of file GlobalLinSysPETScStaticCond.cpp.

         {
             int i, j, n, cnt, gid1, gid2, loc_lda;
             NekDouble sign1, sign2, value;
 
             const int nDirDofs = pLocToGloMap->GetNumGlobalDirBndCoeffs();
 
             DNekScalBlkMatSharedPtr SchurCompl = m_schurCompl;
             DNekScalBlkMatSharedPtr BinvD      = m_BinvD;
             DNekScalBlkMatSharedPtr C          = m_C;
             DNekScalBlkMatSharedPtr invD       = m_invD;
             DNekScalMatSharedPtr    loc_mat;
 
             // Build precon again if we in multi-level static condensation (a
             // bit of a hack)
             if (m_linSysKey.GetGlobalSysSolnType() ==
                     ePETScMultiLevelStaticCond)
             {
                 m_precon = CreatePrecon(m_locToGloMap);
                 m_precon->BuildPreconditioner();
             }
 
             // CALCULATE REORDERING MAPPING
             CalculateReordering(pLocToGloMap->GetGlobalToUniversalBndMap(),
                                 pLocToGloMap->GetGlobalToUniversalBndMapUnique(),
                                 pLocToGloMap);
 
             // SET UP VECTORS AND MATRIX
             SetUpMatVec(pLocToGloMap->GetNumGlobalBndCoeffs(), nDirDofs);
 
             // SET UP SCATTER OBJECTS
             SetUpScatter();
 
             // CONSTRUCT KSP OBJECT
             SetUpSolver(pLocToGloMap->GetIterativeTolerance());
 
             // If we are using the matrix multiplication shell don't try to
             // populate the matrix.
             if (m_matMult == ePETScMatMultShell)
             {
                 return;
             }
 
             // POPULATE MATRIX
             for(n = cnt = 0; n < m_schurCompl->GetNumberOfBlockRows(); ++n)
             {
                 loc_mat = m_schurCompl->GetBlock(n,n);
                 loc_lda = loc_mat->GetRows();
 
                 for(i = 0; i < loc_lda; ++i)
                 {
                     gid1 = pLocToGloMap->GetLocalToGlobalBndMap(cnt + i)-nDirDofs;
                     sign1 = pLocToGloMap->GetLocalToGlobalBndSign(cnt + i);
                     if(gid1 >= 0)
                     {
                         int gid1ro = m_reorderedMap[gid1];
                         for(j = 0; j < loc_lda; ++j)
                         {
                             gid2 = pLocToGloMap->GetLocalToGlobalBndMap(cnt + j)
                                                                     - nDirDofs;
                             sign2 = pLocToGloMap->GetLocalToGlobalBndSign(cnt + j);
                             if(gid2 >= 0)
                             {
                                 int gid2ro = m_reorderedMap[gid2];
                                 value = sign1*sign2*(*loc_mat)(i,j);
                                 MatSetValue(
                                     m_matrix, gid1ro, gid2ro, value, ADD_VALUES);
                             }
                         }
                     }
                 }
                 cnt   += loc_lda;
             }
 
             // ASSEMBLE MATRIX
             MatAssemblyBegin(m_matrix, MAT_FINAL_ASSEMBLY);
             MatAssemblyEnd  (m_matrix, MAT_FINAL_ASSEMBLY);
         }

void Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_BasisInvTransform ( Array< OneD, NekDouble > & pInOut )

protectedvirtual

Reimplemented from Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 314 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::MultiRegions::GlobalLinSysPETSc::m_precon.

         {
             m_precon->DoTransformFromLowEnergy(pInOut);
         }

void Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_BasisTransform	(	Array< OneD, NekDouble > &	pInOut,
		int	offset
	)

protectedvirtual

Reimplemented from Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 307 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::MultiRegions::GlobalLinSysPETSc::m_precon.

         {
             m_precon->DoTransformToLowEnergy(pInOut, offset);
         }

void Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_DoMatrixMultiply	(	const Array< OneD, const NekDouble > &	input,
		Array< OneD, NekDouble > &	output
	)

protectedvirtual

Apply matrix-vector multiplication using local approach and the assembly map.

Parameters

input	Vector input.
output	Result of multiplication.

Todo:: This can possibly be made faster by using the sparse block-matrix multiplication code from the iterative elastic systems.

Implements Nektar::MultiRegions::GlobalLinSysPETSc.

Definition at line 331 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::NekVector< DataType >::GetPtr(), and Nektar::MultiRegions::GlobalLinSysStaticCond::m_locToGloMap.

         {
             int nLocBndDofs = m_locToGloMap->GetNumLocalBndCoeffs();
             int nDirDofs    = m_locToGloMap->GetNumGlobalDirBndCoeffs();
 
             NekVector<NekDouble> in(nLocBndDofs), out(nLocBndDofs);
             m_locToGloMap->GlobalToLocalBnd(input, in.GetPtr(), nDirDofs);
             out = (*m_schurCompl) * in;
             m_locToGloMap->AssembleBnd(out.GetPtr(), output, nDirDofs);
         }

DNekScalBlkMatSharedPtr Nektar::MultiRegions::GlobalLinSysPETScStaticCond::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 from Nektar::MultiRegions::GlobalLinSys.

Definition at line 268 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), Nektar::MultiRegions::GlobalLinSysStaticCond::m_locToGloMap, m_S1Blk, and Nektar::MultiRegions::GlobalLinSysStaticCond::m_schurCompl.

         {
             DNekScalBlkMatSharedPtr schurComplBlock;
             int  scLevel           = m_locToGloMap->GetStaticCondLevel();
             DNekScalBlkMatSharedPtr sc = scLevel == 0 ? m_S1Blk : m_schurCompl;
             DNekScalMatSharedPtr    localMat = sc->GetBlock(n,n);
             unsigned int nbdry    = localMat->GetRows();
             unsigned int nblks    = 1;
             unsigned int esize[1] = {nbdry};
 
             schurComplBlock = MemoryManager<DNekScalBlkMat>
                 ::AllocateSharedPtr(nblks, nblks, esize, esize);
             schurComplBlock->SetBlock(0, 0, localMat);
 
             return schurComplBlock;
         }

void Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_InitObject ( )

protectedvirtual

Reimplemented from Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 138 of file GlobalLinSysPETScStaticCond.cpp.

         {
             m_precon = CreatePrecon(m_locToGloMap);
 
             // Allocate memory for top-level structure
             SetupTopLevel(m_locToGloMap);
 
             // Setup Block Matrix systems
             int n, n_exp = m_expList.lock()->GetNumElmts();
 
             MatrixStorage blkmatStorage = eDIAGONAL;
             const Array<OneD,const unsigned int>& nbdry_size
                     = m_locToGloMap->GetNumLocalBndCoeffsPerPatch();
 
             m_S1Blk = MemoryManager<DNekScalBlkMat>
                 ::AllocateSharedPtr(nbdry_size, nbdry_size, blkmatStorage);
 
             // Preserve original matrix in m_S1Blk
             for (n = 0; n < n_exp; ++n)
             {
                 DNekScalMatSharedPtr mat = m_schurCompl->GetBlock(n, n);
                 m_S1Blk->SetBlock(n, n, mat);
             }
 
             // Build preconditioner
             m_precon->BuildPreconditioner();
 
             // Do transform of Schur complement matrix
             for (n = 0; n < n_exp; ++n)
             {
                 if (m_linSysKey.GetMatrixType() !=
                         StdRegions::eHybridDGHelmBndLam)
                 {
                     DNekScalMatSharedPtr mat = m_S1Blk->GetBlock(n, n);
                     DNekScalMatSharedPtr t = m_precon->TransformedSchurCompl(
                         m_expList.lock()->GetOffset_Elmt_Id(n), mat);
                     m_schurCompl->SetBlock(n, n, t);
                 }
             }
 
             // Construct this level
             Initialise(m_locToGloMap);
         }

DNekScalBlkMatSharedPtr Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_PreSolve	(	int	scLevel,
		NekVector< NekDouble > &	F_GlobBnd
	)

protectedvirtual

Reimplemented from Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 285 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::MultiRegions::GlobalLinSys::CreatePrecon(), Nektar::MultiRegions::GlobalLinSysStaticCond::m_locToGloMap, Nektar::MultiRegions::GlobalLinSysPETSc::m_precon, m_S1Blk, and Nektar::MultiRegions::GlobalLinSysStaticCond::m_schurCompl.

         {
             if (scLevel == 0)
             {
                 // When matrices are supplied to the constructor at the top
                 // level, the preconditioner is never set up.
                 if (!m_precon)
                 {
                     m_precon = CreatePrecon(m_locToGloMap);
                     m_precon->BuildPreconditioner();
                 }
 
                 return m_S1Blk;
             }
             else
             {
                 return m_schurCompl;
             }
         }

GlobalLinSysStaticCondSharedPtr Nektar::MultiRegions::GlobalLinSysPETScStaticCond::v_Recurse	(	const GlobalLinSysKey &	mkey,
		const boost::weak_ptr< ExpList > &	pExpList,
		const DNekScalBlkMatSharedPtr	pSchurCompl,
		const DNekScalBlkMatSharedPtr	pBinvD,
		const DNekScalBlkMatSharedPtr	pC,
		const DNekScalBlkMatSharedPtr	pInvD,
		const boost::shared_ptr< AssemblyMap > &	locToGloMap
	)

protectedvirtual

Implements Nektar::MultiRegions::GlobalLinSysStaticCond.

Definition at line 344 of file GlobalLinSysPETScStaticCond.cpp.

References Nektar::MultiRegions::GlobalLinSysPETSc::m_precon.

         {
             GlobalLinSysPETScStaticCondSharedPtr sys = MemoryManager<
                 GlobalLinSysPETScStaticCond>::AllocateSharedPtr(
                     mkey, pExpList, pSchurCompl, pBinvD, pC, pInvD, l2gMap,
                     m_precon);
             sys->Initialise(l2gMap);
             return sys;
         }

Member Data Documentation

string Nektar::MultiRegions::GlobalLinSysPETScStaticCond::className

static

Initial value:

= GetGlobalLinSysFactory().RegisterCreatorFunction(
                    "PETScStaticCond",
                    GlobalLinSysPETScStaticCond::create,
                    "PETSc static condensation.")

Name of class.

Registers the class with the Factory.

Definition at line 73 of file GlobalLinSysPETScStaticCond.h.

string Nektar::MultiRegions::GlobalLinSysPETScStaticCond::className2

static

Initial value:

= GetGlobalLinSysFactory().RegisterCreatorFunction(
                    "PETScMultiLevelStaticCond",
                    GlobalLinSysPETScStaticCond::create,
                    "PETSc multi-level static condensation.")

Definition at line 74 of file GlobalLinSysPETScStaticCond.h.

DNekScalBlkMatSharedPtr Nektar::MultiRegions::GlobalLinSysPETScStaticCond::m_S1Blk

protected

Definition at line 97 of file GlobalLinSysPETScStaticCond.h.

Referenced by v_GetStaticCondBlock(), v_InitObject(), and v_PreSolve().

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation