Approximating the exponential from a Lie algebra to a Lie group

被引:53
作者
Celledoni, E [1 ]
Iserles, A [1 ]
机构
[1] Univ Cambridge, Dept Appl Math & Theoret Phys, Cambridge CB3 9EW, England
关键词
D O I
10.1090/S0025-5718-00-01223-0
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
Consider a differential equation y' = A(t,y)y, y(0) = y(0) with y(0) is an element of G and A : R+ x G --> g, where a is a Lie algebra of the matricial Lie group G. Every B is an element of g can be mapped to G by the matrix exponential map exp (tB) with t is an element of B. Most numerical methods for solving ordinary differential equations (ODEs) on Lie groups are based on the idea of representing the approximation y(n) of the exact solution y(t(n)) t, is an element of R+, by means of exact exponentials of suitable elements of the Lie algebra, applied to the initial value y(0) This ensures that y(n) is an element of G. When the exponential is difficult to compute exactly, as is the case when the dimension is large, an approximation of exp (tB) 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 exp(tB) 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 B and G are considered, in tandem with general theory. Order conditions are discussed, and a number of numerical experiments conclude the paper.
引用
收藏
页码:1457 / 1480
页数:24
相关论文
共 18 条
[1]   RECURSION FORMULAS FOR THE LIE INTEGRAL [J].
CHACON, RV ;
FOMENKO, AT .
ADVANCES IN MATHEMATICS, 1991, 88 (02) :200-257
[2]   NUMERICAL-INTEGRATION OF ORDINARY DIFFERENTIAL-EQUATIONS ON MANIFOLDS [J].
CROUCH, PE ;
GROSSMAN, R .
JOURNAL OF NONLINEAR SCIENCE, 1993, 3 (01) :1-33
[3]  
FENG K, 1995, NUMER MATH, V71, P451
[4]   6TH-ORDER LIE GROUP INTEGRATORS [J].
FOREST, E .
JOURNAL OF COMPUTATIONAL PHYSICS, 1992, 99 (02) :209-213
[5]   EFFICIENT SOLUTION OF PARABOLIC EQUATIONS BY KRYLOV APPROXIMATION METHODS [J].
GALLOPOULOS, E ;
SAAD, Y .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1992, 13 (05) :1236-1264
[6]   On Krylov subspace approximations to the matrix exponential operator [J].
Hochbruck, M ;
Lubich, C .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1997, 34 (05) :1911-1925
[7]   On the solution of linear differential equations in Lie groups [J].
Iserles, A ;
Norsett, SP .
PHILOSOPHICAL TRANSACTIONS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 1999, 357 (1754) :983-1019
[8]   ON THE NUMERICAL-INTEGRATION OF ORDINARY DIFFERENTIAL-EQUATIONS BY SYMMETRICAL COMPOSITION METHODS [J].
MCLACHLAN, RI .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 1995, 16 (01) :151-168
[9]   Runge-Kutta methods on Lie groups [J].
Munthe-Kaas, H .
BIT, 1998, 38 (01) :92-111
[10]  
SANZSERNA JM, 1994, AMMC 7