Inverting the sum of two singular matrices
Sofia Eriksson, Jonas Nordqvist
TL;DR
This work addresses inverting a sum of a singular matrix and a structured low-rank perturbation, $\widetilde{\mathbf{A}}=\mathbf{A}+\mathbf{e}D\mathbf{f}^*$, under the conditions that $D$ and $\widetilde{\mathbf{A}}$ are invertible and $\operatorname{rank}(\widetilde{\mathbf{A}})=\operatorname{rank}(\mathbf{A})+\operatorname{rank}(\mathbf{e}D\mathbf{f}^*)$. The main contribution is an explicit inverse of $\widetilde{\mathbf{A}}$ in the form $\widetilde{\mathbf{A}}^{-1}=\mathbf{G}+\mathbf{x}D^{-1}\mathbf{y}^*$, with $\mathbf{G}$, $\mathbf{x}$ and $\mathbf{y}$ computed in a way that does not depend on $D$; the result is first established via singular value decomposition and then recast into a non-SVD formulation. The work clarifies connections to generalized inverses, showing partial satisfaction of Penrose conditions and relating $\mathbf{G}$ to the Moore–Penrose inverse $\mathbf{A}^+$ through projections, while also linking to prior results (e.g., Riedel, 1992). Additionally, a singular-matrix determinant lemma is derived, giving $\det(\widetilde{\mathbf{A}})=\det(\mathbf{A}+\mathbf{e}\mathbf{f}^*)\det(D)$ and $\det(\widetilde{\mathbf{A}}^{-1})=\det(\mathbf{G}+\mathbf{x}\mathbf{y}^*)\det(D^{-1})$, thereby extending Woodbury-type identities to singular bases and aiding practical inversion in applications such as SBP-SAT finite difference schemes.
Abstract
Square matrices of the form $\widetilde{\mathbf{A}} =\mathbf{A} + \mathbf{e}D \mathbf{f}^*$ are considered. An explicit expression for the inverse is given, provided $\widetilde{\mathbf{A}}$ and $D$ are invertible with $\text{rank}(\widetilde{\mathbf{A}}) =\text{rank}(\mathbf{A})+\text{rank}(\mathbf{e}D \mathbf{f}^*)$. The inverse is presented in two ways, one that uses singular value decomposition and another that depends directly on the components $\mathbf{A}$, $\mathbf{e}$, $\mathbf{f}$ and $D$. Additionally, a matrix determinant lemma for singular matrices follows from the derivations.
