PublicationsMy Feed Description2017-04-29T00:00:00+01:00Cédric M. Camposhttps://cmcampos.xyz/publicationsPalindromic 3-stage splitting integrators, a roadmaphttps://cmcampos.xyz/publications/s3map2017-04-29T00:00:00+01:002017-04-29T00:00:00+01:00
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 of palindromic, three-stage splitting formulas and identify choices of parameters that may outperform the Strang/Verlet method. One of these choices leads to a method of effective order four suitable to integrate in time some partial differential equations. Other choices may be seen as perturbations of the Strang method that increase efficiency in molecular dynamics simulations and in Hybrid Monte Carlo sampling.
]]>
High order variational integrators in the optimal control of mechanical systemshttps://cmcampos.xyz/publications/hovi-ocms2017-02-01T00:00:00+00:002017-02-01T00:00:00+00:00
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 two high order variational integrators which distinguish themselves in the dimension of the underling space of approximation and we investigate their application to finite-dimensional optimal control problems posed with mechanical systems. The convergence of state and control variables of the approximated problem is shown. Furthermore, by analyzing the adjoint systems of the optimal control problem and its discretized counterpart, we prove that, for these particular integrators, dualization and discretization commute.
]]>
Hamilton-Jacobi theory in Cauchy data spacehttps://cmcampos.xyz/publications/cft-hj2014-11-14T00:00:00+00:002014-11-14T00:00:00+00:00
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.
]]>
Extra Chance Generalized Hybrid Monte Carlohttps://cmcampos.xyz/publications/xhmc2014-07-30T00:00:00+01:002014-07-30T00:00:00+01:00
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 experiments that clearly indicate that the additional work per sample carried out in the extra chance approach clearly pays in terms of the quality of the samples generated.
]]>
Higher Order Variational Integrators: a polynomial approachhttps://cmcampos.xyz/publications/hovi-poly2013-07-23T00:00:00+01:002013-07-23T00:00:00+01:00
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, the evolution of the momentum maps or the energy behaviour. Also they are easily applicable to optimal control problems based on mechanical systems as proposed in Ober-Blöbaum et al. [2011]. Following the same approach, we develop a family of variational integrators to which we refer as symplectic Galerkin schemes in contrast to symplectic partitioned Runge-Kutta. These two families of integrators are, in principle and by construction, different one from the other. Furthermore, the symplectic Galerkin family can as easily be applied in optimal control problems, for which Campos et al. [2012b] is a particular case.
]]>
Higher order variational time discretization of optimal control problemshttps://cmcampos.xyz/publications/hovtd-ocp2012-04-27T00:00:00+01:002012-04-27T00:00:00+01:00
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 the evolution of the momentum maps or the energy behavior. Here, we employ higher order quadrature rules based on polynomial collocation. The resulting variational time discretization decreases the overall computational effort.
]]>
Classical field theories of first order and lagrangian submanifolds of premultisymplectic manifoldshttps://cmcampos.xyz/publications/cft-t32011-10-21T00:00:00+01:002011-10-21T00:00:00+01:00
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 prove that Euler-Lagrange and Hamilton-De Donder-Weyl equations are the local equations defining Lagrangian submanifolds of a premultisymplectic manifold.
]]>
Constrained Variational Calculus for Higher Order Classical Field Theorieshttps://cmcampos.xyz/publications/cft-ho-cons2010-05-12T00:00:00+01:002010-05-12T00:00:00+01:00
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 control theory for partial differential equations.
]]>
Unambiguous Formalism for Higher-Order Lagrangian Field Theorieshttps://cmcampos.xyz/publications/cft-ho-ulf2009-11-16T00:00:00+00:002009-11-16T00:00:00+00:00
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 propose a differential-geometric setting for the dynamics of a higher-order field theory, based on the Skinner and Rusk formalism for mechanics. This approach incorporates aspects of both, the Lagrangian and the Hamiltonian description, since the field equations are formulated using the Lagrangian on a higher-order jet bundle and the canonical multisymplectic form on its dual. As both of these objects are uniquely defined, the Skinner-Rusk approach has the advantage that it does not suffer from the arbitrariness in conventional descriptions. The result is that we obtain a unique and global intrinsic version of the Euler-Lagrange equations for higher-order field theories. Several examples illustrate our construction.
]]>
Functionally Graded Mediahttps://cmcampos.xyz/publications/fgm2007-11-16T00:00:00+00:002007-11-16T00:00:00+00:00
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.
]]>