Block Krylov-type complex moment-based eigensolvers for solving generalized eigenvalue problems

Complex moment-based eigensolvers for solving interior eigenvalue problems have been studied because of their high parallel efficiency. Recently, we proposed the block Arnoldi-type complex moment-based eigensolver without a low-rank approximation. A low-rank approximation plays a very important role in reducing computational cost and stabilizing accuracy in complex moment-based eigensolvers. In this paper, we develop the method and propose block Krylov-type complex moment-based eigensolvers with a low-rank approximation. Numerical experiments indicate that the proposed methods have higher performance than the block SS–RR method, which is one of the most typical complex moment-based eigensolvers.

[1]  Yousef Saad,et al.  A spectrum slicing method for the Kohn-Sham problem , 2012, Comput. Phys. Commun..

[2]  Tetsuya Sakurai,et al.  Efficient Parameter Estimation and Implementation of a Contour Integral-Based Eigensolver , 2013 .

[3]  Eric Polizzi,et al.  A Density Matrix-based Algorithm for Solving Eigenvalue Problems , 2009, ArXiv.

[4]  Tetsuya Sakurai,et al.  Error bounds of Rayleigh–Ritz type contour integral-based eigensolver for solving generalized eigenvalue problems , 2015, Numerical Algorithms.

[5]  T. Sakurai,et al.  CIRR: a Rayleigh-Ritz type method with contour integral for generalized eigenvalue problems , 2007 .

[6]  Tetsuya Sakurai,et al.  Relationships among contour integral-based methods for solving generalized eigenvalue problems , 2016 .

[7]  Ping Tak Peter Tang,et al.  FEAST As A Subspace Iteration Eigensolver Accelerated By Approximate Spectral Projection , 2013, SIAM J. Matrix Anal. Appl..

[8]  Ping Tak Peter Tang,et al.  Zolotarev Quadrature Rules and Load Balancing for the FEAST Eigensolver , 2014, SIAM J. Sci. Comput..

[9]  Tetsuya Sakurai,et al.  CONTOUR INTEGRAL EIGENSOLVER FOR NON-HERMITIAN SYSTEMS: A RAYLEIGH-RITZ-TYPE APPROACH , 2010 .

[10]  Mark Hoemmen,et al.  Communication-avoiding Krylov subspace methods , 2010 .

[11]  T. Sakurai,et al.  A projection method for generalized eigenvalue problems using numerical integration , 2003 .

[12]  Timothy A. Davis,et al.  The university of Florida sparse matrix collection , 2011, TOMS.

[13]  Tetsuya Sakurai,et al.  A filter diagonalization for generalized eigenvalue problems based on the Sakurai-Sugiura projection method , 2008, J. Comput. Appl. Math..

[14]  A. Dubrulle,et al.  Retooling the method of block conjugate gradients. , 2001 .

[15]  Tetsuya Sakurai,et al.  A block Arnoldi-type contour integral spectral projection method for solving generalized eigenvalue problems , 2014, Appl. Math. Lett..