A tensor-based approach to solving linear systems involving Kronecker sum of matrices

Ahmad Y. Al-Dweik,Abdallah Sayyed-Ahmad
2024-10-05
Abstract:A tensor-based formula for solving the linear system involving Kronecker sum is proposed. It is shown that the linear system involving Kronecker sum of matrices is directly related to the matrix and tensor forms of Sylvester equation. The new tensor-based formula indicates that a Sylvester tensor equation has a unique solution if the sum of spectra of the matrices does not contain zero. As an application, we solve the discretized Poisson equation on a rectangular domain and Dirichlet boundary conditions. In 2d, we show that both the new formula and the fast Poisson solver based on diagonalization have the same cost $\mathcal {O}(n^{3/2})$. In 3d, the new formula has the cost $\mathcal {O}(n^{4/3})$ while cyclic reduction and successive over-relaxation methods have a complexity of $\mathcal {O}(n^{3/2})$.
General Mathematics
What problem does this paper attempt to address?