Related papers: Semiexplicit Symplectic Integrators for Non-separa…
Numerical methods for developing port-Hamiltonian representations of general linear time-invariant systems are studied. The approach extends previous port-Hamiltonian characterizations to include the general non-minimal case and the case…
We describe a technique for constructing a symplectic transfer map for a charged particle moving through an accelerator component with arbitrary three-dimensional static magnetic field. The transfer map is constructed by symplectic…
We propose an effective and lightweight learning algorithm, Symplectic Taylor Neural Networks (Taylor-nets), to conduct continuous, long-term predictions of a complex Hamiltonian dynamic system based on sparse, short-term observations. At…
Exponential integrators are a well-known class of time integration methods that have been the subject of many studies and developments in the past two decades. Surprisingly, there have been limited efforts to analyze their stability and…
The construction of modified equations is an important step in the backward error analysis of symplectic integrator for Hamiltonian systems. In the context of partial differential equations, the standard construction leads to modified…
We propose explicit symplectic integrators of molecular dynamics (MD) algorithms for rigid-body molecules in the canonical and isothermal-isobaric ensembles. We also present a symplectic algorithm in the constant normal pressure and lateral…
New classes of Lie-Hamilton systems are obtained from the six-dimensional fundamental representation of the symplectic Lie algebra $\mathfrak{sp}(6,\mathbb{R})$. The ansatz is based on a recently proposed procedure for constructing…
A general high-order fully explicit scheme based on projective integration methods is here presented to solve systems of degenerate parabolic equations in general dimensions. The method is based on a BGK approximation of the…
We present a new method to solve in a semianalytical way the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations at NLO order in the x-space. The method allows to construct an evolution operator expressed in form of a rapidly…
Symplectic integrators have made it possible to study the long-term evolution of planetary systems with direct N-body simulations. In this paper we reassess the accuracy of such simulations by running a convergence test on 20Myr…
The aim of this work is to apply a semi-implicit (SI) strategy within a Rosenbrock-type and IMEX linear multistep (LM) framework to a sequence of 1D time-dependent partial differential equations (PDEs) with high order spatial derivatives.…
On the basis of the previous work by Tang \& Zhang (Appl. Math. Comput. 323, 2018, p. 204--219), in this paper we present a more effective way to construct high-order symplectic integrators for solving second order Hamiltonian equations.…
A variational framework for accelerated optimization was recently introduced on normed vector spaces and Riemannian manifolds in Wibisono et al. (2016) and Duruisseaux and Leok (2021). It was observed that a careful combination of…
In this paper we present a method of constructing a nonlinear accelerator lattice that has an approximate integral of motion that is given upfront. The integral under consideration is a Hamiltonian in normalized (canonical) coordinates that…
We show that symplectic Runge-Kutta methods provide effective symplectic integrators for Hamiltonian systems with index one constraints. These include the Hamiltonian description of variational problems subject to position and velocity…
An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which…
When the singular values of the evolution operator are all smaller or all greater than one, stable integration algorithms are obtained either by explicit or implicit methods. When the singular spectrum mixes greater and smaller than one…
In this paper we consider the numerical solution of the Hamiltonian wave equation in two spatial dimension. We use the Mimetic Finite Difference (MFD) method to approximate the continuous problem combined with a symplectic integration in…
We present a practical algorithm based on symplectic splitting methods to integrate numerically in time the Schr\"odinger equation. When discretized in space, the Schr\"odinger equation can be recast as a classical Hamiltonian system…
This paper is concerned with the strong approximation of a semi-linear stochastic wave equation with strong damping, driven by additive noise. Based on a spatial discretization performed by a spectral Galerkin method, we introduce a kind of…