Continuation of Invariant Subspaces in Large Bifurcation Problems

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.

[1]  G. Stewart.  Error and Perturbation Bounds for Subspaces Associated with Certain Eigenvalue Problems , 1973 .

[2]  C. Chien,et al.  Multiple bifurcations in a reaction-diffusion problem , 1998 .

[3]  Gautam M. Shroff,et al.  Stabilization of unstable procedures: the recursive projection method , 1993 .

[4]  H. Keller,et al.  Homotopy Method for the Large, Sparse, Real Nonsymmetric Eigenvalue Problem , 1997, SIAM J. Matrix Anal. Appl..

[5]  Luca Dieci,et al.  Continuation of invariant subspaces , 2001, Numer. Linear Algebra Appl..

[6]  James Demmel,et al.  LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

[7]  A. Spence,et al.  The numerical analysis of bifurcation problems with application to fluid mechanics , 2000, Acta Numerica.

[8]  Timo Eirola,et al.  On Smooth Decompositions of Matrices , 1999, SIAM J. Matrix Anal. Appl..

[9]  Alan Edelman,et al.  The Geometry of Algorithms with Orthogonality Constraints , 1998, SIAM J. Matrix Anal. Appl..

[10]  G. Dangelmayr.  Degenerate bifurcations near a double eigenvalue in the Brusselator , 1987, The Journal of the Australian Mathematical Society. Series B. Applied Mathematics.

[11]  D. Hu,et al.  Krylov-Subspace Methods for the Sylvester Equation , 2001 .

[12]  Y. Kuznetsov.  Elements of applied bifurcation theory (2nd ed.) , 1998 .

[13]  Jack Dongarra,et al.  Templates for the Solution of Algebraic Eigenvalue Problems , 2000, Software, environments, tools.

[14]  J. Brandts.  The Riccati algorithm for eigenvalues and invariant subspaces of matrices with inexpensive action , 2003 .

[15]  W. Beyn,et al.  Chapter 4 – Numerical Continuation, and Computation of Normal Forms , 2002 .

[16]  James Demmel,et al.  Computing Connecting Orbits via an Improved Algorithm for Continuing Invariant Subspaces , 2000, SIAM J. Sci. Comput..

[17]  I. Prigogine,et al.  Symmetry Breaking Instabilities in Dissipative Systems. II , 1968 .

[18]  P. Ashwin,et al.  A numerical Liapunov-Schmidt method with applications to Hopf bifurcation on a square , 1995 .

[19]  J. A. Kuznecov.  Elements of applied bifurcation theory , 1998 .

[20]  D. Sorensen.  IMPLICITLY RESTARTED ARNOLDI/LANCZOS METHODS FOR LARGE SCALE EIGENVALUE CALCULATIONS , 1996 .

[21]  R. Lehoucq,et al.  Linear stability of flow in a differentially heated cavity via large‐scale eigenvalue calculations , 2004 .

[22]  Willy Govaerts,et al.  Continuation of Homoclinic Orbits in Matlab , 2005, International Conference on Computational Science.

[23]  W. Govaerts.  Computation of Singularities in Large Nonlinear Systems , 1997 .

[24]  Robert E. Mahony,et al.  Cubically Convergent Iterations for Invariant Subspace Computation , 2004, SIAM J. Matrix Anal. Appl..

[25]  W. Govaerts,et al.  MIXED BLOCK ELIMINATION FOR LINEAR-SYSTEMS WITH WIDER BORDERS. , 1993 .

[26]  J. H. Brandts,et al.  Matlab code for sorting real Schur forms , 2002, Numer. Linear Algebra Appl..

[27]  M. Golubitsky,et al.  Singularities and groups in bifurcation theory , 1985 .

[28]  Tosio Kato.  Perturbation theory for linear operators , 1966 .

[29]  Willy Govaerts.  Stable solvers and block elimination for bordered systems , 1991 .

[30]  W. Beyn,et al.  Continuation of Low-Dimensional Invariant Subspaces in Dynamical Systems of Large Dimension , 2001 .

[31]  P. Ashwin,et al.  A Hopf bifurcation with Robin boundary conditions , 1994 .

[32]  H. Shari.  Collocation Methods for Continuation Problems in Nonlinear Elliptic Pdes , .

[33]  J. Varah.  On the Separation of Two Matrices , 1979 .

[34]  Willy Govaerts,et al.  Numerical methods for bifurcations of dynamical equilibria , 1987 .

[35]  N. Higham.  MATRIX NEARNESS PROBLEMS AND APPLICATIONS , 1989 .

[36]  C. Chien,et al.  Numerical Continuation at Double Bifurcation Points of a Reaction–Diffusion Problem , 1998 .

[37]  A. Bunse-Gerstner,et al.  Numerical computation of an analytic singular value decomposition of a matrix valued function , 1991 .

[38]  J. Guckenheimer,et al.  Defining Functions for Multiple Hopf Bifurcations , 1997 .

[39]  James Demmel,et al.  Three methods for refining estimates of invariant subspaces , 1987, Computing.

[40]  G. Stewart,et al.  Matrix Perturbation Theory , 1990 .

[41]  Mark J. Friedman.  Improved Detection of bifurcations in Large nonlinear Systems via the Continuation of Invariant Subspaces Algorithm , 2001, Int. J. Bifurc. Chaos.

[42]  K. Georg.  MATRIX-FREE NUMERICAL CONTINUATION AND BIFURCATION* , 2001 .

[43]  James Demmel,et al.  Bifurcation Analysis of Large Equilibrium Systems in Matlab , 2005, International Conference on Computational Science.

[44]  L. Dieci,et al.  Point-to-Periodic and Periodic-to-Periodic Connections , 2004 .

[45]  Jaroslav Bosek.  Continuation of invariant subspaces in bifurcation problems , 2002 .

[46]  Bernold Fiedler,et al.  Ergodic theory, analysis, and efficient simulation of dynamical systems , 2001 .

[47]  L. Reichel,et al.  Krylov-subspace methods for the Sylvester equation , 1992 .

[48]  Z. Mei.  Numerical Bifurcation Analysis for Reaction-Diffusion Equations , 2000 .

[49]  M. Golubitsky,et al.  Bifurcation analysis near a double eigenvalue of a model chemical reaction , 1981 .

[50]  Matemática,et al.  Society for Industrial and Applied Mathematics , 2010 .

[51]  Luca Dieci,et al.  Continuation of eigendecomposition , 2003, Future Gener. Comput. Syst..

[52]  G. Golub,et al.  A Hessenberg-Schur method for the problem AX + XB= C , 1979 .