Approximate 1-Norm Minimization and Minimum-Rank Structured Sparsity for Various Generalized Inverses via Local Search

Fundamental in matrix algebra and its applications, a generalized inverse of a real matrix $A$ is a matrix $H$ that satisfies the Moore-Penrose (M-P) property $AHA=A$. If $H$ also satisfies the additional M-P property, $HAH=H$, it is called reflexive. Reflexivity is equivalent to minimum rank, so we are particularly interested in reflexive generalized inverses. We consider aspects of symmetry related to the calculation of a sparse reflexive generalized inverse of $A$. As is common, we use (vector) 1-norm minimization for inducing sparsity and for keeping the magnitude of entries under control. When $A$ is symmetric, we may seek a symmetric $H$; while generally such a restriction on $H$ may not lead to a 1-norm minimizing reflexive generalized inverse. We investigate a block construction method to produce a symmetric reflexive generalized inverse that is structured and has guaranteed sparsity, and we provided an approximation algorithm based on it. Another aspect of symmetry that we consider relates to another M-P property: $H$ is ah-symmetric if $AH$ is symmetric. The ah-symmetry property is the key one for solving least-squares problems using $H$. We investigate a column block construction method to produce an ah-symmetric reflexive generalized inverse that is structured and has guaranteed sparsity. Importantly, we provide a theoretically-efficient and practical local-search algorithm to column block construct an approximate 1-norm minimizing ah-symmetric reflexive generalized inverse.

[1]  C. D. Meyer,et al.  Generalized inverses of linear transformations , 1979 .

[2]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[3]  K. S. Banerjee Generalized Inverse of Matrices and Its Applications , 1973 .

[4]  Rémi Gribonval,et al.  Beyond Moore-Penrose Part II: The Sparse Pseudoinverse , 2017, ArXiv.

[5]  Balas K. Natarajan,et al.  Sparse Approximate Solutions to Linear Systems , 1995, SIAM J. Comput..

[6]  R. Penrose A Generalized inverse for matrices , 1955 .

[7]  Nectarios Koziris,et al.  Improving the Performance of the Symmetric Sparse Matrix-Vector Multiplication in Multicore , 2013, 2013 IEEE 27th International Symposium on Parallel and Distributed Processing.

[8]  David P. Williamson,et al.  The Design of Approximation Algorithms , 2011 .

[9]  Rémi Gribonval,et al.  Beyond Moore-Penrose Part I: Generalized Inverses that Minimize Matrix Norms , 2017, ArXiv.

[10]  C. A. Rohde Contributions to the theory, computation and application of generalized inverses. , 1964 .

[11]  Jon Lee,et al.  Diving for Sparse Partially-Reflexive Generalized Inverses , 2019, WCGO.

[12]  Jon Lee,et al.  Sparse pseudoinverses via LP and SDP relaxations of Moore-Penrose , 2016, 1606.06969.

[13]  A. Dresden The fourteenth western meeting of the American Mathematical Society , 1920 .

[14]  A. Bjerhammar Application of calculus of matrices to method of least squares : with special reference to geodetic calculations , 1951 .

[15]  Mihailo Kolundzija,et al.  Beyond Moore-Penrose: Sparse pseudoinverse , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[16]  Jon Lee,et al.  On Sparse Reflexive Generalized Inverses , 2018, Oper. Res. Lett..