Abstract
This article introduces a new map on the symplectic Stiefel manifold. The operation that requires the highest computational cost to compute the novel retraction is a inversion of size 2p-by-2p, which is much less expensive than those required for the available retractions in the literature. Later, with the new retraction, we design a constraint preserving gradient method to minimize smooth functions defined on the symplectic Stiefel manifold. To improve the numerical performance of our approach, we use the non-monotone line-search of Zhang and Hager with an adaptive Barzilai–Borwein type step-size. Our numerical studies show that the proposed procedure is computationally promising and is a very good alternative to solve large-scale optimization problems over the symplectic Stiefel manifold.
Similar content being viewed by others
Data Availability
The data used in this research is freely accessible. This can be downloaded at https://sparse.tamu.edu/.
Notes
The Riemannian gradient methods Cayley and Qgeodesic can be downloaded from https://github.com/opt-gaobin/spopt.
The SuiteSparse Matrix Collection tool-box is available in https://sparse.tamu.edu/.
References
Absil PA, Mahony R, Sepulchre R (2009) Optimization algorithms on matrix manifolds. Princeton University Press, Princeton. https://doi.org/10.1515/9781400830244
Amodio P (2003) A symplectic Lanczos-type algorithm to compute the eigenvalues of positive definite Hamiltonian matrices. In: Peter MAS, David A, Alexander VB, Yuriy EG, Jack JD, Albert YZ (eds) Lecture notes in computer science, vol 2658, pp 139–148. https://doi.org/10.1007/3-540-44862-4_16
Arnol’d VI (2013) Mathematical methods of classical mechanics. Springer, New York. https://doi.org/10.1007/978-1-4757-1693-1
Bendokat T, Zimmermann R (2021) The real symplectic Stiefel and Grassmann manifolds: metrics, geodesics and applications. arXiv:2108.12447. Accessed 30 Jan 2023
Benner P, Fabbender H (1997) An implicitly restarted symplectic Lanczos method for the Hamiltonian eigenvalue problem. Linear Algebra Appl 263:75–111. https://doi.org/10.1016/S0024-3795(96)00524-1
Benner P, Fa\(\beta \)bender H, Watkins DS (1999) SR and SZ algorithms for the symplectic (butterfly) eigenproblem. Linear Algebra Appl 287(1-3):41–76. https://doi.org/10.1016/S0024-3795(98)10090-3
Benner P, Fa\(\beta \)bender H (1998) The symplectic eigenvalue problem, the butterfly form, the SR algorithm, and the Lanczos method. Linear Algebra Appl 275:19-47. https://doi.org/10.1016/S0024-3795(97)10049-0
Davis TA, Hu Y (2011) The university of Florida sparse matrix collection. ACM Trans Math Softw 38(1):1–25. https://doi.org/10.1145/2049662.2049663
Del Buono N, Lopez L, Politi T (2008) Computation of functions of Hamiltonian and skew-symmetric matrices. Math Comput Simul 79(4):1284–1297. https://doi.org/10.1016/j.matcom.2008.03.011
Di Serafino D, Ruggiero V, Toraldo G, Zanni L (2018) On the steplength selection in gradient methods for unconstrained optimization. Appl Math Comput 318:176–195. https://doi.org/10.1016/j.amc.2017.07.037
Dutta B, Mukunda N, Simon R (1995) The real symplectic groups in quantum mechanics and optics. Pramana 45(6):471–497. https://doi.org/10.1007/BF02848172
Eisert J, Tyc T, Rudolph T, Sanders BC (2008) Gaussian quantum marginal problem. Commun Math Phys 280(1):263–280. https://doi.org/10.1007/s00220-008-0442-4
Fa\(\beta \)bender H (2001) The parameterized SR algorithm for symplectic (butterfly) matrices. Math Comput 70(236):1515–1541. https://doi.org/10.1090/S0025-5718-00-01265-5
Fiori S (2016) A Riemannian steepest descent approach over the inhomogeneous symplectic group: application to the averaging of linear optical systems. Appl Math Comput 283:251–264. https://doi.org/10.1016/j.amc.2016.02.018
Gao B, Son NT, Absil PA, Stykel T (2021) Riemannian optimization on the symplectic Stiefel manifold. SIAM J Optim 31(2):1546–1575. https://doi.org/10.1137/20M1348522
Gao B, Son N, Stykel T (2022) Optimization on the symplectic Stiefel manifold: SR decomposition-based retraction and applications. arXiv:2211.09481
Hiroshima T (2006) Additivity and multiplicativity properties of some Gaussian channels for Gaussian inputs. Phys Rev A 73(1):012330. https://doi.org/10.1103/PhysRevA.73.012330
Hofer H, Zehnder E (2012) Symplectic invariants and Hamiltonian dynamics. Birkhäuser. https://doi.org/10.1007/978-3-0348-8540-9
Hu J, Liu X, Wen ZW, Yuan YX (2020) A brief introduction to manifold optimization. J Oper Res Soc China 8(2):199–248. https://doi.org/10.1007/s40305-020-00295-9
Iannazzo B, Porcelli M (2018) The Riemannian Barzilai-Borwein method with nonmonotone line search and the matrix geometric mean computation. IMA J Numer Anal 38(1):495–517. https://doi.org/10.1093/imanum/drx015
Jiang B, Dai YH (2015) A framework of constraint preserving update schemes for optimization on Stiefel manifold. Math Program 153(2):535–575. https://doi.org/10.1007/s10107-014-0816-7
Lei M, Meng G (2016) Symplectic principal component analysis: a noise reduction method for continuous chaotic systems. In: Noor A (ed) Advances in noise analysis, mitigation and control. IntechOpen, p 23. https://doi.org/10.5772/64410
Lei M, Meng G, Zhang W, Wade J, Sarkar N (2016) Symplectic entropy as a novel measure for complex systems. Entropy 18(11):412. https://doi.org/10.3390/e18110412
Nocedal J, Wright SJ (2006) Numerical optimization. Springer, New York. https://doi.org/10.1007/978-0-387-40065-5
Oviedo H (2022) Global convergence of Riemannian line search methods with a Zhang-Hager-type condition. Numer Algor. https://doi.org/10.1007/s11075-022-01298-8
Oviedo H (2022) Implicit steepest descent algorithm for optimization with orthogonality constraints. Optim Lett 16(6):1773–1797. https://doi.org/10.1007/s11590-021-01801-5
Oviedo H, Dalmau O, Herrera R (2021) Two novel gradient methods with optimal step sizes. J Comput Math 39(3):375–391. https://doi.org/10.4208/jcm.2001-m2018-0205
Oviedo H, Dalmau O, Lara H (2021) Two adaptive scaled gradient projection methods for Stiefel manifold constrained optimization. Numer Algor 87(3):1107–1127. https://doi.org/10.1007/s11075-020-01001-9
Oviedo H, Lara H, Dalmau O (2019) A non-monotone linear search algorithm with mixed direction on Stiefel manifold. Optim Methods Softw 34(2):437–457. https://doi.org/10.1080/10556788.2017.1415337
Parra L (1995) Symplectic nonlinear component analysis. Adv Neural Inf Process Syst 1995:437–443
Peng L, Mohseni K (2016) Symplectic model reduction of Hamiltonian systems. SIAM J Sci Comput 38(1):A1–A27. https://doi.org/10.1137/140978922
Raydan M (1997) The Barzilai and Borwein gradient method for the large scale unconstrained minimization problem. SIAM J Optim 7(1):26–33. https://doi.org/10.1137/S1052623494266365
Salam A, Al-Aidarous E (2014) Equivalence between modified symplectic Gram-Schmidt and Householder SR algorithms. BIT Numer Math 54(1):283–302. https://doi.org/10.1007/s10543-013-0441-5
Son NT, Absil PA, Gao B, Stykel T (2021) Computing symplectic eigenpairs of symmetric positive-definite matrices via trace minimization and Riemannian optimization. SIAM J Matrix Anal Appl 42(4):1732–1757. https://doi.org/10.1137/21M1390621
Van Loan C (1984) A symplectic method for approximating all the eigenvalues of a Hamiltonian matrix. Linear Algebra Appl 61:233–251. https://doi.org/10.1016/0024-3795(84)90034-X
Wen Z, Yin W (2013) A feasible method for optimization with orthogonality constraints. Math Program 142(1):397–434. https://doi.org/10.1007/s10107-012-0584-1
Williamson J (1936) On the algebraic problem concerning the normal forms of linear dynamical systems. Am J Math 58(1):141–163. https://doi.org/10.2307/2371062
Wu R, Chakrabarti R, Rabitz H (2008) Optimal control theory for continuous-variable quantum gates. Phys Rev A 77(5):052303. https://doi.org/10.1103/PhysRevA.77.052303
Zhang H, Hager WW (2004) A nonmonotone line search technique and its application to unconstrained optimization. SIAM J Optim 14(4):1043–1056. https://doi.org/10.1137/S1052623403428208
Zhao L (2022) Linear constraint problem of Hermitian unitary symplectic matrices. Linear Multilinear Algebra 70(8):1423–1441. https://doi.org/10.1080/03081087.2020.1762533
Acknowledgements
The first author was financially supported by the Faculty of Engineering and Sciences, Universidad Adolfo Ibáñez, through the FES startup package for scientific research. The second author was financially supported in part by CONACYT (Mexico), Grants 256126.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Communicated by yimin wei.
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Oviedo, H., Herrera, R. A collection of efficient retractions for the symplectic Stiefel manifold. Comp. Appl. Math. 42, 164 (2023). https://doi.org/10.1007/s40314-023-02302-0
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s40314-023-02302-0