Refine
Year of publication
- 2011 (1)
Language
- English (1)
Keywords
- Hamiltonian matrix (1)
- Lagrangian subspace (1)
- graph subspace (1)
- matrix pencil (1)
- optimal control (1)
- structure-preserving doubling algorithm (1)
- symplectic matrix (1)
Project
- C (1) (remove)
Application Area
- C (1)
We derive a new representation of Lagrangian subspaces in the form
%
\[
{\mathrm Im}\Pi^T [I,X]^T,
\]
%
where $\Pi$ is a symplectic matrix which is the product of a permutation matrix and a real orthogonal diagonal matrix, and $X$ satisfies
%
\[
\abs{X_{ij}} \leq \begin{cases}1 & \text{if $i=j$,}\\ \sqrt{2} & \text{if $i\neq j$.} \end{cases}
\]
%
This representation allows to limit element growth in the context of doubling algorithms
for the computation of Lagrangian subspaces and the solution of Riccati equations.
It is shown that a simple doubling algorithm using this representation can reach full machine accuracy on a wide range of problems, obtaining invariant subspaces of the same quality as those computed by the state-of-the-art algorithms based on orthogonal transformations.
The same idea carries over to representations of arbitrary subspaces and can be used
for other types of structured pencils.