ZOLOTAREV QUADRATURE RULES AND LOAD BALANCING FOR THE FEAST EIGENSOLVER - CentraleSupélec Access content directly
Journal Articles SIAM Journal on Scientific Computing Year : 2015

ZOLOTAREV QUADRATURE RULES AND LOAD BALANCING FOR THE FEAST EIGENSOLVER

Abstract

The FEAST method for solving large sparse eigenproblems is equivalent to subspace iteration with an approximate spectral projector and implicit orthogonalization. This relation allows us to characterize the convergence of this method in terms of the error of a certain rational approximant to an indicator function. We propose improved rational approximants leading to FEAST variants with faster convergence, in particular, when using rational approximants based on the work of Zolotarev. Numerical experiments demonstrate the possible computational savings especially for pencils whose eigenvalues are not well separated and when the dimension of the search space is only slightly larger than the number of wanted eigenvalues. The new approach improves both convergence robustness and load balancing when FEAST runs on multiple search intervals in parallel.

Dates and versions

hal-02418265 , version 1 (18-12-2019)

Identifiers

Cite

Stefan Guettel, Eric Polizzi, Ping Tak Peter Tang, Gautier Viaud. ZOLOTAREV QUADRATURE RULES AND LOAD BALANCING FOR THE FEAST EIGENSOLVER. SIAM Journal on Scientific Computing, 2015, 37 (4), pp.A2100-A2122. ⟨10.1137/140980090⟩. ⟨hal-02418265⟩
36 View
0 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More