The extended Krylov subspace method and orthogonal Laurent polynomials
The extended Krylov subspace method and orthogonal Laurent polynomials
复制标题
DOI:
10.1016/j.laa.2009.03.006
复制
发表时间:
2009-07
影响因子:
1.1
通讯作者:
Carl Jagels;L. Reichel
中科院分区:
文献类型:
--
作者:
Carl Jagels;L. Reichel
The need to evaluate expressions of the form f(A)v, where A is a large sparse or structured symmetric matrix, v is a vector, and f is a nonlinear function, arises in many applications. The extended Krylov subspace method can be an attractive scheme for computing approximations of such expressions. This method projects the approximation problem onto an extended Krylov subspace Kℓ,m(A)=span{A-ℓ+1v,…,A-1v,v,Av,…,Am-1v} of fairly small dimension, and then solves the small approximation problem so obtained. We review available results for the extended Krylov subspace method and relate them to properties of Laurent polynomials. The structure of the projected problem receives particular attention. We are concerned with the situations when m=ℓ and m=2ℓ.