Publications

Polyak's Heavy Ball (PHB; Polyak, 1964), a.k.a. Classical Momentum, and Nesterov's Accelerated Gradient (NAG; Nesterov, 1983) are well known examples of momentum-descent methods for optimization. While the latter outperforms the former, solely generalizations of PHB-like methods to nonlinear spaces...

Many of the new developments in machine learning are connected with gradient-based optimization methods. Recently, these methods have been studied using a variational perspective (Betancourt et al., 2018). This has opened up the possibility of introducing variational and symplectic methods using geo...

The implementation of multi-stage splitting integrators is essentially the same as the implementation of the familiar Strang/Verlet method. Therefore multi-stage formulas may be easily incorporated into software that now uses the Strang/Verlet integrator. We study in detail the two-parameter family...

In recent years, much effort in designing numerical methods for the simulation and optimization of mechanical systems has been put into schemes which are structure preserving. One particular class are variational integrators which are momentum preserving and symplectic. In this article, we develop t...

Recently, M. de León el al. ([9]) have developed a geometric Hamilton-Jacobi theory for Classical Field Theories in the setting of multisymplectic geometry. Our purpose in the current paper is to establish the corresponding Hamilton-Jacobi theory in the Cauchy data space, and relate both approaches....

We study a method, Extra Chance Generalized Hybrid Monte Carlo, to avoid rejections in the Hybrid Monte Carlo method and related algorithms. In the spirit of delayed rejection, whenever a rejection would occur, extra work is done to find a fresh proposal that, hopefully, may be accepted. We present...

We reconsider the variational derivation of symplectic partitioned Runge-Kutta schemes. Such type of variational integrators are of great importance since they integrate mechanical systems with high order accuracy while preserving the structural properties of these systems, like the symplectic form,...

We reconsider the variational integration of optimal control problems for mechanical systems based on a direct discretization of the Lagrange-d'Alembert principle. This approach yields discrete dynamical constraints which by construction preserve important structural properties of the system, like t...

A description of classical field theories of first order in terms of Lagrangian submanifolds of premultisymplectic manifolds is presented. For this purpose, a Tulczyjew's triple associated with a fibration is discussed. The triple is adapted to the extended Hamiltonian formalism. Using this triple,...

We develop an intrinsic geometrical setting for higher order constrained field theories. As a main tool we use an appropriate generalization of the classical Skinner-Rusk formalism. Some examples of application are studied, in particular, applications to the geometrical description of optimal contro...

The aim of this paper is to propose an unambiguous intrinsic formalism for higher-order field theories which avoids the arbitrariness in the generalization of the conventional description of field theories, which implies the existence of different Cartan forms and Legendre transformations. We propos...

The notions of uniformity and homogeneity of elastic materials are reviewed in terms of Lie groupoids and frame bundles. This framework is also extended to consider the case Functionally Graded Media, which allows us to obtain some homogeneity conditions.