Numerical computation of bifurcations in large equilibrium systems in MATLAB.

The Continuation of Invariant Subspaces (CIS) algorithm produces a smoothly-varying basis for an invariant subspace R(s) of a parameter-dependent matrix A(s). We have incorporated the CIS algorithm into Cl_matcont, a Matlab package for the study of dynamical systems and their bifurcations. Using subspace reduction, we extend the functionality of Cl_matcont to large-scale computations of bifurcations of equilibria. In this paper, we describe the algorithms and functionality of the resulting Matlab bifurcation package Cl_matcontL. The novel features include: new CIS-based, continuous, well-scaled test functions for codimension 1 and 2 bifurcations; detailed description of locators for large problems; and examples of bifurcation analysis in large sparse problems