English

A new splitting algorithm for dynamical low-rank approximation motivated by the fibre bundle structure of matrix manifolds

Numerical Analysis 2022-03-25 v2 Numerical Analysis

Abstract

In this paper, we propose a new splitting algorithm for dynamical low-rank approximation motivated by the fibre bundle structure of the set of fixed rank matrices. We first introduce a geometric description of the set of fixed rank matrices which relies on a natural parametrization of matrices. More precisely, it is endowed with the structure of analytic principal bundle, with an explicit description of local charts. For matrix differential equations, we introduce a first order numerical integrator working in local coordinates. The resulting algorithm can be interpreted as a particular splitting of the projection operator onto the tangent space of the low-rank matrix manifold. It is proven to be exact in some particular case. Numerical experiments confirm this result and illustrate the robustness of the proposed algorithm.

Keywords

Cite

@article{arxiv.2001.08599,
  title  = {A new splitting algorithm for dynamical low-rank approximation motivated by the fibre bundle structure of matrix manifolds},
  author = {Marie Billaud-Friess and Antonio Falcó and Anthony Nouy},
  journal= {arXiv preprint arXiv:2001.08599},
  year   = {2022}
}