Estimation of the bilinear form y* f(A)x for Hermitian matrices

摘要

For a Hermitian matrix A is an element of C-pxp, given vectors x, y is an element of C-P and for suitable functions f, the bilinear form y* f(A)x is estimated by extending the extrapolation method proposed by C. Brezinski in 1999. Families of one term and two term estimates e(f,v), nu is an element of C and (e) over cap (f,n,k), n, k is an element of Z, respectively, are derived by extrapolation of the moments of the matrix A. For the positive definite case, bounds for the optimal value of v, which leads to an efficient one term estimate in only one matrix vector product, are derived. For f (A) = A(-1) a formula approximating this optimal value of v is specified. Numerical results for several matrix functions and comparisons are provided to demonstrate the effectiveness of the extrapolation method. 2015 Elsevier Inc. All rights reserved.