S3A2

Computation of smooth manifolds via rigorous multi-parameter continuation in infinite dimensions

Marcio Gameiro, Jean-Philippe Lessard and Alessandro Pugliese, 2013.

In this paper, we introduce a constructive rigorous numerical method to compute smooth manifolds implicitly defined by infinite dimensional nonlinear operators. We compute a simplicial triangulation of the manifold using a multi-parameter continuation method on a finite dimensional projection. The triangulation is then used to construct local charts and an atlas of the manifold in the infinite dimensional domain of the operator. The idea behind the construction of the smooth charts is to use the radii polynomials to verify the hypotheses of the uniform contraction principle over a simplex. The construction of the manifold is globalized by proving smoothness along the edge of adjacent simplices. We apply the method to compute portions of a two- dimensional manifold of equilibria of the Cahn-Hilliard equation.

The following links provide the programs sufficient to carry out the main result of the above mentioned work.


The paper in its pdf form can be found here.

Here are files and MATLAB codes used to perform the proof of Theorem 5.6 and Theorem 5.8.
The codes use the interval arithmetic package Intlab.