Continuation of Invariant Subspaces for Large Bifurcation Problems

Authors:
Bindel, David Samuel
Demmel, James
Friedman, Mark
Technical Report Identifier: EECS-2006-13
February 13, 2006

Abstract: We summarize an algorithm 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.