A tensor Alternating Anderson-Richardson method for solving multilinear systems with M-tensors
Jing Niu, Lei Du, Tomohiro Sogabe, Shao-Liang Zhang
TL;DR
The paper tackles solving multilinear systems $\mathcal{A}\boldsymbol{x}^{m-1}=\boldsymbol{b}$ where $\mathcal{A}$ is a nonsingular $\mathcal{M}$-tensor. It extends the Alternating Anderson--Richardson framework to tensors, introducing the Tensor Richardson (TR), Tensor Anderson--Richardson (TAR), and Tensor Alternating Anderson--Richardson (TAAR) methods, and embeds TAR periodically to form TAAR. Through theoretical formulation and extensive experiments, TAAR demonstrates at least an order-of-magnitude acceleration over existing tensor-splitting methods while maintaining accuracy, with favorable computational cost profiles. The results support TAAR as an effective, scalable approach for high-dimensional tensor equations arising in applications like PDE discretizations and data-driven tensor problems.
Abstract
It is well-known that a multilinear system with a nonsingular M-tensor and a positive right-hand side has a unique positive solution. Tensor splitting methods generalizing the classical iterative methods for linear systems have been proposed for finding the unique positive solution. The Alternating Anderson-Richardson (AAR) method is an effective method to accelerate the classical iterative methods. In this study, we apply the idea of AAR for finding the unique positive solution quickly. We first present a tensor Richardson method based on tensor regular splittings, then apply Anderson acceleration to the tensor Richardson method and derive a tensor Anderson-Richardson method, finally, we periodically employ the tensor Anderson-Richardson method within the tensor Richardson method and propose a tensor AAR method. Numerical experiments show that the proposed method is effective in accelerating tensor splitting methods.
