# A Proximal Approach for a Class of Matrix Optimization Problems

arXiv: https://arxiv.org/abs/1801.07452 · HTML: https://mlsanity.com/paper/1801.07452/

## Key facts

- Cited by 8 papers in the map (data through 14 Oct 2026).
- First posted on arXiv 23 Jan 2018.

## Built on by (most cited first, 8 of 8)

- [Learning Graphical Factor Models with Riemannian Optimization](https://mlsanity.com/paper/2210.11950/) (2022)
- [Sparse Graphical Linear Dynamical Systems](https://mlsanity.com/paper/2307.03210/) (2023)
- [Efficient Hierarchical Bayesian Inference for Spatio-temporal Regression Models in Neuroimaging](https://mlsanity.com/paper/2111.01692/) (2021)
- [On Centralized and Distributed Mirror Descent: Convergence Analysis Using Quadratic Constraints](https://mlsanity.com/paper/2105.14385/) (2021)
- [Implicit Differentiation for Hyperparameter Tuning the Weighted Graphical Lasso](https://mlsanity.com/paper/2307.02130/) (2023)
- [Compressive Recovery of Sparse Precision Matrices](https://mlsanity.com/paper/2311.04673/) (2023)
- [Leveraging Low-rank Factorizations of Conditional Correlation Matrices in Graph Learning](https://mlsanity.com/paper/2506.10628/) (2025)
- [Convergence of projected stochastic natural gradient variational inference for various step size and sample or batch size schedules](https://mlsanity.com/paper/2604.00683/) (2026)

## Abstract (arXiv)

In recent years, there has been a growing interest in mathematical models leading to the minimization, in a symmetric matrix space, of a Bregman divergence coupled with a regularization term. We address problems of this type within a general framework where the regularization term is split in two parts, one being a spectral function while the other is arbitrary. A Douglas-Rachford approach is proposed to address such problems and a list of proximity operators is provided allowing us to consider various choices for the fit-to-data functional and for the regularization term. Numerical experiments show the validity of this approach for solving convex optimization problems encountered in the context of sparse covariance matrix estimation. Based on our theoretical results, an algorithm is also proposed for noisy graphical lasso where a precision matrix has to be estimated in the presence of noise. The nonconvexity of the resulting objective function is dealt with a majorization-minimization approach, i.e. by building a sequence of convex surrogates and solving the inner optimization subproblems via the aforementioned Douglas-Rachford procedure. We establish conditions for the convergence of this iterative scheme and we illustrate its good numerical performance with respect to state-of-the-art approaches.