A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix [A projected preconditioned conjugate gradient algorithm for computing a large eigenspace of a Hermitian matrix]
Abstract
Here, we present an iterative algorithm for computing an invariant subspace associated with the algebraically smallest eigenvalues of a large sparse or structured Hermitian matrix A. We are interested in the case in which the dimension of the invariant subspace is large (e.g., over several hundreds or thousands) even though it may still be small relative to the dimension of A. These problems arise from, for example, density functional theory (DFT) based electronic structure calculations for complex materials. The key feature of our algorithm is that it performs fewer Rayleigh–Ritz calculations compared to existing algorithms such as the locally optimal block preconditioned conjugate gradient or the Davidson algorithm. It is a block algorithm, and hence can take advantage of efficient BLAS3 operations and be implemented with multiple levels of concurrency. We discuss a number of practical issues that must be addressed in order to implement the algorithm efficiently on a high performance computer.
- Authors:
-
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Publication Date:
- Research Org.:
- Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Basic Energy Sciences (BES)
- OSTI Identifier:
- 1325863
- Alternate Identifier(s):
- OSTI ID: 1556219
- Report Number(s):
- LLNL-JRNL-695289
Journal ID: ISSN 0021-9991
- Grant/Contract Number:
- AC52-07NA27344; AC02-05CH11231
- Resource Type:
- Journal Article: Accepted Manuscript
- Journal Name:
- Journal of Computational Physics
- Additional Journal Information:
- Journal Volume: 290; Journal Issue: C; Journal ID: ISSN 0021-9991
- Publisher:
- Elsevier
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; 97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; preconditioned eigenvalue solvers; density functional theory based electronic structure calculations
Citation Formats
Vecharynski, Eugene, Yang, Chao, and Pask, John E. A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix [A projected preconditioned conjugate gradient algorithm for computing a large eigenspace of a Hermitian matrix]. United States: N. p., 2015.
Web. doi:10.1016/j.jcp.2015.02.030.
Vecharynski, Eugene, Yang, Chao, & Pask, John E. A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix [A projected preconditioned conjugate gradient algorithm for computing a large eigenspace of a Hermitian matrix]. United States. https://doi.org/10.1016/j.jcp.2015.02.030
Vecharynski, Eugene, Yang, Chao, and Pask, John E. 2015.
"A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix [A projected preconditioned conjugate gradient algorithm for computing a large eigenspace of a Hermitian matrix]". United States. https://doi.org/10.1016/j.jcp.2015.02.030. https://www.osti.gov/servlets/purl/1325863.
@article{osti_1325863,
title = {A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix [A projected preconditioned conjugate gradient algorithm for computing a large eigenspace of a Hermitian matrix]},
author = {Vecharynski, Eugene and Yang, Chao and Pask, John E.},
abstractNote = {Here, we present an iterative algorithm for computing an invariant subspace associated with the algebraically smallest eigenvalues of a large sparse or structured Hermitian matrix A. We are interested in the case in which the dimension of the invariant subspace is large (e.g., over several hundreds or thousands) even though it may still be small relative to the dimension of A. These problems arise from, for example, density functional theory (DFT) based electronic structure calculations for complex materials. The key feature of our algorithm is that it performs fewer Rayleigh–Ritz calculations compared to existing algorithms such as the locally optimal block preconditioned conjugate gradient or the Davidson algorithm. It is a block algorithm, and hence can take advantage of efficient BLAS3 operations and be implemented with multiple levels of concurrency. We discuss a number of practical issues that must be addressed in order to implement the algorithm efficiently on a high performance computer.},
doi = {10.1016/j.jcp.2015.02.030},
url = {https://www.osti.gov/biblio/1325863},
journal = {Journal of Computational Physics},
issn = {0021-9991},
number = C,
volume = 290,
place = {United States},
year = {Wed Feb 25 00:00:00 EST 2015},
month = {Wed Feb 25 00:00:00 EST 2015}
}
Web of Science
Works referenced in this record:
Large-scale ab initio calculations based on three levels of parallelization
journal, April 2008
- Bottin, François; Leroux, Stéphane; Knyazev, Andrew
- Computational Materials Science, Vol. 42, Issue 2
Elemental: A New Framework for Distributed Memory Dense Matrix Computations
journal, February 2013
- Poulson, Jack; Marker, Bryan; van de Geijn, Robert A.
- ACM Transactions on Mathematical Software, Vol. 39, Issue 2
Architecture of Qbox: A scalable first-principles molecular dynamics code
journal, January 2008
- Gygi, F.
- IBM Journal of Research and Development, Vol. 52, Issue 1.2
Some metric inequalities in the space of matrices
journal, January 1955
- Fan, Ky; Hoffman, A. J.
- Proceedings of the American Mathematical Society, Vol. 6, Issue 1
QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials
journal, September 2009
- Giannozzi, Paolo; Baroni, Stefano; Bonini, Nicola
- Journal of Physics: Condensed Matter, Vol. 21, Issue 39, Article No. 395502
A Trace Minimization Algorithm for the Generalized Eigenvalue Problem
journal, December 1982
- Sameh, Ahmed H.; Wisniewski, John A.
- SIAM Journal on Numerical Analysis, Vol. 19, Issue 6
Inhomogeneous Electron Gas
journal, November 1964
- Hohenberg, P.; Kohn, W.
- Physical Review, Vol. 136, Issue 3B, p. B864-B871
The Gradient Projection Method for Nonlinear Programming. Part II. Nonlinear Constraints
journal, December 1961
- Rosen, J. B.
- Journal of the Society for Industrial and Applied Mathematics, Vol. 9, Issue 4
A spectrum slicing method for the Kohn–Sham problem
journal, March 2012
- Schofield, Grady; Chelikowsky, James R.; Saad, Yousef
- Computer Physics Communications, Vol. 183, Issue 3
The Gradient Projection Method for Nonlinear Programming. Part I. Linear Constraints
journal, March 1960
- Rosen, J. B.
- Journal of the Society for Industrial and Applied Mathematics, Vol. 8, Issue 1
ABINIT: First-principles approach to material and nanosystem properties
journal, December 2009
- Gonze, X.; Amadon, B.; Anglade, P. -M.
- Computer Physics Communications, Vol. 180, Issue 12
Simplified method for calculating the energy of weakly interacting fragments
journal, February 1985
- Harris, J.
- Physical Review B, Vol. 31, Issue 4
Block Locally Optimal Preconditioned Eigenvalue Xolvers (BLOPEX) in Hypre and PETSc
journal, January 2007
- Knyazev, A. V.; Argentati, M. E.; Lashuk, I.
- SIAM Journal on Scientific Computing, Vol. 29, Issue 5
Tight-binding models and density-functional theory
journal, June 1989
- Foulkes, W. Matthew C.; Haydock, Roger
- Physical Review B, Vol. 39, Issue 17
Convex programming in Hilbert space
journal, January 1964
- Goldstein, A. A.
- Bulletin of the American Mathematical Society, Vol. 70, Issue 5
The ELPA library: scalable parallel eigenvalue solutions for electronic structure theory and computational science
journal, May 2014
- Marek, A.; Blum, V.; Johanni, R.
- Journal of Physics: Condensed Matter, Vol. 26, Issue 21
Toward the Optimal Preconditioned Eigensolver: Locally Optimal Block Preconditioned Conjugate Gradient Method
journal, January 2001
- Knyazev, Andrew V.
- SIAM Journal on Scientific Computing, Vol. 23, Issue 2
The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices
journal, January 1975
- Davidson, Ernest R.
- Journal of Computational Physics, Vol. 17, Issue 1
A Trace Minimization Algorithm for the Generalized Eigenvalue Problem
journal, December 1982
- Sameh, Ahmed H.; Wisniewski, John A.
- SIAM Journal on Numerical Analysis, Vol. 19, Issue 6
A Simultaneous Iteration Algorithm for Real Matrices
journal, June 1981
- Stewart, William J.; Jennings, Alan
- ACM Transactions on Mathematical Software, Vol. 7, Issue 2
The trace minimization method for the symmetric generalized eigenvalue problem
journal, November 2000
- Sameh, Ahmed; Tong, Zhanye
- Journal of Computational and Applied Mathematics, Vol. 123, Issue 1-2
Self-Consistent Equations Including Exchange and Correlation Effects
journal, November 1965
- Kohn, W.; Sham, L. J.
- Physical Review, Vol. 140, Issue 4A, p. A1133-A1138
A brief introduction to the ABINIT software package
journal, January 2005
- Gonze, Xavier
- Zeitschrift für Kristallographie - Crystalline Materials, Vol. 220, Issue 5/6
Works referencing / citing this record:
Prediction of atomization energy using graph kernel and active learning
journal, January 2019
- Tang, Yu-Hang; de Jong, Wibe A.
- The Journal of Chemical Physics, Vol. 150, Issue 4
Convergence theory for preconditioned eigenvalue solvers in a nutshell
text, January 2014
- Argentati, Merico E.; Knyazev, Andrew V.; Neymeyr, Klaus
- arXiv