High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations

Journal article

Publication Details

Author(s): Wellein G, Alvermann A, Fehske H, Hager G, Kreutzer M, Lang B, Pieper A, Galgon M
Journal: Journal of Computational Physics
Publisher: Elsevier
Publication year: 2016
Volume: 325
Pages range: 226-243
ISSN: 0021-9991


We study Chebyshev filter diagonalization as a tool for the computation of many interior eigenvalues of very large sparse symmetric matrices. In this technique the subspace projection onto the target space of wanted eigenvectors is approximated with filter polynomials obtained from Chebyshev expansions of window functions. After the discussion of the conceptual foundations of Chebyshev filter diagonalization we analyze the impact of the choice of the damping kernel, search space size, and filter polynomial degree on the computational accuracy and effort, before we describe the necessary steps towards a parallel high-performance implementation. Because Chebyshev filter diagonalization avoids the need for matrix inversion it can deal with matrices and problem sizes that are presently not accessible with rational function methods based on direct or iterative linear solvers. To demonstrate the potential of Chebyshev filter diagonalization for large-scale problems of this kind we include as an example the computation of the 10(2) innermost eigenpairs of a topological insulator matrix with dimension 10(9) derived from quantum physics applications. (C) 2016 Elsevier Inc. All rights reserved.

FAU Authors / FAU Editors

Hager, Georg Dr.
Regionales Rechenzentrum Erlangen (RRZE)
Kreutzer, Moritz
Regionales Rechenzentrum Erlangen (RRZE)
Pieper, Andreas
Professur für Höchstleistungsrechnen
Wellein, Gerhard Prof. Dr.
Professur für Höchstleistungsrechnen

How to cite

Wellein, G., Alvermann, A., Fehske, H., Hager, G., Kreutzer, M., Lang, B.,... Galgon, M. (2016). High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations. Journal of Computational Physics, 325, 226-243. https://dx.doi.org/10.1016/j.jcp.2016.08.027

Wellein, Gerhard, et al. "High-performance implementation of Chebyshev filter diagonalization for interior eigenvalue computations." Journal of Computational Physics 325 (2016): 226-243.


Last updated on 2018-04-06 at 19:23