Convergence Rate of the Last Iterate of Stochastic Proximal Algorithms
Kevin Kurian Thomas Vaidyan, Michael P. Friedlander, Ahmet Alacaoglu
TL;DR
This work addresses the lack of theoretical guarantees for last-iterate convergence in proximal stochastic optimization with unbounded gradient variance. By leveraging co-coercivity and a novel last-iterate reduction, it proves a $\widetilde{O}(1/\sqrt{T})$ last-iterate rate for proximal SGD under componentwise smoothness and convexity, and extends the result to randomized incremental proximal methods with a Lipschitz assumption on the component regularizers. It also provides corollaries for projected SGD and stochastic proximal point, and demonstrates the practical viability via BlockProx extensions and numerical experiments showing the last iterate often outperforms averaging. The results apply to graph-guided regularizers common in multi-task and federated learning, offering nonasymptotic, data-independent rates without requiring bounded variance. Overall, the paper advances the understanding of last-iterate behavior in regularized stochastic optimization and broadens the toolkit for large-scale, graph-structured learning tasks.
Abstract
We analyze two classical algorithms for solving additively composite convex optimization problems where the objective is the sum of a smooth term and a nonsmooth regularizer: proximal stochastic gradient method for a single regularizer; and the randomized incremental proximal method, which uses the proximal operator of a randomly selected function when the regularizer is given as the sum of many nonsmooth functions. We focus on relaxing the bounded variance assumption that is common, yet stringent, for getting last iterate convergence rates. We prove the $\widetilde{O}(1/\sqrt{T})$ rate of convergence for the last iterate of both algorithms under componentwise convexity and smoothness, which is optimal up to log terms. Our results apply directly to graph-guided regularizers that arise in multi-task and federated learning, where the regularizer decomposes as a sum over edges of a collaboration graph.
