A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems
Abstract
Polynomial filtering can provide a highly effective means of computing all eigenvalues of a real symmetric (or complex Hermitian) matrix that are located in a given interval, anywhere in the spectrum. This paper describes a technique for tackling this problem by combining a thick-restart version of the Lanczos algorithm with deflation ("locking'') and a new type of polynomial filter obtained from a least-squares technique. Furthermore, the resulting algorithm can be utilized in a “spectrum-slicing” approach whereby a very large number of eigenvalues and associated eigenvectors of the matrix are computed by extracting eigenpairs located in different subintervals independently from one another.
- Authors:
-
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Univ. of Minnesota, Twin Cities, MN (United States)
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Publication Date:
- Research Org.:
- Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
- Sponsoring Org.:
- USDOE
- OSTI Identifier:
- 1438696
- Report Number(s):
- LLNL-JRNL-685856
Journal ID: ISSN 1064-8275
- Grant/Contract Number:
- AC52-07NA27344
- Resource Type:
- Accepted Manuscript
- Journal Name:
- SIAM Journal on Scientific Computing
- Additional Journal Information:
- Journal Volume: 38; Journal Issue: 4; Journal ID: ISSN 1064-8275
- Publisher:
- SIAM
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; Lanczos algorithm; polynomial filtering; Thick-Restart; deflation; spectrum slicing; interior eigenvalue problems
Citation Formats
Li, Ruipeng, Xi, Yuanzhe, Vecharynski, Eugene, Yang, Chao, and Saad, Yousef. A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems. United States: N. p., 2016.
Web. doi:10.1137/15M1054493.
Li, Ruipeng, Xi, Yuanzhe, Vecharynski, Eugene, Yang, Chao, & Saad, Yousef. A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems. United States. https://doi.org/10.1137/15M1054493
Li, Ruipeng, Xi, Yuanzhe, Vecharynski, Eugene, Yang, Chao, and Saad, Yousef. Tue .
"A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems". United States. https://doi.org/10.1137/15M1054493. https://www.osti.gov/servlets/purl/1438696.
@article{osti_1438696,
title = {A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems},
author = {Li, Ruipeng and Xi, Yuanzhe and Vecharynski, Eugene and Yang, Chao and Saad, Yousef},
abstractNote = {Polynomial filtering can provide a highly effective means of computing all eigenvalues of a real symmetric (or complex Hermitian) matrix that are located in a given interval, anywhere in the spectrum. This paper describes a technique for tackling this problem by combining a thick-restart version of the Lanczos algorithm with deflation ("locking'') and a new type of polynomial filter obtained from a least-squares technique. Furthermore, the resulting algorithm can be utilized in a “spectrum-slicing” approach whereby a very large number of eigenvalues and associated eigenvectors of the matrix are computed by extracting eigenpairs located in different subintervals independently from one another.},
doi = {10.1137/15M1054493},
journal = {SIAM Journal on Scientific Computing},
number = 4,
volume = 38,
place = {United States},
year = {Tue Aug 16 00:00:00 EDT 2016},
month = {Tue Aug 16 00:00:00 EDT 2016}
}
Web of Science
Works referenced in this record:
Algorithm 827: irbleigs: A MATLAB program for computing a few eigenpairs of a large sparse Hermitian matrix
journal, September 2003
- Baglama, J.; Calvetti, D.; Reichel, L.
- ACM Transactions on Mathematical Software, Vol. 29, Issue 3
Anasazi software for the numerical solution of large-scale eigenvalue problems
journal, July 2009
- Baker, C. G.; Hetmaniuk, U. L.; Lehoucq, R. B.
- ACM Transactions on Mathematical Software, Vol. 36, Issue 3
JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices
journal, December 2007
- Bollhöfer, Matthias; Notay, Yvan
- Computer Physics Communications, Vol. 177, Issue 12
Strategies for spectrum slicing based on restarted Lanczos methods
journal, March 2012
- Campos, Carmen; Roman, Jose E.
- Numerical Algorithms, Vol. 60, Issue 2
A Filtered Lanczos Procedure for Extreme and Interior Eigenvalue Problems
journal, January 2012
- Fang, Haw-ren; Saad, Yousef
- SIAM Journal on Scientific Computing, Vol. 34, Issue 4
Rounding error analysis of the classical Gram-Schmidt orthogonalization process
journal, May 2005
- Giraud, Luc; Langou, Julien; Rozložník, Miroslav
- Numerische Mathematik, Vol. 101, Issue 1
The loss of orthogonality in the Gram-Schmidt orthogonalization process
journal, October 2005
- Giraud, L.; Langou, J.; Rozloznik, M.
- Computers & Mathematics with Applications, Vol. 50, Issue 7
SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems
journal, September 2005
- Hernandez, Vicente; Roman, Jose E.; Vidal, Vicente
- ACM Transactions on Mathematical Software, Vol. 31, Issue 3
Electronic structure calculations for plane-wave codes without diagonalization
journal, April 1999
- Jay, Laurent O.; Kim, Hanchul; Saad, Yousef
- Computer Physics Communications, Vol. 118, Issue 1
PARSEC – the pseudopotential algorithm for real-space electronic structure calculations: recent advances and novel applications to nano-structures
journal, April 2006
- Kronik, Leeor; Makmal, Adi; Tiago, Murilo L.
- physica status solidi (b), Vol. 243, Issue 5
Deflation Techniques for an Implicitly Restarted Arnoldi Iteration
journal, October 1996
- Lehoucq, R. B.; Sorensen, D. C.
- SIAM Journal on Matrix Analysis and Applications, Vol. 17, Issue 4
Approximating Spectral Densities of Large Matrices
journal, January 2016
- Lin, Lin; Saad, Yousef; Yang, Chao
- SIAM Review, Vol. 58, Issue 1
Computational Variants of the Lanczos Method for the Eigenproblem
journal, January 1972
- Paige, C. C.
- IMA Journal of Applied Mathematics, Vol. 10, Issue 3
Error Analysis of the Lanczos Algorithm for Tridiagonalizing a Symmetric Matrix
journal, January 1976
- Paige, C. C.
- IMA Journal of Applied Mathematics, Vol. 18, Issue 3
Density-matrix-based algorithm for solving eigenvalue problems
journal, March 2009
- Polizzi, Eric
- Physical Review B, Vol. 79, Issue 11
Practical Use of Polynomial Preconditionings for the Conjugate Gradient Method
journal, October 1985
- Saad, Youcef
- SIAM Journal on Scientific and Statistical Computing, Vol. 6, Issue 4
Filtered Conjugate Residual‐type Algorithms with Applications
journal, January 2006
- Saad, Yousef
- SIAM Journal on Matrix Analysis and Applications, Vol. 28, Issue 3
Analysis of Subspace Iteration for Eigenvalue Problems with Evolving Matrices
journal, January 2016
- Saad, Yousef
- SIAM Journal on Matrix Analysis and Applications, Vol. 37, Issue 1
Numerical Methods for Electronic Structure Calculations of Materials
journal, January 2010
- Saad, Yousef; Chelikowsky, James R.; Shontz, Suzanne M.
- SIAM Review, Vol. 52, Issue 1
A projection method for generalized eigenvalue problems using numerical integration
journal, October 2003
- Sakurai, Tetsuya; Sugiura, Hiroshi
- Journal of Computational and Applied Mathematics, Vol. 159, Issue 1
The Lanczos algorithm with partial reorthogonalization
journal, January 1984
- Simon, Horst D.
- Mathematics of Computation, Vol. 42, Issue 165
PRIMME: preconditioned iterative multimethod eigensolver—methods and software description
journal, April 2010
- Stathopoulos, Andreas; McCombs, James R.
- ACM Transactions on Mathematical Software, Vol. 37, Issue 2
Dynamic Thick Restarting of the Davidson, and the Implicitly Restarted Arnoldi Methods
journal, January 1998
- Stathopoulos, Andreas; Saad, Yousef; Wu, Kesheng
- SIAM Journal on Scientific Computing, Vol. 19, Issue 1
A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix
journal, June 2015
- Vecharynski, Eugene; Yang, Chao; Pask, John E.
- Journal of Computational Physics, Vol. 290
Adaptive Projection Subspace Dimension for the Thick-Restart Lanczos Method
journal, September 2010
- Yamazaki, Ichitaro; Bai, Zhaojun; Simon, Horst
- ACM Transactions on Mathematical Software, Vol. 37, Issue 3
Self-consistent-field calculations using Chebyshev-filtered subspace iteration
journal, November 2006
- Zhou, Yunkai; Saad, Yousef; Tiago, Murilo L.
- Journal of Computational Physics, Vol. 219, Issue 1
JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices
dataset, January 2019
- Bollhöfer, Matthias
- Mendeley
Works referencing / citing this record:
Domain decomposition approaches for accelerating contour integration eigenvalue solvers for symmetric eigenvalue problems: Domain decomposition contour integration eigensolvers
journal, February 2018
- Kalantzis, Vassilis; Kestyn, James; Polizzi, Eric
- Numerical Linear Algebra with Applications, Vol. 25, Issue 5