Complexity theory for lie-group solvers

  • Authors:
  • Elena Celledoni;Arieh Iserles;Syvert P. Nørsett;Bojan Orel

  • Affiliations:
  • Institute of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim NTH-7074, Norway;Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Silver Street, Cambridge CB3 9EW, United Kingdom;Institute of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim NTH-7074, Norway;Department of Mathematics and Mechanics, University of Ljubljana, Jadranska 19, 1000 Ljubljana, Slovenia

  • Venue:
  • Journal of Complexity
  • Year:
  • 2002

Quantified Score

Hi-index 0.00

Visualization

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 derivation of the computational cost 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.