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".

Realization of Extremal Spectral Data by Pentadiagonal Matrices

Extension and Related Realization of the Pseudo-Diagonalization Method

A matrix realization of spectral bounds

A general theorem in spectral extremal graph theory

Extremal Eigenvalues of Random Kernel Matrices with Polynomial Scaling

The Projection Methods for Computing the Pseudospectra of Large Scale Matrices

Differential Geometry with Extreme Eigenvalues in the Positive Semidefinite Cone

Fast computation of spectral densities for generalized eigenvalue problems

A Skew-Symmetric Lanczos Bidiagonalization Method for Computing Several Extremal Eigenpairs of a Large Skew-Symmetric Matrix

Inverse eigenvalue problem for pseudo-symmetric Jacobi matrices with two spectra

Extremal inverse eigenvalue problem for matrices described by a connected unicyclic graph

Extremal spectral behavior of weighted random $d$-regular graphs

An Inverse Spectral Problem for Non-Self-Adjoint Jacobi Matrices

Structured Pseudospectra and Structured Sensitivity of Eigenvalues

Sylvester–Kac matrices with quadratic spectra: A comprehensive note

An Inverse Extreme Eigenpair Problem and Its Parallel Iterative Solution for Symmetric Tridiagonal Matrices

A Parallel Iterative Method For Solving Symmetrical Tridiagonal Extreme Eigenpair Problems

Graph Limits and Spectral Extremal Problems for Graphs

Universality of extremal eigenvalues of large random matrices

Partial Proof of a Conjecture with Implications for Spectral Majorization