Numerical implementation of complex orthogonalization, parallel transport on Stiefel bundles, and analyticity

Daniele Avitabile*, Thomas J. Bridges

*Corresponding author for this work

Research output: Contribution to JournalArticleAcademicpeer-review

Abstract

Numerical integration of complex linear systems of ODEs depending analytically on an eigenvalue parameter are considered. Complex orthogonalization, which is required to stabilize the numerical integration, results in non-analytic systems. It is shown that properties of eigenvalues are still efficiently recoverable by extracting information from a non-analytic characteristic function. The orthonormal systems are constructed using the geometry of Stiefel bundles. Different forms of continuous orthogonalization in the literature are shown to correspond to different choices of connection one-form on the Stiefel bundle. For the numerical integration, Gauss-Legendre Runge-Kutta algorithms are the principal choice for preserving orthogonality, and performance results are shown for a range of GLRK methods. The theory and methods are tested by application to example boundary value problems including the Orr-Sommerfeld equation in hydrodynamic stability.

Original languageEnglish
Pages (from-to)1038-1047
Number of pages10
JournalPhysica D: Nonlinear Phenomena
Volume239
Issue number12
DOIs
Publication statusPublished - 15 Jun 2010
Externally publishedYes

Keywords

  • Constraints
  • Geometric integration
  • Newton's method
  • Orthonormalization
  • Poincaré index
  • Stiefel bundles

Fingerprint

Dive into the research topics of 'Numerical implementation of complex orthogonalization, parallel transport on Stiefel bundles, and analyticity'. Together they form a unique fingerprint.

Cite this