Primal-Dual Coarse Correction Optimizer
Implementation & benchmark of arXiv:2608.04848 — Primal-dual multigrid methods for nonsmooth optimization
Source paper: Primal-dual multigrid methods for nonsmooth optimization arXiv:2608.04848 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Add a periodic coarse optimization phase to SGD or Adam that operates on a compressed parameterization and returns a prolongated correction to the full network. Retain nonsmooth constraints or regularizers explicitly through a primal-dual update instead of relying on penalty smoothing. Accept the correction only when it improves a cheap fine-batch merit test, making the method useful even when the coarse objective is only approximately coherent with the fine objective.
Formulas
Mathematical statement
The paper studies composite problems of the form \(\min_{x\in X} F(x)+E(x)+G(Kx)\), where \(F\) and \(G\) are convex, possibly nonsmooth functions, \(E\) is convex and smooth with an \(L\)-Lipschitz gradient, and \(K\in\mathbb{L}(X;Y)\) maps the primal Hilbert space \(X\) to a dual space \(Y\). A coarse primal-dual iterate \(v^{k,m}\) is initialized from the restricted fine pre-iterate \(v^{k,0}:=I_h^H u_+^k\), where \(u_+^k=(x_+^k,y_+^k)\), \(I_h^H\) restricts fine primal and dual variables to the coarse space, and \(I_H^h\) prolongs them back. The resulting correction is \(d^k:=I_H^h(v^{k,m}-v^{k,0})\). We transfer this construction to neural parameters \(\theta\): \(F\) is a nonsmooth constraint or regularizer, \(E\) is the minibatch loss, \(G(K\theta)\) is a structured penalty, \(R\) and \(P\) replace \(I_h^H\) and \(I_H^h\), and \(z\) is a coarse parameter vector. The coarse optimizer uses \(F_H(z)+E_H(z)+G_H(K_Hz)\), with \(K_H:=KP\) and typically \(E_H(z)=E(Pz)\), while the fine correction is \(d_\theta=P(z_m-R\theta)\). The transferable property is that the coarse correction preserves the same composite primal-dual structure, so nonsmooth feasibility and dual regularization survive across levels.
Implementation notes
Integrate this as an optimizer wrapper around a selected parameter block, not as a new forward-network layer. Start with a CNN or small transformer containing a large matrix \(\theta\in\mathbb{R}^n\). Choose a fixed coarse basis \(P\in\mathbb{R}^{n\times r}\), with \(r\ll n\), using block-constant parameters, a 2D Haar/DCT basis for convolution kernels, or rank-r factors for a linear layer. Set the restriction to \(R=P^\top\) after column normalization. Every \(T\) fine optimizer steps, set \(z_0=R\theta\) and run \(m=3\)--\(10\) coarse updates on a fresh minibatch. For each coarse iteration, compute \(q_j=K_Hz_j\); update the dual variable by \(y_{j+1}=\operatorname{prox}_{\sigma G_H^*}(y_j+\sigma q_j)\); then update the primal variable by \(z_{j+1}=\operatorname{prox}_{\tau F_H}(z_j-\tau[\nabla E_H(z_j)+K_H^\top y_{j+1}])\). Use \(\tau\sigma\|K_H\|^2<0.9\), estimating \(\|K_H\|\) with 5--10 power iterations. Form \(d=P(z_m-z_0)\), evaluate the fine loss at \(\theta+d\) on a small validation minibatch, and accept with backtracking \(\alpha\in\{1,1/2,1/4\}\); otherwise skip the correction. For a first nonsmooth test, use \(F_H=\iota_{[-c,c]^r}\), whose proximal map is coordinatewise clipping, or \(G_H(q)=\lambda\|q\|_1\), whose proximal map is soft-thresholding. The paper supplies the composite objective and restriction/prolongation correction structure; minibatch gradients, operator norms, and acceptance statistics are estimated empirically. First experiment: train a 6-layer CNN on CIFAR-10, comparing AdamW with identical fine-step FLOPs against AdamW plus coarse corrections on convolution kernels. Measure training loss versus optimizer steps and wall-clock time, final accuracy, correction acceptance rate, and constraint violations. A successful signal is 10--20% fewer fine steps to reach baseline accuracy without increased wall-clock time, or improved stability at the same learning rate.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.