#include <Perturbation.h>
|
void | dryRun () |
| Dry run to obtain some needed numbers.
|
|
void | run () |
|
◆ 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
-
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. |
References allEigs, deltaMax, errorTol, F, gap, norm, vect, and X.
◆ 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 ) |
◆ checkIdempotencies()
template<typename Treal, typename Tmatrix, typename Tvector>
void per::Perturbation< Treal, Tmatrix, Tvector >::checkIdempotencies |
( |
std::vector< Treal > & | idemErrors | ) |
|
◆ checkMaxSubspaceError()
template<typename Treal, typename Tmatrix, typename Tvector>
void per::Perturbation< Treal, Tmatrix, Tvector >::checkMaxSubspaceError |
( |
Treal & | subsError | ) |
|
◆ dryRun()
template<typename Treal, typename Tmatrix, typename Tvector>
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(), errorTol, gap, mat::Interval< Treal >::length(), mat::Interval< Treal >::low(), mat::Interval< Treal >::midPoint(), nIter, sigma, template_blas_fabs(), threshVal, and mat::Interval< Treal >::upp().
Referenced by perturb().
◆ perturb()
template<typename Treal, typename Tmatrix, typename Tvector>
◆ run()
template<typename Treal, typename Tmatrix, typename Tvector>
◆ allEigs
template<typename Treal, typename Tmatrix, typename Tvector>
◆ deltaMax
template<typename Treal, typename Tmatrix, typename Tvector>
◆ errorTol
template<typename Treal, typename Tmatrix, typename Tvector>
template<typename Treal, typename Tmatrix, typename Tvector>
◆ gap
template<typename Treal, typename Tmatrix, typename Tvector>
◆ nIter
template<typename Treal, typename Tmatrix, typename Tvector>
◆ norm
template<typename Treal, typename Tmatrix, typename Tvector>
◆ sigma
template<typename Treal, typename Tmatrix, typename Tvector>
◆ threshVal
template<typename Treal, typename Tmatrix, typename Tvector>
◆ vect
template<typename Treal, typename Tmatrix, typename Tvector>
template<typename Treal, typename Tmatrix, typename Tvector>
The documentation for this class was generated from the following file: