Numerical determination of partial spectrum of Hermitian matrices using a Lanczos method with selective reorthogonalization
High Energy Physics - Lattice
2015-06-12 v1 Numerical Analysis
Abstract
We introduce a new algorithm for finding the eigenvalues and eigenvectors of Hermitian matrices within a specified region, based upon the LANSO algorithm of Parlett and Scott. It uses selective reorthogonalization to avoid the duplication of eigenpairs in finite-precision arithmetic, but uses a new bound to decide when such reorthogonalization is required, and only reorthogonalizes with respect to eigenpairs within the region of interest. We investigate its performance for the Hermitian Wilson--Dirac operator (\gamma_5D) in lattice quantum chromodynamics, and compare it with previous methods.
Keywords
Cite
@article{arxiv.1211.1180,
title = {Numerical determination of partial spectrum of Hermitian matrices using a Lanczos method with selective reorthogonalization},
author = {Chris Johnson and A. D. Kennedy},
journal= {arXiv preprint arXiv:1211.1180},
year = {2015}
}