Abstract:The distribution of the eigenvalues of a Hermitian matrix (or of a Hermitian matrix pencil) reveals important features of the underlying problem, whether a Hamiltonian system in physics, or a social network in behavioral sciences. However, computing all the eigenvalues explicitly is prohibitively expensive for real-world applications. This paper presents two types of methods to efficiently estimate the spectral density of a matrix pencil $(A, B)$ when both $A$ and $B$ are Hermitian and, in addition, $B$ is positive definite. The first one is based on the Kernel Polynomial Method (KPM) and the second on Gaussian quadrature by the Lanczos procedure. By employing Chebyshev polynomial approximation techniques, we can avoid direct factorizations in both methods, making the resulting algorithms suitable for large matrices. Under some assumptions, we prove bounds that suggest that the Lanczos method converges twice as fast as the KPM method. Numerical examples further indicate that the Lanczos method can provide more accurate spectral densities when the eigenvalue distribution is highly non-uniform. As an application, we show how to use the computed spectral density to partition the spectrum into intervals that contain roughly the same number of eigenvalues. This procedure, which makes it possible to compute the spectrum by parts, is a key ingredient in the new breed of eigensolvers that exploit "spectrum slicing".

Fast parallelizable methods for the Hermitian eigenvalue problem

A Cubically Convergent Parallelizable Method for the Hermitian Eigenvalue Problem

Fast Parallelizable Methods for Computing Invariant Subspaces of Hermitian Matrices

A Quadratically Convergent QR-like Method Without Shifts for the Hermitian Eigenvalue Problem

A Parallelizable Method for Computing the Singular Subspaces of a General Matrix

Accelerating Parallel Jacobi Method for Matrix Eigenvalue Computation in DOA Estimation Algorithm

A Convergence Analysis of the Inexact Rayleigh Quotient Iteration and Simplified Jacobi-Davidson Method for the Large Hermitian Matrix Eigenproblem

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

A Ulm-like Method for Inverse Eigenvalue Problems

A projected preconditioned conjugate gradient algorithm for computing many extreme eigenpairs of a Hermitian matrix

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

The Convergence of Generalized Lanczos Methods for Large Unsymmetric Eigenproblems

Fast computation of spectral densities for generalized eigenvalue problems

Convergence Analysis of Newton-Like Methods for Inverse Eigenvalue Problems with Multiple Eigenvalues

Approximate Cayley Transform Methods for Inverse Eigenvalue Problems and Convergence Analysis

A parallel iterative method for solving symmetric tridiagonal extreme eigenpair problems

On Rayleigh Quotient Iteration for Dual Quaternion Hermitian Eigenvalue Problem

Accelerating Convergence by Augmented Rayleigh-Ritz Projections for Large-Scale Eigenpair Computation

Generalized Jacobi Method for Computing Eigenvalues of Dual Quaternion Hermitian Matrices

The Convergence of Krylov Subspace Methods for Large Unsymmetric Linear Systems