Loading…
Methods for the approximation of the matrix exponential in a Lie‐algebraic setting
Discretization methods for ordinary differential equations based on the use of matrix exponentials have been known for decades. This set of ideas has come off age and acquired greater interest recently, within the context of geometric integration and discretization methods on manifolds based on the...
Saved in:
Published in: | IMA journal of numerical analysis 2001-04, Vol.21 (2), p.463-488 |
---|---|
Main Authors: | , |
Format: | Article |
Language: | English |
Subjects: | |
Citations: | Items that cite this one |
Online Access: | Get full text |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Summary: | Discretization methods for ordinary differential equations based on the use of matrix exponentials have been known for decades. This set of ideas has come off age and acquired greater interest recently, within the context of geometric integration and discretization methods on manifolds based on the use of Lie‐group actions. In the present paper we study the approximation of the matrix exponential in a particular context: given a Lie group G and its Lie algebra g, we seek approximants F(t B) of exp(t B) such that F(t B) ∈ G if B ∈ g. Having fixed a basis V1, …, Vd of g, we write F(t B) as a composition of exponentials of the type exp(αi (t) Vi), where αi for i = 1, 2, …, d are scalar functions. In this manner it becomes possible to increase the order of the approximation without increasing the number of exponentials to evaluate and multiply together. We study order conditions and implementation details and conclude the paper with some numerical experiments. |
---|---|
ISSN: | 0272-4979 1464-3642 |
DOI: | 10.1093/imanum/21.2.463 |