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

Sergio Blanes

Sergio Blanes appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NAastro-ph.EPNumerical Analysisphysics.comp-phcond-mat.othermath-phmath.MPquant-ph

Catalog footprint

What is connected

15works
8topics
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

15 published item(s)

preprint2022arXiv

Positivity-preserving methods for population models

Many important applications are modelled by differential equations with positive solutions. However, it remains an outstanding open problem to develop numerical methods that are both (i) of a high order of accuracy and (ii) capable of preserving positivity. It is known that the two main families of numerical methods, Runge-Kutta methods and multistep methods, face an order barrier: if they preserve positivity, then they are constrained to low accuracy: they cannot be better than first order. We propose novel methods that overcome this barrier: our methods are of second order, and they are guaranteed to preserve positivity. Our methods apply to a large class of differential equations that have a special graph Laplacian structure, which we elucidate. The equations need be neither linear nor autonomous and the graph Laplacian need not be symmetric. This algebraic structure arises naturally in many important applications where positivity is required. We showcase our new methods on applications where standard high order methods fail to preserve positivity, including infectious diseases, Markov processes, master equations and chemical reactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On symmetric-conjugate composition methods in the numerical integration of differential equations

We analyze composition methods with complex coefficients exhibiting the so-called ``symmetry-conjugate'' pattern in their distribution. In particular, we study their behavior with respect to preservation of qualitative properties when projected on the real axis and we compare them with the usual left-right palindromic compositions. New schemes within this family up to order 8 are proposed and their efficiency is tested on several examples. Our analysis shows that higher-order schemes are more efficient even when time step sizes are relatively large.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

New families of symplectic splitting methods for numerical integration in dynamical astronomy

We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We derive in a systematic way an independent set of necessary and sufficient conditions to be satisfied by the coefficients of splitting methods to achieve a prescribed order of accuracy. Splitting methods satisfying such (generalized) order conditions are appropriate in particular for the numerical simulation of the Solar System described in Jacobi coordinates. We show that, when using Poincaré Heliocentric coordinates, the same order of accuracy may be obtained by imposing an additional polynomial equation on the coefficients of the splitting method. We construct several splitting methods appropriate for each of the two sets of coordinates by solving the corresponding systems of polynomial equations and finding the optimal solutions. The experiments reported here indicate that the efficiency of our new schemes is clearly superior to previous integrators when high accuracy is required.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Symplectic integrators for the matrix Hill's equation and its applications to engineering models

We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian function and the fundamental matrix solution is a symplectic matrix. This is a very important property to be preserved by the numerical integrators. In this work we present new sixth-and eighth-order symplectic exponential integrators that are tailored to the Hill's equation. The methods are based on an efficient symplectic approximation to the exponential of high dimensional coupled autonomous harmonic oscillators and yield accurate results for oscillatory problems at a low computational cost. Several numerical examples illustrate the performance of the new methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-order splitting methods for separable non-autonomous parabolic equations

We consider the numerical integration of non-autonomous separable parabolic equations using high order splitting methods with complex coefficients (methods with real coefficients of order greater than two necessarily have negative coefficients). We propose to consider a class of methods in which one set of the coefficients are real and positive numbers, and to split properly the system in the extended phase space where the time is taken as a new coordinate. This allows us to evaluate all time-dependent operators at real values of the time, leading to schemes which are stable and simple to implement. If the system can be considered as the perturbation of an exactly solvable problem and the flow of the dominant part is advanced using the real coefficients, it is possible to build highly efficient methods for these problems. We show the performance of this class of methods on several numerical examples and present some new improved schemes.

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
preprint2014arXiv

The Scaling, Splitting and Squaring Method for the Exponential of Perturbed Matrices

We propose splitting methods for the computation of the exponential of perturbed matrices which can be written as the sum $A=D+\varepsilon B$ of a sparse and efficiently exponentiable matrix $D$ with sparse exponential $e^D$ and a dense matrix $\varepsilon B$ which is of small norm in comparison with $D$. The predominant algorithm is based on scaling the large matrix $A$ by a small number $2^{-s}$, which is then exponentiated by efficient Padé or Taylor methods and finally squared in order to obtain an approximation for the full exponential. In this setting, the main portion of the computational cost arises from dense-matrix multiplications and we present a modified squaring which takes advantage of the smallness of the perturbed matrix $B$ in order to reduce the number of squarings necessary. Theoretical results on local error and error propagation for splitting methods are complemented with numerical experiments and show a clear improvement over existing methods when medium precision is sought.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

High order structure preserving explicit methods for solving linear-quadratic optimal control problems and differential games

We present high order explicit geometric integrators to solve linear-quadratic optimal control problems and $N$-player differential games. These problems are described by a system coupled non-linear differential equations with boundary conditions. We propose first to integrate backward in time the non-autonomous matrix Riccati differential equations and next to integrate forward in time the coupled system of equations for the Riccati and the state vector. This can be achieved by using appropriate splitting methods, which we show they preserve most qualitative properties of the exact solution. Since the coupled system of equations is usually explicitly time dependent, a preliminary analysis has to be considered. We consider the time as two new coordinates, and this allows us to integrate the whole system forward in time using splitting methods while preserving the most relevant qualitative structure of the exact solution. If the system is a perturbation of an exactly solvable problem, the performance of the splitting methods considerably improves. Some numerical examples are also considered which show the performance of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Solving the Schrödinger eigenvalue problem by the imaginary time propagation technique using splitting methods with complex coefficients

The Schrödinger eigenvalue problem is solved with the imaginary time propagation technique. The separability of the Hamiltonian makes the problem suitable for the application of splitting methods. High order fractional time steps of order greater than two necessarily have negative steps and can not be used for this class of diffusive problems. However, there exist methods which use fractional complex time steps with positive real parts which can be used with only a moderate increase in the computational cost. We analyze the performance of this class of schemes and propose new methods which outperform the existing ones in most cases. On the other hand, if the gradient of the potential is available, methods up to fourth order with real and positive coefficients exist. We also explore this case and propose new methods as well as sixth-order methods with complex coefficients. In particular, highly optimized sixth-order schemes for near integrable systems using positive real part complex coefficients with and without modified potentials are presented. A time-stepping variable order algorithm is proposed and numerical results show the enhanced efficiency of the new methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High precision Symplectic Integrators for the Solar System

Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These comparisons are made in Jacobi and Heliocentric coordinates and the implementation of the algorithms is fully detailed for practical use. We conclude that high order integrators should be privileged, with a preference for the new $(10,6,4)$ method of (Blanes et al., 2012)

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Structure preserving integrators for solving linear quadratic optimal control problems with applications to describe the flight of a quadrotor

We present structure preserving integrators for solving linear quadratic optimal control problems. This problem requires the numerical integration of matrix Riccati differential equations whose exact solution is a symmetric positive definite time-dependent matrix which controls the stability of the equation for the state. This property is not preserved, in general, by the numerical methods. We propose second order exponential methods based on the Magnus series expansion which unconditionally preserve positivity for this problem and analyze higher order Magnus integrators. This method can also be used for the integration of nonlinear problems if they are previously linearized. The performance of the algorithms is illustrated with the stabilization of a quadrotor which is an unmanned aerial vehicle.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Error analysis of splitting methods for the time dependent Schrodinger equation

A typical procedure to integrate numerically the time dependent Schrö\-din\-ger equation involves two stages. In the first one carries out a space discretization of the continuous problem. This results in the linear system of differential equations $i du/dt = H u$, where $H$ is a real symmetric matrix, whose solution with initial value $u(0) = u_0 \in \mathbb{C}^N$ is given by $u(t) = \e^{-i t H} u_0$. Usually, this exponential matrix is expensive to evaluate, so that time stepping methods to construct approximations to $u$ from time $t_n$ to $t_{n+1}$ are considered in the second phase of the procedure. Among them, schemes involving multiplications of the matrix $H$ with vectors, such as Lanczos and Chebyshev methods, are particularly efficient. In this work we consider a particular class of splitting methods which also involves only products $Hu$. We carry out an error analysis of these integrators and propose a strategy which allows us to construct different splitting symplectic methods of different order (even of order zero) possessing a large stability interval that can be adapted to different space regularity conditions and different accuracy ranges of the spatial discretization. The validity of the procedure and the performance of the resulting schemes are illustrated on several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fourier methods for the perturbed harmonic oscillator in linear and nonlinear Schrödinger equations

We consider the numerical integration of the Gross-Pitaevskii equation with a potential trap given by a time-dependent harmonic potential or a small perturbation thereof. Splitting methods are frequently used with Fourier techniques since the system can be split into the kinetic and remaining part, and each part can be solved efficiently using Fast Fourier Transforms. To split the system into the quantum harmonic oscillator problem and the remaining part allows to get higher accuracies in many cases, but it requires to change between Hermite basis functions and the coordinate space, and this is not efficient for time-dependent frequencies or strong nonlinearities. We show how to build new methods which combine the advantages of using Fourier methods while solving the timedependent harmonic oscillator exactly (or with a high accuracy by using a Magnus integrator and an appropriate decomposition).

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimized high-order splitting methods for some classes of parabolic equations

We are concerned with the numerical solution obtained by splitting methods of certain parabolic partial differential equations. Splitting schemes of order higher than two with real coefficients necessarily involve negative coefficients. It has been demonstrated that this second-order barrier can be overcome by using splitting methods with complex-valued coefficients (with positive real parts). In this way, methods of orders 3 to 14 by using the Suzuki--Yoshida triple (and quadruple) jump composition procedure have been explicitly built. Here we reconsider this technique and show that it is inherently bounded to order 14 and clearly sub-optimal with respect to error constants. As an alternative, we solve directly the algebraic equations arising from the order conditions and construct methods of orders 6 and 8 that are the most accurate ones available at present time, even when low accuracies are desired. We also show that, in the general case, 14 is not an order barrier for splitting methods with complex coefficients with positive real part by building explicitly a method of order 16 as a composition of methods of order 8.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Splitting methods with complex coefficients

Splitting methods for the numerical integration of differential equations of order greater than two involve necessarily negative coefficients. This order barrier can be overcome by considering complex coefficients with positive real part. In this work we review the composition technique used to construct methods of this class, propose new sixth-order integrators and analyze their main features on a pair of numerical examples, in particular how the errors are propagated along the evolution.

0 citations0 reviewsNo code linkedOpen access