Related papers: A new approach to integrals of discretizations by …
We introduce a method to estimate the size of the domain of definition of the solutions of a meromorphic vector field on a neighborhood of its pole divisor. The corresponding techniques are, in a certain sense, quantitative versions of some…
Many conservative partial differential equations such as the Korteweg-de Vries (KdV) equation, and the nonlinear Schr\"{o}dinger equations, the Klein-Gordon equation have more than one invariant functionals. In this paper, we propose the…
Orthogonal polynomials with respect to the hypergeometric distribution on lattices in polyhedral domains in ${\mathbb R}^d$, which include hexagons in ${\mathbb R}^2$ and truncated tetrahedrons in ${\mathbb R}^3$, are defined and studied.…
A new method is proposed for integrating the equations of motion of an elastic filament. In the standard finite-difference and finite-element formulations the continuum equations of motion are discretized in space and time, but it is then…
We show that the exact integrator for the classical Kepler motion, recently found by Kozlov ({\it J. Phys. A: Math. Theor.\} {\bf 40} (2007) 4529-4539), can be derived in a simple natural way (using well known exact discretization of the…
Compatible discretizations, such as finite element exterior calculus, provide a discretization framework that respect the cohomological structure of the de Rham complex, which can be used to systematically construct stable mixed finite…
This paper studies the expressive and computational power of discrete Ordinary Differential Equations (ODEs), a.k.a. (Ordinary) Difference Equations. It presents a new framework using these equations as a central tool for computation and…
We develop a variational integrator for the shallow-water equations on a rotating sphere. The variational integrator is built around a discretization of the continuous Euler-Poincar\'{e} reduction framework for Eulerian hydrodynamics. We…
Viewing optimization methods as numerical integrators for ordinary differential equations (ODEs) provides a thought-provoking modern framework for studying accelerated first-order optimizers. In this literature, acceleration is often…
A general procedure for constructing conservative numerical integrators for time dependent partial differential equations is presented. In particular, linearly implicit methods preserving a time discretised version of the invariant is…
Using a scheme proposed earlier we set up Hamiltonian path integral quantization for a particle in two dimensions in plane polar coordinates.This scheme uses the classical Hamiltonian, without any $O(\hbar^2)$ terms, in the polar…
In this paper we apply Kahan's nonstandard discretization to three dimensional Lotka-Volterra equations in bi-Hamiltonian form. The periodicity of the solutions and all polynomial and non-polynomial invariants are well preserved in…
Discretization methods for differential-algebraic equations (DAEs) are considered that are based on the integration of an associated inherent ordinary differential equation (ODE). This allows to make use of any discretization scheme…
An interesting family of geometric integrators for Lagrangian systems can be defined using discretizations of the Hamilton's principle of critical action. This family of geometric integrators is called variational integrators. In this…
We present a novel numerical method for solving ODEs while preserving polynomial first integrals. The method is based on introducing multiple quadratic auxiliary variables to reformulate the ODE as an equivalent but higher-dimensional ODE…
Numerical methods that preserve geometric invariants of the system, such as energy, momentum or the symplectic form, are called geometric integrators. Variational integrators are an important class of geometric integrators. The general idea…
Numerical methods that preserves geometric invariants of the system such as energy, momentum and symplectic form, are called geometric integrators. These include variational integrators as an important subclass of geometric integrators. The…
The motion of a rigid body in a quadratic potential is an important example of an integrable Hamiltonian system on a dual to a semidirect product Lie algebra so(n) x Symm(n). We give a Lagrangian derivation of the corresponding equations of…
A variational integrator of arbitrarily high-order on the special orthogonal group $SO(n)$ is constructed using the polar decomposition and the constrained Galerkin method. It has the advantage of avoiding the second-order derivative of the…
We present a remarkable discretization of the classical Kepler problem which preserves its trajectories and all integrals of motion. The points of any discrete orbit belong to an appropriate continuous trajectory.