Abstract:Given two matrices $X,B\in \mathbb{R}^{n\times m}$ and a set $\mathcal{A}\subseteq \mathbb{R}^{n\times n}$, a Procrustes problem consists in finding a matrix $A \in \mathcal{A}$ such that the Frobenius norm of $AX-B$ is minimized. When $\mathcal{A}$ is the set of the matrices whose symmetric part is positive semidefinite, we obtain the so-called non-symmetric positive semidefinite Procrustes (NSPDSP) problem. The NSPDSP problem arises in the estimation of compliance or stiffness matrix in solid and elastic structures. If $X$ has rank $r$, Baghel et al. (Lin. Alg. Appl., 2022) proposed a three-step semi-analytical approach: (1) construct a reduced NSPDSP problem in dimension $r\times r$, (2) solve the reduced problem by means of a fast gradient method with a linear rate of convergence, and (3) post-process the solution of the reduced problem to construct a solution of the larger original NSPDSP problem. In this paper, we revisit this approach of Baghel et al. and identify an unnecessary assumption used by the authors leading to cases where their algorithm cannot attain a minimum and produces solutions with unbounded norm. In fact, revising the post-processing phase of their semi-analytical approach, we show that the infimum of the NSPDSP problem is always attained, and we show how to compute a minimum-norm solution. We also prove that the symmetric part of the computed solution has minimum rank bounded by $r$, and that the skew-symmetric part has rank bounded by $2r$. Several numerical examples show the efficiency of this algorithm, both in terms of computational speed and of finding optimal minimum-norm solutions.

Successive Projection Method for Solving the Unbalanced Procrustes Problem

Conditions for Optimal Solutions of Unbalanced Procrustes Problem on Stiefel Manifold

A SIMPLE SOLUTION OF SYMMETRIC BALANCED PROCRUSTES PROBLEM

An Eigenvalue-Based Method for the Unbalanced Procrustes Problem

Successive Projection for Solving Systems of Nonlinear Equations/Inequalities

A sketch-and-project method for solving the matrix equation AXB = C

Randomized block Kaczmarz method with projection for solving least squares

Proving Inequalities and Solving Global Optimization Problems Via Simplified CAD Projection

A Projection Method Based on Extended Krylov Subspaces for Solving Sylvester Equations

On Descent-Projection Method for Solving the Split Feasibility Problems

Greedy randomized and maximal weighted residual Kaczmarz methods with oblique projection

Matrix iterative solutions to the least squares problem of BXAT = F with some linear constraints

A Non-Krylov subspace Method for Solving Large and Sparse Linear System of Equations

Robust Matrix Completion Via Alternating Projection.

On Kaczmarz method with oblique projection for solving large overdetermined linear systems

An Efficient Gauss--Newton Algorithm for Symmetric Low-Rank Product Matrix Approximations

KRYLOV SUBSPACE METHODS WITH DEFLATION AND BALANCING PRECONDITIONERS FOR LEAST SQUARES PROBLEMS

Minimum-norm solutions of the non-symmetric semidefinite Procrustes problem

Generalized Nonconvex Nonsmooth Low-Rank Matrix Recovery Framework with Feasible Algorithm Designs and Convergence Analysis.

A block Lanczos method for large-scale quadratic minimization problems with orthogonality constraints

Mixed-Projection Conic Optimization: A New Paradigm for Modeling Rank Constraints