We analyze the Lanczos method for matrix function approximation (Lanczos-FA), an iterative algorithm for computing $f(\mathbf{A}) \mathbf{b}$ when $\mathbf{A}$ is a Hermitian matrix and $\mathbf{b}$ is a given mathbftor. Assuming that $f : \mathbb{C} \rightarrow \mathbb{C}$ is piecewise analytic, we give a framework, based on the Cauchy integral formula, which can be used to derive {\em a priori} and \emph{a posteriori} error bounds for Lanczos-FA in terms of the error of Lanczos used to solve linear systems. Unlike many error bounds for Lanczos-FA, these bounds account for fine-grained properties of the spectrum of $\mathbf{A}$, such as clustered or isolated eigenvalues. Our results are derived assuming exact arithmetic, but we show that they are easily extended to finite precision computations using existing theory about the Lanczos algorithm in finite precision. We also provide generalized bounds for the Lanczos method used to approximate quadratic forms $\mathbf{b}^\textsf{H} f(\mathbf{A}) \mathbf{b}$, and demonstrate the effectiveness of our bounds with numerical experiments.


翻译:我们分析矩阵函数近似( Lanczos) 的 Lanczos 方法( Lanczos 方法 ) (Lanczos- FA), 这是计算 $f (\ mathbf{A})\ mathbf{A} 的迭代算法, 当$\ mathbf{A} 是一个 Hermitian 的矩阵和$\ mathb{b} 是一个给定的数学ftor 。 假设$f:\ mathb{C} \ rightrow \ mathb{C} 是一个小分析器, 我们给出了一个框架, 这个框架可以用来计算 $( mathb{A} ) 和\ mapha postioria} 的元值, 当 以 兰czos 用于解决线性体系的错误来计算 时, 计算 。 兰茨- falfb} 错误是用来计算精度的精度 。 我们用精确度 的精确度计算方法来得出 精确度, 我们的精确度的精确度的计算方法也提供了我们现有的精确度。

0
下载
关闭预览

相关内容

专知会员服务
76+阅读 · 2021年3月16日
Linux导论,Introduction to Linux,96页ppt
专知会员服务
77+阅读 · 2020年7月26日
【伯克利-Ke Li】学习优化,74页ppt,Learning to Optimize
专知会员服务
40+阅读 · 2020年7月23日
Fariz Darari简明《博弈论Game Theory》介绍,35页ppt
专知会员服务
109+阅读 · 2020年5月15日
已删除
将门创投
5+阅读 · 2017年10月20日
Arxiv
0+阅读 · 2021年8月19日
Arxiv
0+阅读 · 2021年8月18日
Arxiv
0+阅读 · 2021年8月15日
Arxiv
3+阅读 · 2018年10月18日
VIP会员
相关资讯
已删除
将门创投
5+阅读 · 2017年10月20日
Top
微信扫码咨询专知VIP会员