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

#include <Perturbation.h>

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. More...
 
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

◆ Perturbation()

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

◆ checkCommutators()

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().

◆ checkIdempotencies()

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

Referenced by mainFun().

◆ checkMaxSubspaceError()

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

Referenced by mainFun().

◆ dryRun()

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().

◆ perturb()

template<typename Treal , typename Tmatrix , typename Tvector >
void per::Perturbation< Treal, Tmatrix, Tvector >::perturb ( )
inline

◆ run()

template<typename Treal , typename Tmatrix , typename Tvector >
void per::Perturbation< Treal, Tmatrix, Tvector >::run
protected

Member Data Documentation

◆ allEigs

template<typename Treal , typename Tmatrix , typename Tvector >
mat::Interval<Treal> const& per::Perturbation< Treal, Tmatrix, Tvector >::allEigs
protected

◆ deltaMax

template<typename Treal , typename Tmatrix , typename Tvector >
Treal per::Perturbation< Treal, Tmatrix, Tvector >::deltaMax
protected

◆ errorTol

template<typename Treal , typename Tmatrix , typename Tvector >
Treal per::Perturbation< Treal, Tmatrix, Tvector >::errorTol
protected

◆ F

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

◆ gap

template<typename Treal , typename Tmatrix , typename Tvector >
mat::Interval<Treal> per::Perturbation< Treal, Tmatrix, Tvector >::gap
protected

◆ nIter

template<typename Treal , typename Tmatrix , typename Tvector >
int per::Perturbation< Treal, Tmatrix, Tvector >::nIter
protected

◆ norm

template<typename Treal , typename Tmatrix , typename Tvector >
mat::normType const per::Perturbation< Treal, Tmatrix, Tvector >::norm
protected

◆ sigma

template<typename Treal , typename Tmatrix , typename Tvector >
std::vector<Treal> per::Perturbation< Treal, Tmatrix, Tvector >::sigma
protected

◆ threshVal

template<typename Treal , typename Tmatrix , typename Tvector >
std::vector<Treal> per::Perturbation< Treal, Tmatrix, Tvector >::threshVal
protected

◆ vect

template<typename Treal , typename Tmatrix , typename Tvector >
Tvector& per::Perturbation< Treal, Tmatrix, Tvector >::vect
protected

◆ X

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: