Time-stepping algorithms for semidiscretized linear parabolic PDEs based on rational approximants with distinct real poles

Time dependent problems in Partial Differential Equations (PDEs) are often solved by the Method Of Lines (MOL). For linear parabolic PDEs, the exact solution of the resulting system of first order Ordinary Differential Equations (ODEs) satisfies a recurrence relation involving the matrix exponential function. In this paper, we consider the development of a fourth order rational approximant to the matrix exponential function possessing real and distinct poles which, consequently, readily admits a partial fraction expansion, thereby allowing the distribution of the work in solving the corresponding linear algebraic systems in essentially Backward Euler-like solves on concurrent processors. The resulting parallel algorithm possesses appropriate stability properties, and is implemented on various parabolic PDEs from the literature including the forced heat equation and the advection-diffusion equation.Dedicated to Professor J. Crank on the occasion of his 80th birthday