Error control technique of quadrature-based algorithms for the action of real powers of a Hermitian positive-definite matrix

Abstract

This study considers quadrature-based algorithms to compute $A^{\alpha }\mathbf{b}$, the action of a real power of a Hermitian positive-definite matrix $A$ on a vector $\mathbf{b}$. In these algorithms, the computation of an integral representation of $A^{\alpha }\mathbf{b}$ is reduced to solving several tens or hundreds of shifted linear systems. Current approaches usually analyze the quadrature discretization error, but rarely take into account the additional error introduced by solving these shifted linear systems with iterative solvers. Here, we bound this error with the residual of the approximated solution of these linear systems. This allows the derivation of a stopping criterion for iterative solvers to keep the error of $A^{\alpha }\mathbf{b}$ below a prescribed error tolerance. Numerical results demonstrate that the proposed criterion enables the computation of $A^{\alpha }\mathbf{b}$ within prescribed tolerance limits.