A Short and Unified Convergence Analysis of the SAG, SAGA, and IAG Algorithms
Feng Zhu, Robert W. Heath, Aritra Mitra
TL;DR
The paper delivers a unified, high-probability convergence analysis for SAG, SAGA, and IAG by bounding gradient-staleness due to sub-sampling with Bernstein concentration and by crafting a novel Lyapunov function that accounts for delayed gradient information. This two-step approach yields linear convergence guarantees that hold across stochastic and deterministic variants, and extends to non-convex objectives and Markov sampling. An immediate byproduct is substantially tighter convergence rates for the deterministic IAG method, bringing its performance closer to stochastic VR methods. The framework is modular, simple, and adaptable, offering a foundation for tail bounds in more advanced VR algorithms and settings. Practically, this work clarifies the dynamics of variance-reduced methods and broadens their applicability in large-scale, real-world optimization problems.
Abstract
Stochastic variance-reduced algorithms such as Stochastic Average Gradient (SAG) and SAGA, and their deterministic counterparts like the Incremental Aggregated Gradient (IAG) method, have been extensively studied in large-scale machine learning. Despite their popularity, existing analyses for these algorithms are disparate, relying on different proof techniques tailored to each method. Furthermore, the original proof of SAG is known to be notoriously involved, requiring computer-aided analysis. Focusing on finite-sum optimization with smooth and strongly convex objective functions, our main contribution is to develop a single unified convergence analysis that applies to all three algorithms: SAG, SAGA, and IAG. Our analysis features two key steps: (i) establishing a bound on delays due to stochastic sub-sampling using simple concentration tools, and (ii) carefully designing a novel Lyapunov function that accounts for such delays. The resulting proof is short and modular, providing the first high-probability bounds for SAG and SAGA that can be seamlessly extended to non-convex objectives and Markov sampling. As an immediate byproduct of our new analysis technique, we obtain the best known rates for the IAG algorithm, significantly improving upon prior bounds.
