Continuation of invariant subspaces for large bifurcation problems

Abstract. We summarize an algorithm developed in [17] for computing a smooth orthonormal basis for an invariant subspace of a parameter-dependent matrix, and describe how to extend it for numerical bifurcation analysis. We adapt the continued subspace to track behavior relevant to bifurcations, and use projection methods to deal with large problems. To test our ideas, we have integrated our code into Matcont, a program for numerical continuation and bifurcation analysis. 1. Introduction. Parameter-dependent Jacobia