Initial-value ordinary differential equation solution via variable order Adams method (SIVA/DIVA) package is collection of subroutines for solution of nonstiff ordinary differential equations. There are versions for single-precision and double-precision arithmetic. Requires fewer evaluations of derivatives than other variable-order Adams predictor/ corrector methods. Option for direct integration of second-order equations makes integration of trajectory problems significantly more efficient. Written in FORTRAN 77.

Solving Ordinary Differential Equations

A new class of boundary conditions is described for quantum systems integrable by means of the quantum inverse scattering (R-matrix) method. The method proposed allows the author to treat open quantum chains with appropriate boundary terms in the Hamiltonian. The general considerations are applied to the XXZ and XYZ models, the nonlinear Schrodinger equation and Toda chain.

https://math.berkeley.edu/~reshetik/Sklyanin-boundary.pdf

Boundary conditions for integrable quantum systems

This paper gives a review of integration algorithms for finite dimensional mechanical systems that are based on discrete variational principles. The variational technique gives a unified treatment of many symplectic schemes, including those of higher order, as well as a natural treatment of the discrete Noether theorem. The approach also allows us to include forces, dissipation and constraints in a natural way. Amongst the many specific schemes treated as examples, the Verlet, SHAKE, RATTLE, Newmark, and the symplectic partitioned Runge–Kutta schemes are presented.

/pdf/discrete-mechanics-and-variational-integrators-4qm2t5dk6u.pdf

Discrete mechanics and variational integrators

https://digital.csic.es/bitstream/10261/15682/1/blanes.pdf

The Magnus expansion and some of its applications

Some new similarity reductions of the Boussinesq equation, which arises in several physical applications including shallow water waves and also is of considerable mathematical interest because it is a soliton equation solvable by inverse scattering, are presented. These new similarity reductions, including some new reductions to the first, second, and fourth Painleve equations, cannot be obtained using the standard Lie group method for finding group‐invariant solutions of partial differential equations; they are determined using a new and direct method that involves no group theoretical techniques.

New similarity reductions of the Boussinesq equation

Discrete gradient methods are well-known methods of Geometric Numerical Integration, which preserve the dissipation of gradient systems. In this paper we show that this property of discrete gradient methods can be interesting in the context of variational models for image processing, that is where the processed image is computed as a minimiser of an energy functional. Numerical schemes for computing minimisers of such energies are desired to inherit the dissipative property of the gradient system associated to the energy and consequently guarantee a monotonic decrease of the energy along iterations, avoiding situations in which more computational work might lead to less optimal solutions. Under appropriate smoothness assumptions on the energy functional we prove that discrete gradient methods guarantee a monotonic decrease of the energy towards stationary states, and we promote their use in image processing by exhibiting experiments with convex and non-convex variational models for image deblurring, denoising, and inpainting.

/pdf/discrete-gradient-methods-for-solving-variational-image-4rwjiolb7j.pdf

Discrete gradient methods for solving variational image regularisation models

Volume-preserving integrators have linear error growth

The derivative nonlinear Schrödinger equation and the massive Thirring model

The nonlinear Schrödinger equation and the anisotropic Heisenberg spin chain

The modified Hamiltonian is used to study the nonlinear stability of symplectic integrators, especially for nonlinear oscillators. We give conditions under which an initial condition on a compact energy surface will remain bounded for exponentially long times for sufficiently small time steps. � �

/pdf/on-the-nonlinear-stability-of-symplectic-integrators-4um19w0vtl.pdf

G. R. W. Quispel

Papers

Discrete gradient methods for solving variational image regularisation models

Volume-preserving integrators have linear error growth

The derivative nonlinear Schrödinger equation and the massive Thirring model

The nonlinear Schrödinger equation and the anisotropic Heisenberg spin chain

On the Nonlinear Stability of Symplectic Integrators