Flag manifold

Manifolds.Flag β€” Type
Flag{T,d} <: AbstractDecoratorManifold{ℝ}

Flag manifold of $d$ subspaces of $ℝ^N$ [YWL21]. By default the manifold uses the Stiefel coordinates representation, embedding it in the Stiefel manifold. The other available representation is an embedding in OrthogonalMatrices. It can be utilized using OrthogonalPoint and OrthogonalTangentVector wrappers.

Tangent space is represented in the block-skew-symmetric form.

Constructor

Flag(N, n_1, n_2, ..., n_d; parameter::Symbol=:type)

Generate the manifold $\operatorname{Flag}(n_1, n_2, ..., n_d; N)$ of subspaces

\[𝕍_1 βŠ† 𝕍_2 βŠ† β‹― βŠ† V_d, \quad \operatorname{dim}(𝕍_i) = n_i\]

where $𝕍_i$ for $i ∈ 1, 2, …, d$ are subspaces of $ℝ^N$ of dimension $\operatorname{dim} 𝕍_i = n_i$.

parameter: whether a type parameter should be used to store n. By default size is stored in type. Value can either be :field or :type.

source
Manifolds.OrthogonalPoint β€” Type
OrthogonalPoint <: AbstractManifoldPoint

A type to represent points on a manifold Flag in the orthogonal coordinates representation, i.e. a rotation matrix.

source
Manifolds.ZeroTuple β€” Type
ZeroTuple

Internal structure for representing shape of a Flag manifold. Behaves like a normal tuple, except at index zero returns value 0.

source

The flag manifold represented as points on the Stiefel manifold

ManifoldsBase.check_vector β€” Method
check_vector(M::Flag, p::AbstractMatrix, X::AbstractMatrix; kwargs... )

Check whether X is a tangent vector to point p on the Flag manifold M $\operatorname{Flag}(n_1, n_2, ..., n_d; N)$ in the Stiefel representation, i.e. that X is a matrix of the form

\[X = \begin{bmatrix} 0 & B_{1,2} & β‹― & B_{1,d} \\ -B_{1,2}^\mathrm{T} & 0 & β‹― & B_{2,d} \\ \vdots & \vdots & β‹± & \vdots \\ -B_{1,d}^\mathrm{T} & -B_{2,d}^\mathrm{T} & β‹― & 0 \\ -B_{1,d+1}^\mathrm{T} & -B_{2,d+1}^\mathrm{T} & β‹― & -B_{d,d+1}^\mathrm{T} \end{bmatrix}\]

where $B_{i,j} ∈ ℝ^{(n_i - n_{i-1}) Γ— (n_j - n_{j-1})}$, for $1 ≀ i < j ≀ d+1$.

source
ManifoldsBase.project β€” Method
project(::Flag, p, X)

Project vector X in the Euclidean embedding to the tangent space at point p on Flag manifold. The formula reads [YWL21]:

\[Y_i = X_i - (p_i p_i^{\mathrm{T}}) X_i + \sum_{j \neq i} p_j X_j^{\mathrm{T}} p_i\]

for $i$ from 1 to $d$ where the resulting vector is $Y = [Y_1, Y_2, …, Y_d]$ and $X = [X_1, X_2, …, X_d]$, $p = [p_1, p_2, …, p_d]$ are decompositions into basis vector matrices for consecutive subspaces of the flag.

source
ManifoldsBase.retract β€” Method
retract(M::Flag, p, X, ::PolarRetraction)

Compute the SVD-based retraction PolarRetraction on the Flag M. With $USV = p + X$ the retraction reads

\[\operatorname{retr}_p X = UV^\mathrm{H},\]

where $\cdot^{\mathrm{H}}$ denotes the complex conjugate transposed or Hermitian.

source

The flag manifold represented as orthogonal matrices

ManifoldsBase.check_vector β€” Method
check_vector(M::Flag, p::OrthogonalPoint, X::OrthogonalTangentVector; kwargs... )

Check whether X is a tangent vector to point p on the Flag manifold M $\operatorname{Flag}(n_1, n_2, ..., n_d; N)$ in the orthogonal matrix representation, i.e. that X is block-skew-symmetric with zero diagonal:

\[X = \begin{bmatrix} 0 & B_{1,2} & β‹― & B_{1,d+1} \\ -B_{1,2}^\mathrm{T} & 0 & β‹― & B_{2,d+1} \\ \vdots & \vdots & β‹± & \vdots \\ -B_{1,d+1}^\mathrm{T} & -B_{2,d+1}^\mathrm{T} & β‹― & 0 \end{bmatrix}\]

where $B_{i,j} ∈ ℝ^{(n_i - n_{i-1}) Γ— (n_j - n_{j-1})}$, for $1 ≀ i < j ≀ d+1$.

source
ManifoldsBase.project β€” Method
project(M::Flag, p::OrthogonalPoint, X::OrthogonalTangentVector)

Project vector X to tangent space at point p from Flag manifold M $\operatorname{Flag}(n_1, n_2, ..., n_d; N)$, in the orthogonal matrix representation. It works by first projecting X to the space of SkewHermitianMatrices and then setting diagonal blocks to 0:

\[X = \begin{bmatrix} 0 & B_{1,2} & β‹― & B_{1,d+1} \\ -B_{1,2}^\mathrm{T} & 0 & β‹― & B_{2,d+1} \\ \vdots & \vdots & β‹± & \vdots \\ -B_{1,d+1}^\mathrm{T} & -B_{2,d+1}^\mathrm{T} & β‹― & 0 \end{bmatrix}\]

where $B_{i,j} ∈ ℝ^{(n_i - n_{i-1}) Γ— (n_j - n_{j-1})}$, for $1 ≀ i < j ≀ d+1$.

source