implements SOR preconditioner More...

#include <TSOR.hh>

Inheritance diagram for TSOR:

Public Member Functions
	TSOR (const TSparseMatrix *A, const sor_type_t sor_type, const real omega=1.0, const real damping=1.0)

const TSparseMatrix *	matrix () const
	return internal sparse matrix

sor_type_t	sor_type () const
	return SOR type

real	damping_factor () const
	return damping factor

bool	is_complex () const
	return true, if field type is complex

bool	is_self_adjoint () const
	return true, of operator is self adjoint

virtual void	apply (const TVector x, TVector y, const matop_t op) const

virtual void	apply_add (const real alpha, const TVector x, TVector y, const matop_t op) const

virtual TVector *	domain_vector () const
	return vector in domain space

virtual TVector *	range_vector () const
	return vector in range space

Public Member Functions inherited from TLinearOperator
virtual typeid_t	type () const
	return type ID of matrix

Public Member Functions inherited from TTypeInfo
virtual bool	is_type (const typeid_t t) const
	return true if local object is of given type ID t

virtual std::string	typestr () const
	return string representation of type

Detailed Description

TSOR provides application of a SOR type preconditioner for a given matrix $A = D - L - U$ , with diagonal D, strictly lower triangular matrix L and strictly upper triangular matrix U.

For forward SOR, the applied operator is $\omega ( D - \omega L )^{-1}$ , for backward SOR $\omega ( D - \omega U )^{-1}$ and the combination of both for the symmetric SOR.

The matrix A must be a sparse matrix.

Constructor & Destructor Documentation

TSOR	(	const TSparseMatrix *	A,
		const sor_type_t	sor_type,
		const real	omega = `1.0`,
		const real	damping = `1.0`
	)

constructs SOR preconditioner of type sor_type using coefficients as defined by A

Member Function Documentation

virtual void apply	(	const TVector *	x,
		TVector *	y,
		const matop_t	op
	)		const

virtual

mapping function of linear operator $A$ , e.g. $ y := A(x)$ . Depending on op, either $A$ , $A^T$ or $A^H$ is applied.

Implements TLinearOperator.

virtual void apply_add	(	const real	alpha,
		const TVector *	x,
		TVector *	y,
		const matop_t	op
	)		const

virtual

mapping function with update: $y := y + \alpha A(x)$ . Depending on op, either $A$ , $A^T$ or $A^H$ is applied.

Implements TLinearOperator.

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation