Matrix Inverse Trigonometric and Inverse Hyperbolic Functions: Theory and Algorithms

Research output: Contribution to journalArticle

Abstract

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 the 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 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

Bibliographical metadata

Original languageEnglish
Pages (from-to)1453-1477
Number of pages25
JournalSIAM Journal on Matrix Analysis and Applications
Volume37
Issue number4
DOIs
StatePublished - 8 Nov 2016