Fast symplectic integrator for Nesterov-type acceleration method

06/01/2021
by   Shin-itiro Goto, et al.
0

In this paper, explicit stable integrators based on symplectic and contact geometries are proposed for a non-autonomous ordinarily differential equation (ODE) found in improving convergence rate of Nesterov's accelerated gradient method. Symplectic geometry is known to be suitable for describing Hamiltonian mechanics, and contact geometry is known as an odd-dimensional counterpart of symplectic geometry. Moreover, a procedure, called symplectization, is a known way to construct a symplectic manifold from a contact manifold, yielding Hamiltonian systems from contact ones. It is found in this paper that a previously investigated non-autonomous ODE can be written as a contact Hamiltonian system. Then, by symplectization of a non-autonomous contact Hamiltonian vector field expressing the non-autonomous ODE, novel symplectic integrators are derived. Because the proposed symplectic integrators preserve hidden symplectic and contact structures in the ODE, they should be more stable than the Runge-Kutta method. Numerical experiments demonstrate that, as expected, the second-order symplectic integrator is stable and high convergence rates are achieved.

READ FULL TEXT
research
07/22/2022

Numerical integration of stochastic contact Hamiltonian systems via stochastic Herglotz variational principle

In this work we construct a stochastic contact variational integrator an...
research
05/08/2020

Geometric numerical integration of Lìenard systems via a contact Hamiltonian approach

Starting from a contact Hamiltonian description of Lìenard systems, we i...
research
09/05/2019

Numerical integration in celestial mechanics: a case for contact geometry

Several dynamical systems of interest in celestial mechanics can be writ...
research
02/12/2019

On the information content of the difference from hamiltonian evolution

A dissipative version of hamiltonian mechanics is proposed via a princip...
research
08/05/2020

A variable rest length impedance grasping strategy in the port-Hamiltonian framework

This work is devoted to an impedance grasping strategy for a class of st...
research
02/04/2021

The complex behaviour of Galton rank order statistic

Galton's rank order statistic is one of the oldest statistical tools for...

Please sign up or login with your details

Forgot password? Click here to reset