A Fast Hierarchically Preconditioned Eigensolver Based On Multiresolution Matrix Decomposition

04/10/2018
by   Thomas Y. Hou, et al.
0

In this paper we propose a new iterative method to hierarchically compute a relatively large number of leftmost eigenpairs of a sparse symmetric positive matrix under the multiresolution operator compression framework. We exploit the well-conditioned property of every decomposition components by integrating the multiresolution framework into the Implicitly restarted Lanczos method. We achieve this combination by proposing an extension-refinement iterative scheme, in which the intrinsic idea is to decompose the target spectrum into several segments such that the corresponding eigenproblem in each segment is well-conditioned. Theoretical analysis and numerical illustration are also reported to illustrate the efficiency and effectiveness of this algorithm.

READ FULL TEXT

page 33

page 34

page 35

research
04/14/2022

A New Refinement-Free Preconditioner for the Symmetric Formulation in Electroencephalography

Widely employed for the accurate solution of the electroencephalography ...
research
05/30/2021

An iterative Jacobi-like algorithm to compute a few sparse eigenvalue-eigenvector pairs

In this paper, we describe a new algorithm to compute the extreme eigenv...
research
12/15/2022

A skew-symmetric Lanczos bidiagonalization method for computing several largest eigenpairs of a large skew-symmetric matrix

The spectral decomposition of a real skew-symmetric matrix A can be math...
research
05/17/2021

Full operator preconditioning and the accuracy of solving linear systems

Unless special conditions apply, the attempt to solve ill-conditioned sy...
research
07/24/2023

Robust iterative method for symmetric quantum signal processing in all parameter regimes

This paper addresses the problem of solving nonlinear systems in the con...
research
07/18/2023

Cross-points in the Dirichlet-Neumann method II: a geometrically convergent variant

When considered as a standalone iterative solver for elliptic boundary v...
research
08/04/2020

Well-Conditioned Methods for Ill-Conditioned Systems: Linear Regression with Semi-Random Noise

Classical iterative algorithms for linear system solving and regression ...

Please sign up or login with your details

Forgot password? Click here to reset