Moreau Envelope
   HOME

TheInfoList



OR:

The Moreau envelope (or the Moreau-Yosida regularization) M_f of a proper
lower semi-continuous In mathematical analysis, semicontinuity (or semi-continuity) is a property of extended real-valued functions that is weaker than continuity. An extended real-valued function f is upper (respectively, lower) semicontinuous at a point x_0 if, r ...
convex function In mathematics, a real-valued function is called convex if the line segment between any two distinct points on the graph of a function, graph of the function lies above or on the graph between the two points. Equivalently, a function is conve ...
f is a smoothed version of f. It was proposed by Jean-Jacques Moreau in 1965. The Moreau envelope has important applications in
mathematical optimization Mathematical optimization (alternatively spelled ''optimisation'') or mathematical programming is the selection of a best element, with regard to some criteria, from some set of available alternatives. It is generally divided into two subfiel ...
: minimizing over M_f and minimizing over f are equivalent problems in the sense that the sets of minimizers of f and M_f are the same. However, first-order optimization algorithms can be directly applied to M_f, since f may be non-differentiable while M_f is always continuously differentiable. Indeed, many
proximal gradient method Proximal gradient methods are a generalized form of projection used to solve non-differentiable convex optimization problems. Many interesting problems can be formulated as convex optimization problems of the form \operatorname\limits_ \sum_^ ...
s can be interpreted as a
gradient descent Gradient descent is a method for unconstrained mathematical optimization. It is a first-order iterative algorithm for minimizing a differentiable multivariate function. The idea is to take repeated steps in the opposite direction of the gradi ...
method over M_f.


Definition

The Moreau envelope of a proper lower semi-continuous convex function f from a
Hilbert space In mathematics, a Hilbert space is a real number, real or complex number, complex inner product space that is also a complete metric space with respect to the metric induced by the inner product. It generalizes the notion of Euclidean space. The ...
\mathcal to (-\infty,+\infty] is defined as M_(v) = \inf_ \left(f(x) + \frac \, x - v\, _2^2\right). Given a parameter \lambda \in \mathbb, the Moreau envelope of \lambda f is also called as the Moreau envelope of f with parameter \lambda.


Properties

* The Moreau envelope can also be seen as the infimal convolution between f and (1/2)\, \cdot \, ^2_2. * The
proximal operator In mathematical optimization, the proximal operator is an Operator (mathematics), operator associated with a proper,An Extended real number line, (extended) real-valued function ''f'' on a Hilbert space is said to be ''proper'' if it is not identic ...
of a function is related to the gradient of the Moreau envelope by the following identity: \nabla M_(x) = \frac (x - \mathrm_(x)). By defining the sequence x_ = \mathrm_(x_k) and using the above identity, we can interpret the proximal operator as a gradient descent algorithm over the Moreau envelope. * Using Fenchel's duality theorem, one can derive the following dual formulation of the Moreau envelope: M_(v) = \max_ \left( \langle p, v \rangle - \frac \, p \, ^2 - f^*(p)\right), where f^* denotes the
convex conjugate In mathematics and mathematical optimization, the convex conjugate of a function is a generalization of the Legendre transformation which applies to non-convex functions. It is also known as Legendre–Fenchel transformation, Fenchel transformati ...
of f. Since the subdifferential of a proper, convex, lower semicontinuous function on a Hilbert space is inverse to the subdifferential of its convex conjugate, we can conclude that if p_0 \in \mathcal X is the maximizer of the above expression, then x_0 := v - \lambda p_0 is the minimizer in the primal formulation and vice versa. * By Hopf–Lax formula, the Moreau envelope is a viscosity solution to a
Hamilton–Jacobi equation In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mecha ...
.
Stanley Osher Stanley Osher (born April 24, 1942) is an American mathematician, known for his many contributions in shock capturing, level-set methods, and PDE-based methods in computer vision and image processing. Osher is a professor at the University of ...
and co-authors used this property and
Cole–Hopf transformation The Cole–Hopf transformation is a change of variables that allows to transform a special kind of parabolic partial differential equations (PDEs) with a quadratic nonlinearity into a linear heat equation. In particular, it provides an explicit fo ...
to derive an algorithm to compute approximations to the proximal operator of a function.


See also

*
Proximal operator In mathematical optimization, the proximal operator is an Operator (mathematics), operator associated with a proper,An Extended real number line, (extended) real-valued function ''f'' on a Hilbert space is said to be ''proper'' if it is not identic ...
*
Proximal gradient method Proximal gradient methods are a generalized form of projection used to solve non-differentiable convex optimization problems. Many interesting problems can be formulated as convex optimization problems of the form \operatorname\limits_ \sum_^ ...


References


External links

* {{SpringerEOM, title=Moreau envelope function, oldid=49956
A Hamilton–Jacobi-based Proximal Operator
a
YouTube YouTube is an American social media and online video sharing platform owned by Google. YouTube was founded on February 14, 2005, by Steve Chen, Chad Hurley, and Jawed Karim who were three former employees of PayPal. Headquartered in ...
video explaining an algorithm to approximate the proximal operator Mathematical optimization