English

Rank-one Riemannian Subspace Descent for Nonlinear Matrix Equations

Numerical Analysis 2026-01-22 v1 Numerical Analysis Systems and Control Systems and Control

Abstract

We propose a rank-one Riemannian subspace descent algorithm for computing symmetric positive definite (SPD) solutions to nonlinear matrix equations arising in control theory, dynamic programming, and stochastic filtering. For solution matrices of size n×nn\times n, standard approaches for dense matrix equations typically incur O(n3)\mathcal{O}(n^3) cost per-iteration, while the efficient O(n2)\mathcal{O}(n^2) methods either rely on sparsity or low-rank solutions, or have iteration counts that scale poorly. The proposed method entails updating along the dominant eigen-component of a transformed Riemannian gradient, identified using at most O(log(n))\mathcal{O}(\log(n)) power iterations. The update structure also enables exact step-size selection in many cases at minimal additional cost. For objectives defined as compositions of standard matrix operations, each iteration can be implemented using only matrix--vector products, yielding O(n2)\mathcal{O}(n^2) arithmetic cost. We prove an O(n)\mathcal{O}(n) iteration bound under standard smoothness assumptions, with improved bounds under geodesic strong convexity. Numerical experiments on large-scale CARE, DARE, and other nonlinear matrix equations show that the proposed algorithm solves instances (up to n=10,000n=10{,}000 in our tests) for which the compared solvers, including MATLAB's \texttt{icare}, structure-preserving doubling, and subspace-descent baselines fail to return a solution. These results demonstrate that rank-one manifold updates provide a practical approach for high-dimensional and dense SPD-constrained matrix equations. MATLAB code implementation is publicly available on GitHub : \href{https://github.com/yogeshd-iitk/nonlinear_matrix_equation_R1RSD}{\textcolor{blue}{https://github.com/yogeshd-iitk/nonlinear\_matrix \_equation\_R1RSD}}

Keywords

Cite

@article{arxiv.2601.14933,
  title  = {Rank-one Riemannian Subspace Descent for Nonlinear Matrix Equations},
  author = {Yogesh Darmwal and Ketan Rajawat},
  journal= {arXiv preprint arXiv:2601.14933},
  year   = {2026}
}
R2 v1 2026-07-01T09:14:00.331Z