Search

A new algorithm is developed for computing etAB, where A is an n x n matrix and B is n x n0 with n0 « n. The algorithm works for any A, its computational cost is dominated by the formation of products of A with n x n0 matrices, and the only input parameter is a backward error tolerance. The algorithm can return a single matrix etA B or a sequence Due to image rights restrictions, multiple line equation(s) cannot be graphically displayed. on an equally spaced grid of points tk. It uses the scaling part of the scaling and squaring method together with a truncated Taylor series approximation to the exponential. It determines the amount of scaling and the Taylor degree, using the recent analysis of Al-Mohy and Higham [SIAM J. Matrix Anal. Appl., 31 (2009), pp. 970 989], which provides sharp truncation error bounds expressed in terms of the quantities ||Ak