RWHermEigServer< TypeT > Class Template Reference
[Symmetric Eigenvalue Decomposition]

Abstract base class for the Hermitian eigenvalue server. More...

#include <rw/lapack/heigsrv.h>

Inheritance diagram for RWHermEigServer< TypeT >:

List of all members.

Public Member Functions
virtual RWHermEigDecomp< TypeT >	operator() (const RWHermMat< TypeT > &A)
virtual RWHermEigDecomp< TypeT >	operator() (const RWHermBandMat< TypeT > &A)
virtual RWSymEigDecomp< double >	decompose (const RWTriDiagDecomp< TypeT > &A)=0
virtual bool	computeEigenVectors () const =0

---

Detailed Description

template<class TypeT>
class RWHermEigServer< TypeT >

The class RWHermEigServer<T> is an abstract base class for the Hermitian eigenvalue server. For each precision, the Linear Algebra Module provides several derived types:

• The class RWHermQREigServer<T> uses the implicitly shifted QR method, switching between QR and QL variants in order to handle graded matrices effectively. This is the default server.
• The class RWHermPDQREigServer<T> uses a combination of Cholesky factorization and bidiagonal QR iteration. It applies to positive definite matrices only. The resulting eigenvalue decomposition may be significantly more accurate than those of the other servers.
• The class RWHermRFQREigServer<T> uses a version of the QR algorithm without square roots. It can only compute all the eigenvalues, and no eigenvectors.
• The class RWHermSomeEigServer<T> uses the bisection method to compute some of the eigenvalues and (optionally) eigenvectors. You can use this server to compute the first k largest eigenvalues, the first k smallest eigenvalues, or the i ^th through j ^th largest eigenvalues.
• The class RWHermRangeEigServer<T> uses bisection to compute all the eigenvalues and (optionally) eigenvectors in a real range.

This base class serves two purposes:

• to provide a specification of the commonality between servers, and
• to house implementation code common to several different servers.

The common implementation includes a set of protected routines to modify a given RWHermEigDecomp<T> object. This way, all server classes, including any that you write, can access the internals of the decomposition object without being explicitly listed as friend classes.

Synopsis

 #include <rw/lapack/heigsrv.h>
 
 RWHermEigServer<DComplex> server;
 RWHermEigDecomp<DComplex> eig = server(A);  // A is a RWHermMat<DComplex>

Examples

See the specific server subclasses for examples.

---

Member Function Documentation

template<class TypeT>

virtual bool RWHermEigServer< TypeT >::computeEigenVectors

(

)

const [pure virtual]

Returns true if this server is configured to compute eigenvectors as well as eigenvalues.

Implemented in RWHermQREigServer< TypeT >, RWHermPDQREigServer< TypeT >, RWHermRFQREigServer< TypeT >, RWHermSomeEigServer< TypeT >, and RWHermRangeEigServer< TypeT >.

template<class TypeT>

virtual RWSymEigDecomp<double> RWHermEigServer< TypeT >::decompose

(

const RWTriDiagDecomp< TypeT > &

A

)

[pure virtual]

Computes the eigenvalue decomposition of the tridiagonal matrix inside the tridiagonal decomposition. This pure virtual function is implemented in the specific server classes. It is used to implement the default operator()() functions.

Implemented in RWHermQREigServer< TypeT >, RWHermPDQREigServer< TypeT >, RWHermRFQREigServer< TypeT >, RWHermSomeEigServer< TypeT >, and RWHermRangeEigServer< TypeT >.

template<class TypeT>

virtual RWHermEigDecomp<TypeT> RWHermEigServer< TypeT >::operator()

(

const RWHermBandMat< TypeT > &

A

)

[virtual]

Computes a Hermitian eigenvalue decomposition. The default method implemented in this class first constructs a tridiagonal decomposition from the given matrix, next computes the eigenvalue decomposition for that decomposition, and finally transforms the eigenvectors as required.

If self is an instance of RWHermPDQREigServer<T>, the matrix you pass in must be positive definite.

template<class TypeT>

virtual RWHermEigDecomp<TypeT> RWHermEigServer< TypeT >::operator()

(

const RWHermMat< TypeT > &

A

)

[virtual]

Computes a Hermitian eigenvalue decomposition. The default method implemented in this class first constructs a tridiagonal decomposition from the given matrix, next computes the eigenvalue decomposition for that decomposition, and finally transforms the eigenvectors as required.

If self is an instance of RWHermPDQREigServer<T>, the matrix you pass in must be positive definite.