Geometric integration methods of high order for gradient systems

Abstract

We propose the higher order extension of numerical methods that preserve the energy-diminishing feature of a system with Lyapunov function.