An iterative algorithm for a least squares solution of a matrix equation.

In this article, we propose an iterative algorithm to compute the minimum norm least-squares solution of AXB+CYD=E, based on a matrix form of the algorithm LSQR for solving the least squares problem. We then apply this algorithm to compute the minimum norm least-squares centrosymmetric solution of min X|AXB-E|F. Numerical results are provided to verify the efficiency of the proposed method. [ABSTRACT FROM AUTHOR]