Asymptotic Analysis of Stochastic Splitting Methods for Multivariate Monotone Inclusions
Patrick L. Combettes, Javier I. Madariaga
TL;DR
This work develops a general stochastic framework for convergence of monotone operator splitting methods in Hilbert spaces, accommodating stochastic operator approximations, random coordinate/block updates, and random relaxations. It proves almost surely and in $L^2$ convergence for stochastic proximal point variants and randomized block-iterative saddle projective splitting applied to multivariate inclusions that combine set-valued, cocoercive, and Lipschitzian operators. By leveraging a two-operator abstract model and a robust probabilistic toolkit, the paper unifies and extends deterministic splitting theories to stochastic settings and provides convergence guarantees for complex multivariate minimization problems and randomized Kuhn–Tucker schemes. The results tolerate summable stochastic errors and even allow super-relaxation parameters, broadening applicability to large-scale stochastic optimization and domain-decomposed systems. Overall, the framework significantly enhances the toolkit for solving stochastic monotone inclusions with structured, high-dimensional block architectures.
Abstract
We propose an abstract framework to establish the convergence of the iterates of stochastic versions of a broad range of monotone operator splitting methods in Hilbert spaces. This framework allows for the introduction of stochasticity at several levels: approximation of operators, selection of coordinates and operators in block-iterative implementations, and relaxation parameters. The proposed analysis involves a reduced inclusion model with two operators. At each iteration, stochastic approximations to points in the graphs of these two operators are used to form the update. The results are applied to derive the almost sure and $L^2$ convergence of stochastic versions of the proximal point algorithm, as well as of randomized block-iterative projective splitting methods for solving systems of coupled inclusions involving a mix of set-valued, cocoercive, and Lipschitzian monotone operators combined via various monotonicity-preserving operations.
