Regularization methods for solving hierarchical variational inequalities with complexity guarantees
Daniel Cortild, Meggie Marschner, Mathias Staudigl
TL;DR
The paper addresses solving hierarchical variational inequalities where the lower level is defined by a monotone inclusion, by introducing a double-loop diagonal equilibrium tracking method (DANTE) that combines $T$ikhonov regularization with proximal penalties. The inner loop employs a flexible fixed-point encoding with an inertial Krasnoselskii–Mann iteration to approximate a temporal solution, while the outer loop restarts to update anchors and steer convergence. The authors establish convergence and rate guarantees for the averaged iterates in terms gap functions $\operatorname{Gap}_{\operatorname{opt}}$ and $\operatorname{Gap}_{\operatorname{feas}}$, with explicit complexity bounds and performance across several splitting schemes (two-operator and three-operator). Numerical experiments on equilibrium selection, least-norm least-squares, and image inpainting validate the method’s versatility and show competitive convergence behavior across different problem structures. Overall, the work provides a broad, theory-backed framework for solving hierarchical VI problems with concrete rates and practical splitting strategies that can accommodate infinite-dimensional Hilbert spaces and a wide range of operator structures.
Abstract
We consider hierarchical variational inequality problems, or more generally, variational inequalities defined over the set of zeros of a monotone operator. This framework includes convex optimization over equilibrium constraints and equilibrium selection problems. In a real Hilbert space setting, we combine a Tikhonov regularization and a proximal penalization to develop a flexible double-loop method for which we prove asymptotic convergence and provide rate statements in terms of gap functions. Our method is flexible, and effectively accommodates a large class of structured operator splitting formulations for which fixed-point encodings are available. Finally, we validate our findings numerically on various examples.
