6

Are there high order energy-conserving or symplectic methods for solving $y'=f(y)$?

Maesumi
  • 163
  • 5

1 Answers1

9

There are the geometric integrations written by Ernst Hairer & co:

  • E. Hairer, C. Lubich and G. Wanner (2002): Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential equations. Springer Series in Comput. Math., vol. 31

  • E. Hairer and M. Hairer (2002): GniCodes - Matlab programs for geometric numerical integration

for which you can download the routines in FORTRAN77, Matlab or C++ from Professor Hairer's website. The Implicit Runge-Kutta (IRK) scheme he has implemented has order 8.

GertVdE
  • 6,179
  • 1
  • 21
  • 36