Eigensolver for general Eigenproblems for complex matrices. More...

#include <arpackpp.hh>

Inheritance diagram for eigensolver::ArPackppCmplxGen:

Collaboration diagram for eigensolver::ArPackppCmplxGen:

Public Member Functions
	ArPackppCmplxGen (concepts::Operator< concepts::Cmplx > &OP, concepts::Operator< concepts::Cmplx > &A, concepts::Operator< concepts::Cmplx > &B, int kmax=1, concepts::Cmplx sigma=0.0, Real tol=0.01, int maxiter=300, concepts::Vector< concepts::Cmplx > start=0, concepts::Array< concepts::Cmplx > resid=0, bool schur=false)
	Constructor.
	ArPackppCmplxGen (concepts::Operator< concepts::Cmplx > &OP, concepts::Operator< concepts::Cmplx > &A, concepts::Operator< concepts::Cmplx > &B, int kmax=1, ArPack< Real >::which which=ArPack< Real >::LM, Real tol=0.0, int maxiter=300, concepts::Vector< concepts::Cmplx > start=0, concepts::Array< concepts::Cmplx > resid=0, bool schur=false)
	Constructor.
virtual uint	converged () const
	Returns the number of converged eigenpairs.
virtual concepts::Array < concepts::Vector< Cmplx > * > &	getEF ()
	Getter for the computed eigenvectorsComputes Eigenvalus and Eigenvectors if they were not computed.
virtual const concepts::Array < Cmplx > &	getEV ()
	Returns an array with the eigen values.
concepts::Array< Cmplx >	getRESID ()
	Returns the RESID vector:
virtual uint	iterations () const
	returns the actual number of Arnoldi iterations taken.
virtual	~ArPackppCmplxGen ()
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< Cmplx >	aow_
ARCompGenEig< Real, ArpackOperatorWrapper< Cmplx > , ArpackOperatorWrapper< Cmplx > >	arpackSolver_
bool	computed_
uint	convEigenvalues_
concepts::Array< Cmplx >	eigenValues_
	References into storage for eigenvectors.
concepts::Array < concepts::Vector< Cmplx > * >	eigenVectors_
uint	iter_
ArPack< Real >::modus	modus_
uint	n_
ArPack< Real >::which	which_
char *	whichEV_

Detailed Description

Eigensolver for general Eigenproblems for complex matrices.

Using Cayley mode of arpack++

Definition at line 349 of file arpackpp.hh.

Constructor & Destructor Documentation

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

Constructor.

Parameters:

OP	multiplication operator (A-B*sigma)^-1
A	Stiffness matrix
B	is the Matrix on the right hand side and had to be symmetric and positive definite
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.
maxiter	Maximum number of Arnoldi iterations allowed
target	What sort of eigenvalues to compute
mode	Mode in which ARPACK should be used
sigma	Shift for the shift-invert modes

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

Constructor.

Parameters:

OP	multiplication operator for B^-1
A	is the Matrix on the left hand side of the problem Ax = Bx
B	is the Matrix on the right hand side and had to be positive definite
kmax	Number of eigenpairs to be computed must be lower then the dimension of the matrices
which	- which eigenvalue is searched
tol	Convergence tolerance for the eigenpairs.
maxiter	Maximum number of Arnoldi iterations allowed
target	What sort of eigenvalues to compute
mode	Mode in which ARPACK should be used
sigma	Shift for the shift-invert modes