Stable Computation of Generalized Matrix Functions via Polynomial Interpolation
Stable Computation of Generalized Matrix Functions via Polynomial Interpolation
复制标题
通过多项式插值稳定计算广义矩阵函数
DOI:
10.1137/18m1191786
复制
发表时间:
2019
影响因子:
1.5
通讯作者:
Kalantzis, Vassilis
中科院分区:
文献类型:
--
作者:
Aurentz, Jared L.;Austin, Anthony P.;Benzi, Michele;Kalantzis, Vassilis
Generalized matrix functions (GMFs) extend the concept of a matrix function to rectangular matrices via the singular value decomposition. Several applications involving directed graphs, Hamiltonian dynamical systems, and optimization problems with low-rank constraints require the action of a GMF of a large, sparse matrix on a vector. We present a new method for applying GMFs to vectors based on Chebyshev interpolation. The method is matrix free and requires no orthogonalization and minimal additional storage. Comparisons against existing approaches based on Lanczos bidiagonalization demonstrate the competitiveness of our approach. We prove that our method is backward stable by generalizing the proof of the backward stability of Clenshaw's algorithm to the matrix case.
登录
查看更多内容
影响因子:
1.5
作者:
F. Arrigo;M. Benzi;C. Fenu
通讯作者:
C. Fenu
影响因子:
1.5
作者:
V. Noferini
通讯作者:
V. Noferini
DOI:
10.1016/j.cpc.2017.06.016
发表时间:
2017
期刊:
Comput. Phys. Commun.
影响因子:
--
作者:
J. L. Aurentz;V. Kalantzis;Y. Saad
通讯作者:
Y. Saad
影响因子:
4.6
作者:
N. Buono;L. Lopez;T. Politi
通讯作者:
T. Politi
DOI:
--
发表时间:
2010
期刊:
影响因子:
--
作者:
J. J. Crofts;Ernesto Estrada;D. Higham;Alan Taylor
通讯作者:
Alan Taylor