The Moreau-Yosida Regularization: Difference between revisions

From Optimal Transport Wiki
Jump to navigation Jump to search
(added portmanteau)
Line 1: Line 1:
(to be filled in)
The '''Moreau-Yosida regularization''' is a technique used to approximate lower semicontinuous functions by Lipschitz functions. The main application of this result is to prove Portmanteau's Theorem, which states that integration against a lower semicontinuous and bounded below function is lower semicontinuous with respect to the narrow topology on the space of probability measures.
 
==Motivation==
(to be filled in)
 


==Definitions==
==Definitions==
Line 11: Line 7:


  <math>g_k(x) := \inf\limits_{y \in X} \left[ g(y) + k d(x,y) \right].</math>
  <math>g_k(x) := \inf\limits_{y \in X} \left[ g(y) + k d(x,y) \right].</math>
The distance term <math>d(x,y)</math> may often be raised to a positive exponent. For example, when <math>X</math> is a Hilbert space <ref name="BC"/> <ref name="AGS"/>, <math>g_k</math> is taken to be
<math>g_k(x) := \inf\limits_{y \in X} \left[ g(y) + k \| x - y \|^2 \right].</math>


Note that  
Note that  
Line 27: Line 27:
[[File:Ex 1.png|300px|thumb|Plot of <math>g(x) = x^2</math> and <math>g_k(x)</math> for <math>k = 0, 1, 2, 3</math>.]]
[[File:Ex 1.png|300px|thumb|Plot of <math>g(x) = x^2</math> and <math>g_k(x)</math> for <math>k = 0, 1, 2, 3</math>.]]


==Results==
==Approximating Lower Semicontinuous Functions by Lipschitz Functions==
'''Proposition.''' <ref name="OT"/><ref name="S"/>   
'''Proposition.''' <ref name="OT"/><ref name="S"/>   
* If <math>g</math> is proper and bounded below, so is <math>g_k</math>. Furthermore, <math>g_k</math> is continuous for all <math>k \geq 0</math>.  
* If <math>g</math> is proper and bounded below, so is <math>g_k</math>. Furthermore, <math>g_k</math> is continuous for all <math>k \geq 0</math>.  
Line 53: Line 53:


* By definition, <math>g_k \wedge k \in C_b(X)</math>. Since <math>g_k(x) \nearrow g(x)</math> for all <math>x \in X</math>, <math>g_k(x) \wedge k \nearrow g(x)</math> for all <math>x \in X</math>.
* By definition, <math>g_k \wedge k \in C_b(X)</math>. Since <math>g_k(x) \nearrow g(x)</math> for all <math>x \in X</math>, <math>g_k(x) \wedge k \nearrow g(x)</math> for all <math>x \in X</math>.
==Portmanteau Theorem==
'''Theorem (Portmanteau).''' Let <math>g : X \to (-\infty,+\infty]</math> be lower semicontinuous and bounded below. Then the functional <math>\mu \mapsto \int_X g \, \mathrm{d}\mu</math> is lower semicontinuous with respect to narrow convergence in <math>\mathcal{P}(X)</math>, that is
:<math> \mu_n \to \mu \quad \text{narrowly} \Longrightarrow \liminf\limits_{n \to \infty} \int_X g_n \, \mathrm{d}\mu \geq \int_X g \, \mathrm{d}\mu </math>.
'''Proof.''' By the Moreau-Yosida approximation, for all <math>k \geq 0</math>,
:<math>\liminf\limits_{n \to \infty} \int_X g \, \mathrm{d} \mu_n \geq \liminf\limits_{n \to \infty} \int_X g_k \wedge k \, \mathrm{d}\mu_n  = \int_X g_k \wedge k \, \mathrm{d}\mu </math>.
Taking <math>k \to \infty</math>, Fatou's lemma ensures that
:<math>\liminf\limits_{n \to \infty} \int_X g \, \mathrm{d} \mu_n \geq \liminf\limits_{k \to \infty} \int_X g_k \wedge k \, \mathrm{d} \geq \int_X g \, \mathrm{d}\mu \mu</math>.
==The Mysterious Etymology of Portmanteau==
(spurious stuff, will fill in later)


==References==
==References==
Line 58: Line 74:
<!--Bauschke-Combette Ch 12.<ref name="BC" />; Santambrogio (6)<ref name="S" />; Ambrosio-Gigli-Savare (59-61)<ref name="AGS" />-->
<!--Bauschke-Combette Ch 12.<ref name="BC" />; Santambrogio (6)<ref name="S" />; Ambrosio-Gigli-Savare (59-61)<ref name="AGS" />-->
<references>
<references>
<ref name="AGS">Ambrosio, Luigi, Nicola Gigli, and Giuseppe Savaré. ''Gradient Flows in Metric Spaces and in the Space of Probability Measures.'' Ch. 3.1. Birkhäuser, 2005.</ref>
<ref name="OT">Craig, Katy C. Lower Semicontinuity in the Narrow Topology. Math 260J. Univ. of Ca. at Santa Barbara. Winter 2022.</ref>
<ref name="OT">Craig, Katy C. Lower Semicontinuity in the Narrow Topology. Math 260J. Univ. of Ca. at Santa Barbara. Winter 2022.</ref>
<ref name="BC">Bauschke, Heinz H. and Patrick L. Combettes. ''Convex Analysis and Monotone Operator Theory in Hilbert Spaces, 2nd Ed.'' Ch. 12. Springer, 2017.</ref>
<ref name="S">Santambrogio, Filippo. ''Optimal Transport for Applied Mathematicians: Calculus of Variations, PDEs, and Modeling'' Ch. 1.1. Birkhäuser, 2015.</ref>
<ref name="S">Santambrogio, Filippo. ''Optimal Transport for Applied Mathematicians: Calculus of Variations, PDEs, and Modeling'' Ch. 1.1. Birkhäuser, 2015.</ref>
<!--
<ref name="BC">Bauschke, Heinz H. and Patrick L. Combettes. ''Convex Analysis and Monotone Operator Theory in Hilbert Spaces, 2nd Ed.'' Ch. 12. Springer, 2017.</ref>
<ref name="AGS">Ambrosio, Luigi, Nicola Gigli, and Giuseppe Savaré. ''Gradient Flows in Metric Spaces and in the Space of Probability Measures.'' Ch. 3.1. Birkhäuser, 2005.</ref>
-->


</references>
</references>

Revision as of 19:48, 10 February 2022

The Moreau-Yosida regularization is a technique used to approximate lower semicontinuous functions by Lipschitz functions. The main application of this result is to prove Portmanteau's Theorem, which states that integration against a lower semicontinuous and bounded below function is lower semicontinuous with respect to the narrow topology on the space of probability measures.

Definitions

Let be a metric space. A function is said to be proper if it is not identically equal to , that is, if there exists such that .

For a given function and , its Moreau-Yosida regularization is given by


The distance term may often be raised to a positive exponent. For example, when is a Hilbert space [1] [2], is taken to be


Note that

.

Examples

  • If , then by definition is constant and .
  • If is not proper, then for all .

Take . If is finite-valued and differentiable, we can explicitly write down . Then for a fixed , the map is continuous everywhere and differentiable everywhere except for when , where the derivative does not exist due to the absolute value. Thus we can apply standard optimization techniques from Calculus to solve for : find the critical points of and take the infimum of evaluated at the critical points. One of these values will always be the original function evaluated at , since this corresponds to the critical point for .

  • Let . Then
Plot of and for .

Approximating Lower Semicontinuous Functions by Lipschitz Functions

Proposition. [3][4]

  • If is proper and bounded below, so is . Furthermore, is continuous for all .
  • If, in addition, is lower semicontinuous, then for all .
  • In this case, is continuous and bounded and for all .
Plot of and for .

Proof.

  • Since is proper, there exists such that . Then for any

Thus is proper and bounded below. Next, for a fixed , let . Then as

,

the family is uniformly Lipschitz and hence equicontinuous. Thus is Lipschitz continuous.

  • Suppose that is also lower semicontinuous. Note that for all , . Thus it suffices to show that . This inequality is automatically satisfied when the left hand side is infinite, so without loss of generality assume that . By definition of infimum, for each there exists such that
.

Then

is bounded below by assumption, while the only way is finite in the limit is for to go to zero. Thus converges to in , and by lower semicontinuity of ,

.
  • By definition, . Since for all , for all .

Portmanteau Theorem

Theorem (Portmanteau). Let be lower semicontinuous and bounded below. Then the functional is lower semicontinuous with respect to narrow convergence in , that is

.

Proof. By the Moreau-Yosida approximation, for all ,

.

Taking , Fatou's lemma ensures that

.

The Mysterious Etymology of Portmanteau

(spurious stuff, will fill in later)

References

  1. Bauschke, Heinz H. and Patrick L. Combettes. Convex Analysis and Monotone Operator Theory in Hilbert Spaces, 2nd Ed. Ch. 12. Springer, 2017.
  2. Ambrosio, Luigi, Nicola Gigli, and Giuseppe Savaré. Gradient Flows in Metric Spaces and in the Space of Probability Measures. Ch. 3.1. Birkhäuser, 2005.
  3. Craig, Katy C. Lower Semicontinuity in the Narrow Topology. Math 260J. Univ. of Ca. at Santa Barbara. Winter 2022.
  4. Santambrogio, Filippo. Optimal Transport for Applied Mathematicians: Calculus of Variations, PDEs, and Modeling Ch. 1.1. Birkhäuser, 2015.