Riemannian optimization on the symplectic Stiefel manifold using second-order information

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Jensen, Rasmus, Zimmermann, Ralf
Format: Preprint
Veröffentlicht: 2024
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866917705087451136
author Jensen, Rasmus
Zimmermann, Ralf
author_facet Jensen, Rasmus
Zimmermann, Ralf
contents Riemannian optimization is concerned with problems, where the independent variable lies on a smooth manifold. There is a number of problems from numerical linear algebra that fall into this category, where the manifold is usually specified by special matrix structures, such as orthogonality or definiteness. Following this line of research, we investigate tools for Riemannian optimization on the symplectic Stiefel manifold. We complement the existing set of numerical optimization algorithms with a Riemannian trust region method tailored to the symplectic Stiefel manifold. To this end, we derive a matrix formula for the Riemannian Hessian under a right-invariant metric. Moreover, we propose a novel retraction for approximating the Riemannian geodesics. Finally, we conduct a comparative study in which we juxtapose the performance of the Riemannian variants of the steepest descent, conjugate gradients, and trust region methods on selected matrix optimization problems that feature symplectic constraints.
format Preprint
id arxiv_https___arxiv_org_abs_2404_08463
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Riemannian optimization on the symplectic Stiefel manifold using second-order information
Jensen, Rasmus
Zimmermann, Ralf
Numerical Analysis
53B20, 22E70, 53B25, 65F15, 53Z05, 70G45, 65P10
Riemannian optimization is concerned with problems, where the independent variable lies on a smooth manifold. There is a number of problems from numerical linear algebra that fall into this category, where the manifold is usually specified by special matrix structures, such as orthogonality or definiteness. Following this line of research, we investigate tools for Riemannian optimization on the symplectic Stiefel manifold. We complement the existing set of numerical optimization algorithms with a Riemannian trust region method tailored to the symplectic Stiefel manifold. To this end, we derive a matrix formula for the Riemannian Hessian under a right-invariant metric. Moreover, we propose a novel retraction for approximating the Riemannian geodesics. Finally, we conduct a comparative study in which we juxtapose the performance of the Riemannian variants of the steepest descent, conjugate gradients, and trust region methods on selected matrix optimization problems that feature symplectic constraints.
title Riemannian optimization on the symplectic Stiefel manifold using second-order information
topic Numerical Analysis
53B20, 22E70, 53B25, 65F15, 53Z05, 70G45, 65P10
url https://arxiv.org/abs/2404.08463