A study of rigorous ODE integrators for multi-scale set-oriented computations

Abstract We study the usefulness of two most prominent publicly available rigorous ODE integrators: one provided by the CAPD group ( capd.ii.uj.edu.pl ), the other based on the COSY Infinity project ( cosyinfinity.org ). Both integrators are capable of handling entire sets of initial conditions and provide tight rigorous outer enclosures of the images under a time-T map. We conduct extensive benchmark computations using the well-known Lorenz system, and compare the computation time against the final accuracy achieved. We also discuss the effect of a few technical parameters, such as the order of the numerical integration method, the value of T, and the phase space resolution. We conclude that COSY may provide more precise results due to its ability of avoiding the variable dependency problem. However, the overall cost of computations conducted using CAPD is typically lower, especially when intervals of parameters are involved. Moreover, access to COSY is limited (registration required) and the rigorous ODE integrators are not publicly available, while CAPD is an open source free software project. Therefore, we recommend the latter integrator for this kind of computations. Nevertheless, proper choice of the various integration parameters turns out to be of even greater importance than the choice of the integrator itself.

[1]  Konstantin Mischaikow,et al.  Discrete Morse Theoretic Algorithms for Computing Homology of Complexes and Maps , 2014, Found. Comput. Math..

[2]  M. Berz,et al.  TAYLOR MODELS AND OTHER VALIDATED FUNCTIONAL INCLUSION METHODS , 2003 .

[3]  Martin Berz,et al.  Verified Integration of ODEs and Flows Using Differential Algebraic Methods on High-Order Taylor Models , 1998, Reliab. Comput..

[4]  J. Penrod.  Michigan State University : East Lansing-Michigan , 2007 .

[5]  Piotr Zgliczynski,et al.  C1 Lohner Algorithm , 2002, Found. Comput. Math..

[6]  Kyoko Makino,et al.  Rigorous analysis of nonlinear motion in particle accelerators , 1998 .

[7]  Martin Berz,et al.  SUPPRESSION OF THE WRAPPING EFFECT BY TAYLOR MODEL-BASED VERIFIED INTEGRATORS : LONG-TERM STABILIZATION BY SHRINK WRAPPING , 2006 .

[8]  B. Carreras,et al.  A dynamical model for plasma confinement transitions , 2012 .

[9]  P. Pilarczyk,et al.  Global dynamics in a stage-structured discrete-time population model with harvesting. , 2012, Journal of theoretical biology.

[10]  E. Lorenz.  Deterministic nonperiodic flow , 1963 .

[11]  Pawel Pilarczyk,et al.  Rich Bifurcation Structure in a Two-Patch Vaccination Model , 2015, SIAM J. Appl. Dyn. Syst..

[12]  Daniel Wilczak,et al.  An implicit algorithm for validated enclosures of the solutions to variational equations for ODEs , 2015, Appl. Math. Comput..

[13]  Konstantin Mischaikow,et al.  A Database Schema for the Analysis of Global Dynamics of Multiparameter Systems , 2009, SIAM J. Appl. Dyn. Syst..

[14]  Nathalie Revol,et al.  Taylor models and floating-point arithmetic: proof that arithmetic operations are validated in COSY , 2005, J. Log. Algebraic Methods Program..

[15]  M. Gameiro,et al.  Combinatorial-topological framework for the analysis of global dynamics. , 2012, Chaos.

[16]  Pawel Pilarczyk,et al.  Parallelization Method for a Continuous Property , 2010, Found. Comput. Math..

[17]  Konstantin Mischaikow,et al.  Coarse Dynamics for Coarse Modeling: An Example From Population Biology , 2014, Entropy.

[18]  D. Wilczak,et al.  Cr-Lohner algorithm , 2011 .

[19]  Martin Berz,et al.  SUPPRESSION OF THE WRAPPING EFFECT BY TAYLOR MODEL- BASED VERIFIED INTEGRATORS: LONG-TERM STABILIZATION BY PRECONDITIONING , 2011 .

[20]  M. Mrozek.  The Conley Index on Compact Anr’s Is of Finite Type , 1990 .

[21]  John D. Pryce,et al.  An Effective High-Order Interval Method for Validating Existence and Uniqueness of the Solution of an IVP for an ODE , 2001, Reliab. Comput..

[22]  Tomasz Kapela,et al.  A Lohner-type algorithm for control systems and ordinary differential inclusions , 2007, 0712.0910.

[23]  Konstantin Mischaikow,et al.  Inducing a map on homology from a correspondence , 2014 .

[24]  P. Pilarczyk,et al.  Excision-preserving cubical approach to the algorithmic computation of the discrete Conley index , 2008 .

[25]  Konstantin Mischaikow,et al.  Graph Approach to the Computation of the Homology of Continuous Maps , 2005, Found. Comput. Math..

[26]  D. Wilczak,et al.  C-Lohner algorithm , 2011 .

[27]  C. Conley.  Isolated Invariant Sets and the Morse Index , 1978 .

[28]  A. Szymczak.  The Conley index for discrete semidynamical systems , 1995 .

[29]  Marian Mrozek,et al.  Set arithmetic and the enclosing problem in dynamics , 2000 .