Related papers: Geometric numerical integrators for Hunter-Saxton-…
Variational integrators have traditionally been constructed from the perspective of Lagrangian mechanics, but there have been recent efforts to adopt discrete variational approaches to the symplectic discretization of Hamiltonian mechanics…
Numerical evolution of time-dependent differential equations via explicit Runge-Kutta or Taylor methods typically fails to preserve symmetries of a system. It is known that there exists no numerical integration method that in general…
I introduce an innovative methodology for deriving numerical models of systems of partial differential equations which exhibit the evolution of spatial patterns. The new approach directly produces a discretisation for the evolution of the…
This work develops a symplectic framework for quantum computing to be applied to classical Hamiltonian systems, exploiting the intrinsic geometric compatibility between unitary quantum evolution and symplectic phase-space dynamics in a…
The numerical integration of the Benjamin and Benjamin--Ono equations are considered. They are non-local partial differential equations involving the Hilbert transform, and due to this, so far quite few structure-preserving integrators have…
We expose (without proofs) a unified computational approach to integrable structures (including recursion, Hamiltonian, and symplectic operators) based on geometrical theory of partial differential equations. We adopt a coordinate based…
Several dynamical systems of interest in celestial mechanics can be written in the form of a Newton equation with time-dependent damping, linear in the velocities. For instance, the modified Kepler problem, the spin-orbit model and the…
A supersymmetric extension of the Hunter-Saxton equation is constructed. We present its bi-Hamiltonian structure and show that it arises geometrically as a geodesic equation on the space of superdiffeomorphisms of the circle that leave a…
Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $\MM{u}$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical…
We compare the performances of symplectic and non-symplectic integrators for the computation of normal geodesics and conjugate points in sub-Riemannian geometry at the example of the Martinet case. For this case study we consider first the…
Optimization tasks are crucial in statistical machine learning. Recently, there has been great interest in leveraging tools from dynamical systems to derive accelerated and robust optimization methods via suitable discretizations of…
Geometric discretizations that preserve certain Hamiltonian structures at the discrete level has been proven to enhance the accuracy of numerical schemes. In particular, numerous symplectic and multi-symplectic schemes have been proposed to…
This paper surveys in detail the relations between numerical integration and the Hamiltonian (or hybrid) Monte Carlo method (HMC). Since the computational cost of HMC mainly lies in the numerical integrations, these should be performed as…
This article reviews some integrators particularly suitable for the numerical resolution of differential equations on a large time interval. Symplectic integrators are presented. Their stability on exponentially large time is shown through…
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.…
We give a theoretical framework of stochastic non-canonical Hamiltonian systems as well as their modified symplectic structure which is named stochastic K-symplectic structure. The framework can be applied to the study of the…
We propose a new numerical method for $\alpha$-dissipative solutions of the Hunter-Saxton equation, where $\alpha$ belongs to $W^{1, \infty}(\mathbb{R}, [0, 1))$. The method combines a projection operator with a generalized method of…
The Camassa-Holm (CH) and Hunter-Saxton (HS) equations have an interpretation as geodesic flow equations on the group of diffeomorphisms, preserving the $H^1$ and $\dot{H}^1$ right-invariant metrics correspondingly. There is an analogy to…
The existence of explicit symplectic integrators for general nonseparable Hamiltonian systems is an open and important problem in both numerical analysis and computing in science and engineering, as explicit integrators are usually more…
We present a general framework for constructing structure-preserving numerical integrators for nonholonomically constrained mechanical systems evolving on Lie groups using retraction maps. Retraction maps generalize the exponential map and…