# A geometric approach to subspace updates and orthogonal matrix decompositions under rank-one modifications

Research output: Contribution to journalJournal articleResearchpeer-review

## Abstract

For a matrix $X\in \mathbb{R}^{n\times p}$, we provide an analytic formula that keeps track of an orthonormal basis for the range of $X$ under rank-one modifications. More precisely, we consider rank-one adaptations $X_{new} = X+ab^T$ of a given $X$ with known matrix factorization $X = UW$, where $U\in\mathbb{R}^{n\times p}$ is column-orthogonal and $W\in \R^{p\times p}$ is invertible. Arguably, the most important methods that produce such factorizations are the singular value decomposition (SVD), where $X=UW=U(\Sigma V^T)$, and the QR-decomposition, where $X = UW = QR$.
We give a geometric description of rank-one adaptations and derive a closed-form expression for the geodesic line that travels from the subspace $\mathcal{S}= \colspan(X)$ to the subspace $\mathcal{S}_{new} =\colspan(X_{new}) =\colspan(U_{new}W_{new})$. This leads to update formulas for orthogonal matrix decompositions, where both $U_{new}$ and $W_{new}$ are obtained via elementary rank-one matrix updates in $\mathcal{O}(np)$ time for $n\gg p$.
Moreover, this allows us to determine the subspace distance and the Riemannian midpoint between the subspaces $\mathcal{S}$ and $\mathcal{S}_{new}$ without additional computational effort.
Original language English Mathematics of Computation 90 671-688 0025-5718 https://doi.org/10.1090/mcom/3574 Published - 2021

## Keywords

• singular value decomposition
• QR-decomposition
• subspace tracking
• Grassmann manifold
• rank-one update

## Fingerprint

Dive into the research topics of 'A geometric approach to subspace updates and orthogonal matrix decompositions under rank-one modifications'. Together they form a unique fingerprint.
• ### Geometric subspace updates with applications to online adaptive nonlinear model reduction

Zimmermann, R., Peherstorfer, B. & Willcox, K., 2018, 39, 1, p. 234–261

Research output: Contribution to journalJournal articleResearchpeer-review

Open Access
File