Table of Contents
Fetching ...

Nonconvex Stochastic Bregman Proximal Gradient Method with Application to Deep Learning

Kuangyu Ding, Jingyang Li, Kim-Chuan Toh

TL;DR

This work develops SBPG, a family of nonconvex stochastic optimization methods that replace quadratic gradient approximations with Bregman proximal updates to handle non-Lipschitz gradients. The authors prove convergence to stationary points with optimal $O(\epsilon^{-4})$ sample complexity and introduce MSBPG, which incorporates momentum to relax mini-batch requirements while preserving the same complexity. They apply MSBPG to deep neural network training using a polynomial kernel to ensure smooth adaptivity, showing robustness to stepsize and initialization and competitive generalization against SGD/Adam/AdamW across CNNs and LSTMs. The approach offers practical stability against gradient explosion and demonstrates potential as a universal open-source optimizer for large-scale nonconvex problems. Overall, SBPG/MSBPG provide a theoretically grounded, scalable alternative to traditional stochastic gradient methods in settings where Lipschitz smoothness is violated.

Abstract

Stochastic gradient methods for minimizing nonconvex composite objective functions typically rely on the Lipschitz smoothness of the differentiable part, but this assumption fails in many important problem classes like quadratic inverse problems and neural network training, leading to instability of the algorithms in both theory and practice. To address this, we propose a family of stochastic Bregman proximal gradient (SBPG) methods that only require smooth adaptivity. SBPG replaces the quadratic approximation in SGD with a Bregman proximity measure, offering a better approximation model that handles non-Lipschitz gradients in nonconvex objectives. We establish the convergence properties of vanilla SBPG and show it achieves optimal sample complexity in the nonconvex setting. Experimental results on quadratic inverse problems demonstrate SBPG's robustness in terms of stepsize selection and sensitivity to the initial point. Furthermore, we introduce a momentum-based variant, MSBPG, which enhances convergence by relaxing the mini-batch size requirement while preserving the optimal oracle complexity. We apply MSBPG to the training of deep neural networks, utilizing a polynomial kernel function to ensure smooth adaptivity of the loss function. Experimental results on benchmark datasets confirm the effectiveness and robustness of MSBPG in training neural networks. Given its negligible additional computational cost compared to SGD in large-scale optimization, MSBPG shows promise as a universal open-source optimizer for future applications.

Nonconvex Stochastic Bregman Proximal Gradient Method with Application to Deep Learning

TL;DR

This work develops SBPG, a family of nonconvex stochastic optimization methods that replace quadratic gradient approximations with Bregman proximal updates to handle non-Lipschitz gradients. The authors prove convergence to stationary points with optimal sample complexity and introduce MSBPG, which incorporates momentum to relax mini-batch requirements while preserving the same complexity. They apply MSBPG to deep neural network training using a polynomial kernel to ensure smooth adaptivity, showing robustness to stepsize and initialization and competitive generalization against SGD/Adam/AdamW across CNNs and LSTMs. The approach offers practical stability against gradient explosion and demonstrates potential as a universal open-source optimizer for large-scale nonconvex problems. Overall, SBPG/MSBPG provide a theoretically grounded, scalable alternative to traditional stochastic gradient methods in settings where Lipschitz smoothness is violated.

Abstract

Stochastic gradient methods for minimizing nonconvex composite objective functions typically rely on the Lipschitz smoothness of the differentiable part, but this assumption fails in many important problem classes like quadratic inverse problems and neural network training, leading to instability of the algorithms in both theory and practice. To address this, we propose a family of stochastic Bregman proximal gradient (SBPG) methods that only require smooth adaptivity. SBPG replaces the quadratic approximation in SGD with a Bregman proximity measure, offering a better approximation model that handles non-Lipschitz gradients in nonconvex objectives. We establish the convergence properties of vanilla SBPG and show it achieves optimal sample complexity in the nonconvex setting. Experimental results on quadratic inverse problems demonstrate SBPG's robustness in terms of stepsize selection and sensitivity to the initial point. Furthermore, we introduce a momentum-based variant, MSBPG, which enhances convergence by relaxing the mini-batch size requirement while preserving the optimal oracle complexity. We apply MSBPG to the training of deep neural networks, utilizing a polynomial kernel function to ensure smooth adaptivity of the loss function. Experimental results on benchmark datasets confirm the effectiveness and robustness of MSBPG in training neural networks. Given its negligible additional computational cost compared to SGD in large-scale optimization, MSBPG shows promise as a universal open-source optimizer for future applications.
Paper Structure (19 sections, 21 theorems, 98 equations, 12 figures, 2 tables, 2 algorithms)

This paper contains 19 sections, 21 theorems, 98 equations, 12 figures, 2 tables, 2 algorithms.

Key Result

Lemma 1

robbins1951stochastic Let $\left\{y_k\right\},\left\{u_k\right\},\left\{a_k\right\}$ and $\left\{b_k\right\}$ be non-negative adapted processes with respect to the filtration $\left\{\mathcal{F}_k\right\}$ such that $\sum_{k=0}^\infty a_k<\infty, \sum_{k=0}^\infty b_k<\infty$, and for all $k$, $\mat

Figures (12)

  • Figure 1: For function $F(x)=x^4$, which does not admit a globally Lipschitz continuous gradient. We restrict the feasible set to $[-0.5,2]$. Consider the models \ref{['convention-model']} and \ref{['model-Bregman']} of $F$ at ${x^k}=1$. The Lipschitz smooth constant of $F$ with respect to the kernel $\phi(x)=\frac{1}{2}x^2$ is 48. The smooth adaptivity constant of $F$ with respect to the kernel $\phi(x)=\frac{1}{2}x^2+\frac{1}{4}x^4$ is 4. The figure in (b) is a zoomed-in version of the plot in (a) for the range $[0.6,1]$. The unique minimum of $F(x)$ is at $x=0$.
  • Figure 2: Comparison of SBPG and SPG in terms of their robustness with respect to stepsize and initial point selction. A method is considered non-convergent if it fails to reach an accuracy of $\epsilon_k<10^{-2}$ within 30 seconds or if it collapses. Generally, choosing large stepsize and large radius for the initial point can cause an algorithm to collapse. The safe stepsize threshold is the maximum stepsize (constant schedule) that a method does not collapse. We run 10 tests for each algorithm and report the median of the results.
  • Figure 3: Training loss and test accuracy (%) of VGG16 on CIFAR10 dataset under two frequently used training settings. Here the activation function of VGG16 adopts smoothed ReLU activation function $\sigma_\epsilon$ with different choices of $\epsilon$ ($\epsilon=0$ denotes adopting the original ReLU activation function).
  • Figure 4: Training loss and test accuracy (%) of CNNs on CIFAR10 dataset with learning rate reduced to 0.1 times of the original value at the 150th epoch.
  • Figure 5: Training loss and test accuracy (%) of CNNs on CIFAR100 dataset with learning rate reduced to 0.1 times of the original value at the 150th epoch.
  • ...and 7 more figures

Theorems & Definitions (34)

  • Lemma 1
  • Definition 2
  • Definition 3
  • Lemma 4
  • Lemma 5
  • Definition 6
  • Definition 7
  • Lemma 8
  • Proposition 9
  • Definition 10: Bregman Gradient Mapping
  • ...and 24 more