DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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:
 [1];  [2];  [3];  [3];  [2]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Univ. of Minnesota, Twin Cities, MN (United States)
  3. 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}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 32 works
Citation information provided by
Web of Science

Save / Share:

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
  • DOI: 10.1145/838250.838257

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
  • DOI: 10.1145/1527286.1527287

JADAMILU: a software code for computing selected eigenvalues of large sparse symmetric matrices
journal, December 2007


Strategies for spectrum slicing based on restarted Lanczos methods
journal, March 2012


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
  • DOI: 10.1137/110836535

Rounding error analysis of the classical Gram-Schmidt orthogonalization process
journal, May 2005


The loss of orthogonality in the Gram-Schmidt orthogonalization process
journal, October 2005


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
  • DOI: 10.1145/1089014.1089019

Electronic structure calculations for plane-wave codes without diagonalization
journal, April 1999


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
  • DOI: 10.1002/pssb.200541463

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
  • DOI: 10.1137/S0895479895281484

Approximating Spectral Densities of Large Matrices
journal, January 2016

  • Lin, Lin; Saad, Yousef; Yang, Chao
  • SIAM Review, Vol. 58, Issue 1
  • DOI: 10.1137/130934283

Computational Variants of the Lanczos Method for the Eigenproblem
journal, January 1972


Error Analysis of the Lanczos Algorithm for Tridiagonalizing a Symmetric Matrix
journal, January 1976


Density-matrix-based algorithm for solving eigenvalue problems
journal, March 2009


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
  • DOI: 10.1137/0906059

Filtered Conjugate Residual‐type Algorithms with Applications
journal, January 2006

  • Saad, Yousef
  • SIAM Journal on Matrix Analysis and Applications, Vol. 28, Issue 3
  • DOI: 10.1137/060648945

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
  • DOI: 10.1137/141002037

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
  • DOI: 10.1137/060651653

A projection method for generalized eigenvalue problems using numerical integration
journal, October 2003


The Lanczos algorithm with partial reorthogonalization
journal, January 1984


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
  • DOI: 10.1145/1731022.1731031

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
  • DOI: 10.1137/S1064827596304162

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix
journal, June 2015


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
  • DOI: 10.1145/1824801.1824805

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
  • DOI: 10.1016/j.jcp.2006.03.017

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
  • DOI: 10.1002/nla.2154