Iterative Methods for Symmetric Outer Product Tensor Decompositions

Na Li,Carmeliza Navasca
DOI: https://doi.org/10.48550/arXiv.1312.7629
2013-12-30
Abstract:We study the symmetric outer product decomposition which decomposes a fully (partially) symmetric tensor into a sum of rank-one fully (partially) symmetric tensors. We present iterative algorithms for the third-order partially symmetric tensor and fourth-order fully symmetric tensor. The numerical examples indicate a faster convergence rate for the new algorithms than the standard method of alternating least squares.
Numerical Analysis
What problem does this paper attempt to address?
The main problems that this paper attempts to solve are the inefficiency and slow convergence speed of existing methods in dealing with symmetric tensor decomposition, especially for third - order and fourth - order tensors that are partially symmetric and fully symmetric. Specifically: 1. **Limitations of existing methods**: - **Alternating Least Squares (ALS)**: Although widely used, it has an extremely slow convergence speed when dealing with symmetric tensors, and in some cases, the symmetry of the solution cannot be guaranteed. - **Non - linear sub - problems**: When dealing with partially symmetric tensors, ALS needs to solve three identical non - linear least - squares sub - problems, which leads to an increase in computational complexity. 2. **Proposed new method**: - **Partial Column Least Squares (PCLS)**: The paper proposes a new iterative algorithm - Partial Column Least Squares (PCLS) - for solving Symmetric Outer Product Tensor Decomposition (SOPD). This method improves computational efficiency and convergence speed by optimizing matrixing and solving polynomial roots to avoid non - linear least - squares sub - problems. - **Applicable to tensors of different orders**: PCLS is not only applicable to third - order partially symmetric tensors but also extended to fourth - order fully symmetric tensors, and its superiority has been verified through numerical experiments. 3. **Specific problem description**: - **Third - order partially symmetric tensor**: Given a third - order tensor \( T\in\mathbb{R}^{I\times I\times K} \), where \( t_{ijk} = t_{jik} \), the goal is to find the best sum of rank - one partially symmetric tensors so that the error is minimized. - **Fourth - order fully symmetric tensor**: Given a fourth - order tensor \( T\in\mathbb{R}^{I\times I\times I\times I} \), where \( t_{ijkl}=t_{\sigma(ijkl)} \) for any permutation \( \sigma \), the goal is to find the best sum of rank - one fully symmetric tensors so that the error is minimized. 4. **Results of numerical experiments**: - The paper compares the performance of PCLS and ALS through multiple numerical experiments, and the results show that PCLS is superior to ALS in terms of the number of iterations, CPU time, and convergence stability. In summary, this paper aims to solve the problems of slow convergence and high computational complexity of the existing ALS method in dealing with symmetric tensor decomposition by proposing the PCLS method, thereby providing a more efficient and stable solution.