Table of Contents
Fetching ...

Adaptive Riemannian ADMM for Nonsmooth Optimization: Optimal Complexity without Smoothing

Kangkang Deng, Jiachen Jin, Jiang Hu, Hongxia Wang

TL;DR

This work addresses nonsmooth optimization on compact Riemannian manifolds for a composite objective $f(x) + h(\mathcal A x)$. It introduces Adaptive Riemannian ADMM (ARADMM), which avoids smoothing by adaptively tuning dual stepsizes $\gamma_k$ and penalty parameters $\rho_k$, so that each iteration requires only one Riemannian gradient evaluation and one proximal update. The authors prove an optimal iteration complexity of $O(\epsilon^{-3})$ to obtain an $\epsilon$-approximate KKT point, matching smoothing-based methods while directly handling the original nonsmooth problem. Numerical experiments on sparse PCA and robust subspace recovery demonstrate that ARADMM consistently outperforms state-of-the-art Riemannian ADMM variants in convergence speed and solution quality.

Abstract

We study the problem of minimizing the sum of a smooth function and a nonsmooth convex regularizer over a compact Riemannian submanifold embedded in Euclidean space. By introducing an auxiliary splitting variable, we propose an adaptive Riemannian alternating direction method of multipliers (ARADMM), which, for the first time, achieves convergence without requiring smoothing of the nonsmooth term. Our approach involves only one Riemannian gradient evaluation and one proximal update per iteration. Through careful and adaptive coordination of the stepsizes and penalty parameters, we establish an optimal iteration complexity of order $\mathcal{O}(ε^{-3})$ for finding an $ε$-approximate KKT point, matching the complexity of existing smoothing technique-based Riemannian ADMM methods. Extensive numerical experiments on sparse PCA and robust subspace recovery demonstrate that our ARADMM consistently outperforms state-of-the-art Riemannian ADMM variants in convergence speed and solution quality.

Adaptive Riemannian ADMM for Nonsmooth Optimization: Optimal Complexity without Smoothing

TL;DR

This work addresses nonsmooth optimization on compact Riemannian manifolds for a composite objective . It introduces Adaptive Riemannian ADMM (ARADMM), which avoids smoothing by adaptively tuning dual stepsizes and penalty parameters , so that each iteration requires only one Riemannian gradient evaluation and one proximal update. The authors prove an optimal iteration complexity of to obtain an -approximate KKT point, matching smoothing-based methods while directly handling the original nonsmooth problem. Numerical experiments on sparse PCA and robust subspace recovery demonstrate that ARADMM consistently outperforms state-of-the-art Riemannian ADMM variants in convergence speed and solution quality.

Abstract

We study the problem of minimizing the sum of a smooth function and a nonsmooth convex regularizer over a compact Riemannian submanifold embedded in Euclidean space. By introducing an auxiliary splitting variable, we propose an adaptive Riemannian alternating direction method of multipliers (ARADMM), which, for the first time, achieves convergence without requiring smoothing of the nonsmooth term. Our approach involves only one Riemannian gradient evaluation and one proximal update per iteration. Through careful and adaptive coordination of the stepsizes and penalty parameters, we establish an optimal iteration complexity of order for finding an -approximate KKT point, matching the complexity of existing smoothing technique-based Riemannian ADMM methods. Extensive numerical experiments on sparse PCA and robust subspace recovery demonstrate that our ARADMM consistently outperforms state-of-the-art Riemannian ADMM variants in convergence speed and solution quality.
Paper Structure (17 sections, 9 theorems, 64 equations, 3 figures, 7 tables, 1 algorithm)

This paper contains 17 sections, 9 theorems, 64 equations, 3 figures, 7 tables, 1 algorithm.

Key Result

Proposition 2.1

Let $\mathcal{R}$ be a retraction operator on a compact submanifold $\mathcal{M}$. Then, there exist two positive constants $\alpha, \beta$ such that for all $x\in \mathcal{M}$ and all $u \in T_{x}\mathcal{M}$, we have

Figures (3)

  • Figure 1: Comparison with ADMM-type methods for solving (\ref{['SPCA']}) with different $(n,p)$, $m=n$ and $\mu=0.01$.
  • Figure 2: Comparison with ADMM-type methods for solving (\ref{['DPCP']}) with different $(n,p)$ and $(p_1,p_2)=(150,1000)$.
  • Figure 3: Comparison with ADMM-type methods for solving (\ref{['DPCP']}) with different $(p_1,p_2)$ and $(n,p)=(45,5)$.

Theorems & Definitions (19)

  • Definition 2.1: Retraction, AbsMahSep2008
  • Definition 2.2: Vector transport, AbsMahSep2008
  • Proposition 2.1: grocf
  • Definition 2.3
  • Lemma 2.1
  • Lemma 3.1
  • proof : Proof of Lemma \ref{['lem:bound-lambda']}
  • Theorem 3.1
  • Lemma A.1
  • Definition A.1
  • ...and 9 more