期刊
ELECTRONIC TRANSACTIONS ON NUMERICAL ANALYSIS
卷 55, 期 -, 页码 213-241出版社
KENT STATE UNIVERSITY
DOI: 10.1553/etna_vol55s213
关键词
Schatten p-norm; Monte Carlo estimator; optimal experimental design; Chebyshev polynomials
资金
- National Science Foundation [DMS - 1745654]
This paper presents numerical methods for computing the Schatten p-norm of positive semi-definite matrices. It proposes a matrix-free method using a Monte Carlo estimator and extends the convergence and error analysis to the computation of non-integer and large values of p. The performance of the proposed estimators is demonstrated on test matrices and in an application to optimal experimental design for a model inverse problem.
We present numerical methods for computing the Schatten p-norm of positive semi-definite matrices. Our motivation stems from uncertainty quantification and optimal experimental design for inverse problems, where the Schatten p-norm defines a measure of uncertainty. Computing the Schatten p-norm of high-dimensional matrices is computationally expensive. We propose a matrix-free method to estimate the Schatten p-norm using a Monte Carlo estimator and derive convergence results and error estimates for the estimator. To efficiently compute the Schatten p-norm for non-integer and large values of p, we use an estimator using Chebyshev polynomial approximations and extend our convergence and error analysis to this setting as well. We demonstrate the performance of our proposed estimators on several test matrices and in an application to optimal experimental design for a model inverse problem.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据