Slide 20 — Proof Roadmap for Theorem 2.9
Proof SketchTheorem (Smooth KL Convergence).
Let $f \colon \mathbb{R}^n \to \mathbb{R}$ be $C^1$ with $L$-Lipschitz gradient, $\inf f > -\infty$, and suppose the sublevel set $\{x : f(x) \leq f(x^0)\}$ is bounded. Suppose further that $f$ satisfies the KL property at every critical point. Run gradient descent with step size $\alpha \leq 1/L$: $$x^{k+1} = x^k - \alpha\,\nabla f(x^k).$$ Then the sequence $\{x^k\}$ converges to a critical point of $f$.
Remark. The theorem above only guarantees convergence to a critical point, which could be a local minimum, local maximum, or saddle point. To rule out saddle points, we appeal to a separate result:
Theorem (Lee, Simchowitz, Jordan & Recht, SIAM J. Optim. 2017, Theorem 4.1). Let $f$ be $C^2$ with Lipschitz continuous gradient. Then the set of initializations from which gradient descent converges to a strict saddle point has Lebesgue measure zero.
Combined with the theorem above (see slide o2), gradient descent converges to a local minimizer for almost every initialization.
Closure under composition. If $f \colon \mathbb{R}^m \to \mathbb{R}$ and $g \colon \mathbb{R}^n \to \mathbb{R}^m$ are both semi-algebraic, then $f \circ g$ is semi-algebraic. This follows from the Tarski–Seidenberg theorem (axiom 3 in the o-minimal definition below; see also slide s7): the graph of $f \circ g$ is the projection of $\{(x, z, y) : z = g(x),\, y = f(z)\}$ onto the $(x, y)$ coordinates, and projections of semi-algebraic sets are semi-algebraic.
A full proof via the Curve Selection Lemma and Puiseux series is given in slide 7p of the main deck.
Neural networks: ReLU networks are piecewise polynomial (hence semi-algebraic), so the composition of a smooth loss with a ReLU network is semi-algebraic and satisfies KL. Networks with smooth activations (GELU, sigmoid, tanh) satisfy KL because they are log-exp definable (a larger class that also guarantees KL).
Definition (o-minimal structure). An o-minimal structure on $\mathbb{R}$ is a collection $\mathcal{S} = \{\mathcal{S}_n\}_{n \geq 1}$, where each $\mathcal{S}_n$ is a family of subsets of $\mathbb{R}^n$ (the definable sets), satisfying:
- Each $\mathcal{S}_n$ is a Boolean algebra: closed under finite union, intersection, and complement.
- Closed under Cartesian products: $A \in \mathcal{S}_m$, $B \in \mathcal{S}_n \Rightarrow A \times B \in \mathcal{S}_{m+n}$.
- Closed under projections: if $A \in \mathcal{S}_{n+1}$ then $\pi(A) \in \mathcal{S}_n$ (projection onto first $n$ coordinates). This is the Tarski–Seidenberg property.
- Contains all algebraic sets: $\{x \in \mathbb{R}^n : p(x) = 0\} \in \mathcal{S}_n$ for every polynomial $p$.
- O-minimality: every set in $\mathcal{S}_1$ is a finite union of points and open intervals.
Examples by o-minimal structure:
| Structure | Definition | Examples |
|---|---|---|
| $\mathbb{R}_{\mathrm{alg}}$ (semi-algebraic) | Graph described by finitely many polynomial equalities/inequalities | Polynomials; rational functions (off poles); $|x|^{p/q}$, $\|x\|_p$ ($p/q \in \mathbb{Q}_{>0}$); ReLU, leaky ReLU; indicator of polytope; MSE, $\ell_1$ loss |
| $\mathbb{R}_{\mathrm{an}}$ (globally subanalytic) | Restrictions of real-analytic functions to bounded domains | $\sin(x)$ on $[0, 2\pi]$; Bessel functions on bounded intervals; any real-analytic function restricted to a compact set. (Global $\sin(x)$ on $\mathbb{R}$ is not in this structure: it has infinitely many zeros.) |
| $\mathbb{R}_{\mathrm{an,exp}}$ (log-exp) | Extends $\mathbb{R}_{\mathrm{an}}$ with global $e^x$ and $\log x$; shown o-minimal by van den Dries & Miller (1994) | $e^x$, $\log x$; $x^r$ for any $r \in \mathbb{R}$ (via $x^r = e^{r \log x}$); sigmoid $\sigma(x) = 1/(1+e^{-x})$; $\tanh$; softmax; cross-entropy $-\log \sigma(x)$; GELU approximation $x\cdot\sigma(1.702x)$ |
Mixing structures (e.g., ReLU network + softmax loss). The o-minimal structures form a hierarchy $\mathbb{R}_{\mathrm{alg}} \subseteq \mathbb{R}_{\mathrm{an}} \subseteq \mathbb{R}_{\mathrm{an,exp}}$, so every semi-algebraic function is automatically definable in $\mathbb{R}_{\mathrm{an,exp}}$. Therefore, composing a log-exp definable function with a semi-algebraic function — e.g., a softmax or logistic loss (log-exp) applied to the output of a ReLU network (semi-algebraic) — yields a function definable in $\mathbb{R}_{\mathrm{an,exp}}$, hence KL. You do not need both components to belong to the same smallest structure; it suffices that they share a common o-minimal structure.
ML examples satisfying KL (see slide I1, slide I2, slide I3 for details):
| ML example | Structure | Why KL holds |
|---|---|---|
| MSE loss with smooth network (sigmoid, tanh, GELU, SiLU) — regression, diffusion model, autoencoder (I2) | $\mathbb{R}_{\mathrm{an,exp}}$ | Smooth activations are log-exp definable; MSE is semi-algebraic; composition stays in $\mathbb{R}_{\mathrm{an,exp}}$ |
| MSE or cross-entropy loss with ReLU network — deep nets, ResNets (I1, I2) | $\mathbb{R}_{\mathrm{an,exp}}$ | ReLU is piecewise polynomial (semi-algebraic); cross-entropy involves $\log$ (log-exp); composition in $\mathbb{R}_{\mathrm{an,exp}}$ |
| Matrix factorization / recommender systems (I2) | $\mathbb{R}_{\mathrm{alg}}$ | Loss $\sum_{(i,j)\in\Omega}(R_{ij} - u_i^T v_j)^2 + \lambda\|U\|_F^2$ is a polynomial in $(U,V)$; semi-algebraic |
| BERT / LLM training (GELU, softmax attention, LayerNorm, cross-entropy) (I1, I2) | $\mathbb{R}_{\mathrm{an,exp}}$ | GELU, softmax, LayerNorm ($\sqrt{\cdot}$ is semi-algebraic), cross-entropy all log-exp definable or semi-algebraic; all compositions stay in $\mathbb{R}_{\mathrm{an,exp}}$ |
| LASSO and sparse regularized models (I3) | $\mathbb{R}_{\mathrm{alg}}$ | $\frac{1}{2}\|Ax-b\|^2$ is polynomial; $\lambda\|x\|_1 = \lambda\sum_j|x_j|$ is semi-algebraic; sum is semi-algebraic |
| RL policy gradient (PPO clipped objective, SAC entropy regularization) (I1) | $\mathbb{R}_{\mathrm{alg}}$ / $\mathbb{R}_{\mathrm{an,exp}}$ | PPO's $\min$/$\mathrm{clip}$ operations are piecewise polynomial (semi-algebraic); SAC entropy term $-\log\pi_\theta$ is log-exp definable |
| GANs (I1) | $\mathbb{R}_{\mathrm{an,exp}}$ | Each player's objective (e.g., binary cross-entropy) is log-exp definable, hence KL as a function of its own parameters. Caveat: KL convergence theory applies to single-objective minimization; the min-max structure of GAN training is not directly covered by Theorem 2.9. |
Warm-up: complete proof for the smooth case (gradient descent).
We prove convergence in the simplest setting: $f$ is $C^1$ with $L$-Lipschitz gradient, and the algorithm is gradient descent $x^{k+1} = x^k - \alpha\nabla f(x^k)$ with $\alpha \leq 1/L$. This avoids subdifferentials entirely and makes the role of KL transparent.
Assumption. $f \colon \mathbb{R}^n \to \mathbb{R}$ is $C^1$, $\nabla f$ is $L$-Lipschitz, $\inf f > -\infty$, the sublevel set $\{x : f(x) \leq f(x^0)\}$ is bounded, and $f$ satisfies the KL property at every critical point.
Step 1: Sufficient decrease. The descent lemma (slide r1) gives, using $\|x^{k+1} - x^k\| = \alpha\|\nabla f(x^k)\|$:
Example: $f(x) = \tfrac{1}{2}x^TAx$, $A \succ 0$: $L = \lambda_{\max}(A)$.
Nonconvex example: $f(x) = (x^2-1)^2$: $\nabla^2 f(x) = 12x^2 - 4$, so $L = \sup_{|x|\leq R}|12x^2-4|$ on any bounded region.
Step 2: Preliminary consequences. Summing Step 1 over $k = 0, 1, \ldots$:
where $f^* = \lim_{k\to\infty} f(x^k)$ exists since $f(x^k)$ is nonincreasing and bounded below. In particular, $\|x^{k+1}-x^k\| \to 0$ and $\|\nabla f(x^k)\| \to 0$. Since all iterates lie in the bounded sublevel set $\{f \leq f(x^0)\}$, Bolzano-Weierstrass gives a subsequential cluster point: $x^{k_j} \to \tilde{x}$ for some subsequence. By continuity, $f(\tilde{x}) = f^*$ and $\nabla f(\tilde{x}) = 0$.
Step 3: The KL inequality at $\tilde{x}$. Since $\nabla f(\tilde{x}) = 0$, the KL property applies at $\tilde{x}$: there exist $r > 0$, $\eta > 0$, and a concave $C^1$ function $\varphi\colon[0,\eta)\to[0,\infty)$ with $\varphi(0)=0$, $\varphi'>0$, such that
Step 4: Conditional telescoping bound. Define $\Delta_k = \varphi(f(x^k)-f^*)$. We claim:
Proof. By concavity of $\varphi$ (tangent line lies above the graph):
Now substitute the two key inequalities:
- Sufficient decrease (Step 1): $f(x^k)-f(x^{k+1}) \geq \tfrac{1}{2\alpha}\|x^{k+1}-x^k\|^2$.
- KL inequality (Step 3): $\varphi'(f(x^k)-f^*) \geq \tfrac{1}{\|\nabla f(x^k)\|} = \tfrac{\alpha}{\|x^{k+1}-x^k\|}$.
Step 5: Iterates eventually stay in the KL neighborhood.
Claim. There exists an index $K_0$ such that $x^k \in B(\tilde{x},r)$ and $f^* < f(x^k) < f^*+\eta$ for all $k \geq K_0$.
Proof. Since $x^{k_j}\to\tilde{x}$ and $f(x^k)\to f^*$ and $\varphi$ is continuous with $\varphi(0)=0$, choose $k_j$ large enough that simultaneously:
- $\|x^{k_j}-\tilde{x}\| < r/2$, so $x^{k_j}\in B(\tilde{x},r)$;
- $f^* < f(x^{k_j}) < f^*+\eta$; and
- $2\,\varphi(f(x^{k_j})-f^*) < r/2$.
The Claim asserts two conditions simultaneously: spatial ($x^k\in B(\tilde{x},r)$) and $f$-value ($f^* < f(x^k) < f^*+\eta$). We dismiss the $f$-value exit first:
- Upper bound $f(x^k) < f^*+\eta$: $f(x^k)$ is nonincreasing (Step 1), so $f(x^k)\leq f(x^{k_j}) < f^*+\eta$ for all $k\geq k_j$. The upper bound is never violated.
- Lower bound $f(x^k) \geq f^*$: $f(x^k)$ is nonincreasing and has a subsequence $f(x^{k_j})\to f^*$; hence the full sequence converges to $f^*$ from above, so $f(x^k)\geq f^*$ for all $k$. The value $f(x^k) < f^*$ is impossible. If $f(x^K)=f^*$ for some $K\geq k_j$, then the descent condition gives $f^*+\tfrac{1}{2\alpha}\|x^{K+1}-x^K\|^2\leq f(x^K)=f^*$, forcing $x^{K+1}=x^K$. The sequence is stationary from $K$ on and trivially stays in $B(\tilde{x},r)$; the Claim holds immediately.
So the only way the Claim can fail is through the spatial condition. Suppose for contradiction that some iterate leaves $B(\tilde{x},r)$ after $k_j$. Let $K^*>k_j$ be the first such index. For $k_j \leq k < K^*$, $x^k\in B(\tilde{x},r)$ by definition of $K^*$, and $f^* < f(x^k) < f^*+\eta$ by the two points above. So Step 4 applies at each such $k$, and summing:
But then by the triangle inequality:
a contradiction. So $K^*$ does not exist, and the Claim holds with $K_0 = k_j$. $\square$
Step 6: Convergence and criticality. By the Claim, Step 4 applies for all $k \geq K_0$. Summing and letting $M\to\infty$:
The sequence $\{x^k\}$ has finite total path length, hence is Cauchy in $\mathbb{R}^n$, and converges to some $x^*$. Since $\nabla f$ is continuous and $\|\nabla f(x^k)\|\to 0$, we conclude $\nabla f(x^*)=0$. $\square$
Abstract counterexample. Take $\|x^{k+1}-x^k\| = 1/k$. Then $\sum 1/k^2 = \pi^2/6 < \infty$ but $\sum 1/k = +\infty$. In higher dimensions, such a sequence can spiral indefinitely around a circle of critical points.
Why GD on natural functions never produces $1/k$ steps. Consider $f(x) = \tfrac{1}{2p}x^{2p}$ ($p\geq 1$). GD near $x=0$ gives $|x^{k+1}-x^k|\sim k^{-1-1/(2p-2)}$:
| Function | KL exponent $\theta$ | Step decay | $\ell^1$ summable? |
|---|---|---|---|
| $x^2$ | $\theta=1/2$ | $k^{-2}$ | Yes |
| $x^4$ | $\theta=3/4$ | $k^{-3/2}$ | Yes |
| $x^{10}$ | $\theta=9/10$ | $k^{-9/8}$ | Yes |
| $x^{2p}$, $p\to\infty$ | $\theta\to 1$ | $k^{-1-\epsilon}$, $\epsilon\to 0$ | Yes (barely) |
| hypothetical $\theta=1$ | $\theta=1$ | $k^{-1}$ (harmonic) | No |
General case: proof of Theorem 2.9 (Attouch–Bolte–Svaiter 2013).
The smooth warm-up above is a special case of Theorem 2.9, which covers any algorithm satisfying H1–H3 with a proper l.s.c. KL objective. The full proof extends the warm-up by an index-shift trick (H2 at step $k-1$ bounds the subgradient at $x^k$) and a coupled recurrence ($\star$) that replaces the clean telescoping estimate. A self-contained proof with all six steps is given at slides_KL_nonsmooth.html.