A Generalized Forward-Backward Splitting

This paper introduces a generalized forward-backward splitting algorithm for finding a zero of a sum of maximal monotone operators B + Sigma(n)(i=1) A(i), where B is cocoercive. It involves the computation of B in an explicit (forward) step and the parallel computation of the resolvents of the A(i)'s in a subsequent implicit (backward) step. We prove the algorithm's convergence in infinite dimension and its robustness to summable errors on the computed operators in the explicit and implicit steps. In particular, this allows efficient minimization of the sum of convex functions f + Sigma(n)(i=1) g(i), where f has a Lipschitz-continuous gradient and each g(i) is simple in the sense that its proximity operator is easy to compute. The resulting method makes use of the regularity of f in the forward step, and the proximity operators of the g(i) is are applied in parallel in the backward step. While the forward-backward algorithm cannot deal with more than n = 1 nonsmooth function, we generalize it to the case of arbitrary n. Examples on inverse problems in imaging demonstrate the advantage of the proposed methods in comparison to other splitting algorithms.