Matrix-Free Numerical Torus Bifurcation of Periodic Orbits
E. L. Allgower, U. Garbotz and K. Georg
Colorado State University
We describe a numerical continuation method for tracing branches of periodic solutions of dynamical systems in a matrix-free context i.e., Jacobians are only implemented as actions. This enables us to handle large systems, such as those arising from discretization of PDEs. Of particular interest is the detection and precise numerical approximation of bifurcation points along branches: especially period-doubling and torus bifurcation points. This will also be done in a matrix-free context combining Arnoldi iterations (to obtain coarse information) with the calculation of suitable test functions (for precise approximations). We illustrate the method with the one- and two-dimensional Brusselator.