Related papers: High order explicit symplectic integrators for the…
The usual explicit finite-difference method of solving partial differential equations is limited in stability because it approximates the exact amplification factor by power-series. By adapting the same exponential-splitting method of…
We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators,…
We introduce a novel technique for constructing higher-order variational integrators for Hamiltonian systems of ODEs. In particular, we are concerned with generating globally smooth approximations to solutions of a Hamiltonian system. Our…
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly…
We consider the construction of semi-implicit linear multistep methods which can be applied to time dependent PDEs where the separation of scales in additive form, typically used in implicit-explicit (IMEX) methods, is not possible. As…
We introduce a family of implicit probabilistic integrators for initial value problems (IVPs), taking as a starting point the multistep Adams-Moulton method. The implicit construction allows for dynamic feedback from the forthcoming…
First-order fully implicit as well as implicit--explicit schemes for coupled elliptic-parabolic systems are discussed in [Ern and Meunier, ESAIM: M2AN, 2009] and [Altmann et al., Math.\ Comp., 2021], respectively. The extension of the…
In this paper, we consider the task of efficiently computing the numerical solution of evolutionary complex Ginzburg--Landau equations on Cartesian product domains with homogeneous Dirichlet/Neumann or periodic boundary conditions. To this…
We introduce a new non-resonant low-regularity integrator for the cubic nonlinear Schr\"odinger equation (NLSE) allowing for long-time error estimates which are optimal in the sense of the underlying PDE. The main idea thereby lies in…
We construct a new class of higher-dimensional column-vector loop algebras. Based on it, a method for generating higher-dimensional isospectral-nonisospectral integrable hierarchies is proposed. As an application, we derive a generalized…
We describe and test an easy-to-implement two-step high-order compact (2SHOC) scheme for the Laplacian operator and its implementation into an explicit finite-difference scheme for simulating the nonlinear Schr\"odinger equation (NLSE).…
Dozens of exponential integration formulas have been proposed for the high-accuracy solution of stiff PDEs such as the Allen-Cahn, Korteweg-de Vries and Ginzburg-Landau equations. We report the results of extensive comparisons in MATLAB and…
We introduce a numerical approach to computing the Schr\"odinger map (SM) based on the Hasimoto transform which relates the SM flow to a cubic nonlinear Schr\"odinger (NLS) equation. In exploiting this nonlinear transform we are able to…
We perform a classification of third order integrable systems of evolution equations with respect to higher symmetries. Applying it, we consider polynomial systems that are 0-homogeneous under a suitable weighting of variables with main…
Seismic imaging is a major challenge in geophysics with broad applications. It involves solving wave propagation equations with absorbing boundary conditions (ABC) multiple times. This drives the need for accurate and efficient numerical…
Structure-preserving linearly implicit exponential integrators are constructed for Hamiltonian partial differential equations with linear constant damping. Linearly implicit integrators are derived by polarizing the polynomial terms of the…
In this paper, we propose an efficient exponential integrator finite element method for solving a class of semilinear parabolic equations in rectangular domains. The proposed method first performs the spatial discretization of the model…
We propose new local error estimators for splitting and composition methods. They are based on the construction of lower order schemes obtained at each step as a linear combination of the intermediate stages of the integrator, so that the…
This paper introduces a high-order accurate surface integral equation method for solving 3D electromagnetic scattering for dielectric objects with uniaxially anisotropic permittivity tensors. The N-M\"uller formulation is leveraged…
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…