International Workshop on Validated Numerics for Computer-Assisted Proofs
6-10 July 2026, ICMS, Bayes Centre, Edinburgh
Note: The programme is provisional and subject to changes.
Join us at the ICMS, top floor of the Bayes Centre. Refreshments will be provided.
Jean-Pierre Eckmann - Université de Genève, Switzerland
Although “everybody” might know how this subject evolved, I will try to reconstruct how the early issues and difficulties came up, and how Koch and Lanford (and many others) worked to make this subject a common modern tool in analysis.
Jan Bouwe van den Berg - Vrije Universiteit Amsterdam, Netherlands
An expository talk about CAPs in nonlinear dynamics based on spectral methods.
Hugo Chu - École Polytechnique, France
Rigorously enclosing Lyapunov exponents of stochastic flows is a long-standing and notoriously difficult problem: even determining the sign of the top exponent is often out of reach outside special structures or perturbative regimes. In this paper, we develop a computer-assisted method to obtain rigorous upper and lower bounds on Lyapunov exponents of stochastic flows under mild hypoellipticity assumptions. Our approach starts from the Furstenberg-Khasminskii representation of the top Lyapunov exponent as an ergodic average over the invariant measure of the associated projective process, but crucially avoids any rigorous computation of that invariant measure. Instead, we introduce a numerical adjoint method, amenable to rigorous numerics, and show that an approximate solution of the associated Poisson equation yields a certified enclosure of the Lyapunov exponent via an explicit bound on the residual. This converts a non-rigorous numerical approximation into a rigorous quantitative estimate. The method applies to systems on both compact and non-compact state spaces, does not require special geometric structure, and is not restricted to small-noise or other perturbative settings. We use it to prove positivity of the top Lyapunov exponent for several stochastic systems, including examples exhibiting noise-induced chaos and parameter-dependent sign changes, and we combine it with continuation methods to obtain rigorous bounds over large parameter regions.
Polina Vytnova - University of Surrey, England
(Pending)
Caroline Wormell - The University of Sydney, Australia
(Pending)
Siegfried Rump - Hamburg University of Technology, Germany
The solution of sparse systems is ubiquitous in numerical computations. Despite several efforts, there were no satisfactory method for computing verified bounds for the solution and was one of the “Grand Challenges” in verification methods for decades. In this talk we present two algorithms to compute entrywise error bounds for the solution of general real or complex sparse systems with condition number up to the limit 1016 7 . Our algorithms split into three subalgorithms for symmetric positive definite, symmetric indefinite and general input matrix $A$. A key point is a factorization of $A$ into $L_1L_2$ such that $L_1$ and $L_2$ have identical sets of singular values with the smallest one close to $\sigma_{\mathrm{min}}(A)^{1/2}$. A mathematically correct lower bound on $\sigma_{\mathrm{min}}(L_1) = \sigma_{\mathrm{min}}(L_2)$ is then computed using $L_1^T L_1$. Based on that a second method exploring the inertia of a symmetric/Hermitian matrix is presented. It is often slower but more robust, i.e., it may produce verified inclusions where the first method fails. We show how to compute inclusions with almost maximal accuracy for all entries, i.e., all bounds differ by few bits. That is based on a fast method to compute accurate approximations and bounds for extremely ill-conditioned dot products with a very efficient Matlab implementation. Both approaches are used to compute verified error bounds for the solution of least squares problems and for underdetermined linear systems. Inclusions of the solution of general real or complex systems of nonlinear equations with sparse Jacobi matrix are computed by transforming the problem into a linear system with point matrix and interval right hand side. All algorithms are implemented in pure Matlab/Octave and included from Version 14 of INTLAB, the Matlab/Octave toolbox for verified computations.
Nicolae Mihalache - Univ Paris-Est Creteil, France
(Pending)
Maxime Breden - École Polytechnique, France
(Pending)
George Browne - The Open University, England
(Pending)
Zainab Rahman - University of Portsmouth, England
(Pending)
Benjamin Mestel - Open University (Retired), England
(Pending - On behalf of Farhana Pramy)
Jason Mireles-James - Florida Atlantic University, USA
Some interesting applications to celestial mechanics and geodesic flow, helped by some advances made over the last few years on validated numerics for the DFT (the same technology that helped us do our Feigenbaum fixed point paper).
Maciej Capinski - AGH University Kraków, Poland
We present a computer-assisted methodology based on the propagation of topological discs and discuss its applications to Arnold diffusion and blenders in celestial mechanics.
Matthieu Cadiot - École Polytechnique, France
The Whitham equation is a nonlocal model used to study shallow water waves, capturing complex phenomena such as wave breaking, peaking, and various traveling wave solutions. In this talk, I will present a computer-assisted methodology for studying solitary waves—traveling waves that vanish at infinity. In particular, the method provides a framework for the constructive existence proof of a true solution in the close vicinity of an approximate, numerically computed solution. Furthermore, the tight control obtained over the solution enables us to enclose the spectrum of the linearization at the solitary wave, leading to a rigorous and conclusive stability analysis.
Juan Miranda - Florida Atlantic University, USA
(Pending)
Julia Slipantschuk - University of Bayreuth, Germany
(Pending)
Wednesday has been planned as a half-day at the ICMS to provide participants with opportunities to network and explore the historic city of Edinburgh together in the afternoon. A number of activities are suggested and some will be led by the organisers.
David Kerridges - James Clerk Maxwell Foundation, Scotland
Daniel Wilczak - Jagiellonian University in Kraków, Poland
The CAPD library [1,2] is a C++ library devoted to validated numerics for dynamical systems. It has been developed for over 30 years by a research team at the Jagiellonian University in Kraków, Poland. The library is mostly oriented towards continuous-time systems and implements rigorous integrators for ODEs, higher order variational equations [3,4], differential inclusions [5], and a class of dissipative PDEs [6]. In the talk I’m going to present the basic interface of the library, then give a few one-screen short yet not very trivial case studies. I will present an application of the library to prove the existence of a uniformly hyperbolic chaotic attractor in a 4D non-autonomous ODE [7]. References: [1] https://github.com/CAPDgroup/CAPD [2] T. Kapela, M. Mrozek, D. Wilczak, P. Zgliczyński. CAPD::DynSys: a flexible C++ toolbox for rigorous numerical analysis of dynamical systems. Commun Nonlinear Sci Numer Simul 2021;101:105578. [3] P. Zgliczyński, $C^1$-Lohner algorithm, Found. Comp. Math, 2(4), pp.429-465 [4] D. Wilczak, P. Zgliczyński, $C^r$-Lohner algorithm, Schedae Informaticae 20, 9-46, 2011. [5] T. Kapela, P. Zgliczyński, A Lohner-type algorithm for control systems and ordinary differential inclusions, Disc. Cont. Dyn. Sys-B, 11(2), pp.365-385 [6] D. Wilczak, P. Zgliczyński. A geometric method for infinite-dimensional chaos: Symbolic dynamics for the Kuramoto-Sivashinsky PDE on the line. J Differential Equations 2020;269(10):8509–48 [7] D. Wilczak, Uniformly hyperbolic attractor of the Smale-Williams type for a Poincaré map in the Kuznetsov system, SIADS. 9 (2010), 1263-1283.
Tomasz Kapela - Jagiellonian University in Kraków, Poland
(Pending)
Piotr Zgliczynski - Jagiellonian University in Kraków, Poland
We will describe the proof of the existence of symbolic dynamics for the Kuramoto-Sivashinsky PDE on the line: a semiconjugacy with the Bernoulli shift, the existence of countably many periodic orbits and homo- and heteroclinic connections between them. On this example of computer assisted proof we would like to explain how to perform rigorous numerics for dissipative PDEs with periodic boundary conditions and how to use finite dimensional tools from dynamics in the context of dissipative PDEs.
Jean-Philippe Lessard - McGill University, Canada
(Pending)
Luis Benet Fernández - Universidad Nacional Autónoma de México (UNAM), Mexico
(Pending)
Igors Gorbovickis - Constructor University, Germany
(Pending)
Artem Dudko - Institute of Mathematics of the Polish Academy of Sciences (IM PAN), Poland
In 1990s Bruin, Keller, Nowicki, and van Strien showed that smooth unimodal maps with Fibonacci combinatorics and sufficiently high degree of a critical point have a wild attractor, i.e. their metric and topological attractors do not coincide. However, until recently there were no reasonable estimates on the degree of the critical point needed. In the talk I will present a constructive approach for studying attractors of maps, which are periodic points of a renormalization. Using this approach and rigorous computer estimates, we show that the Fibonacci map of degree d=3.8 does not have a wild attractor, but that for degree d=5.1 the wild attractor exists. The talk is based on a joint work with Denis Gaidashev.
Jeremy Parker - University of Dundee, Scotland
(Pending)
Olivier Henot - National Taiwan University, Taiwan
(Pending)