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.

Convergence, optimization and stability of singular eigenmaps

Continuity of eigenvalues and shape optimisation for Laplace and Steklov problems

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

Convergence Behavior of Gauss-Newton's Method and Extensions of the Smale Point Estimate Theory.

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

On Convergence Lemma and Convergence Stability for Piecewise Analytic Functions

Dirac Eigenvalue Optimisation and Harmonic Maps to Complex Projective Spaces

Relative convergence estimates for the spectral asymptotic in the Large Coupling Limit

Convergence of Laplacian Spectra from Random Samples

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

Convergence of Laplacian Spectra from Point Clouds

Multiple Sets Exponential Concentration and Higher Order Eigenvalues

Variational analysis of unbounded and discontinuous generalized eigenvalue functions with application to topology optimization

Characterizing Parametric and Convergence Stability in Nonconvex and Nonsmooth Optimizations: A Geometric Approach

Optimization on manifolds: A symplectic approach

The Eigenvalue Problem for the Laplacian via Conformal Mapping and the Gohberg--Sigal Theory

Convergence and error control of consistent PINNs for elliptic PDEs

Convergence of Finite Integral Method for Poisson Equation and Corresponding Eigenproblems with Dirichlet Boundary

Convergence of Adaptive Finite Element Approximations for Nonlinear Eigenvalue Problems

Existence of harmonic maps and eigenvalue optimization in higher dimensions

Optimization of Neumann Eigenvalues under convexity and geometric constraints