#include <eigensolve_quda.h>

Inheritance diagram for quda::EigenSolver:

Collaboration diagram for quda::EigenSolver:

[legend]

Public Member Functions
	EigenSolver (QudaEigParam *eig_param, TimeProfile &profile)
	Constructor for base Eigensolver class. More...

virtual	~EigenSolver ()

virtual void	operator() (std::vector< ColorSpinorField *> &kSpace, std::vector< Complex > &evals)=0
	Computes the eigen decomposition for the operator passed to create. More...

void	matVec (const DiracMatrix &mat, ColorSpinorField &out, const ColorSpinorField &in)
	Applies the specified matVec operation: M, Mdag, MMdag, MdagM. More...

void	chebyOp (const DiracMatrix &mat, ColorSpinorField &out, const ColorSpinorField &in)
	Promoted the specified matVec operation: M, Mdag, MMdag, MdagM to a Chebyshev polynomial. More...

Complex	blockOrthogonalize (std::vector< ColorSpinorField > v, std::vector< ColorSpinorField > r, int j)
	Orthogonalise input vector r against vector space v using block-BLAS. More...

void	deflate (std::vector< ColorSpinorField > vec_defl, std::vector< ColorSpinorField > vec, std::vector< ColorSpinorField *> evecs, std::vector< Complex > evals)
	Deflate vector with Eigenvectors. More...

void	deflateSVD (std::vector< ColorSpinorField > vec_defl, std::vector< ColorSpinorField > vec, std::vector< ColorSpinorField *> evecs, std::vector< Complex > evals)
	Deflate vector with both left and Right singular vectors. More...

void	computeSVD (const DiracMatrix &mat, std::vector< ColorSpinorField *> &evecs, std::vector< Complex > &evals)
	Computes Left/Right SVD from pre computed Right/Left. More...

void	computeEvals (const DiracMatrix &mat, std::vector< ColorSpinorField *> &evecs, std::vector< Complex > &evals, int k)
	Compute eigenvalues and their residiua. More...

void	loadFromFile (const DiracMatrix &mat, std::vector< ColorSpinorField *> &eig_vecs, std::vector< Complex > &evals)
	Load and check eigenpairs from file. More...

Static Public Member Functions
static EigenSolver *	create (QudaEigParam *eig_param, const DiracMatrix &mat, TimeProfile &profile)
	Creates the eigensolver using the parameters given and the matrix. More...

static void	loadVectors (std::vector< ColorSpinorField *> &eig_vecs, std::string file)
	Load vectors from file. More...

static void	saveVectors (const std::vector< ColorSpinorField *> &eig_vecs, std::string file)
	Save vectors to file. More...

Protected Attributes
QudaEigParam *	eig_param

TimeProfile &	profile

int	nEv

int	nKr

int	nConv

double	tol

bool	reverse

char	spectrum [3]

bool	converged

int	restart_iter

int	max_restarts

int	check_interval

int	iter

int	iter_converged

int	iter_locked

int	iter_keep

int	num_converged

int	num_locked

int	num_keep

double *	residua

std::vector< ColorSpinorField * >	r

std::vector< ColorSpinorField * >	d_vecs_tmp

ColorSpinorField *	tmp1

ColorSpinorField *	tmp2

Complex *	Qmat

Detailed Description

Definition at line 11 of file eigensolve_quda.h.

Constructor & Destructor Documentation

◆ EigenSolver()

quda::EigenSolver::EigenSolver	(	QudaEigParam *	eig_param,
		TimeProfile &	profile
	)

Constructor for base Eigensolver class.

Parameters

eig_param	MGParam struct that defines all meta data
profile	Timeprofile instance used to profile

Definition at line 27 of file eigensolve_quda.cpp.

References QudaEigParam_s::a_max, QudaEigParam_s::a_min, check_interval, QudaEigParam_s::check_interval, converged, errorQuda, getVerbosity(), iter, iter_converged, iter_keep, iter_locked, max_restarts, QudaEigParam_s::max_restarts, nConv, QudaEigParam_s::nConv, nEv, QudaEigParam_s::nEv, nKr, QudaEigParam_s::nKr, num_converged, num_keep, num_locked, QudaEigParam_s::poly_deg, printfQuda, Qmat, quda::QUDA_PROFILE_INIT, QUDA_SPECTRUM_LI_EIG, QUDA_SPECTRUM_LM_EIG, QUDA_SPECTRUM_LR_EIG, QUDA_SPECTRUM_SI_EIG, QUDA_SPECTRUM_SM_EIG, QUDA_SPECTRUM_SR_EIG, QUDA_VERBOSE, residua, restart_iter, reverse, safe_malloc, spectrum, QudaEigParam_s::spectrum, tol, QudaEigParam_s::tol, and QudaEigParam_s::use_poly_acc.

Here is the call graph for this function:

◆ ~EigenSolver()

quda::EigenSolver::~EigenSolver ( )

virtual

Destructor for EigenSolver class.

Definition at line 490 of file eigensolve_quda.cpp.

References host_free, Qmat, residua, tmp1, and tmp2.

Member Function Documentation

◆ blockOrthogonalize()

Complex quda::EigenSolver::blockOrthogonalize	(	std::vector< ColorSpinorField *>	v,
		std::vector< ColorSpinorField *>	r,
		int	j
	)

Orthogonalise input vector r against vector space v using block-BLAS.

Parameters

[out]	Sum	of inner products
[in]	v	Vector space
[in]	r	Vector to be orthogonalised
[in]	j	Number of vectors in v to orthogonalise against

Definition at line 207 of file eigensolve_quda.cpp.

References quda::blas::caxpy(), quda::blas::cDotProduct(), host_free, quda::s, safe_malloc, and quda::sum().

Referenced by quda::TRLM::lanczosStep().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ chebyOp()

void quda::EigenSolver::chebyOp	(	const DiracMatrix &	mat,
		ColorSpinorField &	out,
		const ColorSpinorField &	in
	)

Promoted the specified matVec operation: M, Mdag, MMdag, MdagM to a Chebyshev polynomial.

Parameters

[in]	mat	Matrix operator
[in]	out	Output spinor
[in]	in	Input spinor

Definition at line 138 of file eigensolve_quda.cpp.

References QudaEigParam_s::a_max, QudaEigParam_s::a_min, quda::blas::caxpby(), quda::blas::caxpbypczw(), quda::blas::copy(), quda::ColorSpinorField::Create(), eig_param, errorQuda, matVec(), QudaEigParam_s::poly_deg, swap(), tmp1, tmp2, and QudaEigParam_s::use_poly_acc.

Referenced by quda::TRLM::lanczosStep().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeEvals()

void quda::EigenSolver::computeEvals	(	const DiracMatrix &	mat,
		std::vector< ColorSpinorField *> &	evecs,
		std::vector< Complex > &	evals,
		int	k
	)

Compute eigenvalues and their residiua.

Parameters

[in]	mat	Matrix operator
[in]	evecs	The eigenvectors
[in]	evals	The eigenvalues
[in]	k	The number to compute

Definition at line 345 of file eigensolve_quda.cpp.

References quda::blas::caxpby(), quda::blas::cDotProduct(), matVec(), quda::blas::norm2(), r, residua, quda::size, and quda::sqrt().

Referenced by loadFromFile(), and quda::TRLM::operator()().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeSVD()

void quda::EigenSolver::computeSVD	(	const DiracMatrix &	mat,
		std::vector< ColorSpinorField *> &	evecs,
		std::vector< Complex > &	evals
	)

Computes Left/Right SVD from pre computed Right/Left.

Parameters

[in]	mat	Matrix operator
[in]	evecs	Computed eigenvectors of NormOp
[in]	evals	Computed eigenvalues of NormOp

Definition at line 263 of file eigensolve_quda.cpp.

References quda::abs(), quda::blas::ax(), quda::blas::cDotProduct(), eig_param, errorQuda, quda::DiracMatrix::Expose(), getVerbosity(), quda::Dirac::M(), nConv, QudaEigParam_s::nConv, quda::norm(), quda::blas::norm2(), printfQuda, QUDA_SUMMARIZE, and quda::sqrt().

Referenced by quda::Solver::constructDeflationSpace().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ create()

EigenSolver * quda::EigenSolver::create	(	QudaEigParam *	eig_param,
		const DiracMatrix &	mat,
		TimeProfile &	profile
	)

static

Creates the eigensolver using the parameters given and the matrix.

Parameters

eig_param	The eigensolver parameters
mat	The operator to solve
profile	Time Profile

Definition at line 109 of file eigensolve_quda.cpp.

References QudaEigParam_s::eig_type, errorQuda, getVerbosity(), printfQuda, QUDA_EIG_IR_ARNOLDI, QUDA_EIG_IR_LANCZOS, QUDA_EIG_TR_LANCZOS, and QUDA_SUMMARIZE.

Referenced by quda::Solver::constructDeflationSpace(), eigensolveQuda(), and quda::MG::generateEigenVectors().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ deflate()

void quda::EigenSolver::deflate	(	std::vector< ColorSpinorField *>	vec_defl,
		std::vector< ColorSpinorField *>	vec,
		std::vector< ColorSpinorField *>	evecs,
		std::vector< Complex >	evals
	)

Deflate vector with Eigenvectors.

Parameters

[in]	vec_defl	The deflated vector
[in]	vec	The input vector
[in]	evecs	The eigenvectors to use in deflation
[in]	evals	The eigenvalues to use in deflation

Definition at line 229 of file eigensolve_quda.cpp.

References quda::blas::caxpy(), quda::blas::cDotProduct(), eig_param, getVerbosity(), host_free, QudaEigParam_s::nConv, printfQuda, QUDA_VERBOSE, quda::s, safe_malloc, and quda::blas::zero().

Referenced by quda::CG::operator()(), and quda::CACG::operator()().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ deflateSVD()

void quda::EigenSolver::deflateSVD	(	std::vector< ColorSpinorField *>	vec_defl,
		std::vector< ColorSpinorField *>	vec,
		std::vector< ColorSpinorField *>	evecs,
		std::vector< Complex >	evals
	)

Deflate vector with both left and Right singular vectors.

Parameters

[in]	vec_defl	The deflated vector
[in]	vec	The input vector
[in]	evecs	The singular vectors to use in deflation
[in]	evals	The singular values to use in deflation