High-Performance Computation of the Exponential of a Large Sparse Matrix
Feng Wu,Kailing Zhang,Li Zhu,Jiayao Hu
DOI: https://doi.org/10.1137/20M1342987
IF: 1.908
2021-12-08
SIAM Journal on Matrix Analysis and Applications
Abstract:SIAM Journal on Matrix Analysis and Applications, Volume 42, Issue 4, Page 1636-1655, January 2021. Computation of the large sparse matrix exponential has been an important topic in many fields, such as network and finite-element analysis. The existing scaling and squaring algorithm (SSA) is not suitable for the computation of the large sparse matrix exponential as it requires more memory and has a higher computational cost than is actually needed. By introducing two novel concepts, i.e., real bandwidth and bandwidth, to measure the sparsity of the matrix, the sparsity of the matrix exponential is analyzed. It is found that for every matrix computed in the squaring phase of the SSA, a corresponding sparse approximate matrix exists. To obtain the sparse approximate matrix, a new filtering technique based on forward error analysis is proposed. Combining the filtering technique with the idea of keeping track of the incremental part, a competitive algorithm is developed for the large sparse matrix exponential. Due to the filtering technique, the proposed method can greatly alleviate the overscaling problem. Three sets of numerical experiments, including one large matrix with a dimension larger than , are conducted. The numerical experiments show that, compared with the expm function in MATLAB, the proposed algorithm can provide higher accuracy at lower computational cost and with less memory.
mathematics, applied