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.

Accurate computations of singular values and linear systems for Polynomial-Vandermonde-type matrices

Accurate singular values of a class of parameterized negative matrices

Extended Newton-type Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

An Extended Ulm-like Method for Inverse Singular Value Problems with Multiple And/or Zero Singular Values

Newton-type Methods for Inverse Singular Value Problems with Multiple Singular Values

Convergence of a Ulm-like Method for Square Inverse Singular Value Problems with Multiple and Zero Singular Values.

Sum-of-squares Chordal Decomposition of Polynomial Matrix Inequalities

Convergence of the Newton-type methods for the square inverse singular value problems with multiple and zero singular values

A Contour Integral-Based Algorithm for Computing Generalized Singular Values

Numerical methods for eigenvalues of singular polynomial eigenvalue problems

Accurate solutions of M -matrix Sylvester equations

An Efficient GIPM Algorithm for Computing the Smallest V-Singular Value of the Partially Symmetric Tensor

High Order Numerical Methods To Approximate The Singular Value Decomposition

Generalizing the SVD of a matrix under non-standard inner product and its applications to linear ill-posed problems

An Inexact Newton-type Method for Inverse Singular Value Problems

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

Mixed-Precision Paterson--Stockmeyer Method for Evaluating Polynomials of Matrices

On the Closed Form Expression of Elementary Symmetric Polynomials and the Inverse of Vandermonde Matrix

Accurate computations with Said-Ball-Vandermonde matrices

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