کد مقاله | کد نشریه | سال انتشار | مقاله انگلیسی | نسخه تمام متن |
---|---|---|---|---|
471757 | 1364829 | 2016 | 10 صفحه PDF | دانلود رایگان |
In this paper we present numerical methods for computing the matrix functions arising in exponential integrators. The matrix functions are the linear combination of the form ∑j=0pφj(A)bj, where φjφj is related to exponential function. By reducing the corresponding systems of linear ODEs with polynomial inhomogeneity based on block Krylov subspace methods, two effective exponentials of low dimension matrix are derived to approximate the objective functions. Some useful a posteriori error estimates are established. The algorithms can be combined with time-stepping strategy to satisfy the accuracy requirements under a moderate dimension of the block Krylov subspace. Several numerical experiments are presented to demonstrate the efficiency of the methods.
Journal: Computers & Mathematics with Applications - Volume 72, Issue 4, August 2016, Pages 846–855