Next Article in Journal
Teaching CAPM for a Pre-Finance Graduate Program at the STEM Undergraduate Level: Linear Algebra Perspective
Next Article in Special Issue
Some Conditions on Trans-Sasakian Manifolds to Be Homothetic to Sasakian Manifolds
Previous Article in Journal
Differential Games for Fractional-Order Systems: Hamilton–Jacobi–Bellman–Isaacs Equation and Optimal Feedback Strategies
Previous Article in Special Issue
Representations of Rectifying Isotropic Curves and Their Centrodes in Complex 3-Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Principal Bundle Structure of Matrix Manifolds

1
Department of Computer Science and Mathematics, Ecole Centrale de Nantes, 1 Rue de la Noë, BP 92101, CEDEX 3, 44321 Nantes, France
2
Departamento de Matemáticas, Física y Ciencias Tecnológicas, Universidad CEU Cardenal Herrera, CEU Universities, San Bartolomé 55, 46115 Alfara del Patriarca, Spain
3
ESI International Chair@CEU-UCH, Universidad Cardenal Herrera-CEU, CEU Universities San Bartolomé 55, 46115 Alfara del Patriarca, Spain
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(14), 1669; https://doi.org/10.3390/math9141669
Submission received: 11 June 2021 / Revised: 13 July 2021 / Accepted: 14 July 2021 / Published: 15 July 2021
(This article belongs to the Special Issue Differential Geometry: Structures on Manifolds and Their Applications)

Abstract

:
In this paper, we introduce a new geometric description of the manifolds of matrices of fixed rank. The starting point is a geometric description of the Grassmann manifold G r ( R k ) of linear subspaces of dimension r < k in R k , which avoids the use of equivalence classes. The set G r ( R k ) is equipped with an atlas, which provides it with the structure of an analytic manifold modeled on R ( k r ) × r . Then, we define an atlas for the set M r ( R k × r ) of full rank matrices and prove that the resulting manifold is an analytic principal bundle with base G r ( R k ) and typical fibre GL r , the general linear group of invertible matrices in R k × k . Finally, we define an atlas for the set M r ( R n × m ) of non-full rank matrices and prove that the resulting manifold is an analytic principal bundle with base G r ( R n ) × G r ( R m ) and typical fibre GL r . The atlas of M r ( R n × m ) is indexed on the manifold itself, which allows a natural definition of a neighbourhood for a given matrix, this neighbourhood being proved to possess the structure of a Lie group. Moreover, the set M r ( R n × m ) equipped with the topology induced by the atlas is proven to be an embedded submanifold of the matrix space R n × m equipped with the subspace topology. The proposed geometric description then results in a description of the matrix space R n × m , seen as the union of manifolds M r ( R n × m ) , as an analytic manifold equipped with a topology for which the matrix rank is a continuous map.

1. Introduction

Low-rank matrices appear in many applications involving high-dimensional data. Low-rank models are commonly used in statistics, machine learning or data analysis (see [1] for a recent survey). Furthermore, low-rank approximation of matrices is the cornerstone of many modern numerical methods for high-dimensional problems in computational science, such as model-order-reduction methods for dynamical systems or parameter-dependent or stochastic equations [2,3,4,5].
These applications yield problems of approximation or optimization in the sets of matrices with fixed rank:
M r ( R n × m ) = { Z R n × m : rank ( Z ) = r } .
Fixed-rank matrices appear also in the theory of characteristics of Partial Differential Equations and Monge-Ampère equations [6]. More precisely, it has been proven [6,7] that Monge-Ampère equations with n independent variables and of Goursat-type are in one-to-one correspondence with the set { Z M r ( R n × n ) : r 2 } . Thus, the parabolic or hyperbolic nature of the Monge-Ampère equation is related to the rank of such matrices.
In [8,9], the authors point out that Algebraic Geometry appears as a natural tool in study of the set M r ( R n × m ) . We wish to mention the papers [10,11,12] that raise the natural question of how large these matrix spaces are.
A usual geometric approach is to endow the set M r ( R n × m ) with the structure of a Riemannian manifold [13,14], which is seen as an embedded submanifold of R n × m equipped with the topology τ R n × m given by matrix norms. Standard algorithms then work in the ambient matrix space R n × m and do not rely on an explicit geometric description of the manifold using local charts (see, e.g., [15,16,17,18]). However, the matrix rank considered as a map is not continuous for the topology τ R n × m , which can yield undesirable numerical issues.
The purpose of this paper is to propose a new geometric description of the sets of matrices with fixed rank, which is amenable for numerical use, and relies on the natural parametrization of matrices in M r ( R n × m ) given by
Z = U G V T ,
where U R n × r and V R m × r are matrices with full rank r < min { n , m } and G R r × r is a non singular matrix. The set M r ( R n × m ) is here endowed with the structure of analytic principal bundle with an explicit description of local charts. This results in a description of the matrix space R n × m as an analytic manifold with a topology induced by local charts that is different from τ R n × m and for which the rank is a continuous map. Note that the representation (1) of a matrix Z is not unique because Z = ( U P ) ( P 1 G P T ) ( V P 1 ) T holds for every invertible matrix P in R r × r . An argument used to dodge this undesirable property is the possibility to uniquely define a tangent space (see for example Section 2.1 in [18]), which is a prerequisite for standard algorithms on differentiable manifolds. The geometric description proposed in this paper avoids this undesirable property. Indeed, the system of local charts for the set M r ( R n × m ) is indexed on the set itself. This allows a natural definition of a neighbourhood for a matrix where all matrices admit a unique representation.
The present work opens the route for new numerical methods for optimization and dynamical low-rank approximation with algorithms working in local coordinates and avoiding the use of a Riemannian structure. In [19], such a framework is introduced for generalising iterative methods in optimization from Euclidean space to manifolds, which ensures that local convergence rates are preserved. Recently, a splitting algorithm relying on the geometric description of the set of fixed rank matrices proposed in this paper has been introduced for dynamical low-rank approximation [20].
The introduction of a principal bundle representation of matrix manifolds is also motivated by the importance of this geometric structure in the concept of gauge potential in physics [21].
Note that the proposed geometric description has a natural extension to the case of fixed-rank operators on infinite dimensional spaces and is consistent with the geometric description of manifolds of tensors with fixed rank proposed by Falcó, Hackbush and Nouy [22] in a tensor Banach space framework.
Before introducing the main results and outline of the paper, we recall some elements of geometry.

1.1. Elements of Geometry

In this paper, we follow the approach of Serge Lang [23] for the definition of a manifold M . In this framework, a set M is equipped with an atlas which gives M the structure of a topological space, with a topology induced by local charts, and the structure of differentiable manifold compatible with this topology. More precisely, the starting point is the definition of a collection of non-empty subsets U α M , with α in a set A, such that { U α } α A is a covering of M . The next step is the explicit construction for any α A of a local chart φ α which is a bijection from U α to an open set X α of the finite dimensional space R N α such that for any pair α , α M such that U α U α , the following properties hold:
(i)
φ α ( U α U α ) and φ α ( U α U α ) are open sets in X α and X α respectively, and
(ii)
the map
φ α φ α 1 : φ α ( U α U α ) φ α ( U α U α )
is a C p differentiable diffeomorphism, with p N { } or p = ω when the map is analytic.
Under the above assumptions, the set A : = { ( U α , φ α ) : α A } is an atlas which endows M with a structure of C p manifold. Then, we can say that ( M , A ) is a C p manifold, or an analytic manifold when p = ω . A consequence of condition ( i i ) is that when U α U α holds for α , α A , then N α = N α . In the particular case where N α = N for all α A , we say that ( M , A ) is a C p manifold modelled on R N . Otherwise, we say that it is a manifold not modelled on a particular finite-dimensional space. A paradigmatic example is the Grassmann manifold G ( R k ) of all linear subspaces of R k , such that
G ( R k ) = 0 r k G r ( R k ) ,
where G 0 ( R k ) = { 0 } and G k ( R k ) = { R k } are trivial manifolds and G r ( R k ) is a manifold modelled on the linear space R ( k r ) × r for 0 < r < k . Consequently, G ( R k ) is a manifold not modelled on a particular finite-dimensional space.
The atlas also endows M with a topology given by
τ A : = φ α 1 ( O ) : α A and O an open set in X α ,
which makes ( M , τ A ) a topological space where each local chart
φ α : ( U α , τ A | U α ) ( X α , τ R N α | X α ) ,
considered as a map between topological spaces is a homeomorphism. (Here ( X , τ ) denotes a topological space, and if X X , then τ | X denotes the subspace topology.)

1.2. Main Results and Outline

Our first remark is that the matrix space R n × m is an analytic manifold modelled on itself, and its geometric structure is fully compatible with the topology τ R n × m induced by a matrix norm. In this paper, we define an atlas on M r ( R n × m ) , which gives this set the structure of an analytic manifold, with a topology induced by the atlas fully compatible with the subspace topology τ R n × m | M r ( R n × m ) . This implies that M r ( R n × m ) is an embedded submanifold of the matrix manifold R n × m modelled on itself. (Note that the set M 0 ( R n × m ) = { 0 } is a trivial manifold, which is trivially embedded in R n × m .) For the topology τ R n × m , the matrix rank considered as a map is not continuous but only lower semi-continuous. However, if R n × m is seen as the disjoint union of sets of matrices with fixed rank,
R n × m = 0 r min { n , m } M r ( R n × m ) ,
then R n × m has the structure of an analytic manifold not modelled on a particular finite-dimensional space equipped with a topology
τ R n × m * = 0 r min { n , m } τ R n × m | M r ( R n × m ) ,
which is not equivalent to τ R n × m , and for which the matrix rank is a continuous map.
Note that in the case where r = n = m , the set M n ( R n × n ) coincides with the general linear group GL n of invertible matrices in R n × n , which is an analytic manifold trivially embedded in R n × n . In all other cases are addressed in this paper, our geometric description of M r ( R n × m ) relies on a geometric description of the Grassmann manifold G r ( R k ) , with k = n or m.
Therefore, we start in Section 2 by introducing a geometric description of G r ( R k ) . A classical approach consists of describing G r ( R k ) as the quotient manifold M r ( R k × r ) / GL r of equivalent classes of full-rank matrices Z in M r ( R k × r ) with the same column space col k , r ( Z ) . Here, we avoid the use of equivalent classes and provide an explicit description of an atlas A k , r = { ( U Z , φ Z ) } Z M r ( R k × r ) for G r ( R k ) , with local chart
φ Z : U Z R ( k r ) × r , φ Z 1 ( X ) = col k , r ( Z + Z X ) ,
where Z R k × ( k r ) is such that Z T Z = 0 (see Remark 1 for a practical choice) and col k , r ( A ) denotes the column space of a matrix A R k × r , and we prove that the neighbourhood U Z has the structure of a Lie group. This parametrization of the Grassmann manifold is introduced in ([24] Section 2), but the authors do not elaborate on it.
Then, in Section 3, we consider the particular case of full-rank matrices. We introduce an atlas B k , r = { ( V Z , ξ Z ) } Z M r ( R k × r ) for the manifold M r ( R k × r ) of matrices with full rank r < k , with local chart
ξ Z : V Z R ( k r ) × r × GL r , ξ Z 1 ( X , G ) = ( Z + Z X ) G ,
and prove that M r ( R k × r ) is an analytic principal bundle with base G r ( R k ) and typical fibre GL r . Moreover, we prove that M r ( R k × r ) is an embedded submanifold of ( R k × r , τ R k × r * ) and that each of the neighbourhoods V Z have the structure of a Lie group.
Finally, in Section 4, we provide an analytic atlas B n , m , r = { ( U Z , θ Z ) } Z M r ( R n × m ) for the set M r ( R n × m ) of matrices Z = U G V T with rank r < min { n , m } , with local chart
θ Z : U Z R ( n r ) × r × R ( m r ) × r × GL r , θ Z 1 ( X , Y , H ) = ( U + U X ) H ( V + V Y ) ,
and we prove that M r ( R n × m ) is an analytic principal bundle with base G r ( R n ) × G r ( R m ) and typical fibre GL r . Then, we prove that M r ( R n × m ) is an embedded submanifold of ( R n × m , τ R n × m * ) and that each of the neighbourhoods U Z have the structure of a Lie group.

2. The Grassmann Manifold G r ( R k )

In this section, we present a geometric description of the Grassmann manifold G r ( R k ) of all subspaces of dimension r in R k , 0 < r < k ,
G r ( R k ) = { V R k : V is a linear subspace with dim ( V ) = r } ,
with an explicit description of local charts. We first introduce the surjective map
col k , r : M r ( R k × r ) G r ( R k ) , Z col k , r ( Z ) ,
where col k , r ( Z ) is the column space of the matrix Z, which is the subspace spanned by the column vectors of Z . Given V G r ( R k ) , there are infinitely many matrices Z such that col k , r ( Z ) = V . Given a matrix Z M r ( R k × r ) , the set of matrices in M r ( R k × r ) with the same column space as Z is
Z GL r : = { Z G : G GL r } .

2.1. An Atlas for G r ( R k )

For a given matrix Z in M r ( R k × r ) , we let Z M k r ( R k × ( k r ) ) be a matrix such that Z T Z = 0 , and we introduce an affine cross section
S Z : = { W M r ( R k × r ) : Z T W = Z T Z } ,
which has the following equivalent characterization.
Lemma 1.
The affine cross section S Z is characterized by
S Z = { Z + Z X : X R ( k r ) × r } ,
and the map
η Z : R ( k r ) × r S Z , X Z + Z X
is bijective.
Proof. 
We first observe that Z T ( Z + Z X Z ) = Z T Z for all X R ( k r ) × r , which implies that { Z + Z X : X R ( k r ) × r } S Z . For the other inclusion, we observe that if W S Z , then Z T W = Z T Z and hence W Z col k , r ( Z ) , the orthogonal subspace to col k , r ( Z ) in R k . Since col k , r ( Z ) = col k , k r ( Z ) , there exists X R ( k r ) × r such that W Z = Z X . Proving that η Z is bijective is straightforward. □
Proposition 1.
For each W M r ( R k × r ) such that det ( Z T W ) 0 , there exists a unique G W GL r such that
W GL r S Z = { W G W 1 }
holds, which means that the set of matrices with the same column space as W intersects S Z at the single point W G W 1 . Furthermore, G W = i d r if and only if W S Z .
Proof. 
By Lemma 1, a matrix A W GL r S Z is such that A = W G W 1 = Z + Z X for a certain G W GL r and a certain X R ( k r ) × r . Then Z T W G W 1 = Z T Z and G W is uniquely defined by G W = ( Z T Z ) 1 ( Z T W ) , which proves that W GL r S Z is the singleton { W G W 1 } , and G W = i d r if and only if W S Z .
Corollary 1.
For each Z M r ( R k × r ) , the map col k , r : S Z G r ( R k ) is injective.
Proof. 
Let us assume the existence of W , W ˜ S Z such that col k , r ( W ) = col k , r ( W ˜ ) . Then W = W ˜ by Proposition 1. □
Lemma 1 and Corollary 1 allow us to construct a system of local charts for G r ( R k ) by defining for each Z M r ( R k × r ) a neighbourhood of col k , r ( Z ) by
U Z : = col k , r ( S Z ) = { col k , r ( W ) : W S Z }
together with the bijective map
φ Z : = ( col k , r η Z ) 1 : U Z R ( k r ) × r
such that
φ Z 1 ( X ) = col k , r ( Z + Z X )
for X R ( k r ) × r . We denote by Z + the Moore–Penrose pseudo-inverse of the full rank matrix Z M r ( R r × k ) , defined by
Z + : = ( Z T Z ) 1 Z T M r ( R r × k ) .
It satisfies Z + Z = i d r and Z + Z = 0 . Moreover, Z Z + R k × k is the projection onto col k , r ( Z ) parallel to col k , r ( Z ) . Finally, we have the following result.
Theorem 1.
The collection A k , r : = { ( U Z , φ Z ) : Z M r ( R k × r ) } is an analytic atlas for G r ( R k ) and hence ( G r ( R k ) , A k , r ) is an analytic r ( k r ) -dimensional manifold modelled on R ( k r ) × r .
Proof. 
Clearly { U Z } Z M r ( R k × r ) is a covering of G r ( R k ) . Now let Z and Z ˜ be such that U Z U Z ˜ . Let V U Z such that V = φ Z 1 ( X ) = col k , r ( Z + Z X ) , with X R k × ( k r ) . We can write Z + Z X = ( Z ˜ + Z ˜ X ˜ ) G with G = Z ˜ + ( Z + Z X ) and X ˜ = Z ˜ + ( Z + Z X ) G 1 . Therefore, V = col k , r ( ( Z ˜ + Z ˜ X ˜ ) G ) = col k , r ( Z ˜ + Z ˜ X ˜ ) = φ Z ˜ 1 ( X ˜ ) U Z ˜ , which implies that U Z = U Z U Z ˜ . Therefore, φ Z ( U Z U Z ˜ ) = φ Z ( U Z ) = R k × ( n k ) is an open set. In the same way, we show that U Z ˜ = U Z U Z ˜ and φ Z ˜ ( U Z ) = R k × ( n k ) is an open set. Finally, the map φ Z ˜ φ Z 1 from R ( k r ) × r to R ( k r ) × r is given by φ Z ˜ φ Z 1 ( X ) = Z ˜ + ( Z + Z X ) G 1 , with G = Z ˜ + ( Z + Z X Z ) , which is clearly an analytic map. □
Remark 1.
A possible choice for Z satisfying Z T Z = 0 is Z = ( i d k Z Z + ) B where B M k r ( R k × ( k r ) ) is such that its column space is a complement of the column space of Z. In practice, we can determine a set of r linear independent rows of Z (see, e.g., [25,26]), with indices I, and then choose B such that ( B ) i , j = δ i , j if i I and 0 if i I , for 1 i k , 1 j k r . For a given X R ( k r ) × r , the computation of Z X does not require Z and has a complexity O ( r 2 k ) .

2.2. Lie Group Structure of Neighbourhoods U Z

Here we prove that each neighbourhood U Z of G r ( R k ) is a Lie group. For that, we first note that a neighbourhood U Z of G r ( R k ) can be identified with the set S Z through the application col k , r : S Z U Z . The next step is to identify S Z with a closed Lie subgroup of GL k , denoted by G Z , with associated Lie algebra g Z isomorphic to R r × ( k r ) , and such that the exponential map exp : g Z G Z is a diffeomorphism. (We recall that the matrix exponential exp : R k × k GL k is defined by exp ( A ) = n = 0 A n n ! . ) To this end, for a given Z M r ( R k × r ) , we introduce the vector space
g Z : = { Z X Z + : X R ( k r ) × r } R k × k .
The following proposition proves that g Z is a commutative subalgebra of R k × k .
Proposition 2.
For all X , X ˜ R ( k r ) × r ,
( Z X Z + ) ( Z X ˜ Z + ) = 0
holds, and g Z is a commutative subalgebra of R k × k . Moreover,
exp ( Z X Z + ) = i d k + Z X Z + ,
exp ( Z X Z + ) Z = Z + Z X ,
and
exp ( Z X Z + ) Z = Z
hold for all X R ( k r ) × r .
Proof. 
Since ( Z X Z + ) ( Z X ˜ Z + ) = 0 holds for all X , X ˜ R ( k r ) × r , the vector space g Z is a closed subalgebra of the matrix unitary algebra R k × k . As a consequence, ( Z X Z + ) p = 0 holds for all X R ( k r ) × r and all p 2 , which proves (6). We directly deduce (7) using Z Z + = i d r and (8) using Z + Z = 0 . □
From Proposition 2 and the definition of S Z , we obtain the following results.
Corollary 2.
The affine cross section S Z satisfies
S Z = { exp ( Z X Z + ) Z : X R ( k r ) × r } ,
and
[ exp ( Z X Z + ) Z | Z ] GL k
for all X R ( k r ) × r , where the brackets · | · are used for matrix concatenation.
Proof. 
From Proposition 2 and (4), we obtain (9) and we can write
[ exp ( Z X Z + ) Z | Z ] = [ exp ( Z X Z + ) Z | exp ( Z X Z + ) Z ] = exp ( Z X Z + ) [ Z | Z ] .
Since exp ( Z X Z + ) , [ Z | Z ] GL k , (10) follows. □
Now we need to introduce the following definition and proposition (see ([27] p. 80)).
Definition 1.
Let ( K , + , · ) be a ring and let ( K , + ) be its additive group. A subset I K is called a two-sided ideal (or simply an ideal) of K if it is an additive subgroup of K such that I · K : = { r · x : r I and x K } I and K · I : = { x · r : r I and x K } I .
Proposition 3.
If g h is a two-sided ideal of the Lie algebra h of a group H , then the subgroup G H generated by exp ( g ) = { exp ( G ) : G g } is normal and closed, with Lie algebra h .
From the above proposition, we deduce the following result.
Lemma 2.
Let Z M r ( R k × r ) and Z M k r ( R k × ( k r ) ) be such that Z T Z = 0 . Then g Z R k × k is a two-sided ideal of the Lie algebra R k × k and hence
G Z : = { exp ( Z X Z + ) : X R ( k r ) × r }
is a closed Lie group with Lie algebra g Z . Furthermore, the map exp : g Z G Z is bijective.
Proof. 
Consider Z X Z + g Z and A R k × k . Noting that Z + Z = i d r and ( Z ) + Z = i d k r , we have that
( Z X Z + ) A = Z ( X Z + A Z ) Z + ,
which proves that g Z · R k × k g Z . Similarly, we have that
A ( Z X Z + ) = Z ( ( Z ) + A Z X ) Z + ,
which proves that R k × k · g Z g Z . This proves that g Z is a two-sided ideal. The map exp is clearly surjective. To prove that it is injective, we assume exp ( Z X Z + ) = exp ( Z X ˜ Z + ) for X , X ˜ R ( k r ) × r . Then, from (6), we obtain Z + Z X = Z + Z X ˜ and hence X = X ˜ , i.e., Z X Z + = Z X ˜ Z + in g Z .
Finally, we can prove the following result.
Theorem 2.
The set S Z together with the group operation × Z defined by
exp ( Z X Z + ) Z × Z exp ( Z X ˜ Z + ) Z = exp ( Z ( X + X ˜ ) Z + ) Z
for X , X ˜ R ( k r ) × r is a Lie group.
Proof. 
To prove that it is a Lie group, we simply note that the multiplication and inversion maps
μ : S Z × S Z S Z , ( W , W ˜ ) exp ( Z ( Z + ( W Z ) + Z + ( W ˜ Z ) ) Z + ) Z
and
δ : S Z S Z , W exp ( Z Z + ( W Z ) Z + ) Z
are analytic. □
It follows that U Z can be identified with a Lie group through the map φ Z .
Theorem 3.
Each neighbourhood U Z of G r ( R k ) together with the group operation Z defined by
V Z V = φ Z 1 ( φ Z ( V ) + φ Z ( V ) )
for V , V U Z , is a Lie group, and the map γ Z : U Z G Z given by
γ Z ( U ) = exp ( Z φ Z ( U ) Z + )
is a Lie group isomorphism.

3. The Non-Compact Stiefel Principal Bundle M r ( R k × r )

In this section, we give a new geometric description of the set M r ( R k × r ) of matrices with full rank r < k , which is based on the geometric description of the Grassmann manifold given in Section 2.

3.1. Principal Bundle Structure of M r ( R k × r )

For Z M r ( R k × r ) , we define a neighbourhood of Z as
V Z : = { W M r ( R k × r ) : det ( Z T W ) 0 } S Z .
From Proposition 1, we know that for a given matrix W V Z , there exists a unique pair of matrices ( X , G ) R ( k r ) × r × GL r such that W = ( Z + Z X ) G . Therefore,
V Z = { ( Z + Z X ) G : X R ( k r ) × r , G GL r } .
It allows us to introduce a parametrisation ξ Z 1 (see Figure 1) defined through the bijection
ξ Z : V Z R ( k r ) × r × GL r ,
such that
ξ Z 1 ( X , G ) = ( Z + Z X ) G
for ( X , G ) R ( k r ) × r × GL r and
ξ Z ( W ) = ( Z + W ( Z + W ) 1 , Z + W )
for W V Z . In particular,
ξ Z 1 ( 0 , i d r ) = Z .
Theorem 4.
The collection B k , r : = { ( V Z , ξ Z ) : Z M r ( R k × r ) } is an analytic atlas for M r ( R k × r ) , and hence ( M r ( R k × r ) , B k , r ) is an analytic k r -dimensional manifold modelled on R ( k r ) × r × R r × r .
Proof. 
{ V Z } Z M r ( R k × r ) is clearly a covering of M r ( R k × r ) . Moreover, since ξ Z is bijective from V Z to R ( k r ) × r × GL r we claim that if V Z V Z ˜ for Z , Z ˜ M r ( R k × r ) , then the following statements hold:
(i)
ξ Z ( V Z V Z ˜ ) and ξ Z ˜ ( V Z V Z ˜ ) are open sets in R ( k r ) × r × GL r and
(ii)
the map ξ Z ˜ ξ Z 1 is analytic from ξ Z ( V Z V Z ˜ ) R ( k r ) × r × GL r to ξ Z ˜ ( V Z V Z ˜ ) R ( k r ) × r × GL r .
In this proof, we equip R k × r with the topology τ R k × r induced by matrix norms. For any Z M r ( R k × r ) , V Z = { W R k × r : det ( Z T W ) 0 } is the inverse image of the open set R { 0 } by the continuous map W det ( Z T W ) from R k × r to R , and therefore, V Z is an open set of R k × r . Since V Z and V Z ˜ are open sets in R k × r , V Z V Z ˜ is also an open set in R k × r and since ξ Z 1 is a continuous map from R ( k r ) × r × GL r to R k × r , the set ξ Z ( V Z V Z ˜ ) , as the inverse image of an open set by a continuous map, is an open set in R ( k r ) × r × GL r . Similarly, ξ Z ˜ ( V Z V Z ˜ ) is an open set. Now let ( X , G ) R ( k r ) × r × GL r such that ξ Z 1 ( X , G ) V Z V Z ˜ . From the expressions of ξ Z 1 and ξ Z ˜ , the map ξ Z ˜ ξ Z 1 is defined by
ξ Z ˜ ξ Z 1 ( X , G ) = ( Z ˜ + ξ Z 1 ( X , G ) ( Z ˜ + ξ Z 1 ( X , G ) ) 1 , Z ˜ + ξ Z 1 ( X , G ) ) ,
with ξ Z 1 ( X , G ) = ( Z + Z X ) G , which is clearly an analytic map. □
Before stating the next result, we recall the definition of a morphism between manifolds and of a fibre bundle. We introduce notions of C p maps and C p manifolds, with p N { } or p = ω . In the latter case, C ω means analytic.
Definition 2.
Let ( M , A ) and ( N , B ) be two C p manifolds. Let F : M N be a map. We say that F is a C p morphism between ( M , A ) and ( N , B ) if given m M , there exists a chart ( U , φ ) A such that m U and a chart ( W , ψ ) B such that F ( m ) W where F ( U ) W , and the map
ψ F φ 1 : φ ( U ) ψ ( W )
is a map of class C p . If it is a C p diffeomorphism, then we say that F is a C p diffeomorphism between manifolds. We say that ψ F φ 1 is a representation of F using a system of local coordinates given by the charts ( U , φ ) and ( W , ψ ) .
Definition 3.
Let B be a C p manifold with atlas A = { ( U b , φ b ) : b B } , and let F be a manifold. A C p fibre bundle E with base B and typical fibre F is a C p manifold which is locally a product manifold; that is, there exists a surjective morphism π : E B such that for each b B there is a C p diffeomorphism between manifolds
χ b : π 1 ( U b ) U b × F ,
such that p b χ b = π where p b : U b × F U b is the projection. For each b B , π 1 ( b ) = E b is called the fibre over b . The C p diffeomorphisms χ b are called fibre bundle charts. If p = 0 , E , B and F are only required to be topological spaces and { U b : b B } an open covering of B . In the case where F is a Lie group, we say that E is a C p principal bundle, and if F is a vector space, we say that it is a C p vector bundle.
Theorem 5.
The set M r ( R k × r ) is an analytic principal bundle with typical fibre GL r and base G r ( R k ) , with a surjective morphism between M r ( R k × r ) and G r ( R k ) given by the map col k , r .
Proof. 
To show that it is an analytic principal bundle, we first observe that
col k , r : ( M r ( R k × r ) , B k , r ) ( G r ( R k ) , A k , r )
is a surjective morphism. Indeed, let Z M r ( R k × r ) and ( V Z , ξ Z ) B k , r and ( U Z , φ Z ) A k , r . Noting that col k , r ( Y G ) = col k , r ( Y ) for all Y S Z , we obtain that col k , r ( V Z ) = U Z . Moreover, a representation of col k , r by using a system of local coordinates given by the charts is
( φ Z col k , r ξ Z 1 ) ( X , G ) = X ,
which is clearly an analytic map from R ( k r ) × r × GL r to R ( k r ) × r such that col k , r 1 ( U Z ) = V Z . Now, a representation of the morphism
χ Z : ( V Z , { ( V Z , ξ Z ) } ) ( U Z , { ( U Z , φ Z ) } ) × ( GL r , { ( GL r , i d R r × r ) } ) , W ( col k , r ( W ) , G )
using the system of local coordinates given by the charts is
( ( φ Z × i d R r × r ) χ Z ξ Z 1 ) : R ( k r ) × r × GL r R ( k r ) × r × GL r ,
defined by
( ( φ Z × i d R r × r ) χ Z ξ Z 1 ) ( X , G ) = ( X , G ) ,
which is clearly an analytic diffeomorphism. To conclude, consider the projection
p Z : U Z × GL r U Z , ( V , G ) V ,
and observe that ( p Z χ Z ) ( W ) = col k , r ( W ) holds for all W V Z .

3.2. M r ( R k × r ) as a Submanifold and Its Tangent Space

Here, we prove that the non-compact Stiefel manifold M r ( R k × r ) equipped with the topology given by the atlas B k , r is an embedded submanifold in R k × r . For that, we have to prove that the standard inclusion map
i : ( M r ( R k × r ) , B k , r ) ( R k × r , { ( R k × r , i d R k × r ) } )
as a morphism is an embedding. To see this, we need to recall some definitions and results.
Definition 4.
Let F : ( M , A ) ( N , B ) be a morphism between C p manifolds and let m M . We say that F is an immersion at mif there exists an open neighbourhood U m of m in M such that the restriction of F to U m induces an isomorphism from U m onto a submanifold of N . We say that F is an immersion if it is an immersion at each point of M .
The next step is to recall the definition of the differential as a morphism which gives a linear map between the tangent spaces of the manifolds (in local coordinates) involved with the morphism. Let us recall that for any m M , we denote by T m M the tangent space of M at m (in local coordinates).
Definition 5.
Let ( M , A ) and ( N , B ) be two C p manifolds. Let F : ( M , A ) ( N , B ) be a morphism of class C p ; i.e., for any m M ,
ψ F φ 1 : φ ( U ) ψ ( W )
is a map of class C p , where ( U , φ ) A is a chart in M containing m and ( W , ψ ) B is a chart in N containing F ( m ) . Then we define
T m F : T m ( M ) T F ( m ) ( N ) , υ D ( ψ F φ 1 ) ( φ ( m ) ) [ υ ] .
For finite dimensional manifolds we have the following criterion for immersions (see Theorem 3.5.7 in [28]).
Proposition 4.
Let ( M , A ) and ( N , B ) be C p manifolds. Let
F : ( M , A ) ( N , B )
be a C p morphism and m M . Then F is an immersion at m if and only if T m F is injective.
A concept related to an immersion between manifolds is given in the following definition.
Definition 6.
Let ( M , A ) and ( N , B ) be C p manifolds and let f : ( M , A ) ( N , B ) be a C p morphism. If f is an injective immersion, then f ( M ) is called an immersed submanifold of N .
Finally, we give the definition of embedding.
Definition 7.
Let ( M , A ) and ( N , B ) be C p manifolds and let f : ( M , A ) ( N , B ) be a C p morphism. If f is an injective immersion, and f : ( M , τ A ) ( f ( M ) , τ B | f ( M ) ) is a topological homeomorphism, then we say that f is an embedding and f ( M ) is called an embedded submanifold of N .
We first note that the representation of the inclusion map i using the system of local coordinates given by the charts ( V Z , ξ Z ) B k , r in M r ( R k × r ) and ( R k × r , i d R k × r ) in R k × r is
( i d R k × r i ξ Z 1 ) = ( i ξ Z 1 ) : R ( k r ) × r × GL r R k × r , ( X , G ) ( Z + Z X ) G .
Then the tangent map T Z i at Z = ξ Z 1 ( 0 , i d r ) , defined by T Z i = D ( i ξ Z 1 ) ( 0 , i d r ) , is
T Z i : R ( k r ) × r × R r × r R k × r , ( X ˙ , G ˙ ) Z X ˙ + Z G ˙ .
Proposition 5.
The tangent map T Z i : R ( k r ) × r × R r × r R k × r at Z M r ( R k × r ) is a linear isomorphism, with inverse ( T Z i ) 1 given by
( T Z i ) 1 ( Z ˙ ) = ( Z + Z ˙ , Z + Z ˙ ) ,
for Z ˙ R k × r . Furthermore, the standard inclusion map i is an embedding from M r ( R k × r ) to R k × r .
Proof. 
Let us assume that T Z i ( X ˙ , G ˙ ) = Z X ˙ + Z G ˙ = 0 . Multiplying this equality by Z + and Z + on the left, we obtain G ˙ = 0 and X ˙ = 0 , respectively, which implies that T Z i is injective. To prove that it is also surjective, we consider a matrix Z ˙ R k × r and observe that X ˙ = Z + Z ˙ R ( k r ) × r and G ˙ = Z + Z ˙ R r × r is such that T Z i ( X ˙ , G ˙ ) = Z ˙ . Since T Z i is injective, the inclusion map i is an immersion.
To prove that it is an embedding, we equip M r ( R k × r ) with the topology τ B k , r given by the atlas and we equip R k × r with the topology τ R k × r induced by matrix norms. We need to check that
i : ( M r ( R k × r ) , τ B k , r ) ( M r ( R k × r ) , τ R k × r | M r ( R k × r ) )
is a topological homeomorphism. Since the topology in ( M r ( R k × r ) , τ B k , r ) has the property that each local chart ξ Z is indeed a homeomorphism from V Z in M r ( R k × r ) to ξ Z ( V Z ) = R ( k r ) × r × GL r (see Section 1.1), we only need to show that the bijection ( i ξ Z 1 ) : R ( k r ) × r × GL r V Z R k × r given by
( i ξ Z 1 ) ( X , G ) = ( Z + Z X ) G
is a topological homeomorphism for all Z M r ( R k × r ) . Observe that D ( i ξ Z 1 ) ( X , G ) L ( R ( k r ) × r × R r × r , R k × r ) is given by
D ( i ξ Z 1 ) ( X , G ) [ ( X ˙ , G ˙ ) ] = Z X ˙ G + ( Z + Z X ) G ˙ .
Assume that Z X ˙ G + ( Z + Z X ) G ˙ = 0 . Multiplying this equality by Z + on the left we obtain G ˙ = 0 , and hence Z X ˙ G = 0 . Multiplying by Z + on the left, we obtain X ˙ G = 0 . Thus, X ˙ = 0 and as a consequence D ( i ξ Z 1 ) ( X , G ) is a linear isomorphism for each ( X , G ) R ( k r ) × r × GL r . The inverse function theorem says us that ( i ξ Z 1 ) is a diffeomorphism, in particular a homeomorphism,, and hence i is an embedding. □
The tangent space to M r ( R k × r ) at Z is the image through T Z i of the tangent space at Z in local coordinates T Z M r ( R k × r ) = R ( k r ) × r × R r × r , i.e.,
T Z M r ( R k × r ) = { Z X ˙ + Z G ˙ : X ˙ R ( k r ) × r , G ˙ R r × r } = R k × r ,
and can be decomposed into a vertical tangent space
T Z V M r ( R k × r ) = { Z G ˙ : G ˙ R r × r } ,
and a horizontal tangent space
T Z H M r ( R k × r ) = { Z X ˙ : X ˙ R ( k r ) × r } .

3.3. Lie Group Structure of Neighbourhoods V Z

We here prove that each neighbourhood V Z of M r ( R k × r ) has the structure of a Lie group. For that, we first note that V Z can be identified with S Z × GL r , with S Z given by (9). Noting that S Z can be identified with the Lie group G Z defined in (11), we then have that V Z can be identified with a product of two Lie groups G Z × GL r , which is a Lie group with the group operation Z given by
( exp ( Z X Z + ) , G ) Z ( exp ( Z X Z + ) , G ) = ( exp ( Z ( X + X ) Z + ) , G G ) ,
for X , X R ( k r ) × r and G , G GL r . This allows us to define a group operation Z over V Z defined for W = ξ Z 1 ( X , G ) and W = ξ Z 1 ( X , G ) by
W Z W = ξ Z 1 ( X + X , G G ) ,
and to state the following result.
Theorem 6.
The set V Z together with the group operation Z defined by (15) is a Lie group and the map η Z : V Z G Z × GL r given by
η Z ( ξ Z 1 ( X , G ) ) = ( exp ( Z X Z + ) , G , )
is a Lie group isomorphism.

4. The Principal Bundle M r ( R n × m ) for 0 < r < min ( m , n )

In this section, we give a geometric description of the set of matrices M r ( R n × m ) with rank r < min ( m , n ) .

4.1. M r ( R n × m ) as a Principal Bundle

For Z M r ( R n × m ) , there exists U M r ( R n × r ) , V M r ( R m × r ) , and G GL r such that
Z = U G V T ,
where the column space of Z is col n , r ( U ) and the row space of Z is col m , r ( V ) .
Let us first introduce the surjective map
ϱ r : M r ( R n × m ) G r ( R n ) × G r ( R m ) , U G V T ( col n , r ( U ) , col m , r ( V ) ) .
The set
ϱ r 1 ( col n , r ( U ) , col m , r ( V ) ) = { U H V T : H GL r }
can be identified with GL r . Let us consider U M n r ( R n × ( n r ) ) such that U T U = 0 and V M m r ( R m × ( m r ) ) such that V T V = 0 (see Remark 1 for a practical definition). Then we define a neighbourhood of U G V T in the set M r ( R n × m ) by
U Z : = ϱ r 1 ( U U × U V ) ,
where U U and U V are the neighbourhoods of col n , r ( U ) and col m , r ( V ) , respectively (see Section 2.2). Noting that U U = φ U 1 ( R ( n r ) × r ) = col n , r ( S U ) and U V = φ V 1 ( R ( m r ) × r ) = col m , r ( S V ) , where S U and S V are the affine cross sections of U and V, respectively (defined by (4)), the neighbourhood of U G V T can be written
U Z = { ( U + U X ) H ( V + V Y ) T : ( X , Y , H ) R ( n r ) × r × R ( m r ) × r × GL r } .
We can associate with U Z the parametrisation θ Z 1 given by the chart (see Figure 2)
θ Z : U Z R ( n r ) × r × R ( m r ) × r × GL r
defined by
θ Z 1 ( X , Y , H ) = ( U + U X ) H ( V + V Y ) T
for ( X , Y , H ) R ( n r ) × r × R ( m r ) × r × GL r , and
θ Z ( A ) = ( U + A ( V + ) T ( U + A ( V + ) T ) 1 , V + A T ( U + ) T ( V + A T ( U + ) T ) 1 , U + A ( V + ) T )
for A U Z . In particular, we have θ Z 1 ( 0 , 0 , G ) = Z . We point out that U Z = U Z and θ Z = θ Z for every Z = U G V T with G G .
Theorem 7.
The collection B n , m , r : = { ( U Z , θ Z ) : Z M r ( R n × m ) } is an analytic atlas for M r ( R n × m ) , and hence ( M r ( R n × m ) , B n , m , r ) is an analytic r ( n + m r ) -dimensional manifold modelled on R ( n r ) × r × R ( m r ) × r × R r × r .
Proof. 
{ U Z } Z M r ( R n × m ) is clearly a covering of M r ( R n × m ) . Moreover, since θ Z is bijective from U Z to R ( n r ) × r × R ( m r ) × r × GL r , we claim that if U Z U Z ˜ for Z = U G V T and Z ˜ = U ˜ G ˜ V ˜ T M r ( R n × m ) , then the following statements hold:
(i)
θ Z ( U Z U Z ˜ ) and θ Z ˜ ( U Z U Z ˜ ) are open sets in R ( n r ) × r × R ( m r ) × r × GL r and
(ii)
the map θ Z ˜ θ Z 1 is analytic from θ Z ( U Z U Z ˜ ) R ( n r ) × r × R ( m r ) × r × GL r to θ Z ˜ ( U Z U Z ˜ ) R ( n r ) × r × R ( m r ) × r × GL r .
In this proof, we equip R n × m with the topology τ R n × m induced by matrix norms. We first observe that the set U Z = { A M r ( R n × m ) : det ( U T A V ) 0 } = O Z M r ( R n × m ) , where O Z = { A R n × m : det ( U T A V ) 0 } , as the inverse image of the open set R { 0 } through the continuous map A det ( U T A V ) from R n × m to R , is an open set in R n × m . In the same way, we have that U Z ˜ = O Z ˜ M r ( R n × m ) , with U Z ˜ as an open set in R n × m . Since U Z U Z ˜ = O Z O Z ˜ M r ( R n × m ) , and since the image of θ Z 1 is in M r ( R n × m ) , we have
θ Z ( U Z U Z ˜ ) = ( θ Z 1 ) 1 ( U Z U Z ˜ ) = ( θ Z 1 ) 1 ( O Z O Z ˜ ) ,
the inverse image through θ Z 1 of the open set O Z O Z ˜ in R n × m . Since θ Z 1 is a continuous map from R ( n r ) × r × R ( m r ) × r × GL r to R n × m , we deduce that θ Z ( U Z U Z ˜ ) is an open set in R ( n r ) × r × R ( m r ) × r × GL r . Similarly, θ Z ˜ ( U Z U Z ˜ ) is an open set in R ( n r ) × r × R ( m r ) × r × GL r . Now, let ( X , Y , H ) R ( n r ) × r × R ( m r ) × r × GL r such that θ Z 1 ( X , Y , H ) U Z U Z ˜ . From the expressions of θ Z 1 and θ Z ˜ , the map θ Z ˜ θ Z 1 is defined by
θ Z ˜ θ Z 1 ( X , Y , H ) = ( U ˜ + θ Z 1 ( X , Y , H ) ( V ˜ + ) T ( U ˜ + θ Z 1 ( X , Y , H ) ( V ˜ + ) T ) 1 , V ˜ + θ Z 1 ( X , Y , H ) T ( U ˜ + ) T ( V ˜ + θ Z 1 ( X , Y , H ) T ( U ˜ + ) T ) 1 , U ˜ + θ Z 1 ( X , Y , H ) ( V ˜ + ) T ) ,
with θ Z 1 ( X , Y , H ) = ( U + U X ) H ( V + V Y ) T , which is clearly an analytic map. □
Theorem 8.
The set M r ( R n × m ) is an analytic principal bundle with typical fibre GL r and base G r ( R n ) × G r ( R m ) with surjective morphism ϱ r between M r ( R n × m ) and G r ( R n ) × G r ( R m ) given by ϱ r .
Proof. 
To prove that it is an analytic principal bundle, we consider the surjective map
ϱ r : M r ( R n × m ) G r ( R n ) × G r ( R m ) , U G V T ( col n , r ( U ) , col m , r ( V ) ) ,
the atlas A n , r : = { ( U U , φ U ) : U M r ( R n × r ) } of G r ( R n ) and the atlas A m , r : = { ( U V , φ V ) : V M r ( R m × r ) } of G r ( R m ) . Recall that
U Z = { col k , r ( Z + Z X ) : X R ( k r ) × r } ,
with k = n if Z = U or k = m if Z = V , and hence
ϱ r 1 ( U U , U V ) = ( U + U X ) H ( V + V Y ) T : X R ( n r ) × r , Y R ( m r ) × r , H GL r .
Observe that for each fixed G GL r , we have that ϱ r 1 ( U U , U V ) = U Z , where Z = U G V T . Since U Z = U Z holds for Z = U G V T , where G GL r , the map
χ Z : U Z U U × U V × GL r
defined by
χ Z ( U H ( V ) T ) : = ( col n , r ( U ) , col m , r ( V ) , H ) ,
is independent of the choice of Z = U G V T , where G GL r . Now, the representation of χ Z in local coordinates is the map
( ( φ U × φ V × i d R r × r ) χ Z θ Z 1 ) : R ( n r ) × r × R ( m r ) × r × GL r R ( n r ) × r × R ( m r ) × r × GL r
given by ( ( φ U × φ V × i d R r × r ) χ Z θ Z 1 ) ( X , Y , H ) = ( X , Y , H ) , which is an analytic diffeomorphism. Moreover, let p Z : U U × U V × GL r U U × U V be the projection over the first two components. Then
( p Z χ Z ) ( U H V T ) = ( col n , r ( U ) , col m , r ( V ) ) = ϱ r ( U H V T )
and the theorem follows. □

4.2. M r ( R n × m ) as a Submanifold and Its Tangent Space

Here, we prove that M r ( R n × m ) equipped with the topology given by the atlas B n , m , r is an embedded submanifold in R n × m . For that, we have to prove that the standard inclusion map i : M r ( R n × m ) R n × m is an embedding. Noting that the inclusion map restricted to the neighbourhood U Z of Z = U G V T is identified with
( i θ Z 1 ) : R ( n r ) × r × R ( m r ) × r × GL r R n × m , ( X , Y , H ) ( U + U X ) H ( V + V Y ) T ,
the tangent map T Z i at Z = θ Z 1 ( 0 , 0 , G ) , defined by T Z i = D ( i θ Z 1 ) ( 0 , 0 , G ) , is
T Z i : R ( n r ) × r × R ( m r ) × r × R r × r R n × m , ( X ˙ , Y ˙ , H ˙ ) U X ˙ G V T + U G ( V Y ˙ ) T + U H ˙ V T .
Proposition 6.
The tangent map T Z i : R ( n r ) × r × R ( m r ) × r × R r × r R n × m at Z = U G V T M r ( R n × m ) is a linear isomorphism with inverse ( T Z i ) 1 given by
( T Z i ) 1 ( Z ˙ ) = ( U + Z ˙ ( V + ) T G 1 , V + Z ˙ T ( U + ) T G T , U + Z ˙ ( V + ) T )
for Z ˙ R n × m . Furthermore, the standard inclusion map i is an embedding from M r ( R n × m ) to R n × m .
Proof. 
Let us suppose that T Z i ( X ˙ , Y ˙ , H ˙ ) = 0 . Multiplying this equality by ( U ) + and U + on the left leads to
X ˙ G V T = 0 and G ( V Y ˙ ) T + H ˙ V T = 0 ,
respectively. By multiplying the first equation by ( V + ) T on the right, we obtain X ˙ = 0 . By multiplying the second equation on the right by ( V + ) T and ( V + ) T , we respectively obtain H ˙ = 0 and Y ˙ = 0 . Then, T Z i is injective and then i is an immersion. For Z ˙ R n × m , we note that X ˙ = U + Z ˙ ( V + ) T G 1 R n × r , Y ˙ = V + Z ˙ T ( U + ) T G T R m × r , and G ˙ = U + Z ˙ ( V + ) T R r × r is such that T Z i ( X ˙ , Y ˙ , G ˙ ) = Z ˙ , and T Z i is also surjective. Let us now equip M r ( R n × m ) with the topology τ B n , m , r given by the atlas and R n × m with the topology τ R n × m induced by matrix norms. We have to prove that
i : ( M r ( R n × m ) , τ B n , m , r ) ( M r ( R n × m ) , τ R n × m | M r ( R n × m ) )
is a topological isomorphism. The topology in ( M r ( R n × m ) , τ B n , m , r ) is such that a local chart θ Z is a homeomorphism from U Z M r ( R n × m ) to θ Z ( U Z ) = R ( n r ) × r × R ( m r ) × r × GL r (see Section 1.1). Then, to prove that the map i is an embedding, we need to show that the bijection
( i θ Z 1 ) : R ( n r ) × r × R ( m r ) × r × GL r U Z R n × m
is a topological homeomorphism. For that, observe that its differential
D ( i θ Z 1 ) ( X , Y , H ) L ( R ( n r ) × r × R ( m r ) × r × R r × r , R n × m )
at ( X , Y , H ) R ( n r ) × r × R ( m r ) × r × GL r is given by
D ( i θ Z 1 ) ( X , Y , H ) [ ( X ˙ , Y ˙ , H ˙ ) ] = ( U X ˙ ) H ( V + V Y ) T + ( U + U X ) H ( V Y ˙ ) T + ( U + U X ) H ˙ ( V + V Y ) T .
Assume that
( U X ˙ ) H ( V + V Y ) T + ( U + U X ) H ( V Y ˙ ) T + ( U + U X ) H ˙ ( V + V Y ) T = 0 .
Multiplying on the left by U + and on the right by ( V + ) T , we obtain H ˙ = 0 . Multiplying on the left by U + and on the right by ( V + ) T we deduce that X ˙ H = 0 , that is, X ˙ = 0 . Finally, multiplying on the left by U + and on the right by ( V + ) T , we obtain H Y ˙ T = 0 , and hence Y ˙ = 0 . Thus, D ( i θ Z 1 ) ( X , Y , H ) is a linear isomorphism from R ( n r ) × r × R ( m r ) × r × R r × r to D ( i θ Z 1 ) ( X , Y , H ) [ R ( n r ) × r × R ( m r ) × r × R r × r ] for each ( X , Y , H ) R ( n r ) × r × R ( m r ) × r × GL r . The inverse function theorem tells us that ( i θ Z 1 ) is a diffeomorphism from R ( n r ) × r × R ( m r ) × r × GL r to U Z = ( i θ Z 1 ) ( R ( n r ) × r × R ( m r ) × r × GL r ) and, in particular, that it is a topological homeomorphism. In consequence, the map i is an embedding. □
The tangent space to M r ( R n × m ) at Z = U G V T , which is the image through T Z i of the tangent space in local coordinates T Z M r ( R n × m ) = R ( n r ) × r × R ( m r ) × r × R r × r , is
T Z M r ( R n × m ) = { U X ˙ G V T + U G ( V Y ˙ ) T + U G ˙ V T : X ˙ R ( n r ) × r , Y ˙ R ( m r ) × r , G ˙ R r × r } ,
and can be decomposed into a vertical tangent space
T Z V M r ( R n × m ) = { U G ˙ V T : G ˙ R r × r } ,
and a horizontal tangent space
T Z H M r ( R n × m ) = { U X ˙ G V T + U G ( V Y ˙ ) T : X ˙ R ( n r ) × r , Y ˙ R ( m r ) × r } .

4.3. Lie Group Structure of Neighbourhoods U Z

We here prove that M r ( R n × m ) locally has the structure of a Lie group by proving that the neighbourhoods U Z can be identified with Lie groups.
Let Z = U G V T M r ( R n × m ) . We first note that U Z can be identified with S U × S V × GL r , with S U and S V defined by (9). Noting that S U and S V can be identified with Lie groups G U and G V defined in (11), we then have that U Z can be identified with a product of three Lie groups, which is a Lie group with the group operation Z given by
( exp ( U X U + ) , exp ( V Y V + ) , G ) Z ( exp ( U X U + ) , exp ( V Y V + ) , G ) = ( exp ( U ( X + X ) U + ) , exp ( V ( Y + Y ) V + ) , G G ) .
This allows us to define a group operation Z over U Z defined for W = θ Z 1 ( X , Y , G ) and W = θ Z 1 ( X , Y , G ) by
W Z W = θ Z 1 ( X + X , Y + Y , G G ) ,
and to state the following result.
Theorem 9.
Let Z = U G V T M r ( R n × m ) . Then the set U Z together with the group operation Z defined by (17) is a Lie group with identity element U V T , and the map η Z : U Z G U × G V × GL r given by
η Z ( θ Z 1 ( X , Y , H ) ) = ( exp ( U X U + ) , exp ( V Y V + ) , H )
is a Lie group isomorphism.

Author Contributions

M.B.-F., A.F. and A.N. equally contributed. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the RTI2018-093521-B-C32 grant from the Ministerio de Ciencia, Innovación y Universidades and by the grant number INDI20/13 from Universidad CEU Cardenal Herrera.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Zhou, X.; Yang, C.; Zhao, H.; Yu, W. Low-Rank Modeling and Its Applications in Image Analysis. ACM Comput. Surv. 2014, 47, 1–33. [Google Scholar] [CrossRef] [Green Version]
  2. Antoulas, A.; Sorensen, D.; Gugercin, S. A survey of model reduction methods for large-scale systems. Contemp. Math. 2001, 280, 193–220. [Google Scholar]
  3. Benner, P.; Gugercin, S.; Willcox, K. A survey of projection-based model reduction methods for parametric dynamical systems. SIAM Rev. 2015, 57, 483–531. [Google Scholar] [CrossRef]
  4. Nouy, A. Low-rank tensor methods for model order reduction. In Handbook of Uncertainty Quantification; Ghanem, R., Higdon, D., Owhadi, H., Eds.; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar]
  5. Benner, P.; Cohen, A.; Ohlberger, M.; Willcox, K. (Eds.) Model Reduction and Approximation: Theory and Algorithms; SIAM: Philadelphia, PA, USA, 2017. [Google Scholar]
  6. Kushner, A.; Lychagin, V.; Rubtsov, V. Contact Geometry and Non-Linear Differential Equations; Encyclopedia of Mathematics and Its Applications 101; Cambridge University Press: Cambridge, UK, 2007. [Google Scholar]
  7. Goursat, E. Sur les équations du second ordre à n variables analogues à l’équation de Monge-Ampère. Bull. Soc. Math. Fr. 1899, 27, 1–34. [Google Scholar] [CrossRef] [Green Version]
  8. Sylvester, J. On the dimension of spaces of linear transformations satisfying rank conditions. Linear Algebra Its Appl. 1986, 78, 1–10. [Google Scholar] [CrossRef] [Green Version]
  9. Eisenbud, D.; Harris, J. Vector spaces of matrices of low rank. Adv. Math. 1988, 70, 135–155. [Google Scholar] [CrossRef] [Green Version]
  10. Westwick, R. Spaces of matrices of fixed rank. Linear Multilinear Algebra 1987, 20, 171–174. [Google Scholar] [CrossRef]
  11. Westwick, R. Spaces of matrices of fixed rank II. Linear Algebra Its Appl. 1996, 235, 163–169. [Google Scholar] [CrossRef] [Green Version]
  12. Ellia, P.; Menegatti, P. Spaces of matrices of constant rank and uniform vector bundles. Linear Algebra Its Appl. 2016, 507, 474–485. [Google Scholar] [CrossRef] [Green Version]
  13. Smith, S.T. Optimization techniques on Riemannian manifolds. Fields Inst. Commun. 1994, 3, 113–135. [Google Scholar]
  14. Absil, P.-A.; Mahony, R.; Sepulchre, R. Optimization Algorithms on Matrix Manifolds; Princeton University Press: Princeton, NJ, USA, 2008. [Google Scholar]
  15. Vandereycken, B. Low-Rank Matrix Complet. Riemannian Optim. SIAM J. Optim. 2013, 23, 1214–1236. [Google Scholar] [CrossRef] [Green Version]
  16. Mishra, B.; Meyer, G.; Bach, F.; Sepulchre, R. Low-rank optimization with trace norm penalty. SIAM J. Optim. 2013, 23, 2124–2149. [Google Scholar] [CrossRef] [Green Version]
  17. Mishra, B.; Meyer, G.; Bonnabel, S.; Sepulchre, R. Fixed-rank matrix factorizations and Riemannian low-rank optimization. Comput. Stat. 2014, 29, 591–621. [Google Scholar] [CrossRef] [Green Version]
  18. Koch, O.; Lubich, C. Dynamical low-rank approximation. SIAM J. Matrix Anal. Appl. 2007, 29, 434–454. [Google Scholar] [CrossRef]
  19. Manton, J.H. A framework for generalising the Newton method and other iterative methods from Euclidean Space to manifolds. Numer. Math. 2015, 129, 91–125. [Google Scholar] [CrossRef] [Green Version]
  20. Billaud-Friess, M.; Falcó, A.; Nouy, A. A New Splitting Algorithm for Dynamical Low-Rank Approximation Motivated by the Fibre Bundle Structure of Matrix Manifolds. BIT Numerical Mathematics, Accepted. Available online: https://arxiv.org/pdf/2001.08599 (accessed on 11 July 2021).
  21. Michor, P.W. Gauge theory for the diffeomorphism group. In Proceedings of the Conference Differential Geometric Methods in Theoretical Physics, Como, Italy, 24–29 August 1987; Bleuler, K., Werner, M., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1987; pp. 345–371. [Google Scholar]
  22. Falcó, A.; Hackbusch, A.W.; Nouy, A. On the Dirac-Frenkel variational principle on tensor Banach spaces. Found. Comput. Math. 2019, 19, 159–204. [Google Scholar] [CrossRef] [Green Version]
  23. Lang, S. Differential and Riemannian Manifolds, 3rd ed.; Graduate Texts in Mathematics; Springer: New York, NY, USA, 1995; Volume 160. [Google Scholar]
  24. Absil, P.-A.; Mahony, R.; Sepulchre, R. Riemannian Geometry of Grassmann Manifolds with a View on Algorithmic Computation. Acta Appl. Math. 2004, 80, 199–220. [Google Scholar] [CrossRef]
  25. Hoog, F.R.D.; Mattheij, R.M.M. Subset selection for matrices. Linear Algebra Its Appl. 2007, 422, 349–359. [Google Scholar] [CrossRef] [Green Version]
  26. Goreinov, S.A.; Oseledets, I.V.; Savostyanov, D.V.; Tyrtyshnikov, E.E.; Zamarashkin, N.L. How to find a good submatrix. In Matrix Methods: Theory, Algorithms, Applications; World Scientific: Hackensack, NY, USA, 2010; pp. 247–256. [Google Scholar]
  27. Procesi, C. Lie Groups: An Approach through Invariants and Representations, 1st ed.; Springer: New York, NY, USA, 2007. [Google Scholar]
  28. Abraham, R.; Marsden, J.E.; Ratiu, T. Manifolds, Tensor Analysis, and Applications, 2nd ed.; Springer: New York, NY, USA, 1988; Volume 75. [Google Scholar]
Figure 1. Illustration of the chart ξ Z which associates with W = ( Z + Z X ) G V Z M r ( R k × r ) , the parameters ( X , G ) in R ( k r ) × r × GL r .
Figure 1. Illustration of the chart ξ Z which associates with W = ( Z + Z X ) G V Z M r ( R k × r ) , the parameters ( X , G ) in R ( k r ) × r × GL r .
Mathematics 09 01669 g001
Figure 2. Illustration of the chart θ Z which associates with W = ( U + U X ) H ( V + V Y ) T U Z M r ( R n × m ) , the parameters ( X , Y , G ) in R ( n r ) × r × R ( m r ) × r × GL r .
Figure 2. Illustration of the chart θ Z which associates with W = ( U + U X ) H ( V + V Y ) T U Z M r ( R n × m ) , the parameters ( X , Y , G ) in R ( n r ) × r × R ( m r ) × r × GL r .
Mathematics 09 01669 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Billaud-Friess, M.; Falcó, A.; Nouy, A. Principal Bundle Structure of Matrix Manifolds. Mathematics 2021, 9, 1669. https://doi.org/10.3390/math9141669

AMA Style

Billaud-Friess M, Falcó A, Nouy A. Principal Bundle Structure of Matrix Manifolds. Mathematics. 2021; 9(14):1669. https://doi.org/10.3390/math9141669

Chicago/Turabian Style

Billaud-Friess, Marie, Antonio Falcó, and Anthony Nouy. 2021. "Principal Bundle Structure of Matrix Manifolds" Mathematics 9, no. 14: 1669. https://doi.org/10.3390/math9141669

APA Style

Billaud-Friess, M., Falcó, A., & Nouy, A. (2021). Principal Bundle Structure of Matrix Manifolds. Mathematics, 9(14), 1669. https://doi.org/10.3390/math9141669

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop