Boosting the computation of the matrix exponential

Abstract This paper presents new Taylor algorithms for the computation of the matrix exponential based on recent new matrix polynomial evaluation methods. Those methods are more efficient than the well known Paterson–Stockmeyer method. The cost of the proposed algorithms is reduced with respect to previous algorithms based on Taylor approximations. Tests have been performed to compare the MATLAB implementations of the new algorithms to a state-of-the-art Pade algorithm for the computation of the matrix exponential, providing higher accuracy and cost performances.

[1]  Emilio Defez,et al.  Two algorithms for computing the matrix cosine function , 2017, Appl. Math. Comput..

[2]  J. Sastre.  Efficient evaluation of matrix polynomials , 2018 .

[3]  Emilio Defez,et al.  A new efficient and accurate spline algorithm for the matrix exponential computation , 2018, J. Comput. Appl. Math..

[4]  Emilio Defez,et al.  High performance computing of the matrix exponential , 2016, J. Comput. Appl. Math..

[5]  Emilio Defez,et al.  Accurate matrix exponential computation to solve coupled differential models in engineering , 2011, Math. Comput. Model..

[6]  Nicholas J. Higham,et al.  The Scaling and Squaring Method for the Matrix Exponential Revisited , 2005, SIAM J. Matrix Anal. Appl..

[7]  Nicholas J. Higham,et al.  Functions of matrices - theory and computation , 2008 .

[8]  Awad H. Al-Mohy,et al.  A New Scaling and Squaring Algorithm for the Matrix Exponential , 2009, SIAM J. Matrix Anal. Appl..

[9]  Jorge J. Moré,et al.  Digital Object Identifier (DOI) 10.1007/s101070100263 , 2001 .

[10]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[11]  Emilio Defez,et al.  Accurate and efficient matrix exponential computation , 2014, Int. J. Comput. Math..

[12]  Larry J. Stockmeyer,et al.  On the Number of Nonscalar Multiplications Necessary to Evaluate Polynomials , 1973, SIAM J. Comput..

[13]  Nicholas J. Higham,et al.  FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation , 1988, TOMS.