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 of the cosine, sine, hyperbolic cosine, and hyperbolic sine of a matrix on a matrix , without first computing . The algorithm can compute and simultaneously, and likewise for and , and it uses only real arithmetic when is real. The algorithm exploits an existing algorithm \texttt{expmv} of Al-Mohy and Higham for 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 .
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}
}