1. A tensor-based approach to solving linear systems involving Kronecker sum of matrices
- Author
-
Al-Dweik, Ahmad Y. and Sayyed-Ahmad, Abdallah
- Subjects
Mathematics - General Mathematics - 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})$.
- Published
- 2024