Title
The Scaling, Splitting, and Squaring Method for the Exponential of Perturbed Matrices.
Abstract
We propose splitting methods for the computation of the exponential of perturbed matrices which can be written as the sum A = D+epsilon B of a sparse and efficiently exponentiable matrix D with sparse exponential e(D) and a dense matrix epsilon B which is of small norm in comparison with D. The predominant algorithm is based on scaling the large matrix A by a small number 2(-s), which is then exponentiated by efficient Pade or Taylor methods and finally squared in order to obtain an approximation for the full exponential. In this setting, the main portion of the computational cost arises from dense-matrix multiplications and we present a modified squaring which takes advantage of the smallness of the perturbed matrix B in order to reduce the number of squarings necessary. Theoretical results on local error and error propagation for splitting methods are complemented with numerical experiments and show a clear improvement over existing methods when medium precision is sought.
Year
DOI
Venue
2015
10.1137/14098003X
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS
Keywords
Field
DocType
matrix exponential,scaling and squaring method,splitting method,Pade approximation,backward error analysis
Mathematical optimization,Square (algebra),Exponential function,Padé approximant,Matrix (mathematics),Mathematical analysis,Matrix exponential,Scaling,Sparse matrix,Matrix splitting,Mathematics
Journal
Volume
Issue
ISSN
36
2
0895-4798
Citations 
PageRank 
References 
2
0.39
4
Authors
3
Name
Order
Citations
PageRank
Philipp Bader1163.20
S. Blanes24210.47
Muaz Seydaoglu351.96