Partager cette page :

# Linear algebraic aspects of exponential integrators

le 12 décembre 2012

14H - Groupe de travail "Applications des Mathématiques"

ENS Rennes Bâtiment Sauvy, Salle 5 (rdc)

Séminaire de Antti Koskela (University of Innsbruck) au groupe de travail "Applications des mathématiques"

We condsider the time integration of large stiff systems of ordinary differential equations, which result from semidiscretization of partial differential equations and which are of the form $u'(t) = Au(t)+g(u(t))$. We introduce so called exponential Taylor methods, which can be obtained by Taylor series expansion of the nonlinear part at each numerical approximation. When implementing exponential integrators, the size of the operator plays however a crucial role, as one needs to evaluate the action of a matrix exponential (and closely related functions) on a vector. We discuss the use of Krylov subspace methods for these problems. We introduce a so called moment-matching Arnoldi iteration for computing series of the form $\exp(hA)u_n+ \sum_{k=1}^p h^k \varphi_k(hA)w_k$, where the matrix functions $\varphi_k$ are related to the exponential function. Using Cauchy's integral formula we give a representation for the error of the moment-matching approximation and derive a priori error bounds. In addition an efficient a posteriori estimate is derived.

Thématique(s)
Recherche - Valorisation
Contact
Erwan Faou et Yannick Privat

Mise à jour le 4 décembre 2012