ergo
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

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:
FVector with matrices (input).
DVector with matrices (output).
gapBand gap.
allEigsInterval containing all eigenvalues of X0 + delta*X1 + delta^2*X2 + ... for all delta in [0, deltaMax] for initial X.
deltaMaxLargest allowed delta.
errorTolError tolerance.
normNorm for truncation etc.
vectVector.

References per::Perturbation< Treal, Tmatrix, Tvector >::allEigs, per::Perturbation< Treal, Tmatrix, Tvector >::F, per::Perturbation< Treal, Tmatrix, Tvector >::gap, and per::Perturbation< Treal, Tmatrix, Tvector >::X.


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 
)

Referenced by mainFun().

template<typename Treal , typename Tmatrix , typename Tvector >
void per::Perturbation< Treal, Tmatrix, Tvector >::checkIdempotencies ( std::vector< Treal > &  idemErrors)

Referenced by mainFun().

template<typename Treal , typename Tmatrix , typename Tvector >
void per::Perturbation< Treal, Tmatrix, Tvector >::checkMaxSubspaceError ( Treal &  subsError)

Referenced by mainFun().

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.

References mat::Interval< Treal >::empty(), mat::Interval< Treal >::length(), mat::Interval< Treal >::low(), mat::Interval< Treal >::midPoint(), template_blas_fabs(), and mat::Interval< Treal >::upp().

Referenced by per::Perturbation< Treal, Tmatrix, Tvector >::perturb().

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: