Mathematics > Numerical Analysis
[Submitted on 9 Sep 2024 (v1), last revised 24 Sep 2024 (this version, v2)]
Title:Uniform Approximation of Eigenproblems of a Large-Scale Parameter-Dependent Hermitian Matrix
View PDFAbstract: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.
Submission history
From: Emre Mengi [view email][v1] Mon, 9 Sep 2024 16:51:09 UTC (138 KB)
[v2] Tue, 24 Sep 2024 16:19:54 UTC (144 KB)
Current browse context:
math.NA
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Connected Papers (What is Connected Papers?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
alphaXiv (What is alphaXiv?)
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Hugging Face (What is Huggingface?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.