Laurent Condat, Abdurakhmon Sadiev, Peter Richtárik
10 min
This paper by Laurent Condat, Abdurakhmon Sadiev, and Peter Richtárik introduces APAPC (Accelerated Proximal Alternating Predictor–Corrector), a breakthrough algorithm that brings Nesterov acceleration to primal-dual splitting methods for solving structured convex optimization problems. These problems, common in signal processing, machine learning, and inverse problems, take the form:
Here, is smooth and convex (like a data fidelity term), and are nonsmooth convex (regularizers or indicators), and is a linear operator (e.g., blurring or subsampling matrix). Traditional primal-dual methods like PDHG or Condat-Vũ split the problem using proximal operators and gradients, but accelerating them—especially w.r.t. the smooth —is tricky due to 'rotational dynamics' in the primal-dual space that destabilize momentum.
Nesterov acceleration boosts first-order methods from to convergence by adding momentum: instead of stepping from the current point , you predict an extrapolated point and step from there. In proximal gradient (), this works beautifully. But in primal-dual settings with , the coupling via and creates oscillations that break naive momentum. The paper's key insight: exploit dual strong convexity (induced by ) to stabilize accelerated primal updates, allowing seamless integration into a forward-backward primal-dual scheme.
APAPC is an accelerated version of the Proximal Alternating Predictor–Corrector (PAPC) method. It alternates primal and dual updates with prediction-correction steps and Nesterov momentum on the primal side:
Stepsizes satisfy and , offering larger ranges than competitors like Condat-Vũ (CV). This fully splits the problem: only needs , , (or ), , . When , it recovers accelerated three-operator splitting (TOS).
Using a unified Lyapunov framework, the authors prove:
The Lyapunov function combines objective errors, squared distances to optima, and a novel correction term to handle momentum stability.
Primal-dual iterates converge weakly to a saddle-point satisfying KKT conditions , . This leverages recent accelerated gradient descent results. The work generalizes PD3O/PDDY algorithms, improves stepsize flexibility over CV/PAPC, and characterizes primal-dual condition number scaling, with acceleration reducing dependence from linear to square-root. Applications span imaging (TV denoising), ML (structured SVMs), and control.
In sum, APAPC makes Nesterov acceleration practical for a wide class of nonsmooth problems, closing a long-standing gap in proximal splitting methods.
We investigate the integration of Nesterov-type acceleration into primal-dual methods for structured convex optimization. While proximal splitting algorithms efficiently handle composite problems of the form $\min_x f(x)+g(x)+h(Kx)$, accelerating their convergence with respect to the smooth term $f$ is notoriously challenging due to the rotational dynamics in the primal-dual space. In this paper, we overcome this barrier by proposing the Accelerated Proximal Alternating Predictor-Corrector algorithm (APAPC), focusing on the setting where $g(x)=\frac{μ_g}{2}\|x\|^2$. Our analysis reveals that Nesterov momentum can be seamlessly integrated into a primal-dual forward-backward scheme by exploiting the strong convexity of the dual problem to stabilize the accelerated primal updates. Using a unified Lyapunov framework, we establish optimal $O(1/t^2)$ sublinear convergence rates, as well as accelerated linear convergence when $μ_g > 0$, across three regimes of dual strong convexity: (i) when $h$ is smooth, (ii) when the linear operator $K^*$ is bounded below, and (iii) for linearly constrained optimization. Furthermore, leveraging recent results on accelerated gradient descent, we characterize the weak convergence of the primal-dual iterates to a saddle-point solution.
Alex: Wait, so the dual bowl counters the wobble from primal lunges? But how do they prove it doesn't still diverge?
Sam: They craft a single score, called a Lyapunov function, that bundles distances from the current points to the ideal solution, plus a term for primal-dual imbalance—like a dashboard tracking speed, position error, and steering alignment all shrinking together. Each step shows this score drops reliably, under three dual-strength cases: smooth composite terms, bounded-below ones, or linear constraints. This yields optimal progress—about twice as fast long-term as standard methods—and, newly, proves the points weakly converge to balance. For image denoising, where you smooth data but enforce edges sharply, it cuts iterations meaningfully.
Alex: Huh, so that dashboard glues it all, turning wobbles into steady shrink. Makes the speed-up credible for real tasks like cleaning noisy photos.
Alex: That Lyapunov dashboard is clever—it bundles everything to show steady progress. But in the APGD building block you mentioned, how does it nail those faster rates when there's strong convexity on the nonsmooth side?
Sam: In APGD, strong convexity means the nonsmooth cost has a reliable bowl shape pulling toward the best point, like a ball in a valley that always rolls back if nudged. When that's present alongside smooth parts, the method switches to linear convergence: after initial speedup, error to the optimum shrinks by a fixed fraction each step. They prove this via the Lyapunov, showing distances and function gaps contract at an optimal rate for first-order methods; to hit tiny error takes far fewer steps than sublinear.
Alex: So it halves error with steps squared long-term, but linearly when bowl-shaped—like compound interest in reverse. Does it actually reach the exact solution, or just get close in some averaged way?
Sam: Historically, accelerated methods left a gap: function values approached optimum, but the points themselves might wobble without settling. This paper fills it for APGD, proving weak convergence—the sequence has one cluster point that's the minimizer, so iterates approach it gently. Bounded auxiliary points ensure averages settle; practically, choose momentum ramping to blend sublinear burst then linear lock-in.
Alex: Huh, so points don't just circle—they cluster at the goal. How does this feed into the full primal-dual APAPC without losing that?
Sam: APAPC slots APGD's predictor into PAPC's predict-correct loop for primal and dual. It predicts primal with current dual, updates dual using that guess, corrects primal with new dual, then averages for estimates—like lookahead passes in a team relay stabilizing the handoff. Dual strength counters primal swings; the Lyapunov extends, yielding same optimal sublinear rates and first weak convergence for fully-split accelerated primal-dual, meaningful for denoising where slow standards drag.
Alex: Right—that relay analogy clicks; lookahead plus correction, bowl-stabilized, proves the points converge weakly too. Solid bridge to real speed in mixed problems.
Alex: So that relay stabilizes via the dual bowl. But walk me through how the Lyapunov inequality actually captures one step's progress—like, what pieces does it track to show everything shrinks?
Sam: The Lyapunov bundles three main distances: how far the primal trial point is from ideal, the dual trial from ideal, and a coupling term measuring mismatch—like a tug-of-war tension between teams. It also includes the Lagrangian gap at the averaged estimates, which gauges how well current guesses balance the full cost. Each iteration proves this bundle at the next step is smaller than before, minus positive progress terms—ensuring nonnegative under stepsize rules. When dual strong convexity holds—meaning the dual cost curves upward reliably—the extra push shrinks dual error firmly.
Alex: Okay, so the gap extends the primal gap to dual balance, and the bundle contracts reliably in those three dual cases. Does smooth dual cost give the cleanest speedup?
Sam: Yes, when the dual cost is smooth—meaning its cost changes steadily without jumps—it has strong convexity like a deep bowl always pulling back. They define Lyapunov weighting the bundle by momentum parameters growing roughly with steps, proving it shrinks by a factor under one. Choosing those parameters yields gaps shrinking about twice as fast asymptotically as standard methods. If primal also strongly convex, it locks into linear convergence after warmup.
Alex: Huh, so smooth dual bowls steadily for that speedup, linear if primal bowls too. What about the bounded coupling case—does it work without smooth dual?
Sam: There, bounded coupling below means dual error translates to primal via a minimum stretch—like no slack in the link, forcing strong pull-back. Lyapunov adjusts weights, proving similar contraction. Yields same rates on a refined gap, first for fully-split accelerated primal-dual. Linear if primal strongly convex too.
Alex: Right, coupling strength or smoothness bowls the dual either way, unlocking acceleration without divergence. Meaningful for denoising, where dual enforces edges sharply but smoothly.
Sam: Precisely. The unified Lyapunov proves single-loop APAPC preserves optimal rates across cases, with weak convergence to saddle-point—points cluster at solution—extending recent primal results. Practical edge over lumped smooth methods: decoupled stepsizes boost constants and robustness.
Alex: One last bit—beyond gaps, do the points themselves converge, not just averages?
Sam: In sublinear regime with dual strong but primal flat, they prove weak point convergence of primal iterates to a saddle under bounded trials—like the path clusters at solution without circling forever. Strong dual pulls dual home; full strong convergence if both sides bowl-shaped via linear rates. Ties the analysis complete for fully-split acceleration. The paper establishes optimal rates on the Lagrangian gap and dual error, with linear convergence when both sides have bowl shapes, across smooth dual, bounded coupling, or linear constraints. It also proves weak convergence of primal iterates to the minimizer—first for fully-split accelerated primal-dual. For denoising, this means fewer steps to balance smooth fitting and sharp penalties. The analysis assumes a structured primal nonsmooth part—like a quadratic bowl. General nonsmooth needs more complex splitting, left as open work.
Alex: Makes sense—a meaningful step for those cases, with clear paths ahead. Thanks, Sam, for breaking down how they stabilized acceleration without divergence.
Sam: Pleasure, Alex. This work by Condat, Sadiev, and Richtárik pushes primal-dual methods toward reliable speed. Thanks for listening to ResearchPod.