You are here: MIMS > EPrints
MIMS EPrints

2016.4: Matrix Inverse Trigonometric and Inverse Hyperbolic Functions: Theory and Algorithms

2016.4: Mary Aprahamian and Nicholas J. Higham (2016) Matrix Inverse Trigonometric and Inverse Hyperbolic Functions: Theory and Algorithms. SIAM Journal on Matrix Analysis and Applications, 37 (4). pp. 1453-1477. ISSN 1095-7162

This is the latest version of this eprint.

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
459 Kb

DOI: 10.1137/080716426


Theoretical and computational aspects of matrix inverse trigonometric and inverse hyperbolic functions are studied. Conditions for existence are given, all possible values are characterized, and principal values $\acos$, $\asin$, $\acosh$, and $\asinh$ are defined and shown to be unique primary matrix functions. Various functional identities are derived, some of which are new even in the scalar case, with care taken to specify precisely the choices of signs and branches. New results include a ``round trip'' formula that relates $\mathrm{acos}(\cos A)$ to $A$ and similar formulas for the other inverse functions. Key tools used in the derivations are the matrix unwinding function and the matrix sign function. A new inverse scaling and squaring type algorithm employing a Schur decomposition and variable-degree Pad\'e approximation is derived for computing acos, and it is shown how it can also be used to compute asin, acosh, and asinh. In numerical experiments the algorithm is found to behave in a forward stable fashion and to be superior to computing these functions via logarithmic formulas.

Item Type:Article
Uncontrolled Keywords:matrix function, inverse trigonometric functions, inverse hyperbolic functions, matrix inverse sine, matrix inverse cosine, matrix inverse hyperbolic sine, matrix inverse hyperbolic cosine, matrix exponential, matrix logarithm, matrix sign function, rational approximation, Pad\'{e} approximation, MATLAB, GNU Octave, Fr\'echet derivative, condition number
Subjects:MSC 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2016.4
Deposited By:Nick Higham
Deposited On:19 October 2016

Available Versions of this Item

Download Statistics: last 4 weeks
Repository Staff Only: edit this item