ergo
|
Class template for building Krylov subspaces with Lanczos. More...
#include <Lanczos.h>
Public Member Functions | |
Lanczos (Tmatrix const &AA, Tvector const &startVec, int maxIt=100, int cap=100) | |
void | restart (Tvector const &startVec) |
virtual void | run () |
void | copyTridiag (MatrixTridiagSymmetric< Treal > &Tricopy) |
virtual | ~Lanczos () |
Protected Member Functions | |
void | increaseCapacity (int const newCapacity) |
void | step () |
void | getEigVector (Tvector &eigVec, Treal const *const eVecTri) const |
virtual void | update ()=0 |
virtual bool | converged () const =0 |
Protected Attributes | |
Tmatrix const & | A |
Tvector * | v |
Tvector | r |
Vectors spanning Krylov subspace. More... | |
MatrixTridiagSymmetric< Treal > | Tri |
Residual vector. More... | |
int | capacity |
int | j |
int | maxIter |
Current step. More... | |
Private Attributes | |
Treal | alpha |
Treal | beta |
Class template for building Krylov subspaces with Lanczos.
Build up Krylov subspace for symmetric matrix with a Lanczos process.
Treal: Type for real numbers
Tmatrix: The matrix class
Tvector: Vector class
|
inline |
|
inlinevirtual |
References mat::arn::Lanczos< Treal, Tmatrix, Tvector >::v.
|
protectedpure virtual |
|
inline |
|
protected |
|
protected |
|
inline |
References mat::arn::Lanczos< Treal, Tmatrix, Tvector >::alpha, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::beta, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::capacity, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::j, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::r, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::Tri, and mat::arn::Lanczos< Treal, Tmatrix, Tvector >::v.
Referenced by mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::rerun().
|
inlinevirtual |
Reimplemented in mat::arn::LanczosLargestMagnitudeEigIfSmall< Treal, Tmatrix, Tvector >, and mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >.
References mat::arn::Lanczos< Treal, Tmatrix, Tvector >::converged(), mat::arn::Lanczos< Treal, Tmatrix, Tvector >::j, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::maxIter, mat::arn::Lanczos< Treal, Tmatrix, Tvector >::step(), and mat::arn::Lanczos< Treal, Tmatrix, Tvector >::update().
Referenced by mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::rerun(), mat::arn::LanczosLargestMagnitudeEig< Treal, Tmatrix, Tvector >::run(), and mat::arn::LanczosLargestMagnitudeEigIfSmall< Treal, Tmatrix, Tvector >::run().
|
protected |
References A, and mat::transpose().
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::run().
|
protectedpure virtual |
|
protected |
|
private |
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::restart().
|
private |
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::restart().
|
protected |
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::restart().
|
protected |
|
protected |
Current step.
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::run().
|
protected |
Vectors spanning Krylov subspace.
In step j: Vectors 0 : j-2 is on file Vectors j-1 : j is in memory
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::Lanczos(), and mat::arn::Lanczos< Treal, Tmatrix, Tvector >::restart().
|
protected |
Residual vector.
Referenced by mat::arn::Lanczos< Treal, Tmatrix, Tvector >::copyTridiag(), and mat::arn::Lanczos< Treal, Tmatrix, Tvector >::restart().
|
protected |