Search

MiS Preprint Repository

We have decided to discontinue the publication of preprints on our preprint server as of 1 March 2024. The publication culture within mathematics has changed so much due to the rise of repositories such as ArXiV (www.arxiv.org) that we are encouraging all institute members to make their preprints available there. An institute's repository in its previous form is, therefore, unnecessary. The preprints published to date will remain available here, but we will not add any new preprints here.

MiS Preprint
48/2018

Pencil-based algorithms for tensor rank decomposition are not stable

Carlos Beltrán, Paul Breiding and Nick Vannieuwenhoven

Abstract

We prove the existence of an open set of $n_1\times n_2 \times n_3$ tensors of rank $r$ on which a popular and efficient class of algorithms for computing tensor rank decompositions based on a reduction to a linear matrix pencil, typically followed by a generalized eigendecomposition, is arbitrarily numerically forward unstable. Our analysis shows that this problem is caused by the fact that the condition number of the tensor rank decomposition can be much larger for $n_1 \times n_2 \times 2$ tensors than for the $n_1\times n_2 \times n_3$ input tensor. Moreover, we present a lower bound for the limiting distribution of the condition number of random tensor rank decompositions of third-order tensors. The numerical experiments illustrate that for random tensor rank decompositions one should anticipate a loss of precision of a few digits.

Received:
12.07.18
Published:
17.07.18

Related publications

inJournal
2019 Repository Open Access
Carlos Beltrán, Paul Breiding and Nick Vannieuwenhoven

Pencil-based algorithms for tensor rank decomposition are not stable

In: SIAM journal on matrix analysis and applications, 40 (2019) 2, pp. 739-773