Slide 20 — Proof of KL Convergence (Smooth Case)

Extracted from the full slide deck. All cross-references link back to the main slides.

Slide 20 — Proof Roadmap for Theorem 2.9

Proof Sketch

Theorem (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.

Recall: what is the KL property? (See slide s5 for the full definition.) A function $f$ satisfies the Kurdyka-Łojasiewicz (KL) property at a critical point $x^*$ if there exist a neighborhood $U$ of $x^*$, a constant $\eta > 0$, and a continuous concave function $\varphi \colon [0, \eta) \to [0, \infty)$ with $\varphi(0) = 0$, $\varphi' > 0$ on $(0, \eta)$, such that for all $x \in U$ with $f(x^*) < f(x) < f(x^*) + \eta$: $$\varphi'\bigl(f(x) - f(x^*)\bigr) \cdot \|\nabla f(x)\| \;\geq\; 1.$$ In words: near a critical point, the gradient cannot be too small relative to the function value gap, after reparameterization by $\varphi$. The function $\varphi$ is called the desingularizing function (slide s4). When $\varphi(s) = c\,s^{1-\theta}$, the exponent $\theta \in [0,1)$ controls the convergence rate (slide 22).
When does KL hold? By the theorem of Bolte-Sabach-Teboulle (slide s7), every proper l.s.c. semi-algebraic function satisfies KL.

Definition (semi-algebraic function). A function $f \colon \mathbb{R}^n \to \mathbb{R}$ is semi-algebraic if its graph $\{(x,y) : y = f(x)\}$ can be described by finitely many polynomial equalities and inequalities.

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.
Theorem 3 (Bolte–Sabach–Teboulle, Math. Programming 2014). Every proper l.s.c. semi-algebraic function $f \colon \mathbb{R}^n \to \mathbb{R} \cup \{+\infty\}$ is a KL function: it satisfies the KL property at every point of $\mathrm{dom}\,\partial f$.
A full proof via the Curve Selection Lemma and Puiseux series is given in slide 7p of the main deck.
Example: All polynomial functions are semi-algebraic, so KL holds. This includes $f(x) = (x^2-1)^2$, $f(x) = x^4$, and the MSE loss $f(\theta) = \frac{1}{n}\sum_i \|\Phi(a_i;\theta) - b_i\|^2$ when $\Phi$ is a polynomial-activation network.
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:
  1. Each $\mathcal{S}_n$ is a Boolean algebra: closed under finite union, intersection, and complement.
  2. Closed under Cartesian products: $A \in \mathcal{S}_m$, $B \in \mathcal{S}_n \Rightarrow A \times B \in \mathcal{S}_{m+n}$.
  3. 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.
  4. Contains all algebraic sets: $\{x \in \mathbb{R}^n : p(x) = 0\} \in \mathcal{S}_n$ for every polynomial $p$.
  5. O-minimality: every set in $\mathcal{S}_1$ is a finite union of points and open intervals.
A function $f \colon \mathbb{R}^n \to \mathbb{R}$ is definable in $\mathcal{S}$ if its graph $\{(x, f(x))\} \in \mathcal{S}_{n+1}$. The o-minimality axiom (5) is the key tameness condition: it forbids infinite oscillation (e.g., $\sin(1/x)$ near 0 is not definable in any o-minimal structure). Every structure satisfying (1)–(5) is automatically closed under composition, sum, product, and partial derivatives of smooth definable functions. By Kurdyka (1998) [ref 5, References], every function definable in an o-minimal structure satisfies the KL property.

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)$
Composition within o-minimal structures. If $f \colon \mathbb{R}^m \to \mathbb{R}$ and $g \colon \mathbb{R}^n \to \mathbb{R}^m$ are both definable in the same o-minimal structure (e.g., both semi-algebraic, or both log-exp definable), then $f \circ g$ is also definable in that structure and hence KL. This follows because o-minimal structures are closed under composition — a defining property of the framework; see Kurdyka (1998) [ref 5, References]. Calculus rules for the KL exponent of a composition (e.g., computing the exponent $\theta$ of $f \circ g$ from those of $f$ and $g$) are given in Li–Pong (2018) [ref 16, References].

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.
Caution: composition of abstract KL functions. If $f$ and $g$ satisfy the KL gradient inequality but are not assumed definable in any o-minimal structure, the composition $f \circ g$ need not be KL. The closure property above is specific to the o-minimal/semi-algebraic setting and does not hold for KL functions defined purely by the abstract inequality.

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)\|$:

$$f(x^{k+1}) \;\leq\; f(x^k) - \frac{1}{2\alpha}\|x^{k+1} - x^k\|^2.$$
When does the descent lemma hold? It requires $\|\nabla f(x) - \nabla f(y)\| \leq L\|x-y\|$ for all $x,y$, equivalently $f(y) \leq f(x) + \langle \nabla f(x), y-x\rangle + \frac{L}{2}\|y-x\|^2$, with step size $\alpha \leq 1/L$.
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$:

$$\sum_{k=0}^{\infty}\|x^{k+1} - x^k\|^2 \;\leq\; 2\alpha\bigl(f(x^0) - f^*\bigr) \;<\; \infty,$$

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$.

A cluster point is not the same as a limit. At this stage we only know some subsequence converges to $\tilde{x}$. The full sequence could drift between multiple cluster points — even with $\|x^{k+1}-x^k\|\to 0$, iterates can walk along a circle with shrinking steps and never settle. Steps 3–5 use the KL inequality to rule this out and prove the full sequence converges.
When is the sublevel set bounded? Whenever $f$ is coercive: $f(x)\to+\infty$ as $\|x\|\to\infty$. Any $f$ with a weight-decay term $+\tfrac{\lambda}{2}\|x\|^2$ is coercive. For $f(x)=(x^2-1)^2$: coercive since $f(x)\to\infty$ as $|x|\to\infty$.

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

$$\varphi'\bigl(f(x)-f^*\bigr)\cdot\|\nabla f(x)\| \;\geq\; 1 \qquad \text{for all } x\in B(\tilde{x},r) \text{ with } f^* < f(x) < f^*+\eta.$$
Why $f(x^k) > f^*$? Since $f(x^k)$ is nonincreasing with limit $f^*$, we have $f(x^k) \geq f^*$. If equality holds at some finite step $K$, then Step 1 gives $\|x^{K+1}-x^K\|^2 \leq 2\alpha(f(x^K)-f^*) = 0$, so the sequence is stationary from step $K$ onward — convergence is trivial. In the interesting case, $f(x^k) > f^*$ for all $k$.

Step 4: Conditional telescoping bound. Define $\Delta_k = \varphi(f(x^k)-f^*)$. We claim:

$$\text{if } x^k \in B(\tilde{x},r) \text{ and } f^*< f(x^k) < f^*+\eta, \quad\text{then}\quad \|x^{k+1}-x^k\| \;\leq\; 2(\Delta_k - \Delta_{k+1}).$$

Proof. By concavity of $\varphi$ (tangent line lies above the graph):

$$\Delta_k - \Delta_{k+1} \;\geq\; \varphi'(f(x^k)-f^*)\cdot\bigl(f(x^k)-f(x^{k+1})\bigr).$$
Why does concavity give this? For concave $\varphi$ and $s_1 \geq s_2 \geq 0$: $\varphi(s_2) \leq \varphi(s_1) + \varphi'(s_1)(s_2-s_1)$, so $\varphi(s_1)-\varphi(s_2) \geq \varphi'(s_1)(s_1-s_2)$. Set $s_1 = f(x^k)-f^*$, $s_2 = f(x^{k+1})-f^*$.

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\|}$.
$$\Delta_k - \Delta_{k+1} \;\geq\; \frac{\alpha}{\|x^{k+1}-x^k\|}\cdot\frac{1}{2\alpha}\|x^{k+1}-x^k\|^2 \;=\; \frac{1}{2}\|x^{k+1}-x^k\|. \quad\square$$

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:

$$\sum_{k=k_j}^{K^*-1}\|x^{k+1}-x^k\| \;\leq\; 2\Delta_{k_j} \;=\; 2\,\varphi(f(x^{k_j})-f^*) \;<\; r/2.$$

But then by the triangle inequality:

$$r \;\leq\; \|x^{K^*}-\tilde{x}\| \;\leq\; \underbrace{\|x^{K^*}-x^{k_j}\|}_{\leq\,\sum_{m=k_j}^{K^*-1}\|x^{m+1}-x^m\|\,<\,r/2} + \underbrace{\|x^{k_j}-\tilde{x}\|}_{<\,r/2} \;<\; r,$$

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$:

$$\sum_{k=K_0}^{\infty}\|x^{k+1}-x^k\| \;\leq\; 2\,\varphi(f(x^{K_0})-f^*) \;<\; \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$

Concrete example: $f(x) = (x^2-1)^2$. Critical points: $x=0$ (local max), $x=\pm 1$ (global minima). All assumptions hold: $\nabla f$ is Lipschitz on bounded sets, $f$ is coercive, and $f$ is semi-algebraic hence KL. Near $x=1$: $f\approx 4u^2$, $|f'|\approx 8|u|$ ($u=x-1$), so KL holds with $\theta=1/2$ (slide o2), giving linear convergence. The proof guarantees GD converges to $x=1$ or $x=-1$ for almost all initializations.
Where exactly did KL help? Step 2 already gives $\sum\|x^{k+1}-x^k\|^2<\infty$ ($\ell^2$ summability), but this does not imply $\sum\|x^{k+1}-x^k\|<\infty$ ($\ell^1$ summability), and without $\ell^1$ summability the sequence need not be Cauchy. The KL inequality in Step 4 supplies the missing factor: it bounds $\varphi'(f(x^k)-f^*)$ from below by $\alpha/\|x^{k+1}-x^k\|$, which cancels one power of $\|x^{k+1}-x^k\|$ from the squared term, upgrading $\ell^2$ to $\ell^1$.

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
KL with $\theta<1$ guarantees the step decay exponent stays strictly above $1$, keeping $\sum\|x^{k+1}-x^k\|$ convergent. The harmonic $1/k$ threshold is exactly the boundary KL prevents.

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.