Pencil-Based Algorithms for Tensor Rank Decomposition are not Stable

Carlos Beltrán,Paul Breiding,Nick Vannieuwenhoven
DOI: https://doi.org/10.1137/18M1200531
IF: 1.908
2019-10-07
SIAM Journal on Matrix Analysis and Applications
Abstract:We prove the existence of an open set of $n_1\times n_2 \times n_3$ tensors of rank $r$ for which popular and efficient algorithms for computing tensor rank decompositions based on a reduction to a linear matrix pencil, typically followed by a generalized eigendecomposition, are arbitrarily numerically forward unstable. Our analysis shows that this problem is caused by the fact that the condition number of 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.
mathematics, applied
What problem does this paper attempt to address?