SSOR-preconditioned CG-based preconditioner. More...

#include <preconditioners.h>

Inheritance diagram for SSORCGPreconditioner:

Collaboration diagram for SSORCGPreconditioner:

Public Member Functions
	SSORCGPreconditioner (Parallel const &par, InputFile const &inp, std::vector< std::pair< int, int >> const &ll, Bspline const &bspline, CommandLine const &cmd)

virtual RadialIntegrals const &	rad () const
	Get radial integrals. More...

virtual void	rhs (cArrayView chi, int ienergy, int instate) const
	Return the right-hand side. More...

virtual void	multiply (const cArrayView p, cArrayView q) const
	Multiply by the matrix equation. More...

virtual void	precondition (const cArrayView r, cArrayView z) const
	Precondition the equation. More...

virtual void	setup ()
	Initialize the preconditioner. More...

virtual void	update (double E)
	Update the preconditioner for the next energy. More...

virtual void	CG_prec (int iblock, const cArrayView r, cArrayView z) const

Public Member Functions inherited from CGPreconditioner
	CGPreconditioner (Parallel const &par, InputFile const &inp, std::vector< std::pair< int, int >> const &ll, Bspline const &bspline, CommandLine const &cmd)

virtual void	CG_mmul (int iblock, const cArrayView p, cArrayView q) const

Public Member Functions inherited from NoPreconditioner
	NoPreconditioner (Parallel const &par, InputFile const &inp, std::vector< std::pair< int, int >> const &ll, Bspline const &bspline, CommandLine const &cmd)

Static Public Attributes
static const std::string	name = "SSOR"

static const std::string	description = "Block inversion using SSOR-preconditioned conjugate gradients."

Static Public Attributes inherited from CGPreconditioner
static const std::string	name = "cg"

static const std::string	description = "Block inversion using plain conjugate gradients."

Static Public Attributes inherited from NoPreconditioner
static const std::string	name = "none"

static const std::string	description = "Preconditioning by the identity matrix."

Protected Attributes
std::vector< SymDiaMatrix >	SSOR_

Protected Attributes inherited from NoPreconditioner
CommandLine const &	cmd_

Parallel const &	par_

InputFile const &	inp_

std::vector< std::pair< int, int > > const &	l1_l2_

std::vector< SymDiaMatrix >	dia_blocks_

Bspline	s_bspline_

RadialIntegrals	s_rad_

SymDiaMatrix	S_kron_S_

SymDiaMatrix	S_kron_Mm1_tr_

SymDiaMatrix	S_kron_Mm2_

SymDiaMatrix	Mm1_tr_kron_S_

SymDiaMatrix	Mm2_kron_S_

SymDiaMatrix	half_D_minus_Mm1_tr_

SymDiaMatrix	half_D_minus_Mm1_tr_kron_S_

SymDiaMatrix	S_kron_half_D_minus_Mm1_tr_

Detailed Description

Enhances CGPreconditioner conjugate gradients solver by SSOR preconditioning. This is done by redefining virtual function CG_prec.

Constructor & Destructor Documentation

SSORCGPreconditioner::SSORCGPreconditioner	(	Parallel const &	par,
		InputFile const &	inp,
		std::vector< std::pair< int, int >> const &	ll,
		Bspline const &	bspline,
		CommandLine const &	cmd
	)

inline

Member Function Documentation

void SSORCGPreconditioner::CG_prec	(	int	iblock,
		const cArrayView	r,
		cArrayView	z
	)		const

virtual

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.

virtual void SSORCGPreconditioner::multiply	(	const cArrayView	p,
		cArrayView	q
	)		const

inlinevirtual

This function implements matrix multiplication by the matrix of the set of equations that is to be solved.

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.

virtual void SSORCGPreconditioner::precondition	(	const cArrayView	r,
		cArrayView	z
	)		const

inlinevirtual

This function preconditions the equation, solving the preconditioner equation

\[ \mathbf{M}\mathbf{z} = \mathbf{r} \ . \]

It may use the MPI environment.

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.

virtual RadialIntegrals const& SSORCGPreconditioner::rad ( ) const

inlinevirtual

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.

virtual void SSORCGPreconditioner::rhs	(	cArrayView	chi,
		int	ienergy,
		int	instate
	)		const

inlinevirtual

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.

void SSORCGPreconditioner::setup ( )

virtual

This function contains all computation intensive preparations for the preconditioner, e.g. computation of radial integrals. It may use only SMP environment.

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.

void SSORCGPreconditioner::update ( double E )

virtual

This function updates the preconditioner for another right hand side. It may use the MPI environment.

Reimplemented from CGPreconditioner.

Reimplemented in TwoLevelPreconditioner.