# Proximal point methods in mathematical programming

The proximal point method for finding a zero of a maximal monotone operator generates a sequence , starting with any , whose iteration formula is given by (a1)

where and is a bounded sequence of positive real numbers. The origin of the method can be traced back to the study of regularization of ill-posed problems (see [a3]; Regularization method; Ill-posed problems) and is closely related to Moreau–Yoshida convolution (see [a5]). A survey on the method can be found in [a4]. The basic properties of the method, as established in [a7], are as follows.

i) The sequence is well defined (in the sense that has a unique zero), converges to a zero of if such a zero exists, and is unbounded otherwise.

ii) If is strongly monotone, then the convergence rate of is linear, and super-linear when . When is the subdifferential of a convex polyhedral function (i.e. the epigraph of is a polyhedron), then convergence is finite.

iii) Convergence properties are preserved under inexact computation of the zero of , in the following sense: (a1) is replaced by , where and .

When (the subdifferential of a convex function ), then (a1) is equivalent to (a2)

and the zeros of are the minimizers of , so that converges to a minimizer of . If is the subdifferential of the dual objective of a convex optimization problem, then the sequence coincides with the dual sequence generated by the augmented Lagrangian method (see [a6]).

All these properties hold also for maximal monotone operators in a Hilbert space, in which case convergence is understood in the sense of the weak topology.

The proximal point method can be applied to problems with convex constraints, e.g. the variational inequality problem , for a closed and convex set , which consists of finding a such that there exists an satisfying for all . When with convex , reduces to minimizing subject to . For applying the proximal point method to it suffices to replace by , where is the normal cone operator of . In such a case, (a1) becomes equivalent to stating that solves and (a2) becomes There is another way to deal with constrained problems in the context of the proximal point method, without introducing the normal cone of , for the case when has non-empty interior and when a boundary-coercive Bregman function with zone is available. In this case the method starts with (the interior of ) and the iterative formula is still (a1) but now . In the optimization case, i.e. , this formula is equivalent to where is the Bregman distance associated with . This generalization of the proximal point method is an interior point algorithm, in the sense that the sequence lies in the interior of (cf. also Interior-point methods in mathematical programming). Under certain technical hypotheses on , satisfied, e.g., when , property i) above still holds, i.e. the sequence converges to a solution of when a solution exists and is unbounded otherwise, [a8]. A property of this generalization, not shared by the original proximal point method (a1), is that for certain problems, including linear programming and quadratic programming, the limit of the sequence is the solution of the problem which minimizes the Bregman distance to the initial iterate (see [a2], [a10]). This procedure also generates augmented Lagrangian-type methods when is the subdifferential of the dual objective of a convex optimization problem. For instance, when is the non-negative orthant of and , the sequence coincides with the dual sequence of the Bertsekas exponential multipliers method (see [a1], [a11]).

Finally, some results are available for the proximal point method with Bregman distances under inexact solutions of the subproblems, in the spirit of property iii) above, [a9]. Convergence is guaranteed for the sequence generated by , with , where and, for any monotone operator and any , is defined as  