English

Computing the Action of Trigonometric and Hyperbolic Matrix Functions

Numerical Analysis 2017-05-02 v2

Abstract

We derive a new algorithm for computing the action f(A)Vf(A)V of the cosine, sine, hyperbolic cosine, and hyperbolic sine of a matrix AA on a matrix VV, without first computing f(A)f(A). The algorithm can compute cos(A)V\cos(A)V and sin(A)V\sin(A)V simultaneously, and likewise for cosh(A)V\cosh(A)V and sinh(A)V\sinh(A)V, and it uses only real arithmetic when AA is real. The algorithm exploits an existing algorithm \texttt{expmv} of Al-Mohy and Higham for eAV\mathrm{e}^AV and its underlying backward error analysis. Our experiments show that the new algorithm performs in a forward stable manner and is generally significantly faster than alternatives based on multiple invocations of \texttt{expmv} through formulas such as cos(A)V=(eiAV+eiAV)/2\cos(A)V = (\mathrm{e}^{\mathrm{i}A}V + \mathrm{e}^{\mathrm{-i}A}V)/2.

Cite

@article{arxiv.1607.04012,
  title  = {Computing the Action of Trigonometric and Hyperbolic Matrix Functions},
  author = {Nicholas J. Higham and Peter Kandolf},
  journal= {arXiv preprint arXiv:1607.04012},
  year   = {2017}
}
R2 v1 2026-06-22T14:54:21.596Z