Complexity theory for Lie-group solvers


E. Celledoni, A. Iserles, S.P. Nørsett and B. Orel

submitted to J. Complexity

Abstract:
Commencing with a brief survey of Lie-group theory and differential equations evolving on Lie groups, we describe a number of numerical algorithms designed to respect Lie-group structure: Runge--Kutta--Munthe-Kaas schemes, Fer and Magnus expansions. This is followed by complexity analysis of Fer and Magnus expansions, whose conclusion is that for order four, six and eight an appropriately discretized Magnus method is always cheaper than a Fer method of the same order.
Each Lie-group method of the kind surveyed in this paper requires the computation of a matrix exponential. Classical methods, e.g.\ Krylov-subspace and rational approximants, may fail to map elements in a Lie algebra to a Lie group. Therefore we survey a number of approximants based on the splitting approach and demonstrate that their cost is compatible (and often superior) to classical methods.
Submitted by ai@damtp.cam.ac.uk Fri, 15 Sep 2000

Email of authors:
elenac@math.ntnu.no
ai@damtp.cam.ac.uk
norsett@math.ntnu.no
bojan@rosika.fri.uni-lj.si

URL of authors:
http://www.math.ntnu.no/~elenac/
http://www.damtp.cam.ac.uk/user/na/people/Arieh/
http://www.math.ntnu.no/~norsett/

Download:
Postscript-format Compressed-postscript-format
1999-006