{"title":"大规模参数相关赫米矩阵特征问题的统一逼近","authors":"Mattia Manucci, Emre Mengi, Nicola Guglielmi","doi":"arxiv-2409.05791","DOIUrl":null,"url":null,"abstract":"We consider the approximation of the smallest eigenvalue of a large\nparameter-dependent Hermitian matrix over a continuum compact domain. Our\napproach is based on approximating the smallest eigenvalue by the one obtained\nby projecting the large matrix onto a suitable small subspace, a practice\nwidely employed in the literature. The projection subspaces are constructed\niteratively (to reduce the error of the approximation where it is large) with\nthe addition of the eigenvectors of the parameter-dependent matrix at the\nparameter values where a surrogate error is maximal. The surrogate error is the\ngap between the approximation and a lower bound for the smallest eigenvalue\nproposed in [Sirkovic and Kressner, SIAM J. Matrix Anal. Appl., 37(2), 2016].\nUnlike the classical approaches, such as the successive constraint method, that\nmaximize such surrogate errors over a discrete and finite set, we maximize the\nsurrogate error over the continuum of all permissible parameter values\nglobally. We put particular attention to the lower bound, which enables us to\nformally prove the global convergence of our framework both in\nfinite-dimensional and infinite-dimensional settings. In the second part, we\nfocus on the approximation of the smallest singular value of a large\nparameter-dependent matrix, in case it is non-Hermitian, and propose another\nsubspace framework to construct a small parameter-dependent non-Hermitian\nmatrix whose smallest singular value approximates the original large-scale\nsmallest singular value. We perform numerical experiments on synthetic\nexamples, as well as on real examples arising from parametric PDEs. The\nnumerical experiments show that the proposed techniques are able to drastically\nreduce the size of the large parameter-dependent matrix, while ensuring an\napproximation error for the smallest eigenvalue/singular value below the\nprescribed tolerance.","PeriodicalId":501162,"journal":{"name":"arXiv - MATH - Numerical Analysis","volume":"26 1","pages":""},"PeriodicalIF":0.0000,"publicationDate":"2024-09-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Uniform Approximation of Eigenproblems of a Large-Scale Parameter-Dependent Hermitian Matrix\",\"authors\":\"Mattia Manucci, Emre Mengi, Nicola Guglielmi\",\"doi\":\"arxiv-2409.05791\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"We consider the approximation of the smallest eigenvalue of a large\\nparameter-dependent Hermitian matrix over a continuum compact domain. Our\\napproach is based on approximating the smallest eigenvalue by the one obtained\\nby projecting the large matrix onto a suitable small subspace, a practice\\nwidely employed in the literature. The projection subspaces are constructed\\niteratively (to reduce the error of the approximation where it is large) with\\nthe addition of the eigenvectors of the parameter-dependent matrix at the\\nparameter values where a surrogate error is maximal. The surrogate error is the\\ngap between the approximation and a lower bound for the smallest eigenvalue\\nproposed in [Sirkovic and Kressner, SIAM J. Matrix Anal. Appl., 37(2), 2016].\\nUnlike the classical approaches, such as the successive constraint method, that\\nmaximize such surrogate errors over a discrete and finite set, we maximize the\\nsurrogate error over the continuum of all permissible parameter values\\nglobally. We put particular attention to the lower bound, which enables us to\\nformally prove the global convergence of our framework both in\\nfinite-dimensional and infinite-dimensional settings. In the second part, we\\nfocus on the approximation of the smallest singular value of a large\\nparameter-dependent matrix, in case it is non-Hermitian, and propose another\\nsubspace framework to construct a small parameter-dependent non-Hermitian\\nmatrix whose smallest singular value approximates the original large-scale\\nsmallest singular value. We perform numerical experiments on synthetic\\nexamples, as well as on real examples arising from parametric PDEs. The\\nnumerical experiments show that the proposed techniques are able to drastically\\nreduce the size of the large parameter-dependent matrix, while ensuring an\\napproximation error for the smallest eigenvalue/singular value below the\\nprescribed tolerance.\",\"PeriodicalId\":501162,\"journal\":{\"name\":\"arXiv - MATH - Numerical Analysis\",\"volume\":\"26 1\",\"pages\":\"\"},\"PeriodicalIF\":0.0000,\"publicationDate\":\"2024-09-09\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"arXiv - MATH - Numerical Analysis\",\"FirstCategoryId\":\"1085\",\"ListUrlMain\":\"https://doi.org/arxiv-2409.05791\",\"RegionNum\":0,\"RegionCategory\":null,\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"\",\"JCRName\":\"\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"arXiv - MATH - Numerical Analysis","FirstCategoryId":"1085","ListUrlMain":"https://doi.org/arxiv-2409.05791","RegionNum":0,"RegionCategory":null,"ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"","JCRName":"","Score":null,"Total":0}
Uniform Approximation of Eigenproblems of a Large-Scale Parameter-Dependent Hermitian Matrix
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.