Fast Taylor polynomial evaluation for the computation of the matrix cosine

Abstract In this work we introduce a new method to compute the matrix cosine. It is based on recent new matrix polynomial evaluation methods for the Taylor approximation and a mixed forward and backward error analysis. The matrix polynomial evaluation methods allow to evaluate the Taylor polynomial approximation of the matrix cosine function more efficiently than using Paterson–Stockmeyer method. A sequential Matlab implementation of the new algorithm is provided, giving better efficiency and accuracy than state-of-the-art algorithms. Moreover, we provide an implementation in Matlab that can use NVIDIA GPUs easily and efficiently.

[1]  Nicholas J. Higham,et al.  Efficient algorithms for the matrix cosine and sine , 2005, Numerical Algorithms.

[2]  Emilio Defez,et al.  Efficient and accurate algorithms for computing matrix trigonometric functions , 2017, J. Comput. Appl. Math..

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

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

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

[6]  Awad H. Al-Mohy,et al.  New Algorithms for Computing the Matrix Sine and Cosine Separately or Simultaneously , 2015, SIAM J. Sci. Comput..

[7]  Pedro Alonso,et al.  Computing matrix trigonometric functions with GPUs through Matlab , 2018, The Journal of Supercomputing.

[8]  Emilio Defez,et al.  Computing matrix functions arising in engineering models with orthogonal matrix polynomials , 2013, Math. Comput. Model..

[9]  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.

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

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

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

[13]  Emilio Defez,et al.  Efficient computation of the matrix cosine , 2013, Appl. Math. Comput..

[14]  Emilio Defez,et al.  Boosting the computation of the matrix exponential , 2019, Appl. Math. Comput..