In this paper, we propose a quasi-Newton method for solving smooth and monotone nonlinear equations, including unconstrained minimization and minimax optimization as special cases. For the strongly monotone setting, we establish two global convergence bounds: (i) a linear convergence rate that matches the rate of the celebrated extragradient method, and (ii) an explicit global superlinear convergence rate that provably surpasses the linear convergence rate after at most
$$\mathcal {O}(d)$$
O
(
d
)
iterations, where
d
is the problem’s dimension. In addition, for the case where the operator is only monotone, we prove a global convergence rate of
$$\mathcal {O}(\min \{\frac{1}{k},\frac{\sqrt{d}}{k^{1.25}}\})$$
O
(
min
{
1
k
,
d
k
1.25
}
)
in terms of the duality gap. This matches the rate of the extragradient method when
$$k = \mathcal {O}(d^2)$$
k
=
O
(
d
2
)
and is faster when
$$k = \varOmega (d^2)$$
k
=
Ω
(
d
2
)
. These results are the first global convergence results to demonstrate a provable advantage of a quasi-Newton method over the extragradient method, without querying the Jacobian of the operator. Unlike classical quasi-Newton methods, we achieve this by using the hybrid proximal extragradient framework and a novel online learning approach for updating the Jacobian approximation matrices. Specifically, guided by the convergence analysis, we formulate the Jacobian approximation update as an online convex optimization problem over non-symmetric matrices, relating the regret of the online problem to the convergence rate of our method. To facilitate efficient implementation, we further develop a tailored online learning algorithm based on an approximate separation oracle, which preserves structures such as symmetry and sparsity in the Jacobian matrices.
We study stochastic simple bilevel optimization with smooth, possibly nonconvex upper- and lower-level objectives accessed only through stochastic gradient oracles. A key challenge is that the dual multiplier induced by the lower-level constraint may become unbounded near lower-level stationary points, invalidating bounded-dual analyses and destabilizing stochastic gradient estimates. To address this, we propose \emph{Stochastic Dynamic Barrier Perturbed Gradient} (SDBPG), a single-loop method that adaptively perturbs the dual formulation to regularize this degeneracy. The perturbation stabilizes the multiplier and yields controlled bias and variance even near the lower-level stationarity region. Under a rare-visit assumption governed by a parameter $\delta \in (0, \tfrac{1}{2}]$, SDBPG finds an $(\epsilon, \epsilon)$-stationary point in $\mathcal{O}(\epsilon^{-1/\delta})$ iterations, with sample gradient complexities $\mathcal{O}(\epsilon^{-2/\delta})$ and $\mathcal{O}(\epsilon^{-3/\delta})$ for the upper- and lower-level objectives, where larger $\delta$ corresponds to rarer visits to the bad region describing the negative alignment between the two objectives when the lower-level gradient is small. We further develop PR-SDBPG, a penalty-regularized variant that eliminates the rare-visit assumption, and VR-PR-SDBPG, which improves the resulting sample complexities entirely through variance reduction. To our knowledge, these are the first explicit $(\epsilon_f,\epsilon_g)$-stationarity guarantees for stochastic nonconvex-nonconvex simple bilevel optimization.
Mohammad Mahdi Ahmadi, Jincheng Cao, Aryan Mokhtari et al.· 0 citations