Eigensolver for general Eigenproblems for symmetric matrices. More...

#include <arpackpp.hh>

Inheritance diagram for eigensolver::ArPackppSymGen:

Collaboration diagram for eigensolver::ArPackppSymGen:

Public Member Functions
	ArPackppSymGen (concepts::Operator< Real > &OP, concepts::Operator< Real > &A, concepts::Operator< Real > &B, int kmax=1, Real sigma=0.1, Real tol=0.0, int maxiter=300, concepts::Vector< Real > start=0, concepts::Array< Real > resid=0, bool schur=false)
	Constructor for SHIFTandINVERT mode, Cayley mode.
	ArPackppSymGen (char invertType, concepts::Operator< Real > &OP, concepts::Operator< Real > &A, int kmax=1, Real sigma=0.1, Real tol=0.0, int maxiter=300, concepts::Vector< Real > start=0, concepts::Array< Real > resid=0, bool schur=false)
	Constructor for SHIFTandINVERT mode and Buckling mode.
virtual uint	converged () const
	Returns the number of converged eigenpairs.
virtual concepts::Array < concepts::Vector< Real > * > &	getEF ()
	Getter for the computed eigenvectorsComputes Eigenvalus and Eigenvectors if they were not computed.
virtual const concepts::Array < Real > &	getEV ()
	Getter for the computed eigenvalues.
concepts::Array< Real >	getRESID ()
	Returns the RESID vector:
virtual uint	iterations () const
	Returns the actual number of Arnoldi iterations taken.
virtual	~ArPackppSymGen ()
	Constructor for SHIFTandINVERT mode and Buckling mode.
Protected Member Functions
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
Private Member Functions
void	compute_ ()
Private Attributes
ArpackOperatorWrapper < concepts::Real, concepts::Real, concepts::Real >	aow_
ARSymGenEig< Real, ArpackOperatorWrapper < concepts::Real, concepts::Real, concepts::Real > , ArpackOperatorWrapper < concepts::Real, concepts::Real, concepts::Real > >	arpackSolver_
bool	computed_
uint	convEigenvalues_
concepts::Array< Real >	eigenValues_
	References into storage for eigenvectors.
concepts::Array < concepts::Vector< Real > * >	eigenVectors_
uint	iter_
ArPack< Real >::modus	modus_
uint	n_
char *	whichEV_

Detailed Description

Eigensolver for general Eigenproblems for symmetric matrices.

Using Cayley mode of Arpack.

Definition at line 303 of file arpackpp.hh.

Constructor & Destructor Documentation

eigensolver::ArPackppSymGen::ArPackppSymGen	(	concepts::Operator< Real > &	OP,
		concepts::Operator< Real > &	A,
		concepts::Operator< Real > &	B,
		int	kmax = `1`,
		Real	sigma = `0.1`,
		Real	tol = `0.0`,
		int	maxiter = `300`,
		concepts::Vector< Real > *	start = `0`,
		concepts::Array< Real > *	resid = `0`,
		bool	schur = `false`
	)

Constructor for SHIFTandINVERT mode, Cayley mode.

Warning:: B have to be a symmetric positive definite matrix

Parameters:

OP	is the multiplication operator (A-sigma)^(-1)
A	Stiffness matrix
B	is the Matrix on the right hand side
kmax	Number of eigenpairs to be computed, must be lower then the dimension of the matrices -1

Warning:: kmax have to be lower than the dimension of the matrix -1!

Parameters:

tol	Convergence tolerance for the eigenpairs. The default value 0.0 is replaced by `DLAMCH`('EPS') from LAPACK.
maxiter	Maximum number of Arnoldi iterations allowed
target	What sort of eigenvalues to compute
sigma	Shift for the shift-invert modes

Warning:: In the cayley mode sigma should be a real number that is NOT ZERO!

Parameters:

start	Initial vector for iteration
schur	Calculate Schur vector basis instead of eigenvectors

eigensolver::ArPackppSymGen::ArPackppSymGen	(	char	invertType,
		concepts::Operator< Real > &	OP,
		concepts::Operator< Real > &	A,
		int	kmax = `1`,
		Real	sigma = `0.1`,
		Real	tol = `0.0`,
		int	maxiter = `300`,
		concepts::Vector< Real > *	start = `0`,
		concepts::Array< Real > *	resid = `0`,
		bool	schur = `false`
	)

Constructor for SHIFTandINVERT mode and Buckling mode.

Warning:: B have to be a symmetric positive definite matrix

Parameters:

OP	is the multiplication operator (A-sigma)^(-1)
A	is the matrix of the left hand side in Buckling Mode (invertType = 'B')
A	is the Matrix on the right hand side in standard mode (invertType = 'S')
kmax	Number of eigenpairs to be computed, must be lower then the dimension of the matrices -1

Warning:: kmax have to be lower than the dimension of the matrix -1!

Parameters:

tol	Convergence tolerance for the eigenpairs. The default value 0.0 is replaced by `DLAMCH`('EPS') from LAPACK.
maxiter	Maximum number of Arnoldi iterations allowed
target	What sort of eigenvalues to compute
sigma	Shift for the shift-invert modes

Warning:: In the buckling mode sigma should be a real number that is NOT ZERO!

Parameters:

start	Initial vector for iteration
schur	Calculate Schur vector basis instead of eigenvectors

virtual eigensolver::ArPackppSymGen::~ArPackppSymGen ( ) [inline, virtual]

Constructor for SHIFTandINVERT mode and Buckling mode.

Warning:: B have to be a symmetric positive definite matrix

Parameters:

OP	is the multiplication operator B^-1
A	Stiffness matrix
B	is the Matrix on the right hand side
kmax	Number of eigenpairs to be computed, must be lower then the dimension of the matrices -1

Warning:: kmax have to be lower than the dimension of the matrix -1!

Parameters:

tol	Convergence tolerance for the eigenpairs. The default value 0.0 is replaced by `DLAMCH`('EPS') from LAPACK.
maxiter	Maximum number of Arnoldi iterations allowed
target	What sort of eigenvalues to compute
sigma	Shift for the shift-invert modes