Grid 0.7.0
ImplicitlyRestartedLanczos< Field > Class Template Reference

#include <ImplicitlyRestartedLanczos.h>

Public Member Functions

 ImplicitlyRestartedLanczos (LinearFunction< Field > &PolyOp, LinearFunction< Field > &HermOp, ImplicitlyRestartedLanczosTester< Field > &Tester, int _Nstop, int _Nk, int _Nm, RealD _eresid, int _MaxIter, RealD _betastp=0.0, int _MinRestart=0, int _orth_period=1, IRLdiagonalisation _diagonalisation=IRLdiagonaliseWithEigen)
 ImplicitlyRestartedLanczos (LinearFunction< Field > &PolyOp, LinearFunction< Field > &HermOp, int _Nstop, int _Nk, int _Nm, RealD _eresid, int _MaxIter, RealD _betastp=0.0, int _MinRestart=0, int _orth_period=1, IRLdiagonalisation _diagonalisation=IRLdiagonaliseWithEigen)
void orthogonalize (Field &w, std::vector< Field > &evec, int k)
void calc (std::vector< RealD > &eval, std::vector< Field > &evec, const Field &src, int &Nconv, bool reverse=false)

Static Public Member Functions

template<typename T>
static RealD normalise (T &v)

Private Member Functions

void step (std::vector< RealD > &lmd, std::vector< RealD > &lme, std::vector< Field > &evec, Field &w, int Nm, int k)
void diagonalize_Eigen (std::vector< RealD > &lmd, std::vector< RealD > &lme, int Nk, int Nm, Eigen::MatrixXd &Qt, GridBase *grid)
void QR_decomp (std::vector< RealD > &lmd, std::vector< RealD > &lme, int Nk, int Nm, Eigen::MatrixXd &Qt, RealD Dsh, int kmin, int kmax)
void diagonalize (std::vector< RealD > &lmd, std::vector< RealD > &lme, int Nk, int Nm, Eigen::MatrixXd &Qt, GridBase *grid)
void diagonalize_lapack (std::vector< RealD > &lmd, std::vector< RealD > &lme, int Nk, int Nm, Eigen::MatrixXd &Qt, GridBase *grid)
void diagonalize_QR (std::vector< RealD > &lmd, std::vector< RealD > &lme, int Nk, int Nm, Eigen::MatrixXd &Qt, GridBase *grid)

Private Attributes

const RealD small
int MaxIter
int MinRestart
int Nstop
int Nk
int Nm
IRLdiagonalisation diagonalisation
int orth_period
RealD OrthoTime
RealD eresid
RealD betastp
LinearFunction< Field > & _PolyOp
LinearFunction< Field > & _HermOp
ImplicitlyRestartedLanczosTester< Field > & _Tester
ImplicitlyRestartedLanczosHermOpTester< Field > SimpleTester

Detailed Description

template<class Field>
class ImplicitlyRestartedLanczos< Field >

Definition at line 97 of file ImplicitlyRestartedLanczos.h.

Constructor & Destructor Documentation

◆ ImplicitlyRestartedLanczos() [1/2]

template<class Field>
ImplicitlyRestartedLanczos< Field >::ImplicitlyRestartedLanczos ( LinearFunction< Field > & PolyOp,
LinearFunction< Field > & HermOp,
ImplicitlyRestartedLanczosTester< Field > & Tester,
int _Nstop,
int _Nk,
int _Nm,
RealD _eresid,
int _MaxIter,
RealD _betastp = 0.0,
int _MinRestart = 0,
int _orth_period = 1,
IRLdiagonalisation _diagonalisation = IRLdiagonaliseWithEigen )
inline

Definition at line 142 of file ImplicitlyRestartedLanczos.h.

References _HermOp, _PolyOp, _Tester, betastp, diagonalisation, eresid, IRLdiagonaliseWithEigen, MaxIter, MinRestart, Nk, Nm, Nstop, orth_period, and SimpleTester.

◆ ImplicitlyRestartedLanczos() [2/2]

template<class Field>
ImplicitlyRestartedLanczos< Field >::ImplicitlyRestartedLanczos ( LinearFunction< Field > & PolyOp,
LinearFunction< Field > & HermOp,
int _Nstop,
int _Nk,
int _Nm,
RealD _eresid,
int _MaxIter,
RealD _betastp = 0.0,
int _MinRestart = 0,
int _orth_period = 1,
IRLdiagonalisation _diagonalisation = IRLdiagonaliseWithEigen )
inline

Definition at line 159 of file ImplicitlyRestartedLanczos.h.

References _HermOp, _PolyOp, _Tester, betastp, diagonalisation, eresid, IRLdiagonaliseWithEigen, MaxIter, MinRestart, Nk, Nm, Nstop, orth_period, and SimpleTester.

Member Function Documentation

◆ normalise()

template<class Field>
template<typename T>
RealD ImplicitlyRestartedLanczos< Field >::normalise ( T & v)
inlinestatic

Definition at line 178 of file ImplicitlyRestartedLanczos.h.

References norm2().

Referenced by calc(), orthogonalize(), and step().

◆ orthogonalize()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::orthogonalize ( Field & w,
std::vector< Field > & evec,
int k )
inline

Definition at line 186 of file ImplicitlyRestartedLanczos.h.

References basisOrthogonalize(), normalise(), OrthoTime, and usecond().

Referenced by step().

◆ calc()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::calc ( std::vector< RealD > & eval,
std::vector< Field > & evec,
const Field & src,
int & Nconv,
bool reverse = false )
inline

big constraint on the basis

Definition at line 211 of file ImplicitlyRestartedLanczos.h.

References _HermOp, _Tester, B, basisRotate(), basisRotateJ(), basisSortInPlace(), diagonalisation, diagonalize(), eresid, eval(), GridLogError, GridLogIRL, GridLogMessage, innerProduct(), IRLdiagonaliseWithDSTEGR, IRLdiagonaliseWithEigen, IRLdiagonaliseWithQR, MaxIter, MinRestart, Nk, Nm, norm2(), normalise(), Nstop, OrthoTime, QR_decomp(), real(), and step().

Referenced by LocalCoherenceLanczos< Fobj, CComplex, nbasis >::calcCoarse(), and LocalCoherenceLanczos< Fobj, CComplex, nbasis >::calcFine().

◆ step()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::step ( std::vector< RealD > & lmd,
std::vector< RealD > & lme,
std::vector< Field > & evec,
Field & w,
int Nm,
int k )
inlineprivate

Definition at line 459 of file ImplicitlyRestartedLanczos.h.

References _PolyOp, GridLogDebug, GridLogIRL, innerProduct(), Nm, normalise(), orth_period, orthogonalize(), and real().

Referenced by calc().

◆ diagonalize_Eigen()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::diagonalize_Eigen ( std::vector< RealD > & lmd,
std::vector< RealD > & lme,
int Nk,
int Nm,
Eigen::MatrixXd & Qt,
GridBase * grid )
inlineprivate

Definition at line 501 of file ImplicitlyRestartedLanczos.h.

References Nk, and Nm.

Referenced by diagonalize().

◆ QR_decomp()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::QR_decomp ( std::vector< RealD > & lmd,
std::vector< RealD > & lme,
int Nk,
int Nm,
Eigen::MatrixXd & Qt,
RealD Dsh,
int kmin,
int kmax )
inlineprivate

Definition at line 527 of file ImplicitlyRestartedLanczos.h.

References Nk, and Nm.

Referenced by calc(), and diagonalize_QR().

◆ diagonalize()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::diagonalize ( std::vector< RealD > & lmd,
std::vector< RealD > & lme,
int Nk,
int Nm,
Eigen::MatrixXd & Qt,
GridBase * grid )
inlineprivate

Definition at line 587 of file ImplicitlyRestartedLanczos.h.

References diagonalisation, diagonalize_Eigen(), diagonalize_lapack(), diagonalize_QR(), IRLdiagonaliseWithDSTEGR, IRLdiagonaliseWithEigen, IRLdiagonaliseWithQR, Nk, and Nm.

Referenced by calc().

◆ diagonalize_lapack()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::diagonalize_lapack ( std::vector< RealD > & lmd,
std::vector< RealD > & lme,
int Nk,
int Nm,
Eigen::MatrixXd & Qt,
GridBase * grid )
inlineprivate

Definition at line 612 of file ImplicitlyRestartedLanczos.h.

References CartesianCommunicator::_Nprocessors, CartesianCommunicator::_processor, CartesianCommunicator::GlobalSumVector(), Nk, and Nm.

Referenced by diagonalize().

◆ diagonalize_QR()

template<class Field>
void ImplicitlyRestartedLanczos< Field >::diagonalize_QR ( std::vector< RealD > & lmd,
std::vector< RealD > & lme,
int Nk,
int Nm,
Eigen::MatrixXd & Qt,
GridBase * grid )
inlineprivate

Definition at line 694 of file ImplicitlyRestartedLanczos.h.

References GridLogError, Nk, Nm, and QR_decomp().

Referenced by diagonalize().

Field Documentation

◆ small

template<class Field>
const RealD ImplicitlyRestartedLanczos< Field >::small
private

Definition at line 99 of file ImplicitlyRestartedLanczos.h.

◆ MaxIter

template<class Field>
int ImplicitlyRestartedLanczos< Field >::MaxIter
private

Definition at line 100 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ MinRestart

template<class Field>
int ImplicitlyRestartedLanczos< Field >::MinRestart
private

Definition at line 101 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ Nstop

template<class Field>
int ImplicitlyRestartedLanczos< Field >::Nstop
private

Definition at line 102 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ Nk

template<class Field>
int ImplicitlyRestartedLanczos< Field >::Nk
private

Definition at line 103 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), diagonalize(), diagonalize_Eigen(), diagonalize_lapack(), diagonalize_QR(), ImplicitlyRestartedLanczos(), ImplicitlyRestartedLanczos(), and QR_decomp().

◆ Nm

template<class Field>
int ImplicitlyRestartedLanczos< Field >::Nm
private

Definition at line 105 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), diagonalize(), diagonalize_Eigen(), diagonalize_lapack(), diagonalize_QR(), ImplicitlyRestartedLanczos(), ImplicitlyRestartedLanczos(), QR_decomp(), and step().

◆ diagonalisation

template<class Field>
IRLdiagonalisation ImplicitlyRestartedLanczos< Field >::diagonalisation
private

Definition at line 106 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), diagonalize(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ orth_period

template<class Field>
int ImplicitlyRestartedLanczos< Field >::orth_period
private

Definition at line 107 of file ImplicitlyRestartedLanczos.h.

Referenced by ImplicitlyRestartedLanczos(), ImplicitlyRestartedLanczos(), and step().

◆ OrthoTime

template<class Field>
RealD ImplicitlyRestartedLanczos< Field >::OrthoTime
private

Definition at line 109 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), and orthogonalize().

◆ eresid

template<class Field>
RealD ImplicitlyRestartedLanczos< Field >::eresid
private

Definition at line 110 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ betastp

template<class Field>
RealD ImplicitlyRestartedLanczos< Field >::betastp
private

Definition at line 110 of file ImplicitlyRestartedLanczos.h.

Referenced by ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ _PolyOp

template<class Field>
LinearFunction<Field>& ImplicitlyRestartedLanczos< Field >::_PolyOp
private

Definition at line 114 of file ImplicitlyRestartedLanczos.h.

Referenced by ImplicitlyRestartedLanczos(), ImplicitlyRestartedLanczos(), and step().

◆ _HermOp

template<class Field>
LinearFunction<Field>& ImplicitlyRestartedLanczos< Field >::_HermOp
private

Definition at line 115 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ _Tester

template<class Field>
ImplicitlyRestartedLanczosTester<Field>& ImplicitlyRestartedLanczos< Field >::_Tester
private

Definition at line 116 of file ImplicitlyRestartedLanczos.h.

Referenced by calc(), ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

◆ SimpleTester

template<class Field>
ImplicitlyRestartedLanczosHermOpTester<Field> ImplicitlyRestartedLanczos< Field >::SimpleTester
private

Definition at line 118 of file ImplicitlyRestartedLanczos.h.

Referenced by ImplicitlyRestartedLanczos(), and ImplicitlyRestartedLanczos().

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