Approximating the exponential from a Lie algebra to a Lie group

Authors:
Elena Celledoni and Arieh Iserles

Journal:
Math. Comp. **69** (2000), 1457-1480

MSC (1991):
Primary 65D15; Secondary 22E99, 65F30

Published electronically:
March 15, 2000

MathSciNet review:
1709149

Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Consider a differential equation with and , where is a Lie algebra of the matricial Lie group . Every can be mapped to by the matrix exponential map with .

Most numerical methods for solving ordinary differential equations (ODEs) on Lie groups are based on the idea of representing the approximation of the exact solution , , by means of exact exponentials of suitable elements of the Lie algebra, applied to the initial value . This ensures that .

When the exponential is difficult to compute exactly, as is the case when the dimension is large, an approximation of plays an important role in the numerical solution of ODEs on Lie groups. In some cases rational or polynomial approximants are unsuitable and we consider alternative techniques, whereby is approximated by a product of simpler exponentials.

In this paper we present some ideas based on the use of the Strang splitting for the approximation of matrix exponentials. Several cases of and are considered, in tandem with general theory. Order conditions are discussed, and a number of numerical experiments conclude the paper.

**1.**R. V. Chacon and A. T. Fomenko,*Recursion formulas for the Lie integral*, Adv. Math.**88**(1991), no. 2, 200–257. MR**1120613**, 10.1016/0001-8708(91)90008-U**2.**P. E. Crouch and R. Grossman,*Numerical integration of ordinary differential equations on manifolds*, J. Nonlinear Sci.**3**(1993), no. 1, 1–33. MR**1216986**, 10.1007/BF02429858**3.**Feng Kang and Zai Jiu Shang,*Volume-preserving algorithms for source-free dynamical systems*, Numer. Math.**71**(1995), no. 4, 451–463. MR**1355050**, 10.1007/s002110050153**4.**Etienne Forest,*Sixth-order Lie group integrators*, J. Comput. Phys.**99**(1992), no. 2, 209–213. MR**1158207**, 10.1016/0021-9991(92)90203-B**5.**E. Gallopoulos and Y. Saad,*Efficient solution of parabolic equations by Krylov approximation methods*, SIAM J. Sci. Statist. Comput.**13**(1992), no. 5, 1236–1264. MR**1177807**, 10.1137/0913071**6.**Marlis Hochbruck and Christian Lubich,*On Krylov subspace approximations to the matrix exponential operator*, SIAM J. Numer. Anal.**34**(1997), no. 5, 1911–1925. MR**1472203**, 10.1137/S0036142995280572**7.**A. Iserles and S. P. Nørsett.

On the Solution of Linear Differential Equations in Lie Groups. Roy. Soc. London Philos. Trans. Ser. A Math. Phys. Eng. Sci.**357**(1999), 983-1019. CMP**99:13****8.**Robert I. McLachlan,*On the numerical integration of ordinary differential equations by symmetric composition methods*, SIAM J. Sci. Comput.**16**(1995), no. 1, 151–168. MR**1311683**, 10.1137/0916010**9.**Hans Munthe-Kaas,*Runge-Kutta methods on Lie groups*, BIT**38**(1998), no. 1, 92–111. MR**1621080**, 10.1007/BF02510919**10.**Hans Munthe-Kaas and Antonella Zanna,*Numerical integration of differential equations on homogeneous manifolds*, Foundations of computational mathematics (Rio de Janeiro, 1997) Springer, Berlin, 1997, pp. 305–315. MR**1661989****11.**J. M. Sanz-Serna and M. P. Calvo,*Numerical Hamiltonian problems*, Applied Mathematics and Mathematical Computation, vol. 7, Chapman & Hall, London, 1994. MR**1270017****12.**M. Schatzman.

Higher order alternate directions methods.

Technical report, Université de Lyon, Laboratoire d'analyse numérique, 1997.**13.**Q. Sheng.*Solving partial differential equations by exponential splitting*.

PhD thesis, Cambridge University, 1989.**14.**Gilbert Strang,*Accurate partial difference methods. II. Non-linear problems*, Numer. Math.**6**(1964), 37–46. MR**0166942****15.**Masuo Suzuki,*General theory of fractal path integrals with applications to many-body theories and statistical physics*, J. Math. Phys.**32**(1991), no. 2, 400–407. MR**1088360**, 10.1063/1.529425**16.**V. S. Varadarajan,*Lie groups, Lie algebras, and their representations*, Graduate Texts in Mathematics, vol. 102, Springer-Verlag, New York, 1984. Reprint of the 1974 edition. MR**746308****17.**Haruo Yoshida,*Construction of higher order symplectic integrators*, Phys. Lett. A**150**(1990), no. 5-7, 262–268. MR**1078768**, 10.1016/0375-9601(90)90092-3**18.**A. Zanna.

The method of iterated commutators for ordinary differential equations on Lie groups.

Technical Report 1996/NA12, Department of Applied Mathematics and Theoretical Physics, University of Cambridge, England, 1996.

Retrieve articles in *Mathematics of Computation*
with MSC (1991):
65D15,
22E99,
65F30

Retrieve articles in all journals with MSC (1991): 65D15, 22E99, 65F30

Additional Information

**Elena Celledoni**

Affiliation:
DAMTP, Cambridge University, Silver Street, England CB3 9EW

Address at time of publication:
Department of Mathematical Sciences, NTNU 7491 Trondheim, Norway

Email:
elenac@math.ntnu.no

**Arieh Iserles**

Affiliation:
DAMTP, Cambridge University, Silver Street, England CB3 9EW

Email:
A.Iserles@damtp.cam.ac.uk

DOI:
http://dx.doi.org/10.1090/S0025-5718-00-01223-0

Received by editor(s):
January 1, 1998

Received by editor(s) in revised form:
October 27, 1998

Published electronically:
March 15, 2000

Article copyright:
© Copyright 2000
American Mathematical Society