English

A Riemannian Trust Region Method for the Canonical Tensor Rank Approximation Problem

Numerical Analysis 2022-09-02 v2

Abstract

The canonical tensor rank approximation problem (TAP) consists of approximating a real-valued tensor by one of low canonical rank, which is a challenging non-linear, non-convex, constrained optimization problem, where the constraint set forms a non-smooth semi-algebraic set. We introduce a Riemannian Gauss-Newton method with trust region for solving small-scale, dense TAPs. The novelty of our approach is threefold. First, we parametrize the constraint set as the Cartesian product of Segre manifolds, hereby formulating the TAP as a Riemannian optimization problem, and we argue why this parametrization is among the theoretically best possible. Second, an original ST-HOSVD-based retraction operator is proposed. Third, we introduce a hot restart mechanism that efficiently detects when the optimization process is tending to an ill-conditioned tensor rank decomposition and which often yields a quick escape path from such spurious decompositions. Numerical experiments show improvements of up to three orders of magnitude in terms of the expected time to compute a successful solution over existing state-of-the-art methods.

Keywords

Cite

@article{arxiv.1709.00033,
  title  = {A Riemannian Trust Region Method for the Canonical Tensor Rank Approximation Problem},
  author = {Paul Breiding and Nick Vannieuwenhoven},
  journal= {arXiv preprint arXiv:1709.00033},
  year   = {2022}
}
R2 v1 2026-06-22T21:29:38.463Z