Related papers: Higher Order Variational Integrators: a polynomial…

Explicit symplectic integrators have been important tools for accurate and efficient approximations of mechanical systems with separable Hamiltonians. For the first time, the article proposes for arbitrary Hamiltonians similar integrators,…

We consider Hamiltonian systems driven by multi-dimensional Gaussian processes in rough path sense, which include fractional Brownian motions with Hurst parameter $H\in(1/4,1/2]$. We indicate that the phase flow preserves the symplectic…

Geometric integration of non-autonomous classical engineering problems, such as rotor dynamics, is investigated. It is shown, both numerically and by backward error analysis, that geometric (structure preserving) integration algorithms are…

We study the problem of efficient integration of variational equations in multi-dimensional Hamiltonian systems. For this purpose, we consider a Runge-Kutta-type integrator, a Taylor series expansion method and the so-called `Tangent Map'…

Based on reasonable testing model problems, we study the preservation by symplectic Runge-Kutta method (SRK) and symplectic partitioned Runge-Kutta method (SPRK) of structures for fixed points of linear Hamiltonian systems. The…

We consider a Lagrangian system $L(q,\dot q) = \sum_{l=1}^{N}L^{\{l\}}(q,\dot q)$, where the $q$-variable is treated by a Generalized Additive Runge--Kutta (GARK) method. Applying the technique of discrete variations, we show how to…

The design of numerical integrators for solving stochastic dynamics with high weak order relies on tedious calculations and is subject to a high number of order conditions. The original approaches from the literature consider strong…

Runge-Kutta methods are affine equivariant: applying a method before or after an affine change of variables yields the same numerical trajectory. However, for some applications, one would like to perform numerical integration after a…

The paper proposes a scheme by combining the Runge-Kutta discontinuous Galerkin method with a {\delta}-mapping algorithm for solving hyperbolic conservation laws with discontinuous fluxes. This hybrid scheme is particularly applied to…

Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent…

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.…

Retractions maps are used to define a discretization of the tangent bundle of the configuration manifold as two copies of the configuration manifold where the dynamics take place. Such discretization maps can be conveniently lifted to a…

In this work we propose a new numerical approach to distinguish between regular and chaotic orbits in Hamiltonian systems, based on the simultaneous integration of both the orbit and the deviation vectors using a symplectic scheme, hereby…

Discrete Hamiltonian variational integrators are derived from Type II and Type III generating functions for symplectic maps, and in this paper we establish a variational error analysis result that relates the order of accuracy of the…

Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These…

Computer simulations in QCD are based on the discretization of the theory on a Euclidean lattice. To compute the mean value of an observable, usually the Hybrid Monte Carlo method is applied. Here equations of motion, derived from an…

Two families of symplectic methods specially designed for second-order time-dependent linear systems are presented. Both are obtained from the Magnus expansion of the corresponding first-order equation, but otherwise they differ in…

We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian…

Multi-symplectic integrators are typically regarded as a discretization of the Hamiltonian partial differential equations. This is due to the fact that, for generic finite-dimensional Hamiltonian systems, there exists only one independent…

Modified Hamiltonian Monte Carlo (MHMC) methods combine the ideas behind two popular sampling approaches: Hamiltonian Monte Carlo (HMC) and importance sampling. As in the HMC case, the bulk of the computational cost of MHMC algorithms lies…