The submatrix constraint problem of matrix equation \(AXB+CYD=E\) (Q1044427)

This paper deals with the real matrix equation \(AXB+CYD=E\), where \(X\) of order \(n\) is symmetric centrosymmetric \(-\) that is: \(x_{ij}=x_{ji}=x_{n-j+1,n-i+1}\), and \(Y\) of order \(m\) is symmetric. It is presented an algorithm to minimize the \(\| AXB+CYD-E\|\) in the Frobenious norm, using the conjugate gradient square method. Then, a second algorithm for the optimal approximation of the solution set is derived. Also the stability of the algorithms and numerical behaviors are studied.