Projects per year
Abstract / Description of output
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.
Original language | English |
---|---|
Pages (from-to) | 689–697 |
Journal | Computer Physics Communications |
Volume | 184 |
Issue number | 3 |
DOIs | |
Publication status | Published - 6 Nov 2012 |
Fingerprint
Dive into the research topics of 'Numerical determination of partial spectrum of Hermitian matrices using a Lanczos method with selective reorthogonalization'. Together they form a unique fingerprint.Projects
- 2 Finished
-
Science and Innovation: Numerical Algorithms and Intelligent Software for the Evolving HPC Platform
1/08/09 → 31/07/14
Project: Research
-
The Standard Model and beyond
Kenway, R., Ball, R., Berera, A., Binoth, T., Boyle, P., Del Debbio, L., Kennedy, A., Pendleton, B. & Plehn, T.
1/10/08 → 30/09/11
Project: Research