Application of Weeks method for the numerical inversion of the Laplace transform to the matrix exponential. (English) Zbl 1092.65115

Summary: W. T. Week’s method [J. Assoc. Comput. Math. 13, 419–429 (1966; Zbl 0141.33401)] is a well established algorithm for the numerical inversion of scalar Laplace space functions. In this paper, we extend the method to the inversion of matrix functions of a single time variable and assess the qualities of this approach. To illustrate and quantify our discussion, we compute the matrix exponential by means of an fast Fourier transform based algorithm. Particular attention is paid to a comparison of algorithms for the automated selection of two tuning parameters. In addition to selection algorithms from the literature, we introduce a pseudospectra based approach for the particular case of the matrix exponential. Finally, applications involving both pathological matrices and the numerical solution of differential equations highlight the utility of the method.


65R10 Numerical methods for integral transforms
44A10 Laplace transform
35Q60 PDEs in connection with optics and electromagnetic theory
35K15 Initial value problems for second-order parabolic equations


Zbl 0141.33401
Full Text: DOI