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
51/2016
Perturbation of higher-order singular values
Wolfgang Hackbusch, Daniel Kressner and André Uschmajew
Abstract
The higher-order singular values for a tensor of order $d$ are defined as the singular values of the $d$ different matricizations associated with the multilinear rank. When $d\geq3$, the singular values are generally different for different matricizations but not completely independent. Characterizing the set of feasible singular values turns out to be difficult. In this work, we contribute to this question by investigating which first-order perturbations of the singular values for a given tensor are possible. We prove that, except for trivial restrictions, any perturbation of the singular values can be achieved for almost every tensor with identical mode sizes. This settles a conjecture from [Hackbusch and Uschmajew, 2016] for the case of identical mode sizes. Our theoretical results are used to develop and analyze a variant of the Newton method for constructing a tensor with specified higher-order singular values or, more generally, with specified Gramians for the matricizations. We establish local quadratic convergence and demonstrate the robust convergence behavior with numerical experiments.
