Abstract:We consider the approximation of the smallest eigenvalue of a large parameter-dependent Hermitian matrix over a continuum compact domain. Our approach is based on approximating the smallest eigenvalue by the one obtained by projecting the large matrix onto a suitable small subspace, a practice widely employed in the literature. The projection subspaces are constructed iteratively (to reduce the error of the approximation where it is large) with the addition of the eigenvectors of the parameter-dependent matrix at the parameter values where a surrogate error is maximal. The surrogate error is the gap between the approximation and a lower bound for the smallest eigenvalue proposed in [Sirkovic and Kressner, SIAM J. Matrix Anal. Appl., 37(2), 2016]. Unlike the classical approaches, such as the successive constraint method, that maximize such surrogate errors over a discrete and finite set, we maximize the surrogate error over the continuum of all permissible parameter values globally. We put particular attention to the lower bound, which enables us to formally prove the global convergence of our framework both in finite-dimensional and infinite-dimensional settings. In the second part, we focus on the approximation of the smallest singular value of a large parameter-dependent matrix, in case it is non-Hermitian, and propose another subspace framework to construct a small parameter-dependent non-Hermitian matrix whose smallest singular value approximates the original large-scale smallest singular value. We perform numerical experiments on synthetic examples, as well as on real examples arising from parametric PDEs. The numerical experiments show that the proposed techniques are able to drastically reduce the size of the large parameter-dependent matrix, while ensuring an approximation error for the smallest eigenvalue/singular value below the prescribed tolerance.

Matrix-less spectral approximation for large structured matrices

Matrix‐less methods for the spectral approximation of large non‐Hermitian Toeplitz matrices: A concise theoretical analysis and a numerical study

Fast Non-Hermitian Toeplitz Eigenvalue Computations, Joining Matrixless Algorithms and FDE Approximation Matrices

Fast Toeplitz eigenvalue computations, joining interpolation-extrapolation matrix-less algorithms and simple-loop theory

A note on eigenvalues and singular values of variable Toeplitz matrices and matrix-sequences, with application to variable two-step BDF approximations to parabolic equations

Eigenvalue superposition for Toeplitz matrix-sequences with matrix order dependent symbols

Block structured matrix-sequences and their spectral and singular value canonical distributions: a general theory

Approximating matrix eigenvalues by subspace iteration with repeated random sparsification

Asymptotic Spectra of Large (grid) Graphs with a Uniform Local Structure (part II): Numerical Applications

Exact Eigenvalues and Eigenvectors for Some n-Dimensional Matrices

Uniform Approximation of Eigenproblems of a Large-Scale Parameter-Dependent Hermitian Matrix

Sublinear Time Eigenvalue Approximation via Random Sampling

A smoothing spline algorithm to interpolate and predict the eigenvalues of matrices extracted from the sequence of preconditioned banded symmetric Toeplitz matrices

A Note on Generalized Locally Toeplitz Operators

A new iterative method for simultaneous computation of several eigenpair derivatives of a large matrix

Optimal Structured Matrix Approximation for Robustness to Incomplete Biosequence Data

Diagonal Matrix Sequences and their Spectral Symbols

A numerical projection technique for large-scale eigenvalue problems

Sublinear Time Spectral Density Estimation

A Quadratically Convergent Algorithm for Structured Low-Rank Approximation

Fine-grained Analysis and Faster Algorithms for Iteratively Solving Linear Systems