UNIFORM APPROXIMATION OF phi-FUNCTIONS IN EXPONENTIAL INTEGRATORS BY A RATIONAL KRYLOV SUBSPACE METHOD WITH SIMPLE POLES

摘要

We consider the approximation of the matrix phi-functions that appear in exponential integrators for stiff systems of differential equations. For stiff systems, the field-of-values of the occurring matrices is large and lies somewhere in the left complex half-plane. In order to obtain an efficient method uniformly for all matrices with a field-of-values in the left complex half-plane, we consider the approximation by a rational Krylov subspace method with equidistant poles of order one on the line Re z = gamma > 0. We present error bounds that predict a faster convergence rate as for the resolvent Krylov subspace approximation using a single repeated pole at gamma > 0. Poles of order one allow moreover for a parallel implementation of the corresponding rational Krylov subspace decomposition. We analyze the convergence of the proposed rational Krylov subspace method and present numerical experiments that illustrate our results.