---
title: Cayley-Based Retractions
url: https://www.emergentmind.com/topics/cayley-based-retractions
type: topic
---

# Cayley-Based Retractions

A Cayley-based retraction is a specific class of retraction map used for locally parametrizing points on matrix manifolds, most prominently the real symplectic Stiefel manifold, in computational geometric and optimization contexts. Retractions approximate exponential maps, enabling efficient first-order embeddings between a tangent space and the manifold, and are pivotal in Riemannian algorithms where movements off the manifold must be projected back. On the symplectic Stiefel manifold, the Cayley retraction stands out due to its closed-form inverse and computational efficiency, balancing accuracy and stability while preserving symplectic structure [2605.06399].

## 1. Mathematical Setting and Definition

The Cayley-based retraction is defined on the real symplectic Stiefel manifold
$$
SpSt(2n,2p) = \{ Y \in \mathbb{R}^{2n \times 2p} \mid Y^T J_n Y = J_p \}, \;\; J_n = \begin{pmatrix}0 & I_n \\ -I_n & 0\end{pmatrix},
$$
where $J_n$ is the standard symplectic matrix and $J_p$ its $2p$-dimensional analog. For a base point $Y \in SpSt(2n,2p)$ and tangent vector $S \in T_Y SpSt$, define the $2p \times 2p$ Hamiltonian matrix
$$
A := Y^+S = J_p^T Y^T J_n S \in sp(2p),
$$
where the Moore–Penrose pseudoinverse $Y^+$ is used. $sp(2p)$ denotes the Lie algebra of Hamiltonian matrices satisfying $(J_pA)^T = J_pA$.

The Cayley map for Hamiltonian matrices is
$$
\mathrm{Cay}(A) = (I_{2p} - \frac{1}{2}A)^{-1}(I_{2p} + \frac{1}{2}A) \in Sp(2p).
$$
The Cayley retraction from the tangent bundle to the manifold is then
$$
R_Y^{\mathrm{Cay}}(S) = Y\,\mathrm{Cay}(A) + (I_{2n} - Y Y^+)S \in SpSt(2n,2p), \quad A = Y^+S.
$$
This mapping provides a computationally tractable surrogate to the Riemannian exponential map with first-order accuracy [2605.06399].

## 2. Closed-Form Inverse and Domain

One of the distinguishing features of the Cayley retraction on $SpSt$ is its closed-form inverse. Given $\widetilde Y \in SpSt(2n,2p)$ in a neighborhood of $Y$ such that $W := Y^+\widetilde{Y} \in Sp(2p)$ and $\det(W + I_{2p}) \neq 0$, the inverse, sometimes called the Cayley logarithm, is
$$
A = \mathrm{Cay}^{-1}(W) = 2(W - I_{2p})(W + I_{2p})^{-1} \in sp(2p).
$$
The tangent vector is then recovered as
$$
S = R_Y^{\mathrm{Cay}^{-1}}(\widetilde Y) = Y A + (I_{2n} - Y Y^+)\widetilde{Y}.
$$
Domain conditions for invertibility require
$$
\operatorname{spec}(Y^+\widetilde{Y}) \cap \{-1\} = \emptyset, \qquad \operatorname{spec}(Y^+S) \cap \{\pm 2\} = \emptyset,
$$
ensuring the relevant linear systems are solvable. This feature enables practical implementations requiring mapping between ambient and tangent representations [2605.06399].

## 3. Structural and Accuracy Properties

The Cayley-based retraction is characterized by several important mathematical properties:
- **First-order accuracy:** The expansion
  $$
  R_Y^{\mathrm{Cay}}(tS) = Y + tS + O(t^2)
  $$
  shows that it correctly recovers the tangent step to first order, matching the exponential map up to $O(t^2)$.
- **Symplectic preservation:** For $\Delta = R_Y^{\mathrm{Cay}}(S)$,
  $$
  \Delta^T J_n \Delta = (\mathrm{Cay}(A))^T J_p \mathrm{Cay}(A) = J_p,
  $$
  so symplectic structure is maintained. The cross-term vanishes by the algebraic properties of $Y$ and $S$.
- **Hamiltonicity and orthogonality:** The construction guarantees $A \in sp(2p)$, with $J_pA$ symmetric, $Y^+Y = I_{2p}$, and $(I - Y Y^+)^2 = (I - Y Y^+)$, ensuring the relevant invariances and projection properties.

These properties make the Cayley retraction particularly suitable for geometric optimization and integration schemes where structure preservation is critical [2605.06399].

## 4. Computational Complexity and Numerical Stability

The cost profile of the Cayley retraction is determined by several dominating operations:
- Computing $A = Y^+S$: $O(n p^2)$.
- Forming $\mathrm{Cay}(A)$: $O(p^3)$, due to matrix inversion or solving two linear systems.
- Multiplying $Y \mathrm{Cay}(A)$: $O(n p^2)$.
- Computing the orthogonal complement $(I - Y Y^+)S$: $O(n p^2)$.

Thus, overall complexity per evaluation is
$$
O(p^3 + n p^2)
$$
for both the forward and inverse maps. Stability requires maintaining distance from singularities (eigenvalue $-1$ of $W$ or $\pm2$ of $A$), which in practice is managed by monitoring $\|A\|$ or the minimal singular value of $(W + I)$. The Cayley map is efficient for moderate step sizes due to its avoidance of computationally intensive decompositions [2605.06399].

## 5. Comparison with Polar-Factor Retraction

Zimmermann (2024) introduced a contrasting retraction, the symplectic polar-factor (or "polar-light") retraction, defined as
$$
R^f_U(D) = \left( U \exp(U^+ D) + (I - U U^+)D \right) G_f,
$$
where $G_f \in Sh(2p)$ is a skew-Hamiltonian symplectifier obtained from a Schur decomposition and matrix square root. Comparing the Cayley and polar-factor retractions:

| Retraction     | Cost per Step               | Residual $\|Y^+_{\text{out}} Y_{\text{out}} - I\|$      |
|----------------|----------------------------|-----------------------------------------|
| Cayley         | $O(p^3 + n p^2)$           | $\approx 10^{-13}$                      |
| Polar-factor   | $O(p^3 + n p^2)$ (larger constant) | $\approx 10^{-12}$                      |

In practical scenarios ($n \gg p$), Cayley retraction is about three times faster and delivers slightly better symplectic residuals for moderate step sizes. The polar-factor’s use of the exact matrix exponential confers better robustness for larger steps, at the expense of higher per-iteration cost [2605.06399].

## 6. Implementation Details

MATLAB-style pseudocode for the Cayley retraction and its inverse, aligning with Oviedo–Herrera (2023) and Gao et al. (2020), is as follows:

```matlab
function Ynew = CayleyRetraction(Y, S, Jn, Jp)
  A = Jp'*Y'*Jn * S;                 % 2p×2p Hamiltonian
  M = (I2p - 0.5*A)\(I2p + 0.5*A);   % Cayley(A)
  Ynew = Y*M + (eye(2n)-Y*Y^+)*S;
end
```
and for the inverse:
```matlab
function S = CayleyInverseRetraction(Y, Ytil, Jn, Jp)
  W = Jp'*Y'*Jn * Ytil;              % 2p×2p symplectic
  A = 2*(W - I2p)/(W + I2p);         % Cay^-1(W)
  S = Y*A + (eye(2n)-Y*Y^+) * Ytil;
end
```

A plausible implication is that these routines are immediately applicable in manifold optimization libraries where efficient, structure-preserving updates are required.

## 7. References and Further Directions

Relevant contributions to the theory and implementation of Cayley retractions on $SpSt$ are provided by Gao-Son-Absil-Stykel (SIAM J. Optim. 2021), Bendokat–Zimmermann [arXiv:2108.12447], and Oviedo–Herrera (2023). Recent work by Zimmermann has advanced alternatives, particularly the polar-factor retraction, and benchmarks these schemes for cost and accuracy [2605.06399]. Further exploration is warranted in the direction of large-step stability, hybrid schemes, and the incorporation of such retractions into large-scale optimization and numerical integration routines.

Source: https://www.emergentmind.com/topics/cayley-based-retractions