#include <SOR.hpp>

Inheritance diagram for lama::SOR:

Data Structures
struct	SORRuntime
Public Member Functions
	SOR (const std::string &id)
	SOR (const std::string &id, const Scalar omega)
	SOR (const std::string &id, LoggerPtr logger)
	SOR (const std::string &id, const Scalar omega, LoggerPtr logger)
	SOR (const SOR &other)
	Copy constructor that copies the status independent solver information.
virtual	~SOR ()
void	initialize (const Matrix &coefficients)
	This abstract method is used by derived solvers to initialize a omega solver.
void	iterate ()
	This abstract method is used by derived solvers to represent a solver iteration.
virtual SORRuntime &	getRuntime ()
	Returns the complete configuration of the derived class.
virtual const SORRuntime &	getConstRuntime () const
	Returns the complete const configuration of the derived class.
virtual SolverPtr	copy ()
	Copies the status independent solver informations to create a new instance of the same type.
void	setOmega (const Scalar omega)
	Sets the omega parameter of this.
Scalar	getOmega () const
	Returns omega.
virtual void	solveImpl ()
	Solves the equation system.
void	setStoppingCriterion (const CriterionPtr criterion)
	set a new StoppingCriterion to the solver.
void	setPreconditioner (SolverPtr const conditioner)
	Sets the preconditioner of this solver.
const SolverPtr	getPreconditioner () const
	returns the preconditioner of this solver
int	getIterationCount () const
	returns the number of iterations, this solver has done so far.
virtual void	solve (Vector &solution, const Vector &rhs)
	Solves the equation system based on the given rhs.
virtual void	solveInit (Vector &solution, const Vector &rhs)
	Initializes the solver with rhs and solution.
virtual void	solveFinalize ()
	Finalizes the solving process.
const std::string &	getId () const
	Returns the ID of this solver.
const Vector &	getResidual () const
	Contingently calculates the current residual based on the coefficients, rhs and solution currently associated with this.
const Matrix &	getCoefficients () const
	Gets the matrix A from A*u=f.
void	setLogger (LoggerPtr logger)
	Redefines mLogger.
void	setLogLevel (LogLevel::LogLevel level)
	Switches the loglevel of mLogger.
virtual void	setContext (ContextPtr context)
	Sets the context where this solver should be executed.
Protected Member Functions
bool	criteriaAreSatisfied () const
	Checks if all of the stopping criteria are satisfied.
void	logStartSolve ()
	Logging methods to maintain code-readability.
void	logEndSolve ()
void	logIterationEndAndResidual ()
void	logIterationStart ()
virtual void	writeAt (std::ostream &stream) const
	Overrides Printable::writeAt with more useful stuff.
Protected Attributes
SORRuntime	mSORRuntime
Scalar	mOmega
SolverPtr	mPreconditioner
	The preconditioner of this solver.
CriterionPtr	mCriterionRootComponent
	The root stopping criterion evaluated every iteration in the solve method.
std::string	mId
	The ID of this solver.
LoggerPtr	mLogger
	The solver logger.
ContextPtr	mContext
	For forcing the context solver dependent.
Private Member Functions
template<typename T >
void	iterateImpl ()
	LAMA_LOG_DECL_STATIC_LOGGER (logger)
Private Attributes
std::auto_ptr< const Matrix >	mIterationMatrix

Constructor & Destructor Documentation

lama::SOR::SOR ( const std::string & id )

Referenced by copy().

lama::SOR::SOR	(	const std::string &	id,
		const Scalar	omega
	)

lama::SOR::SOR	(	const std::string &	id,
		LoggerPtr	logger
	)

lama::SOR::SOR	(	const std::string &	id,
		const Scalar	omega,
		LoggerPtr	logger
	)

lama::SOR::SOR ( const SOR & other )

Copy constructor that copies the status independent solver information.

lama::SOR::~SOR ( ) [virtual]

Member Function Documentation

SolverPtr lama::SOR::copy ( ) [virtual]

Copies the status independent solver informations to create a new instance of the same type.

Returns:: shared pointer of the copied solver

Implements lama::OmegaSolver.

References SOR().

bool lama::IterativeSolver::criteriaAreSatisfied ( ) const [protected, inherited]

Checks if all of the stopping criteria are satisfied.

return Whether the criteria are satisfied or not

References lama::Solver::getId(), LAMA_THROWEXCEPTION, and lama::IterativeSolver::mCriterionRootComponent.

Referenced by lama::IterativeSolver::solveImpl().

const Matrix & lama::Solver::getCoefficients ( ) const [inherited]

Gets the matrix A from A*u=f.

Returns:: The coefficient matrix A.

References lama::Solver::getConstRuntime(), LAMA_ASSERT_DEBUG, and lama::Solver::SolverRuntime::mCoefficients.

Referenced by lama::InverseSolver::invert(), and lama::LUSolver::piamax_own().

const SOR::SORRuntime & lama::SOR::getConstRuntime ( ) const [virtual]

Returns the complete const configuration of the derived class.

Implements lama::OmegaSolver.

References mSORRuntime.

const std::string & lama::Solver::getId ( ) const [inherited]

Returns the ID of this solver.

Returns:: The ID of this solver.

References lama::Solver::mId.

Referenced by lama::IterativeSolver::criteriaAreSatisfied(), lama::IterativeSolver::solveImpl(), and lama::SolverConfigGrammar::SolverConfigGrammar().

int lama::IterativeSolver::getIterationCount ( ) const [inherited]

returns the number of iterations, this solver has done so far.

References lama::IterativeSolver::getConstRuntime(), and lama::IterativeSolver::IterativeSolverRuntime::mIterations.

Referenced by lama::IterationCount::isSatisfied(), lama::ResidualThreshold::isSatisfied(), lama::CG::iterate(), and lama::GMRES::iterate().

Scalar lama::OmegaSolver::getOmega ( ) const [inherited]

Returns omega.

Returns:: Omega.

References lama::OmegaSolver::mOmega.

const SolverPtr lama::IterativeSolver::getPreconditioner ( ) const [inherited]

returns the preconditioner of this solver

Returns:: the preconditioner

References lama::IterativeSolver::mPreconditioner.

Referenced by lama::MetaSolver::initializePreconditioner(), and lama::IterativeSolver::IterativeSolver().

const Vector & lama::Solver::getResidual ( ) const [inherited]

Contingently calculates the current residual based on the coefficients, rhs and solution currently associated with this.

Should be used internally only, because the three vectors mentioned above have to be initialized.

Returns:: The current residual

References lama::LogLevel::completeInformation, lama::Vector::create(), lama::SolutionProxy::getConstReference(), lama::Solver::getConstRuntime(), lama::SolutionProxy::isDirty(), LAMA_ASSERT_DEBUG, lama::Solver::SolverRuntime::mCoefficients, lama::Solver::mLogger, lama::Solver::SolverRuntime::mResidual, lama::Solver::SolverRuntime::mRhs, lama::Solver::SolverRuntime::mSolution, and lama::SolutionProxy::setDirty().

Referenced by lama::ResidualStagnation::isSatisfied(), lama::ResidualThreshold::isSatisfied(), lama::CG::iterate(), lama::GMRES::iterate(), and lama::Logger::logResidual().

SOR::SORRuntime & lama::SOR::getRuntime ( ) [virtual]

Returns the complete configuration of the derived class.

Implements lama::OmegaSolver.

References mSORRuntime.

Referenced by iterate(), and iterateImpl().

void lama::SOR::initialize ( const Matrix & coefficients ) [virtual]

This abstract method is used by derived solvers to initialize a omega solver.

Reimplemented from lama::OmegaSolver.

References lama::CSR, lama::File::DOUBLE, lama::File::FLOAT, lama::_MatrixStorage::getFormat(), lama::_SparseMatrix::getHaloStorage(), lama::_SparseMatrix::getLocalStorage(), lama::Matrix::getValueType(), LAMA_REGION, LAMA_THROWEXCEPTION, and mIterationMatrix.

void lama::SOR::iterate ( ) [virtual]

This abstract method is used by derived solvers to represent a solver iteration.

Implements lama::OmegaSolver.

References lama::Scalar::DOUBLE, lama::Scalar::FLOAT, getRuntime(), LAMA_REGION, and LAMA_THROWEXCEPTION.

template<typename T >

void lama::SOR::iterateImpl ( ) [private]

References lama::SparseMatrix< T >::getHalo(), lama::SparseMatrix< T >::getHaloStorage(), lama::DenseVector< T >::getHaloValues(), lama::CSRStorage< T >::getIA(), lama::CSRStorage< T >::getJA(), lama::SparseMatrix< T >::getLocalStorage(), lama::DenseVector< T >::getLocalValues(), lama::_MatrixStorage::getNumRows(), getRuntime(), lama::Scalar::getValue(), lama::CSRStorage< T >::getValues(), lama::_MatrixStorage::hasDiagonalProperty(), LAMA_ASSERT_ERROR, LAMA_REGION, mIterationMatrix, lama::OmegaSolver::mOmega, and lama::DenseVector< T >::updateHalo().

lama::SOR::LAMA_LOG_DECL_STATIC_LOGGER ( logger ) [private]

Reimplemented from lama::OmegaSolver.

void lama::IterativeSolver::logEndSolve ( ) [protected, inherited]

References lama::IterativeSolver::getRuntime(), lama::Solver::mLogger, and lama::LogLevel::solverInformation.

Referenced by lama::IterativeSolver::solveImpl().

void lama::IterativeSolver::logIterationEndAndResidual ( ) [protected, inherited]

References lama::LogLevel::advancedInformation, lama::LogLevel::convergenceHistory, lama::IterativeSolver::getRuntime(), lama::l2Norm(), lama::IterativeSolver::IterativeSolverRuntime::mIterations, and lama::Solver::mLogger.

Referenced by lama::IterativeSolver::solveImpl().

void lama::IterativeSolver::logIterationStart ( ) [protected, inherited]

References lama::LogLevel::completeInformation, lama::IterativeSolver::getRuntime(), and lama::Solver::mLogger.

Referenced by lama::IterativeSolver::solveImpl().

void lama::IterativeSolver::logStartSolve ( ) [protected, inherited]

Logging methods to maintain code-readability.

References lama::LogLevel::convergenceHistory, lama::l2Norm(), lama::Solver::mLogger, and lama::LogLevel::solverInformation.

Referenced by lama::IterativeSolver::solveImpl().

void lama::Solver::setContext ( ContextPtr context ) [virtual, inherited]

Sets the context where this solver should be executed.

Sets the context where this solver should be executed. Caution: This overrides the context of the coefficients matrix A from A * u = f used to inializ this solver.

Parameters:

[in] context the context where this solver should be executed.

Reimplemented in lama::InverseSolver.

References lama::Solver::mContext.

Referenced by lama::GMRES::initialize().

void lama::Solver::setLogger ( LoggerPtr logger ) [inherited]

Redefines mLogger.

References lama::Solver::mLogger.

void lama::Solver::setLogLevel ( LogLevel::LogLevel level ) [inherited]

Switches the loglevel of mLogger.

Returns:

References lama::Solver::mLogger.

void lama::OmegaSolver::setOmega ( const Scalar omega ) [inherited]

Sets the omega parameter of this.

Parameters:

[in] omega The omega parameter of the omega solver.

References lama::OmegaSolver::mOmega.

void lama::IterativeSolver::setPreconditioner ( SolverPtr const conditioner ) [inherited]

Sets the preconditioner of this solver.

Preconditioner should be set before initializing the this solver, otherwise the preconditioner initialization should be executed manually

Parameters:

conditioner The preconditioner

References lama::IterativeSolver::mPreconditioner.

void lama::IterativeSolver::setStoppingCriterion ( const CriterionPtr criterion ) [inherited]

set a new StoppingCriterion to the solver.

Parameters:

[in] criterion the new criterion.

References LAMA_ASSERT_ERROR, and lama::IterativeSolver::mCriterionRootComponent.

Referenced by lama::SingleGridSetup::initialize().

void lama::Solver::solve	(	Vector &	solution,
		const Vector &	rhs
	)		`[virtual, inherited]`

Solves the equation system based on the given rhs.

The solver needs to be initialized first with the matrix from the equation to solve, e.g. A from A*u=f (call solver::initialize(A) for example) This method is abstract. It has to be implemented by a class which inherits from this class.

Parameters:

rhs	The right hand side of A*u=f.
solution	The solution from A*u=f. Mostly used as starting solution for an IterativeSolver.

Reimplemented in lama::LUSolver, lama::DefaultJacobi, and lama::SpecializedJacobi.

References lama::Solver::getConstRuntime(), LAMA_REGION, lama::Solver::solveFinalize(), lama::Solver::solveImpl(), and lama::Solver::solveInit().

Referenced by lama::SimpleAMG::cycle().

void lama::Solver::solveFinalize ( ) [virtual, inherited]

Finalizes the solving process.

Reimplemented in lama::DefaultJacobi, and lama::SpecializedJacobi.

Referenced by lama::Solver::solve().

void lama::IterativeSolver::solveImpl ( ) [virtual, inherited]

Solves the equation system.

Rhs and starting solution have to be initialized first! (call solveInit( rhs, solution ) ). The solver needs to be initialized first with the matrix from the equation to solve, e.g. A from A*u=f (call solver::initialize(A) for example)

This method solves the equation system by using the given rhs and solution. For most iterative solvers the solution-vector is used as a starting solution for the solve process. This class does not take responsibility for deleting the vectors after the solver! Make sure you do not delete the vectors during the solver process.

Implements lama::Solver.

References lama::IterativeSolver::criteriaAreSatisfied(), lama::IterativeSolver::getConstRuntime(), lama::Solver::getId(), lama::IterativeSolver::getRuntime(), lama::IterativeSolver::iterate(), LAMA_THROWEXCEPTION, lama::IterativeSolver::logEndSolve(), lama::IterativeSolver::logIterationEndAndResidual(), lama::IterativeSolver::logIterationStart(), lama::IterativeSolver::logStartSolve(), and lama::IterativeSolver::IterativeSolverRuntime::mIterations.