“…where A ∈ R n×n , B ∈ R m×m and X, C ∈ R n×m . This matrix equation plays an important role in many practical applications such as surfaces fitting in computer aided geometric design (CAGD), signal and image processing, photogrammetry, etc., see for example [14,19,20,30,31] and a large literature therein. Given the existence of many and important applications, the theory and algorithms for matrix equations have intrigued the researchers for decades, see [7,9,15,21,28,29,30,31,34,37].…”