An Arbitrary Precision Scaling and Squaring Algorithm for the Matrix Exponential

Fasi, Massimiliano and Higham, Nicholas J. (2018) An Arbitrary Precision Scaling and Squaring Algorithm for the Matrix Exponential. [MIMS Preprint]

WarningThere is a more recent version of this item available.
[img] Text

Download (569kB)


The most popular algorithms for computing the matrix exponential are those based on the scaling and squaring technique. For optimal efficiency these are usually tuned to a particular precision of floating-point arithmetic. We design a new scaling and squaring algorithm that takes the unit roundoff of the arithmetic as input and chooses the algorithmic parameters in order to keep the forward error in the underlying Padé approximation below the unit roundoff. To do so, we derive an explicit expression for all the coefficients in an error expansion for Padé approximants to the exponential and use it to obtain a new bound for the truncation error. We also derive a new technique for selecting the internal parameters used by the algorithm, which at each step decides whether to scale or to increase the degree of the approximant. The algorithm can employ diagonal Padé approximants or Taylor approximants and can be used with a Schur decomposition or in transformation-free form. Our numerical experiments show that the new algorithm performs in a forward stable way for a wide range of precisions and that the most accurate of our implementations, the Taylor-based transformation-free variant, is superior to existing alternatives.

Item Type: MIMS Preprint
Subjects: MSC 2010, the AMS's Mathematics Subject Classification > 15 Linear and multilinear algebra; matrix theory
MSC 2010, the AMS's Mathematics Subject Classification > 65 Numerical analysis
Depositing User: Massimiliano Fasi
Date Deposited: 25 Nov 2018 09:30
Last Modified: 25 Nov 2018 09:30

Available Versions of this Item

Actions (login required)

View Item View Item