arXiv:1612.07142v1 [math.AT] 21 Dec 2016 and agn spaces. tangent where ai ftepicplfibration principal the of basis stesto matrices of set the is nti aew have we case this In with T oecnetosadnttosi s nti ae n state and paper this in use in notations and conventions some osrcinwsgnrlzdb .Goe-aoadtefirst the G´omez-Tato and A. by A generalized was construction that h agn space tangent the h quaternions the orsod oteotooa ru O( group Sp( orthogonal group the to corresponds hsmpstse h equality the satisfies map This sdfie by defined is aeoy Optimisation. category; A ∈ nti work, this In e 0 Let eoeby Denote Let 2010 Date e od n phrases. and words Key G( G( x AA n NIU MAC ENRIQUE I ac 5 2018. 25, March : ι = ) ahmtc ujc Classification. Subject Mathematics = n ( utri-cnrlanctgr n piiainproble optimisation and category Lusternik-Schnirelmann Abstract. ∈ B ∗ samap a as ) ≤ n ρ = = ) K { ( ). k ALYTASOMO TEE MANIFOLDS STIEFEL ON TRANSFORM CAYLEY X A n I × n ≤ c ,w eoeby denote we ), K  ∈ I n where , ( n eteiett arx h lsia alytransform, Cayley classical The matrix. identity the be X ecntutCye rnfrso tee manifolds Stiefel on transforms Cayley construct we H h ler fete h elnumbers real the either of algebra the K B 0 h opc tee aiodO manifold Stiefel compact The . edfieaCye rnfr nSiflmnfls Applicati manifolds. Stiefel on transform Cayley a define We e G( Let . n ( = ) T × ´ c IAS-VIRG c I c A n I A − 0 G( x A k Ω( : 1 | A Ω(  ( I ∈ utrincSiflmnfl;Cye rnfr;Lustern transform; Cayley manifold; Stiefel Quaternionic A X = n ∗ − A and ∗ = ) K = ( = ) A X c n = ) A X n ) A O( = ) × S MAR OS, ∗ G( ¯ ´ + → )( ρ t k n hr sadffoopimbtentetnetspace tangent the between diffeomorphism a is there and { I ρ  c X I n stecnuaetasoe eedn on Depending transpose. conjugate the is X stepoeto notelast the onto projection the is uhthat such ∗ I 2 X Ω( − A + ∗ − d Moreover id. = A : A A ∈ ∈ n, X Introduction k T ∗ ∗ 0 = ´ ) AJOS IA K X rmr 38;Scnay5C9 55M30. 58C99; Secondary 43A85; Primary A ) endby defined ) K K − −→ G( n )( eteLegopo matrices of group Lie the be ) n 1 × ι } n × A n n ihΩ( with ,teuiaygopU( group unitary the ), n n Ω( and x G( 1 ) | PEREIRA-S E ∗ + ´ ( x | → A n X X = ) + ) T −→ − I + A I x ρ X c k 1 = ) O ∗ I hsmnfl per loa the as also appears manifold This . X = ) ) n,k n,k O nue iemrhs between diffeomorphism a induces − ∩ ∗ 1 c n,k  I G( forthonormal of saepresented. are ms exists 0 = E,ADDNE TANR DANIEL AND AEZ, ( h a nue ewe the between induced map the ´ X R A , n h ope numbers complex the , ∗ ∈ ). X } K

) u anresults. main our n Ω( and A k n . w uhr 4 oany to [4] authors two × ∗ oun.If columns. n , n | rtesymplectic the or ) ( efis specify first We . I I c n othe to ons ) + A I k ∩ Ω( : fae in -frames X ik-Schnirelmann ∈ K G( ) K − I hsgroup this A n 1 ) n .This ). × → exists ∈ E ´ n G( C Ω( such K (1) n I or

n ) ) . 2 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´

The next statement contains the existence and the main properties of a Cayley trans- form in Stiefel manifolds. T Theorem A. Let 0 ≤ k ≤ n and x = ∈ O with P ∈ Kk×k. We choose  P  n,k α T A = ∈ G(n). Then there exists a map  β P  A γ : TxOn,k → On,k A such that γ ◦ ρ∗A = ρ ◦ cA. Moreover we have the following properties. 1) The map γA is injective on the open subset

x X −1 Γ = v = A ∈ TxO | (βX + P ) exists .   Y  n,k  This subset Γx does not depend on the choice of A such that ρ(A)= x. Furthermore, A x if γ is injective on an open subset U ⊂ TxOn,k then we have U ⊂ Γ . A x 2) The map γ induces a diffeomorphism between Γ ⊂ TxOn,k and the open subset τ Ωx = ∈ O | (π + P ∗)−1 exists .  π n,k  An explicit formula for γA is given in Definition 1.3. Also, the expression of the inverse A −1 map (γ |Γx ) appears in Equations (12) and (13). ∗ As we said before for the group G(n), the Cayley transform cA : TAG(n) → Ω(A ) ∩ G(n) is a diffeomorphism. Therefore the Cayley open subset Ω(A∗)∩G(n) is contractible. This property cannot be extended as it stands in the case of a Stiefel manifold. However the image of the injectivity domain of a Cayley transform in On,k is contractible in On,k. x Theorem B. For every x ∈ On,k the open subset Ω is contractible in On,k. This property is a consequence of the existence of a local section (see Proposition 3.3) A x s :Ω → G(n) of the projection ρ: G(n) → On,k and the contractibility of the Cayley open subsets Ω(A∗) ∩ G(n). The contents of the paper are as follows. Section 1 contains the construction of the A Cayley transform γ : TxOn,k → On,k. The study of the injectivity of its derivative is done in Section 2 and the proofs of Theorems A and B occupy Section 3. Finally, Section 4 is devoted to applications of this construction to Lusternik-Schnirelmann cat- egory of the quaternionic Stiefel manifolds and to optimisation problems on real Stiefel manifolds.

1. Construction Let Kn be either the real vector space Rn, the complex vector space Cn or the quater- nionic vector space Hn (with the structure of a right H-vector space) endowed with the ∗ inner product hu, vi = u v. Let 0 ≤ k ≤ n. The compact Stiefel manifold On,k of n n×k ∗ orthonormal k-frames in K is the set of matrices x ∈ K such that x x = Ik. It is standard to denote On,k by Vn,k in the real case, Wn,k in the complex case and Xn,k in the quaternionic case. CAYLEYTRANSFORMONSTIEFELMANIFOLDS 3

T Usually we write x = ∈ O , with T ∈ K(n−k)×k and P ∈ Kk×k. The linear left  P  n,k 0 action of G(n) on On,k is transitive and the isotropy group of x0 = is isomorphic  Ik to G(n−k). Therefore On,k is diffeomorphic to G(n)/G(n−k) and we have the principal ι ρ fibration G(n − k) −→ G(n) −→ On,k. Let x ∈ On,k. We complete x to a matrix A ∈ G(n) such that ρ(A) = x. The tangent space TI G(n) of the group G(n) at the identity is the set of skew-Hermitian (skew-symmetric in the real case) matrices and the tangent space TAG(n) at A equals A · TI G(n). On the other side, recall that Kn×k ∗ ∗ TxOn,k = {v ∈ | v x + x v = 0} = A · Tx0 On,k, 0 where x0 = = ρ(In). So each tangent vector v ∈ TxOn,k can be written as  Ik X v = A , with X ∈ K(n−k)×k, Y ∈ Kk×k and Y + Y ∗ = 0. Grants to the principal  Y  fibration defining the Stiefel manifold, the tangent space TxOn,k can be identified to the ⊥ orthogonal (TAG(n − k)) of the image of the inclusion of G(n − k) in G(n),

⊥ 0 X ∗ (TAG(n − k)) = A | Y + Y = 0 =∼ TxO .   −X∗ Y   n,k

X ∗ With this identification, the tangent space TxO = A | Y + Y = 0 is consid- n,k   Y   ered as a subspace of TAG(n) and we may apply the Cayley map cA of G(n) on it. From Equation (1) we have

cA = RA∗ ◦ cI ◦ LA∗ , (2) where LA∗ and RA∗ denote as usual the left and right multiplications in a Lie group. ⊥ Thus, we have first to determine cI on the elements of (TI G(n − k)) . Lemma 1.1. Let X ∈ K(n−k)×k and Y ∈ Kk×k such that Y + Y ∗ = 0. Then the matrix ∗ Ik + X X + Y is invertible. 0 X Proof. The skew-symmetric matrix M = cannot have real eigenvalues,  −X∗ Y  In−k X then In + M = ∗ is invertible. In the following product,  −X Ik + Y 

In−k 0 In−k X ∗ · (In + M)= ∗ , (3)  X Ik  0 X X + Ik + Y  the two factors on the left-hand side admit an inverse. So, the matrix on the right-hand side admits an inverse, and the results follows. 

We denote ∗ −1 b = (Ik + X X + Y ) . (4) 4 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´

0 X Proposition 1.2. Let M = ∈ TI G(n), then we have  −X∗ Y  ∗ In−k − 2XbX −2Xb cI (M)= ∗ . (5)  2bX −Ik + 2b Proof. From Equation (3) we have −1 −1 In−k X In−k 0 In−k −Xb In−k 0 (In + M) = −1 ∗ = ∗ .  0 b   X Ik  0 b  X Ik

By applying the definition of cI : TI G(n) → G(n), we get: −1 cI (M) = (In − M) (In + M) ∗ In−k −X −1 In−k − 2XbX −2Xb  = ∗ (In + M) = ∗ .  X Ik − Y   2bX −Ik + 2b A computation from (2) and (5) gives directly:

∗ ∗ ∗ ∗ ∗ ∗ (In−k − 2XbX )α − 2XbT (In−k − 2XbX )β − 2XbP cA (AM)= ∗ ∗ ∗ ∗ ∗ ∗ .  2bX α + (−Ik + 2b)T 2bX β + (−Ik + 2b)P 

The Cayley transform is now obtained by projecting this expression on On,k. T α T Definition 1.3. Let x = ∈ O . We choose A = ∈ G(n), and consider  P  n,k  β P  X A v = A ∈ TxO . The Cayley transform on the Stiefel manifold, γ : TxO →  Y  n,k n,k On,k, is defined by ∗ ∗ ∗ ∗ A (In−k − 2XbX )β − 2XbP −Xb ∗ β γ (v)= ∗ ∗ ∗ = 2 (βX + P ) + ∗ , (6)  2bX β + (−Ik + 2b)P   b  −P  where b is given in Equation (4). Remark 1.4. The map γA depends on the choice of A such that ρ(A) = x. With the previous notation, the elements of G(n) that are sent on x are the matrices A•E := E 0 αE T X E∗X A = with E ∈ G(n−k). We observe v = A = (A•E)  0 Ik  βE P   Y   Y  ∗ ∗ ∗ ∗ ∗ and Ik + (E X) (E X)+ Y = Ik + X X + Y . Thus, in (6), if we replace X by E X, β by βE and keep unchanged b and Y , we get E∗ 0 γA•E(v)= γA(v).  0 Ik

We end this section by noticing that the behavior of γA is different from that of the T Cayley transform cA in G(n). For instance, when n − k ≥ k, if we choose x =  0 0 β∗ and v = A , we have γA(v) = , which does not depend on Y . Thus γA is  Y   0  not injective on the tangent space TxOn,k. We address the determination of a domain of injectivity for γA in Section 3 but, before that, we study the differential of γA. CAYLEYTRANSFORMONSTIEFELMANIFOLDS 5

2. Differential The results of this section are used in the study of the domain of injectivity of the Cayley transform γA. T α T X Let x = ∈ O and A = ∈ G(n), as before. Let v = ∈  P  n,k  β P  0  Y  A A Kn×k Tx0 On,k. The differential of γ as a map γ : TxOn,k → at the point v = Av0 ∈ TxOn,k is denoted A n×k (γ )∗v : TvTxOn,k =∼ TxOn,k → K .

A M ∗ We compute (γ ) v(w) for any w = A ∈ TxO , that is, N + N = 0. Since, with ∗  N  n,k ⊥ the identification TxOn,k =∼ (TAG(n − k)) , we have A γ = ρ ◦ RA∗ ◦ cI ◦ LA∗ , (7) A the differential of γ is determined by that of cI . Therefore, we first consider M d X M (cI )∗v0 = cI + t .  N  dt |t=0 Y  N  Let −1 ∗ bt := Ik + (X + tM) (X + tM)+ Y + tN. d −1 Its derivative dt |t=0 bt is denoted ξ and equals ξ := X∗M + M ∗X + N. (8) −1 −1 −1 ′ Moreover b0 = b . From btbt = I, we deduce b0 = −bξb. Then, a direct computation from Equation (5) gives M −2MbX∗ + 2XbξbX∗ − 2XbM ∗ −2Mb + 2Xbξb (cI ) v0 = . (9) ∗  N   −2bξbX∗ + 2bM ∗ −2bξb  Proposition 2.1. With the previous notations, the differential of the Cayley transform of the Stiefel manifold On,k is given by ∗ ∗ ∗ ∗ ∗ A (−2MbX + 2XbξbX − 2XbM )β + (−2Mb + 2Xbξb)P (γ ) v(w)= . (10) ∗  (−2bξbX∗ + 2bM ∗)β∗ − 2bξbP ∗  A ∗ Proof. The equality (7) gives by the chain rule (γ )∗v(w) = ρ ((cI )∗v0 (w0) · A ), where ∗ ∗ M ∗ α β w = and A = , because the projection ρ and the translations RA∗ 0  N   T ∗ P ∗ and LA∗ are linear maps. Then formula (9) gives the value (10).  A Proposition 2.2. The differential (γ )∗v is injective if and only if the matrix βX + P is invertible.

A M Proof. According to (10), the kernel of (γ ) v is the space of solutions of the ∗  N  system (i) 2Mb(βX + P )∗ = 2Xbξb(βX + P )∗ − 2XbM ∗β∗,  (ii) ξb(βX + P )∗ = M ∗β∗, 6 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´ where we have used that the matrix b is invertible. Then we get Mb(X∗β∗ + P ∗) = 0, so the first system is equivalent to (iii) Mb(βX + P )∗ = 0,  (iv) ξb(βX + P )∗ = (βM)∗. • If we suppose the matrix βX + P invertible, then the equation (iii) gives M = 0 and the equation (iv) gives ξ = 0. Finally, from the definition of ξ in Equation (8) we have M = N = 0. • Conversely, we suppose the kernel of βX + P not reduced to 0 and we look for an A element in the kernel of (γ )∗v of the particular type M = 0. In this case, the equation (iii) is trivially satisfied and the equation (iv) may be reduced to Nb(βX + P )∗ = 0. We consider the singular value decomposition of b(βX + P )∗ ∈ Kk×k (for the quaternionic case see [13]): Z 0 b(βX + P )∗ = µ 1 ν∗, µ,ν ∈ G(k),  0 Z2 r×r where Z1 is diagonal without zero value on it and Z2 = 0 ∈ K . As βX + P is not invertible, we have r > 0. The existence of solutions in the equation Nb(βX + P )∗ = 0 is then equivalent to the existence of solutions in Z 0 Z 0 Nµ 1 = 0 ⇐⇒ µ∗Nµ 1 = 0.  0 Z2  0 Z2 The fact that r > 0 allows the choice of a non-zero, skew-symmetric matrix µ∗Nµ A 0 satisfying the last equation. Thus N 6= 0 is skew-symmetric and (γ ) v = 0.  ∗  N

3. Properties This section consists of the proof of Theorems A and B. Recall the notations x = T α T ∈ O with T ∈ K(n−k)×k, P ∈ Kk×k and the choice of A = ∈ G(n).  P  n,k  β P  Proof of Theorem A. 1) First, we look at the independence of Γx on the choice of A. With the notations of αE T Remark 1.4, any matrix projecting on x can be written as A•E = with E ∈  βE P  X E∗X G(n − k). An element v ∈ TxO may be expressed as v = A = (A•E) . n,k  Y   Y  The fact that Γx does not depend on the choice of A comes from (βE)(E∗X)+ P = βX + P . X1 X2 x As for the injectivity, let v1 = A and v2 = A be two vectors of Γ . From  Y1   Y2  A A Definition 1.3, the equality γ (v1)= γ (v2) is equivalent to the system ∗ ∗ (i) X1b1(βX1 + P ) = X2b2(βX2 + P ) , ∗ ∗  (ii) b1(βX1 + P ) = b2(βX2 + P ) , CAYLEYTRANSFORMONSTIEFELMANIFOLDS 7 from which we deduce ∗ (X1 − X2)b1(βX1 + P ) = 0. (11) x As v1 ∈ Γ means that βX1 + P is invertible, we get from (11) the equality X1 = X2. Then the equation (ii) implies b1 = b2, from which and (4) we deduce that Y1 = Y2 and the injectivity of γA on Γx. Conversely, suppose γA invertible on an open subset U. This implies the injectivity A x of the differential (γ )∗v for any v ∈ U and Proposition 2.1 gives the inclusion U ⊂ Γ . τ τ 2) The values of τ and π such that γA(v) = are given by (6). Let such  π  π that π + P ∗ is invertible. We are looking for matrices X ∈ K(n−k)×k and Y ∈ Kk×k, with Y skew-symmetric, such that βX + P is invertible and the following system, which is equivalent to Definition 1.3, is satisfied: (a) τ − β∗ = −2Xb(βX + P )∗,  (b) π + P ∗ = 2b(βX + P )∗. In particular, we have from (b) that the matrix π +P ∗ is invertible if and only if βX +P is so. From (a) we get the value of X, X = −(τ − β∗)(π + P ∗)−1. (12) Also from (b) we obtain 1 b = (π + P ∗) [(βX + P )∗]−1 2 and the expression of Y follows from the fact that Y is the skew-symmetric part of b−1, that is, 2Y = b−1 − (b−1)∗. (13) We need not the explicit expression. If we replace those values in (6), we get γA(v) = τ , so we have proved the existence of a right inverse to the map γA : Γx → Ωx. Since  π γA is injective we obtain the desired result.  Remark 3.1. For any Stiefel manifold it is possible to prove that the domain of injectivity x Γ is not the whole vector space TxOn,k. T Definition 3.2. For any x = ∈ O , the open subset  P  n,k τ Ωx = ∈ O | π + P ∗ invertible  π n,k  is called a Cayley open subset of the Stiefel manifold. We continue with an explicit trivialization of the fibration ρ over each Cayley open set. T Proposition 3.3. Let x = ∈ O and let Ωx be the open subset of the elements  P  n,k τ ∈ O such that π + P ∗ is invertible. Then the projection ρ: G(n) → O admits  π n,k n,k a local section sA :Ωx → G(n). 8 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´

⊥ A Proof. With the identification TxOn,k =∼ (TAG(n − k)) , and from the definition of γ A x we can write γ = ρ ◦ cA on Γ . Moreover, we have proved in Theorem A that the A x x A −1 restriction γ |Γx : Γ → Ω is a diffeomorphism whose inverse is denoted (γ |Γx ) . We set A A −1 x s = cA ◦ (γ |Γx ) :Ω → G(n) and verify A A −1 A A −1 x ρ ◦ s = ρ ◦ cA ◦ (γ |Γx ) = γ |Γx ◦ (γ |Γx ) = idΩ . Notice that sA(Ωx) ⊂ Ω(A∗) ∩ G(n).  A A −1 An explicit formula for s could be obtained from those of cA and (γ |Γx ) . α T Proof of Theorem B. We choose A = ∈ G(n). Let sA :Ωx → G(n) be the  β P  local section of Proposition 3.3. With the notations of the statement, we consider the x application H :Ω × [0, 1] → On,k defined by A H(y,t)= ρ(cA(tcA∗ (s (y)))). ∗ ∗ β A This map verifies H(y, 0) = ρ(cA(0)) = ρ(A ) = = γ (0) and H(y, 1) =  P ∗ ρ(sA(y)) = y. Therefore, it is a contraction of Ωx on the point γA(0). 

4. Some applications 4.1. Lusternik-Schnirelmann category of some quaternionic Stiefel manifolds. Let K = H be the algebra of quaternions, G(n) = Sp(n) the symplectic group and On,k = Xn,k the quaternionic Stiefel manifold. With the notations of the proof of Theorem B, we observe that, in general, the point γA(0) does not belong to Ωx. Therefore, our proof x does not imply the contractibility of Ω but only its contractibility in Xn,k, as stated. This property suffices for our first application. Recall that the Lusternik-Schnirelmann category (henceforth LS-category) of a topo- logical space X is the least integer m ≥ 0 such that X admits a covering by (m + 1) open sets which are contractible in X, see [2] for more details. We denote it cat X. The LS-category has applications in a wide range of fields coming from dynamical systems to homotopy theory, but it has also proven to be difficult to determine. For instance, a longstanding problem is the computation of the LS-category of Lie groups. In the case of unitary groups, W. Singhof ([11]) proved that cat U(n) = n by using an argument based on the eigenvalues. Nevertheless, this method cannot be carried out for the quaternionic group Sp(n), see [7]. However, some results have been obtained for small n as cat Sp(2) = 3 ([10]) and cat Sp(3) = 5 ([3]) and for the determination of n+1 some bounds as cat Sp(n) ≥ n + 2 when n ≥ 3 ([5]) and cat Sp(n) ≤ 2 ([8]). The quaternionic Stiefel manifolds Xn,k are more accessible in certain ranges. For
instance, we know that cat Xn,k = k when n ≥ 2k. For proving that in [9] T. Nishimoto uses the number of eigenvalues of a complex matrix in a way similar to Singhof’s approach. This has also been established in [6] by H. Kadzisa and M. Mimura from the Morse-Bott functions defined on Xn,k. In the next proposition, we give a short proof of this result with the Cayley open subsets of Definition 3.2. CAYLEYTRANSFORMONSTIEFELMANIFOLDS 9

Proposition 4.1 ([9], [6]). If n ≥ 2k, we have cat Xn,k ≥ 2k.

Tθ Proof. Let θ ∈]0, π/2[ and take xθ = ∈ Xn,k, with Pθ = (cos θ)Ik and Tθ =  Pθ

0 xθ . We know from Theorem B that Ω is contractible in Xn,k.  (sin θ)Ik We choose (k + 1) numbers θi such that 0 < θ0 < θ2 < · · · < θk < π/2 and observe Hk×k that an element π ∈ such that π + Pθi is not invertible for all i should have (k + 1) xθi distinct real eigenvalues. This is impossible and the family (Ω )0≤i≤k is an open cover of Xn,k by subsets contractible in Xn,k. 

4.2. Optimisation theory. Let G(n) = O(n) be the orthogonal group and On,k = Vn,k the real Stiefel manifold. In optimisation theory, the problems with orthogonality constraints are widely known and have concrete applications in many different areas (see [1] for instance). A typical example is looking for k orthogonal n-vectors that are optimal with respect to some parameter f like cost or likehood. This kind of problems can be seen as optimisation problems on a real Stiefel manifold. The most popular method for this study is the gradient descent method which can be summarized as follows. Let x = x0 be an initial trial point in the Stiefel manifold Vn,k and let F be the negative gradient of f at x. Then a curve α(t) must be found on the manifold such that α(0) = x and α′(0) = F . By fixing a step size τ small enough, the next iterate is obtained by curvilinear search, that is, putting x1 = α(τ). Under certain conditions the sequence x0,x1,... converges to a local minimun of the function f. Most existing methods either use matrix factorizations (such as the SVD decomposi- tion) or require the determination of geodesic curves, which is computationally expensive. A different algorithm has been proposed in [12], where the curve is not a geodesic but is constructed from the Cayley transform in the orthogonal group. Specifically, one con- siders the skew-symmetric matrix A = Fx∗ − xF ∗ and computes the Cayley transform Q(t) = cI (tA) on the group O(n). Since the group acts on the Stiefel manifold, the desired curve can be given by α(t)= Q(t)x. Our construction of a Cayley transform is intrinsic to the Stiefel manifold and should lead to more efficient methods.

Acknowledgements This research was supported by Xunta de Galicia 2015-PG006, by ESF ACAT Grant 6442 and by the ANR-11-LABX-0007-01 “CEMPI” which gave us the opportunity to work together. The three authors are also partially supported by the MINECO and FEDER research project MTM2016-78647-P.

References [1] P.-A. Absil, R. Mahony, and R. Sepulchre, Optimization algorithms on matrix manifolds, Princeton University Press, Princeton, NJ, 2008, With a foreword by Paul Van Dooren. MR 2364186 [2] Octav Cornea, Gregory Lupton, John Oprea, and Daniel Tanr´e, Lusternik-Schnirelmann category, Mathematical Surveys and Monographs, vol. 103, American Mathematical Society, Providence, RI, 2003. MR 1990857 (2004e:55001) 10 E.MAC´IAS-VIRGOS,´ M.-J. PEREIRA-SAEZ,´ AND D. TANRE´

[3] Luc´ıa Fern´andez-Su´arez, Antonio G´omez-Tato, Jeffrey Strom, and Daniel Tanr´e, The Lusternik- Schnirelmann category of Sp(3), Proc. Amer. Math. Soc. 132 (2004), no. 2, 587–595. MR 2022385 (2004m:55005) [4] A. G´omez-Tato, E. Mac´ıas-Virg´os, and M. J. Pereira-S´aez, Trace map, Cayley transform and LS cat- egory of Lie groups, Ann. Global Anal. Geom. 39 (2011), no. 3, 325–335. MR 2769302 (2012d:22011) [5] Norio Iwase and Mamoru Mimura, L-S categories of simply-connected compact simple Lie groups of low rank, Categorical decomposition techniques in algebraic topology (Isle of Skye, 2001), Progr. Math., vol. 215, Birkh¨auser, Basel, 2004, pp. 199–212. MR 2039767 (2005a:55001) [6] Hiroyuki Kadzisa and Mamoru Mimura, Morse-Bott functions and the Lusternik-Schnirelmann cat- egory, J. Fixed Point Theory Appl. 10 (2011), no. 1, 63–85. MR 2825740 (2012k:58016) [7] E. Mac´ıas-Virg´os and M. J. Pereira-S´aez, Symplectic matrices with predetermined left eigenvalues, Linear Algebra Appl. 432 (2010), no. 1, 347–350. MR 2566482 (2010m:15020) [8] , An upper bound for the Lusternik-Schnirelmann category of the symplectic group, Math. Proc. Cambridge Philos. Soc. 155 (2013), no. 2, 271–276. MR 3091519 [9] Tetsu Nishimoto, On the Lusternik-Schnirelmann category of Stiefel manifolds, Topology Appl. 154 (2007), no. 9, 1956–1960. MR 2319267 (2008a:55003) [10] Paul A. Schweitzer, Secondary cohomology operations induced by the diagonal mapping, Topology 3 (1965), 337–355. MR 0182969 (32 #451) [11] Wilhelm Singhof, On the Lusternik-Schnirelmann category of Lie groups, Math. Z. 145 (1975), no. 2, 111–116. MR 0391075 (52 #11897) [12] Zaiwen Wen and Wotao Yin, A feasible method for optimization with orthogonality constraints, Math. Program. 142 (2013), no. 1-2, Ser. A, 397–434. MR 3127080 [13] Fuzhen Zhang, Quaternions and matrices of quaternions, Linear Algebra Appl. 251 (1997), 21–57. MR 1421264

Departamento de Matematicas,´ Universidade de Santiago de Compostela, 15782 Spain E-mail address: [email protected]

Facultade de Econom´ıa e Empresa, Universidade da Coruna,˜ 15071 Spain. E-mail address: [email protected]

Departement´ de Mathematiques,´ Universite´ de Lille 1, 59655 Villeneuve d’Ascq Cedex, France E-mail address: [email protected]