per::Perturbation< Treal, Tmatrix, Tvector > Class Template Reference

#include <Perturbation.h>

List of all members.

Public Member Functions
	Perturbation (std::vector< Tmatrix * > const &F, std::vector< Tmatrix * > &D, mat::Interval< Treal > const &gap, mat::Interval< Treal > const &allEigs, Treal const deltaMax, Treal const errorTol, mat::normType const norm, Tvector &vect)
void	perturb ()
void	checkIdempotencies (std::vector< Treal > &idemErrors)
template<typename TmatNoSymm>
void	checkCommutators (std::vector< Treal > &commErrors, TmatNoSymm const &dummyMat)
void	checkMaxSubspaceError (Treal &subsError)
Protected Member Functions
void	dryRun ()
	Dry run to obtain some needed numbers.
void	run ()
Protected Attributes
std::vector< Tmatrix * > const &	F
std::vector< Tmatrix * > &	X
mat::Interval< Treal >	gap
mat::Interval< Treal > const &	allEigs
Treal	deltaMax
Treal	errorTol
mat::normType const	norm
Tvector &	vect
int	nIter
std::vector< Treal >	threshVal
std::vector< Treal >	sigma

template<typename Treal, typename Tmatrix, typename Tvector>
class per::Perturbation< Treal, Tmatrix, Tvector >

---

Constructor & Destructor Documentation

template<typename Treal, typename Tmatrix, typename Tvector>

per::Perturbation< Treal, Tmatrix, Tvector >::Perturbation	(	std::vector< Tmatrix * > const &	F,
		std::vector< Tmatrix * > &	D,
		mat::Interval< Treal > const &	gap,
		mat::Interval< Treal > const &	allEigs,
		Treal const	deltaMax,
		Treal const	errorTol,
		mat::normType const	norm,
		Tvector &	vect
	)

Parameters:

F	Vector with matrices (input).
D	Vector with matrices (output).
gap	Band gap.
allEigs	Interval containing all eigenvalues of X0 + delta*X1 + delta^2*X2 + ... for all delta in [0, deltaMax] for initial X.
deltaMax	Largest allowed delta.
errorTol	Error tolerance.
norm	Norm for truncation etc.
vect	Vector.

---

Member Function Documentation

template<typename Treal, typename Tmatrix, typename Tvector>

template<typename TmatNoSymm>

void per::Perturbation< Treal, Tmatrix, Tvector >::checkCommutators	(	std::vector< Treal > &	commErrors,
		TmatNoSymm const &	dummyMat
	)

template<typename Treal, typename Tmatrix, typename Tvector>

void per::Perturbation< Treal, Tmatrix, Tvector >::checkIdempotencies

(

std::vector< Treal > &

idemErrors

)

template<typename Treal, typename Tmatrix, typename Tvector>

void per::Perturbation< Treal, Tmatrix, Tvector >::checkMaxSubspaceError

(

Treal &

subsError

)

template<typename Treal, typename Tmatrix, typename Tvector>

void per::Perturbation< Treal, Tmatrix, Tvector >::dryRun

(

)

[protected]

Dry run to obtain some needed numbers.

After call to this function we know:

• number of iterations (nIter),
• threshold values (threshVal), and
• polyunomials to choose (sigma = -1 | = 1)

If requested accuracy is too high or gap too small, an exception is thrown.

template<typename Treal, typename Tmatrix, typename Tvector>

void per::Perturbation< Treal, Tmatrix, Tvector >::perturb

(

)

[inline]

template<typename Treal, typename Tmatrix, typename Tvector>

void per::Perturbation< Treal, Tmatrix, Tvector >::run

(

)

[protected]

---

Member Data Documentation

template<typename Treal, typename Tmatrix, typename Tvector>

mat::Interval<Treal> const& per::Perturbation< Treal, Tmatrix, Tvector >::allEigs [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

Treal per::Perturbation< Treal, Tmatrix, Tvector >::deltaMax [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

Treal per::Perturbation< Treal, Tmatrix, Tvector >::errorTol [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

std::vector<Tmatrix *> const& per::Perturbation< Treal, Tmatrix, Tvector >::F [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

mat::Interval<Treal> per::Perturbation< Treal, Tmatrix, Tvector >::gap [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

int per::Perturbation< Treal, Tmatrix, Tvector >::nIter [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

mat::normType const per::Perturbation< Treal, Tmatrix, Tvector >::norm [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

std::vector<Treal> per::Perturbation< Treal, Tmatrix, Tvector >::sigma [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

std::vector<Treal> per::Perturbation< Treal, Tmatrix, Tvector >::threshVal [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

Tvector& per::Perturbation< Treal, Tmatrix, Tvector >::vect [protected]

template<typename Treal, typename Tmatrix, typename Tvector>

std::vector<Tmatrix *>& per::Perturbation< Treal, Tmatrix, Tvector >::X [protected]

---

The documentation for this class was generated from the following file:

• source/matrix/Perturbation.h

---