#include <ImplicitlyRestartedBlockLanczos.h>

Public Member Functions
	ImplicitlyRestartedBlockLanczos (LinearOperatorBase< Field > &Linop, LinearOperatorBase< Field > &SLinop, GridRedBlackCartesian FrbGrid, GridRedBlackCartesian SFrbGrid, int _mrhs, OperatorFunction< Field > &poly, int _Nstop, int _Nconv_test_interval, int _Nu, int _Nk, int _Nm, RealD _eresid, int _MaxIter, IRBLdiagonalisation _diagonalisation=IRBLdiagonaliseWithEigen)
void	orthogonalize (Field &w, std::vector< Field > &evec, int k, int if_print=0)
void	reorthogonalize (Field &w, std::vector< Field > &evec, int k)
void	orthogonalize (std::vector< Field > &w, int _Nu, std::vector< Field > &evec, int k, int if_print=0)
void	orthogonalize_blas (std::vector< Field > &w, std::vector< Field > &evec, int R, int do_print=0)
void	orthogonalize_blockhead (Field &w, std::vector< Field > &evec, int k, int Nu)
void	calc (std::vector< RealD > &eval, std::vector< Field > &evec, const std::vector< Field > &src, int &Nconv, LanczosType Impl)
void	calc_irbl (std::vector< RealD > &eval, std::vector< Field > &evec, const std::vector< Field > &src, int &Nconv)
void	calc_rbl (std::vector< RealD > &eval, std::vector< Field > &evec, const std::vector< Field > &src, int &Nconv)

Static Public Member Functions
static RealD	normalize (Field &v, int if_print=0)

Data Fields
int	split_test

Private Member Functions
void	blockwiseStep (std::vector< std::vector< ComplexD > > &lmd, std::vector< std::vector< ComplexD > > &lme, std::vector< Field > &evec, std::vector< Field > &w, std::vector< Field > &w_copy, int b)
void	diagonalize_Eigen (std::vector< RealD > &eval, std::vector< std::vector< ComplexD > > &lmd, std::vector< std::vector< ComplexD > > &lme, int Nu, int Nk, int Nm, Eigen::MatrixXcd &Qt, GridBase *grid)
void	diagonalize (std::vector< RealD > &eval, std::vector< std::vector< ComplexD > > &lmd, std::vector< std::vector< ComplexD > > &lme, int Nu, int Nk, int Nm, Eigen::MatrixXcd &Qt, GridBase *grid)
void	unpackHermitBlockTriDiagMatToEigen (std::vector< std::vector< ComplexD > > &lmd, std::vector< std::vector< ComplexD > > &lme, int Nu, int Nb, int Nk, int Nm, Eigen::MatrixXcd &M)
void	packHermitBlockTriDiagMatfromEigen (std::vector< std::vector< ComplexD > > &lmd, std::vector< std::vector< ComplexD > > &lme, int Nu, int Nb, int Nk, int Nm, Eigen::MatrixXcd &M)
void	shiftedQRDecompEigen (Eigen::MatrixXcd &M, int Nu, int Nm, RealD Dsh, Eigen::MatrixXcd &Qprod)
void	exampleQRDecompEigen (void)

Private Attributes
std::string	cname
int	MaxIter
int	Nstop
int	Nu
int	Nk
int	Nm
int	Nblock_k
int	Nblock_m
int	Nconv_test_interval
RealD	eresid
IRBLdiagonalisation	diagonalisation
SortEigen< Field >	_sort
LinearOperatorBase< Field > &	_Linop
LinearOperatorBase< Field > &	_SLinop
OperatorFunction< Field > &	_poly
GridRedBlackCartesian *	f_grid
GridRedBlackCartesian *	sf_grid
int	mrhs
int	Nevec_acc

Detailed Description

template<class Field>
class Grid::ImplicitlyRestartedBlockLanczos< Field >

Definition at line 117 of file ImplicitlyRestartedBlockLanczos.h.

Constructor & Destructor Documentation

◆ ImplicitlyRestartedBlockLanczos()

template<class Field>

Grid::ImplicitlyRestartedBlockLanczos< Field >::ImplicitlyRestartedBlockLanczos	(	LinearOperatorBase< Field > &	Linop,
		LinearOperatorBase< Field > &	SLinop,
		GridRedBlackCartesian *	FrbGrid,
		GridRedBlackCartesian *	SFrbGrid,
		int	_mrhs,
		OperatorFunction< Field > &	poly,
		int	_Nstop,
		int	_Nconv_test_interval,
		int	_Nu,
		int	_Nk,
		int	_Nm,
		RealD	_eresid,
		int	_MaxIter,
		IRBLdiagonalisation	_diagonalisation = IRBLdiagonaliseWithEigen )

inline

Definition at line 156 of file ImplicitlyRestartedBlockLanczos.h.

References _Linop, _poly, _SLinop, diagonalisation, eresid, f_grid, Grid::IRBLdiagonaliseWithEigen, MaxIter, mrhs, Nblock_k, Nblock_m, Nconv_test_interval, Nevec_acc, Nk, Nm, Nstop, Nu, sf_grid, and split_test.

Member Function Documentation

◆ normalize()

template<class Field>

RealD Grid::ImplicitlyRestartedBlockLanczos< Field >::normalize	(	Field &	v,
		int	if_print = 0 )

inlinestatic

Definition at line 183 of file ImplicitlyRestartedBlockLanczos.h.

References norm2(), and sqrt().

Referenced by orthogonalize(), orthogonalize(), orthogonalize_blas(), and orthogonalize_blockhead().

◆ orthogonalize() [1/2]

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::orthogonalize	(	Field &	w,
		std::vector< Field > &	evec,
		int	k,
		int	if_print = 0 )

inline

Definition at line 191 of file ImplicitlyRestartedBlockLanczos.h.

References Glog, innerProduct(), norm2(), and normalize().

Referenced by blockwiseStep(), calc_irbl(), calc_rbl(), and reorthogonalize().

◆ reorthogonalize()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::reorthogonalize	(	Field &	w,
		std::vector< Field > &	evec,
		int	k )

inline

Definition at line 211 of file ImplicitlyRestartedBlockLanczos.h.

References orthogonalize().

◆ orthogonalize() [2/2]

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::orthogonalize	(	std::vector< Field > &	w,
		int	_Nu,
		std::vector< Field > &	evec,
		int	k,
		int	if_print = 0 )

inline

Definition at line 216 of file ImplicitlyRestartedBlockLanczos.h.

References innerProduct(), and normalize().

◆ orthogonalize_blas()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::orthogonalize_blas	(	std::vector< Field > &	w,
		std::vector< Field > &	evec,
		int	R,
		int	do_print = 0 )

inline

Definition at line 232 of file ImplicitlyRestartedBlockLanczos.h.

References AcceleratorWrite, autoView, CUDA_COMPLEX, CUDA_FLOAT, CUDA_GEMM, CartesianCommunicator::GlobalSumVector(), Glog, GridBase::lSites(), MAKE_CUDA_COMPLEX, Nevec_acc, normalize(), and Nu.

Referenced by blockwiseStep().

◆ orthogonalize_blockhead()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::orthogonalize_blockhead	(	Field &	w,
		std::vector< Field > &	evec,
		int	k,
		int	Nu )

inline

Definition at line 321 of file ImplicitlyRestartedBlockLanczos.h.

References innerProduct(), normalize(), and Nu.

◆ calc()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::calc	(	std::vector< RealD > &	eval,
		std::vector< Field > &	evec,
		const std::vector< Field > &	src,
		int &	Nconv,
		LanczosType	Impl )

inline

Definition at line 333 of file ImplicitlyRestartedBlockLanczos.h.

References calc_irbl(), calc_rbl(), CUDA_COMPLEX, eval(), Grid::irbl, GridBase::lSites(), Nevec_acc, Nu, Grid::rbl, and GridBase::show_decomposition().

◆ calc_irbl()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::calc_irbl	(	std::vector< RealD > &	eval,
		std::vector< Field > &	evec,
		const std::vector< Field > &	src,
		int &	Nconv )

inline

Definition at line 371 of file ImplicitlyRestartedBlockLanczos.h.

References _Linop, _sort, B, blockwiseStep(), cname, diagonalisation, diagonalize(), eresid, eval(), Glog, innerProduct(), Grid::IRBLdiagonaliseWithEigen, MaxIter, Nblock_k, Nblock_m, Nk, Nm, norm2(), Nstop, Nu, orthogonalize(), packHermitBlockTriDiagMatfromEigen(), real(), shiftedQRDecompEigen(), and unpackHermitBlockTriDiagMatToEigen().

Referenced by calc().

◆ calc_rbl()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::calc_rbl	(	std::vector< RealD > &	eval,
		std::vector< Field > &	evec,
		const std::vector< Field > &	src,
		int &	Nconv )

inline

Definition at line 576 of file ImplicitlyRestartedBlockLanczos.h.

References _Linop, _sort, AdviseInfrequentUse, B, blockwiseStep(), cname, diagonalisation, diagonalize(), eresid, eval(), Glog, innerProduct(), Grid::IRBLdiagonaliseWithEigen, MaxIter, Nblock_k, Nconv_test_interval, Nk, Nm, norm2(), Nstop, Nu, orthogonalize(), and real().

Referenced by calc().

◆ blockwiseStep()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::blockwiseStep	(	std::vector< std::vector< ComplexD > > &	lmd,
		std::vector< std::vector< ComplexD > > &	lme,
		std::vector< Field > &	evec,
		std::vector< Field > &	w,
		std::vector< Field > &	w_copy,
		int	b )

inlineprivate

Definition at line 777 of file ImplicitlyRestartedBlockLanczos.h.

References _Linop, _poly, _SLinop, conjugate(), f_grid, Glog, Grid_split(), Grid_unsplit(), innerProduct(), mrhs, Nm, norm2(), Nu, orthogonalize(), orthogonalize_blas(), real(), sf_grid, and split_test.

Referenced by calc_irbl(), and calc_rbl().

◆ diagonalize_Eigen()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::diagonalize_Eigen	(	std::vector< RealD > &	eval,
		std::vector< std::vector< ComplexD > > &	lmd,
		std::vector< std::vector< ComplexD > > &	lme,
		int	Nu,
		int	Nk,
		int	Nm,
		Eigen::MatrixXcd &	Qt,
		GridBase *	grid )

inlineprivate

Definition at line 925 of file ImplicitlyRestartedBlockLanczos.h.

References conjugate(), eval(), Nk, Nm, and Nu.

Referenced by diagonalize().

◆ diagonalize()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::diagonalize	(	std::vector< RealD > &	eval,
		std::vector< std::vector< ComplexD > > &	lmd,
		std::vector< std::vector< ComplexD > > &	lme,
		int	Nu,
		int	Nk,
		int	Nm,
		Eigen::MatrixXcd &	Qt,
		GridBase *	grid )

inlineprivate

Definition at line 1107 of file ImplicitlyRestartedBlockLanczos.h.

References diagonalisation, diagonalize_Eigen(), eval(), Grid::IRBLdiagonaliseWithEigen, Nk, Nm, and Nu.

Referenced by calc_irbl(), and calc_rbl().

◆ unpackHermitBlockTriDiagMatToEigen()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::unpackHermitBlockTriDiagMatToEigen	(	std::vector< std::vector< ComplexD > > &	lmd,
		std::vector< std::vector< ComplexD > > &	lme,
		int	Nu,
		int	Nb,
		int	Nk,
		int	Nm,
		Eigen::MatrixXcd &	M )

inlineprivate

Definition at line 1127 of file ImplicitlyRestartedBlockLanczos.h.

References conjugate(), Nk, Nm, and Nu.

Referenced by calc_irbl().

◆ packHermitBlockTriDiagMatfromEigen()

template<class Field>

void Grid::ImplicitlyRestartedBlockLanczos< Field >::packHermitBlockTriDiagMatfromEigen	(	std::vector< std::vector< ComplexD > > &	lmd,
		std::vector< std::vector< ComplexD > > &	lme,
		int	Nu,
		int	Nb,
		int	Nk,
		int	Nm,
		Eigen::MatrixXcd &	M )