Nektar::LibUtilities::NekSys Class Reference

#include <NekSys.h>

Inheritance diagram for Nektar::LibUtilities::NekSys:

Public Member Functions
	NekSys (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::CommSharedPtr &vComm, const int nDimen, const NekSysKey &pKey)
void	InitObject ()
virtual	~NekSys ()
void	SetSysOperators (const NekSysOperators &in)
const NekSysOperators &	GetSysOperators ()
int	SolveSystem (const int nGlobal, const Array< OneD, const NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const int nDir, const NekDouble tol=1.0E-7, const NekDouble factor=1.0)
bool	ConvergenceCheck (const int nIteration, const Array< OneD, const NekDouble > &Residual, const NekDouble tol=1.0E-7)
virtual void	v_NekSysInitialGuess (const Array< OneD, const NekDouble > &pInput, Array< OneD, NekDouble > &pguess)
void	SetFlagWarnings (bool in)

Static Public Member Functions
static NekSysSharedPtr	CreateInstance (const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::CommSharedPtr &vComm, const int nDimen, const NekSysKey &pKey)

Protected Member Functions
virtual void	v_InitObject ()
virtual int	v_SolveSystem (const int nGlobal, const Array< OneD, const NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const int nDir, const NekDouble tol, const NekDouble factor)
virtual bool	v_ConvergenceCheck (const int nIteration, const Array< OneD, const NekDouble > &Residual, const NekDouble tol)

Protected Attributes
int	m_maxiter
	Maximum iterations. More...
NekDouble	m_tolerance
	Tolerance of iterative solver. More...
LibUtilities::CommSharedPtr	m_Comm
	Communicate. More...
bool	m_converged
	Whether the iteration has been converged. More...
bool	m_root
	Root if parallel. More...
bool	m_verbose
	Verbose. More...
bool	m_FlagWarnings
NekSysOperators	m_operator
	Operators. More...
int	m_SysDimen
	The dimension of the system. More...

Friends
class	MemoryManager< NekSys >
	Support creation through MemoryManager. More...

Detailed Description

Solves a nonlinear or linear system.

Definition at line 220 of file NekSys.h.

Constructor & Destructor Documentation

NekSys()

Nektar::LibUtilities::NekSys::NekSys	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const LibUtilities::CommSharedPtr &	vComm,
		const int	nDimen,
		const NekSysKey &	pKey
	)

Definition at line 50 of file NekSys.cpp.

{
    m_tolerance = pKey.m_Tolerance;
    m_verbose   = false;
    m_root      = false;
    m_Comm      = vComm;
 
    m_FlagWarnings = true;
 
    if (0 == m_Comm->GetRank())
    {
        m_root = true;
    }
    m_verbose = pSession->DefinesCmdLineArgument("verbose");
 
    m_converged = false;
 
    m_SysDimen = nDimen;
}

References Nektar::LibUtilities::NekSysKey::m_Tolerance.

~NekSys()

Nektar::LibUtilities::NekSys::~NekSys ( )

virtual

Definition at line 72 of file NekSys.cpp.

{
}

Member Function Documentation

ConvergenceCheck()

bool Nektar::LibUtilities::NekSys::ConvergenceCheck ( const int  nIteration, const Array< OneD, const NekDouble > &  Residual, const NekDouble  tol = 1.0E-7  )

inline

Definition at line 263 of file NekSys.h.

    {
        return v_ConvergenceCheck(nIteration, Residual, tol);
    }

References v_ConvergenceCheck().

Referenced by Nektar::LibUtilities::NekLinSysIterFixedpointJacobi::v_SolveSystem().

CreateInstance()

static NekSysSharedPtr Nektar::LibUtilities::NekSys::CreateInstance ( const LibUtilities::SessionReaderSharedPtr &  pSession, const LibUtilities::CommSharedPtr &  vComm, const int  nDimen, const NekSysKey &  pKey  )

inlinestatic

Definition at line 226 of file NekSys.h.

    {
        NekSysSharedPtr p = MemoryManager<NekSys>::
            AllocateSharedPtr(pSession, vComm, nDimen, pKey);
        return p;
    }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and CellMLToNektar.cellml_metadata::p.

GetSysOperators()

const NekSysOperators& Nektar::LibUtilities::NekSys::GetSysOperators ( )

inline

Definition at line 250 of file NekSys.h.

    {
        return m_operator;
    }

References m_operator.

InitObject()

void Nektar::LibUtilities::NekSys::InitObject ( )

inline

Definition at line 239 of file NekSys.h.

    {
        v_InitObject();
    }

References v_InitObject().

SetFlagWarnings()

void Nektar::LibUtilities::NekSys::SetFlagWarnings ( bool  in )

inline

Definition at line 274 of file NekSys.h.

    {
        m_FlagWarnings = in;
    }

References m_FlagWarnings.

SetSysOperators()

void Nektar::LibUtilities::NekSys::SetSysOperators ( const NekSysOperators &  in )

inline

Definition at line 245 of file NekSys.h.

    {
        m_operator = in;
    }

References m_operator.

SolveSystem()

int Nektar::LibUtilities::NekSys::SolveSystem ( const int  nGlobal, const Array< OneD, const NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput, const int  nDir, const NekDouble  tol = 1.0E-7, const NekDouble  factor = 1.0  )

inline

Definition at line 255 of file NekSys.h.

    {
        return v_SolveSystem(nGlobal, pInput, pOutput, nDir, tol, factor);
    }

References v_SolveSystem().

v_ConvergenceCheck()

bool Nektar::LibUtilities::NekSys::v_ConvergenceCheck ( const int  nIteration, const Array< OneD, const NekDouble > &  Residual, const NekDouble  tol  )

protectedvirtual

Reimplemented in Nektar::LibUtilities::NekNonlinSysNewton.

Definition at line 76 of file NekSys.cpp.

{
    bool converged = false;
    int ntotal     = Residual.size();
    boost::ignore_unused(nIteration);
 
    NekDouble SysResNorm = Vmath::Dot(ntotal, Residual, Residual);
    m_Comm->AllReduce(SysResNorm, Nektar::LibUtilities::ReduceSum);
 
    if (SysResNorm < tol)
    {
        converged = true;
    }
    return converged;
}

References Vmath::Dot(), and Nektar::LibUtilities::ReduceSum.

Referenced by ConvergenceCheck().

v_InitObject()

virtual void Nektar::LibUtilities::NekSys::v_InitObject ( )

inlineprotectedvirtual

Reimplemented in Nektar::LibUtilities::NekNonlinSysNewton, Nektar::LibUtilities::NekNonlinSys, Nektar::LibUtilities::NekLinSysIterGMRES, Nektar::LibUtilities::NekLinSysIterFixedpointJacobi, Nektar::LibUtilities::NekLinSysIterCG, and Nektar::LibUtilities::NekLinSysIter.

Definition at line 298 of file NekSys.h.

    {
    }

Referenced by InitObject(), Nektar::LibUtilities::NekLinSysIter::v_InitObject(), Nektar::LibUtilities::NekNonlinSys::v_InitObject(), and Nektar::LibUtilities::NekNonlinSysNewton::v_InitObject().

v_NekSysInitialGuess()

void Nektar::LibUtilities::NekSys::v_NekSysInitialGuess ( const Array< OneD, const NekDouble > &  pInput, Array< OneD, NekDouble > &  pguess  )

virtual

Natural guess

Definition at line 97 of file NekSys.cpp.

{
    size_t ndim = pInput.size();
    if(pguess.size() != ndim)
    {
        pguess = Array<OneD, NekDouble> {ndim};
    }
 
    Vmath::Vcopy(ndim, pInput, 1, pguess, 1);
}

References Vmath::Vcopy().

v_SolveSystem()

virtual int Nektar::LibUtilities::NekSys::v_SolveSystem ( const int  nGlobal, const Array< OneD, const NekDouble > &  pInput, Array< OneD, NekDouble > &  pOutput, const int  nDir, const NekDouble  tol, const NekDouble  factor  )

inlineprotectedvirtual

Reimplemented in Nektar::LibUtilities::NekNonlinSysNewton, Nektar::LibUtilities::NekLinSysIterGMRES, Nektar::LibUtilities::NekLinSysIterFixedpointJacobi, and Nektar::LibUtilities::NekLinSysIterCG.

Definition at line 302 of file NekSys.h.

    {
        boost::ignore_unused(nGlobal, pInput, pOutput, nDir, tol, factor);
        ASSERTL0(false, "LinSysIterSolver NOT CORRECT.");
        return 0;
    }

References ASSERTL0.

Referenced by SolveSystem().

Friends And Related Function Documentation

MemoryManager< NekSys >

friend class MemoryManager< NekSys >

friend

Support creation through MemoryManager.

Definition at line 212 of file NekSys.h.

Member Data Documentation

m_Comm

LibUtilities::CommSharedPtr Nektar::LibUtilities::NekSys::m_Comm

protected

Communicate.

Definition at line 285 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterGMRES::DoArnoldi(), Nektar::LibUtilities::NekLinSysIterCG::DoConjugateGradient(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGmresRestart(), Nektar::LibUtilities::NekNonlinSys::NekNonlinSys(), Nektar::LibUtilities::NekLinSysIter::Set_Rhs_Magnitude(), and Nektar::LibUtilities::NekNonlinSysNewton::v_ConvergenceCheck().

m_converged

bool Nektar::LibUtilities::NekSys::m_converged

protected

Whether the iteration has been converged.

Definition at line 287 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGmresRestart(), Nektar::LibUtilities::NekLinSysIterFixedpointJacobi::v_SolveSystem(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SolveSystem().

m_FlagWarnings

bool Nektar::LibUtilities::NekSys::m_FlagWarnings

protected

Definition at line 292 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), SetFlagWarnings(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SolveSystem().

m_maxiter

int Nektar::LibUtilities::NekSys::m_maxiter

protected

Maximum iterations.

Definition at line 281 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterCG::DoConjugateGradient(), Nektar::LibUtilities::NekLinSysIter::NekLinSysIter(), Nektar::LibUtilities::NekLinSysIterGMRES::NekLinSysIterGMRES(), Nektar::LibUtilities::NekNonlinSys::NekNonlinSys(), Nektar::LibUtilities::NekLinSysIter::SetNekLinSysMaxIterations(), Nektar::LibUtilities::NekNonlinSys::SetNekNonlinSysMaxIterations(), Nektar::LibUtilities::NekLinSysIterFixedpointJacobi::v_SolveSystem(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SolveSystem().

m_operator

NekSysOperators Nektar::LibUtilities::NekSys::m_operator

protected

Operators.

Definition at line 294 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterGMRES::DoArnoldi(), Nektar::LibUtilities::NekLinSysIterCG::DoConjugateGradient(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGmresRestart(), GetSysOperators(), SetSysOperators(), Nektar::LibUtilities::NekNonlinSysNewton::v_SetupNekNonlinSystem(), Nektar::LibUtilities::NekLinSysIterFixedpointJacobi::v_SolveSystem(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SolveSystem().

m_root

bool Nektar::LibUtilities::NekSys::m_root

protected

Root if parallel.

Definition at line 289 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterCG::DoConjugateGradient(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SolveSystem().

m_SysDimen

int Nektar::LibUtilities::NekSys::m_SysDimen

protected

The dimension of the system.

Definition at line 296 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekNonlinSys::NekNonlinSys(), Nektar::LibUtilities::NekNonlinSys::SetRefResidual(), Nektar::LibUtilities::NekLinSysIter::setUniversalUniqueMap(), Nektar::LibUtilities::NekNonlinSysNewton::v_InitObject(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SetupNekNonlinSystem().

m_tolerance

NekDouble Nektar::LibUtilities::NekSys::m_tolerance

protected

Tolerance of iterative solver.

Definition at line 283 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterCG::DoConjugateGradient(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGmresRestart(), Nektar::LibUtilities::NekLinSysIter::NekLinSysIter(), Nektar::LibUtilities::NekNonlinSys::NekNonlinSys(), Nektar::LibUtilities::NekLinSysIter::SetNekLinSysTolerance(), Nektar::LibUtilities::NekNonlinSys::SetNekNonlinSysTolerance(), Nektar::LibUtilities::NekLinSysIterCG::v_SolveSystem(), Nektar::LibUtilities::NekLinSysIterFixedpointJacobi::v_SolveSystem(), and Nektar::LibUtilities::NekLinSysIterGMRES::v_SolveSystem().

m_verbose

bool Nektar::LibUtilities::NekSys::m_verbose

protected

Verbose.

Definition at line 291 of file NekSys.h.

Referenced by Nektar::LibUtilities::NekLinSysIterCG::DoConjugateGradient(), Nektar::LibUtilities::NekLinSysIterGMRES::DoGMRES(), and Nektar::LibUtilities::NekNonlinSysNewton::v_SolveSystem().