BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

J. M. Sanz-Serna

J. M. Sanz-Serna appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysismath.DSmath.OCmath.PRphysics.comp-phComputationMachine Learningmath.CO

Catalog footprint

What is connected

10works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2021arXiv

The connections between Lyapunov functions for some optimization algorithms and differential equations

In this manuscript, we study the properties of a family of second-order differential equations with damping, its discretizations and their connections with accelerated optimization algorithms for $m$-strongly convex and $L$-smooth functions. In particular, using the Linear Matrix Inequality LMI framework developed by \emph{Fazlyab et. al. $(2018)$}, we derive analytically a (discrete) Lyapunov function for a two-parameter family of Nesterov optimization methods, which allows for the complete characterization of their convergence rate. In the appropriate limit, this family of methods may be seen as a discretization of a family of second-order ordinary differential equations for which we construct(continuous) Lyapunov functions by means of the LMI framework. The continuous Lyapunov functions may alternatively, be obtained by studying the limiting behaviour of their discrete counterparts. Finally, we show that the majority of typical discretizations of the family of ODEs, such as the Heavy ball method, do not possess Lyapunov functions with properties similar to those of the Lyapunov function constructed here for the Nesterov method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A technique for studying strong and weak local errors of splitting stochastic integrators

We present a technique, based on so-called word series, to write down in a systematic way expansions of the strong and weak local errors of splitting algorithms for the integration of Stratonovich stochastic differential equations. Those expansions immediately lead to the corresponding order conditions. Word series are similar to, but simpler than, the B-series used to analyze Runge-Kutta and other one-step integrators. The suggested approach makes it unnecessary to use the Baker-Campbell-Hausdorff formula. As an application, we compare two splitting algorithms recently considered by Leimkuhler and Matthews to integrate the Langevin equations. The word series method bears out clearly reasons for the advantages of one algorithm over the other.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Adaptive multi-stage integrators for optimal energy conservation in molecular simulations

We introduce a new Adaptive Integration Approach (AIA) to be used in a wide range of molecular simulations. Given a simulation problem and a step size, the method automatically chooses the optimal scheme out of an available family of numerical integrators. Although we focus on two-stage splitting integrators, the idea may be used with more general families. In each instance, the system-specific integrating scheme identified by our approach is optimal in the sense that it provides the best conservation of energy for harmonic forces. The AIA method has been implemented in the BCAM-modified GROMACS software package. Numerical tests in molecular dynamics and hybrid Monte Carlo simulations of constrained and unconstrained physical systems show that the method successfully realises the fail-safe strategy. In all experiments, and for each of the criteria employed, the AIA is at least as good as, and often significantly outperforms the standard Verlet scheme, as well as fixed parameter, optimized two-stage integrators. In particular, the sampling efficiency found in simulations using the AIA is up to 5 times better than the one achieved with other tested schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Vibrational resonance: a study with high-order word-series averaging

We study a model problem describing vibrational resonance by means of a high-order averaging technique based on so-called word series. With the tech- nique applied here, the tasks of constructing the averaged system and the associ- ated change of variables are divided into two parts. It is first necessary to build recursively a set of so-called word basis functions and, after that, all the required manipulations involve only scalar coefficients that are computed by means of sim- ple recursions. As distinct from the situation with other approaches, with word- series, high-order averaged systems may be derived without having to compute the associated change of variables. In the system considered here, the construction of high-order averaged systems makes it possible to obtain very precise approxima- tions to the true dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computing normal forms and formal invariants of dynamical systems by means of word series

We show how to use extended word series in the reduction of continuous and discrete dynamical systems to normal form and in the computation of formal invariants of motion in Hamiltonian systems. The manipulations required involve complex numbers rather than vector fields or diffeomorphisms. More precisely we construct a group G and a Lie algebra g in such a way that the elements of G and g are families of complex numbers; the operations to be performed involve the multiplication F in G and the bracket of g and result in universal coefficients that are then applied to write the normal form or the invariants of motion of the specific problem under consideration.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symplectic Runge-Kutta schemes for adjoint equations, automatic differentiation, optimal control and more

It is well known that symplectic Runge-Kutta and Partitioned Runge-Kutta methods exactly preserve {\em quadratic} first integrals (invariants of motion) of the system being integrated. While this property is often seen as a mere curiosity (it does not hold for arbitrary first integrals), it plays an important role in the computation of numerical sensitivities, optimal control theory and Lagrangian mechanics, as described in this paper, which, together with some new material, presents in a unified way a number of results now scattered or implicit in the literature. Some widely used procedures, such as the direct method in optimal control theory and the computation of sensitivities via reverse accumulation imply "hidden" integrations with symplectic Partitioned Runge-Kutta schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Word series for dynamical systems and their numerical integrators

We study word series and extended word series, classes of formal series for the analysis of some dynamical systems and their discretizations. These series are similar to but more compact than B-series. They may be composed among themselves by means of a simple rule. While word series have appeared before in the literature, extended word series are introduced in this paper. We exemplify the use of extended word series by studying the reduction to normal form and averaging of some perturbed integrable problems. We also provide a detailed analysis of the behaviour of splitting numerical methods for those problems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Extra Chance Generalized Hybrid Monte Carlo

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical integrators for the Hybrid Monte Carlo method

We construct numerical integrators for Hamiltonian problems that may advantageously replace the standard Verlet time-stepper within Hybrid Monte Carlo and related simulations. Past attempts have often aimed at boosting the order of accuracy of the integrator and/or reducing the size of its error constants; order and error constant are relevant concepts in the limit of vanishing step-length. We propose an alternative methodology based on the performance of the integrator when sampling from Gaussian distributions with not necessarily small step-lengths. We construct new splitting formulae that require two, three or four force evaluations per time-step. Limited, proof-of-concept numerical experiments suggest that the new integrators may provide an improvement on the efficiency of the standard Verlet method, especially in problems with high dimensionality.

0 citations0 reviewsNo code linkedOpen access