Hugo Jiménez-Pérez

Hugo Jiménez-Pérez, Jacques Laskar

We introduce a time-parallel algorithm for solving numerically almost integrable Hamiltonian systems in action-angle coordinates. This algorithm is a refinement of that introduced by Saha, Stadel and Tremaine in 1997 (SST97) for the same type of problems. Our refined algorithm has a better convergence obtained from the use of derivatives of the perturbing term not considered in the original SST97 algorithm. An advantage of this algorithm is its independence of the step-size for the parallelized procedures which can be consider as a particular case of the parareal scheme.

Hugo Jiménez-Pérez

In this article we introduce a new method for constructing implicit symplectic maps using special symplectic manifolds and Liouvillian forms. This method extends, in a natural way, the method of generating functions to 1-forms which are globally defined on the symplectic manifold. The maps constructed by this method, are related to the symplectic Cayley's transformation and belong to a continuous space of dimension n(2n+1). Applying the implicit map to the discrete Hamilton equations we obtain the generalized symplectic Euler scheme. We show the relations of the elements of this family with other discrete symplectic mapping, in particular 1) with the mappings obtained by generating functions of type I, II, and III and IV; 2) with the symplectic Euler methods A and B; and 3) with the mid-point rule. Moreover, we show the corresponding symplectic diffeomorphisms and their Liouvillian forms on the product symplectic manifold. We illustrate the details of the method in constructing two different families of implicit symplectic maps for $n=1$. This is a geometrical method which overcomes the difficulties of the Hamilton-Jacobi theory and generating functions.

Hugo Jiménez-Pérez

In this note we give simple symplecticity conditions for implicit schemes in the linear case. We consider implicit maps on generic symplectic manifold and we introduce the concept of consistent implicit maps, to generalize the symplecticity conditions to symplectic manifolds. Additionally, we give a preliminary geometrical interpretation of those conditions.

Hugo Jiménez-Pérez

Symplectic integrators constructed from Hamiltonian and Lie formalisms are obtained as symplectic maps whose flow follows the exact solution of a "sourrounded" Hamiltonian K = H + h^k H_1. Those modified Hamiltonians depends virtually on the time by h. When the numerical integration of a Hamiltonian system involves more than one symplectic scheme as in the parallel-in-time algorithms, there are not a simple way to control the dynamical behavior of the error Hamiltonian. The interplay of to different symplectic integrators can degenerate their behavior if both have different dynamical properties, reflected in the number of iterations to approximate the sequential solution. Considered as flows of time-dependent Hamiltonians we use the Hofer's geometry to search for the optimal coupling of symplectic schemes. As a result we obtain the constraints in the Parareal method to have a good behavior for Hamiltonian dynamics.

Hugo Jiménez-Pérez, Jean-Pierre Vilotte, Barbara Romanowicz

We implement and investigate the numerical properties of a new family of integrators connecting both variants of the symplectic Euler schemes, and including an alternative to the classical symplectic mid-point scheme, with some additional terms. This family is derived from a new method, introduced in a previous study, for generating symplectic integrators based on the concept of special symplectic manifold. The use of symplectic rotations and a particular type of projection keeps the whole procedure within the symplectic framework. We show that it is possible to define a set of parameters that control the additional terms providing a way of "tuning" these new symplectic schemes. We test the "tuned" symplectic integrators with the perturbed pendulum and we compare its behavior with an explicit scheme for perturbed systems. Remarkably, for the given examples, the error in the energy integral can be reduced considerably. There is a natural geometrical explanation, sketched at the end of this paper. This is the subject of a parallel article where a finer analysis is performed. Numerical results obtained in this paper open a new point of view on symplectic integrators and Hamiltonian error.