A framework for bilevel optimization that enables stochastic and global variance reduction algorithms
Mathieu Dagréou, Pierre Ablin, Samuel Vaiter, Thomas Moreau
TL;DR
The paper tackles stochastic bilevel optimization by introducing a joint dynamic framework where the inner minimizer and a linear-system solution evolve together with the outer variable. It defines unbiased directions $D_z$, $D_v$, and $D_x$ as simple sample-based sums, enabling single-loop updates and easy integration of variance-reduction techniques. The authors present SOBA, a single-sample SGD-like method, and SABA, a SAGA-like variance-reduction algorithm, and provide convergence guarantees: SOBA achieves $O(T^{-1/2})$ (or $O(rac{ ext{log}T}{\sqrt{T}})$ with decreasing steps) for the gradient norm, while SABA attains $O(N^{2/3}T^{-1})$ in general, and linear convergence under a PL inequality. Empirical results on hyperparameter selection and data cleaning show that SABA often outperforms existing stochastic bilevel methods, highlighting the practical impact of global variance reduction in bilevel settings.
Abstract
Bilevel optimization, the problem of minimizing a value function which involves the arg-minimum of another function, appears in many areas of machine learning. In a large scale empirical risk minimization setting where the number of samples is huge, it is crucial to develop stochastic methods, which only use a few samples at a time to progress. However, computing the gradient of the value function involves solving a linear system, which makes it difficult to derive unbiased stochastic estimates. To overcome this problem we introduce a novel framework, in which the solution of the inner problem, the solution of the linear system, and the main variable evolve at the same time. These directions are written as a sum, making it straightforward to derive unbiased estimates. The simplicity of our approach allows us to develop global variance reduction algorithms, where the dynamics of all variables is subject to variance reduction. We demonstrate that SABA, an adaptation of the celebrated SAGA algorithm in our framework, has $O(\frac1T)$ convergence rate, and that it achieves linear convergence under Polyak-Lojasciewicz assumption. This is the first stochastic algorithm for bilevel optimization that verifies either of these properties. Numerical experiments validate the usefulness of our method.
