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). In the case when A(s) is the Jacobian matrix for a system that comes from a spatial discretization of a partial differential equation, it will typically be large and sparse. Cl matcont is a user-friendly matlab package for the study of dynamical systems and their bifurcations. We incorporate the CIS algorithm into cl matcont to extend its functionality to large scale bifurcation computations via subspace reduction.
David Bindel, James Demmel, Mark J. Friedman, Will