Rational spectral filters with optimal convergence rate

01/13/2020
by   Konrad Kollnig, et al.
0

In recent years, contour-based eigensolvers have emerged as a standard approach for the solution of large and sparse eigenvalue problems. Building upon recent performance improvements through non-linear least square optimization of so-called rational filters, we introduce a systematic method to design these filters by minimizing the worst-case convergence ratio and eliminate the parametric dependence on weight functions. Further, we provide an efficient way to deal with the box-constraints which play a central role for the use of iterative linear solvers in contour-based eigensolvers. Indeed, these parameter-free filters consistently minimize the number of iterations and the number of FLOPs to reach convergence in the eigensolver. As a byproduct, our rational filters allow for a simple solution to load balancing when the solution of an interior eigenproblem is approached by the slicing of the sought after spectral interval.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
04/11/2017

Non-Linear Least-Squares Optimization of Rational Filters for the Solution of Interior Eigenvalue Problems

Rational filter functions can be used to improve convergence of contour-...
research
07/06/2020

An Iterative Method for Contour-Based Nonlinear Eigensolvers

Contour integration techniques have become a popular choice for solving ...
research
07/06/2020

An Iterative Method for Contour-Based Nonlinear Egensolvers

Contour integration techniques have become a popular choice for solving ...
research
01/31/2017

Interior Eigensolver for Sparse Hermitian Definite Matrices Based on Zolotarev's Functions

This paper proposes an efficient method for computing selected generaliz...
research
01/12/2022

Majorization-type cluster robust bounds for block filters and eigensolvers

Convergence analysis of block iterative solvers for Hermitian eigenvalue...
research
11/08/2020

Fourier-based and Rational Graph Filters for Spectral Processing

Data are represented as graphs in a wide range of applications, such as ...
research
10/19/2020

Twice is enough for dangerous eigenvalues

We analyze the stability of a class of eigensolvers that target interior...

Please sign up or login with your details

Forgot password? Click here to reset