A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems

Li, Ruipeng; Xi, Yuanzhe; Vecharynski, Eugene; Yang, Chao; Saad, Yousef

doi:10.1137/15M1054493

A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems

Journal Article · Tue Aug 16 00:00:00 EDT 2016 · SIAM Journal on Scientific Computing

DOI:https://doi.org/10.1137/15M1054493· OSTI ID:1438696

Li, Ruipeng ^[1]; Xi, Yuanzhe ^[2]; Vecharynski, Eugene ^[3]; Yang, Chao ^[3]; Saad, Yousef ^[2]

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)

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.

View Accepted Manuscript (DOE)

Research Organization:: Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)

Sponsoring Organization:: USDOE

Grant/Contract Number:: AC52-07NA27344

OSTI ID:: 1438696

Report Number(s):: LLNL-JRNL--685856

Journal Information:: SIAM Journal on Scientific Computing, Journal Name: SIAM Journal on Scientific Computing Journal Issue: 4 Vol. 38; ISSN 1064-8275

Publisher:: SIAMCopyright Statement

Country of Publication:: United States

Language:: English

References (28)

PARSEC – the pseudopotential algorithm for real-space electronic structure calculations: recent advances and novel applications to nano-structures Kronik, Leeor; Makmal, Adi; Tiago, Murilo L. physica status solidi (b), Vol. 243, Issue 5 https://doi.org/10.1002/pssb.200541463	journal	April 2006
Rounding error analysis of the classical Gram-Schmidt orthogonalization process Giraud, Luc; Langou, Julien; Rozložník, Miroslav Numerische Mathematik, Vol. 101, Issue 1 https://doi.org/10.1007/s00211-005-0615-4	journal	May 2005
Strategies for spectrum slicing based on restarted Lanczos methods Campos, Carmen; Roman, Jose E. Numerical Algorithms, Vol. 60, Issue 2 https://doi.org/10.1007/s11075-012-9564-z	journal	March 2012
Electronic structure calculations for plane-wave codes without diagonalization Jay, Laurent O.; Kim, Hanchul; Saad, Yousef Computer Physics Communications, Vol. 118, Issue 1 https://doi.org/10.1016/S0010-4655(98)00192-1	journal	April 1999
A projection method for generalized eigenvalue problems using numerical integration Sakurai, Tetsuya; Sugiura, Hiroshi Journal of Computational and Applied Mathematics, Vol. 159, Issue 1 https://doi.org/10.1016/S0377-0427(03)00565-X	journal	October 2003
The loss of orthogonality in the Gram-Schmidt orthogonalization process Giraud, L.; Langou, J.; Rozloznik, M. Computers & Mathematics with Applications, Vol. 50, Issue 7 https://doi.org/10.1016/j.camwa.2005.08.009	journal	October 2005
JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices Bollhöfer, Matthias; Notay, Yvan Computer Physics Communications, Vol. 177, Issue 12 https://doi.org/10.1016/j.cpc.2007.08.004	journal	December 2007
Self-consistent-field calculations using Chebyshev-filtered subspace iteration Zhou, Yunkai; Saad, Yousef; Tiago, Murilo L. Journal of Computational Physics, Vol. 219, Issue 1 https://doi.org/10.1016/j.jcp.2006.03.017	journal	November 2006
A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix Vecharynski, Eugene; Yang, Chao; Pask, John E. Journal of Computational Physics, Vol. 290 https://doi.org/10.1016/j.jcp.2015.02.030	journal	June 2015
The Lanczos algorithm with partial reorthogonalization Simon, Horst D. Mathematics of Computation, Vol. 42, Issue 165 https://doi.org/10.1090/S0025-5718-1984-0725988-X	journal	January 1984
Computational Variants of the Lanczos Method for the Eigenproblem Paige, C. C. IMA Journal of Applied Mathematics, Vol. 10, Issue 3 https://doi.org/10.1093/imamat/10.3.373	journal	January 1972
Error Analysis of the Lanczos Algorithm for Tridiagonalizing a Symmetric Matrix Paige, C. C. IMA Journal of Applied Mathematics, Vol. 18, Issue 3 https://doi.org/10.1093/imamat/18.3.341	journal	January 1976
Density-matrix-based algorithm for solving eigenvalue problems Polizzi, Eric Physical Review B, Vol. 79, Issue 11 https://doi.org/10.1103/PhysRevB.79.115112	journal	March 2009
Filtered Conjugate Residual‐type Algorithms with Applications Saad, Yousef SIAM Journal on Matrix Analysis and Applications, Vol. 28, Issue 3 https://doi.org/10.1137/060648945	journal	January 2006
Numerical Methods for Electronic Structure Calculations of Materials Saad, Yousef; Chelikowsky, James R.; Shontz, Suzanne M. SIAM Review, Vol. 52, Issue 1 https://doi.org/10.1137/060651653	journal	January 2010
Practical Use of Polynomial Preconditionings for the Conjugate Gradient Method Saad, Youcef SIAM Journal on Scientific and Statistical Computing, Vol. 6, Issue 4 https://doi.org/10.1137/0906059	journal	October 1985
A Filtered Lanczos Procedure for Extreme and Interior Eigenvalue Problems Fang, Haw-ren; Saad, Yousef SIAM Journal on Scientific Computing, Vol. 34, Issue 4 https://doi.org/10.1137/110836535	journal	January 2012
Approximating Spectral Densities of Large Matrices Lin, Lin; Saad, Yousef; Yang, Chao SIAM Review, Vol. 58, Issue 1 https://doi.org/10.1137/130934283	journal	January 2016
Analysis of Subspace Iteration for Eigenvalue Problems with Evolving Matrices Saad, Yousef SIAM Journal on Matrix Analysis and Applications, Vol. 37, Issue 1 https://doi.org/10.1137/141002037	journal	January 2016
Deflation Techniques for an Implicitly Restarted Arnoldi Iteration Lehoucq, R. B.; Sorensen, D. C. SIAM Journal on Matrix Analysis and Applications, Vol. 17, Issue 4 https://doi.org/10.1137/S0895479895281484	journal	October 1996
Thick-Restart Lanczos Method for Large Symmetric Eigenvalue Problems Wu, Kesheng; Simon, Horst SIAM Journal on Matrix Analysis and Applications, Vol. 22, Issue 2 https://doi.org/10.1137/S0895479898334605	journal	January 2000
Dynamic Thick Restarting of the Davidson, and the Implicitly Restarted Arnoldi Methods Stathopoulos, Andreas; Saad, Yousef; Wu, Kesheng SIAM Journal on Scientific Computing, Vol. 19, Issue 1 https://doi.org/10.1137/S1064827596304162	journal	January 1998
SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems Hernandez, Vicente; Roman, Jose E.; Vidal, Vicente ACM Transactions on Mathematical Software, Vol. 31, Issue 3 https://doi.org/10.1145/1089014.1089019	journal	September 2005
Anasazi software for the numerical solution of large-scale eigenvalue problems Baker, C. G.; Hetmaniuk, U. L.; Lehoucq, R. B. ACM Transactions on Mathematical Software, Vol. 36, Issue 3 https://doi.org/10.1145/1527286.1527287	journal	July 2009
PRIMME: preconditioned iterative multimethod eigensolver—methods and software description Stathopoulos, Andreas; McCombs, James R. ACM Transactions on Mathematical Software, Vol. 37, Issue 2 https://doi.org/10.1145/1731022.1731031	journal	April 2010
Adaptive Projection Subspace Dimension for the Thick-Restart Lanczos Method Yamazaki, Ichitaro; Bai, Zhaojun; Simon, Horst ACM Transactions on Mathematical Software, Vol. 37, Issue 3 https://doi.org/10.1145/1824801.1824805	journal	September 2010
Algorithm 827: irbleigs: A MATLAB program for computing a few eigenpairs of a large sparse Hermitian matrix Baglama, J.; Calvetti, D.; Reichel, L. ACM Transactions on Mathematical Software, Vol. 29, Issue 3 https://doi.org/10.1145/838250.838257	journal	September 2003
JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices Bollhöfer, Matthias Mendeley https://doi.org/10.17632/dc7p5gjmjf	dataset	January 2019

Cited By (1)

Domain decomposition approaches for accelerating contour integration eigenvalue solvers for symmetric eigenvalue problems: Domain decomposition contour integration eigensolvers Kalantzis, Vassilis; Kestyn, James; Polizzi, Eric Numerical Linear Algebra with Applications, Vol. 25, Issue 5 https://doi.org/10.1002/nla.2154	journal	February 2018

Similar Records

Extending the eigCG algorithm to nonsymmetric Lanczos for linear systems with multiple right-hand sides

Journal Article · Fri Aug 01 00:00:00 EDT 2014 · Numerical Linear Algebra with Applications · OSTI ID:1156889

Projective block Lanczos algorithm for dense, Hermitian eigensystems

Journal Article · Thu Feb 29 23:00:00 EST 1996 · Journal of Computational Physics · OSTI ID:274230

Thick-restart Lanczos method for electronic structure calculations

Journal Article · Wed Sep 01 00:00:00 EDT 1999 · Journal of Computational Physics · OSTI ID:20000631

Related Subjects

97 MATHEMATICS AND COMPUTING
Lanczos algorithm
Thick-Restart
deflation
interior eigenvalue problems
polynomial filtering
spectrum slicing

A Thick-Restart Lanczos Algorithm with Polynomial Filtering for Hermitian Eigenvalue Problems

Citation Formats

References (28)

Cited By (1)

Similar Records

Related Subjects