High performance computing of the matrix exponential

Abstract

This work presents a new algorithm for matrix exponential computation that significantly simplifies a Taylor scaling and squaring algorithm presented previously by the authors, preserving accuracy. A Matlab version of the new simplified algorithm has been compared with the original algorithm, providing similar results in terms of accuracy, but reducing processing time. It has also been compared with two state-of-the-art implementations based on Fade approximations, one commercial and the other implemented in Matlab, getting better accuracy and processing time results in the majority of cases. (C) 2015 Elsevier B.V. All rights reserved.