Empirical risk minimization on massive datasets naturally exhibits a nested double finite-sum structure, where
N=nm total samples are logically or physically partitioned into
n blocks of size
m (e.g., in pooled data silos, out-of-core learning, or deliberate stratification). While variance-reduced methods achieve optimal oracle complexities for nonconvex objectives, they suffer from severe scaling bottlenecks in this centralized regime. Recursive estimators, such as PAGE, require periodic global full-gradient refreshes over all
nm samples, which are computationally expensive. Conversely, single-loop methods, such as SILVER, avoid such refreshes but require an impractical
O(nm) memory footprint to store a control variate for every sample. In this paper, we propose SILAGE, a variance-reduced algorithm that addresses this trade-off. By actively exploiting the double-sum structure, SILAGE eliminates periodic global full-gradient refreshes over all
nm components (evaluating at most one local group gradient per iteration) while requiring only
O(n) memory. Furthermore, we provide a tight convergence analysis that avoids pessimistic worst-case Lipschitz constants. Instead, SILAGE's complexity natively adapts to the underlying data geometry via nested functional similarities: across-group (
δ1) and within-group (
δ2) heterogeneity. Our results improve existing state-of-the-art bounds in several practically relevant regimes.