Computing \(exp(-\tau A) b\) with Laguerre polynomials
Electronic transactions on numerical analysis, Tome 37 (2010), pp. 147-165
This paper discusses a method based on Laguerre polynomials combined with a Filtered Conjugate Residual (FCR) framework to compute the product of the exponential of a matrix by a vector. The method implicitly uses an expansion of the exponential function in a series of orthogonal Laguerre polynomials, much like existing methods based on Chebyshev polynomials do. Owing to the fact that orthogonal polynomials satisfy a three-term recurrence, what these series expansion methods offer over other approaches such as Krylov subspace methods lies in the elimination of inner products and the economy in storage since there is no need to compute and keep a set of basis vectors. Compared with Chebyshev polynomials that are orthogonal within a restricted interval and need estimates of the outermost eigenvalues, Laguerre polynomials offer the added feature that they are orthogonal on the half real line, alleviating therefore the need to estimate eigenvalues.
Classification :
65F50, 65L05, 41A10
Keywords: conjugate residual, filtered conjugate residual, polynomial filtering, exponential propagation, orthogonal polynomials
Keywords: conjugate residual, filtered conjugate residual, polynomial filtering, exponential propagation, orthogonal polynomials
@article{ETNA_2010__37__a17,
author = {Sheehan, Bernard N. and Saad, Yousef and Sidje, Roger B.},
title = {Computing \(exp(-\tau {A)} b\) with {Laguerre} polynomials},
journal = {Electronic transactions on numerical analysis},
pages = {147--165},
year = {2010},
volume = {37},
zbl = {1205.65165},
language = {en},
url = {http://geodesic.mathdoc.fr/item/ETNA_2010__37__a17/}
}
TY - JOUR AU - Sheehan, Bernard N. AU - Saad, Yousef AU - Sidje, Roger B. TI - Computing \(exp(-\tau A) b\) with Laguerre polynomials JO - Electronic transactions on numerical analysis PY - 2010 SP - 147 EP - 165 VL - 37 UR - http://geodesic.mathdoc.fr/item/ETNA_2010__37__a17/ LA - en ID - ETNA_2010__37__a17 ER -
Sheehan, Bernard N.; Saad, Yousef; Sidje, Roger B. Computing \(exp(-\tau A) b\) with Laguerre polynomials. Electronic transactions on numerical analysis, Tome 37 (2010), pp. 147-165. http://geodesic.mathdoc.fr/item/ETNA_2010__37__a17/