A standard sufficient form of Brenier theorem assumes have finite second moments and , that is, absolute continuity of measures with respect to Lebesgue measure. For the cost , there exists a unique optimal transport plan, and it is induced by a transport map:
Here is a proper convex function, which may be chosen sequentially lower semicontinuous, and its gradient exists -almost everywhere. The map is unique -almost everywhere and is the unique minimizer of the Monge optimal transport problem; its cost equals the Kantorovich optimal transport problem minimum. Equivalently, it is the unique gradient of a convex function transporting to .
The uniqueness claim concerns the map and the transport plan, not a globally unique potential. The potential may be shifted by a constant, and additional nonuniqueness away from the source can occur. No density assumption is required on . A primary reference is Brenier's Polar factorization and monotone rearrangement of vector-valued functions.
Define ; the printed quotient is undefined at the origin, but this continuous extension changes no transport cost because . In polar coordinates, the two probability density functions give
Both angular distributions are uniform, with independence of angle and radius. The proposed transport map preserves the angle and sends to . Therefore
and preservation of the angle proves . As a local check using the Jacobian determinant, its radial and tangential derivatives for are and , respectively, so and .
Now use the convex function
It is convex because the Euclidean norm is convex and is increasing and convex on . It is differentiable, including at zero, and
Its graph transport plan lies in the graph of , so the Knott–Smith optimality criterion proves quadratic optimality. Since that plan is induced by a map, part 1(b) proves optimality for the Monge optimal transport problem as well. The minimum provides a useful independent check:
The empirical risk is
For ,
so is a convex function. Each signed margin is an affine function of , and composition with an affine function preserves convexity. Equivalently,
is a positive semidefinite matrix. Therefore is convex.
For and , convexity gives
Taking the maximum over and then bounding each term by the corresponding endpoint maxima yields
Thus the finite pointwise maximum of convex functions is convex. Summing the original inequalities over proves that is convex as well.
Taking gives , so every matrix is -invertible. The smallest admissible value is
The map inside the norm is affine in , and a norm composed with an affine map is a convex function. The space of matrices is a convex feasible set, so this is a convex optimization problem.
A -minimizing solution is an exact solution satisfying
The source condition in variational regularization says that some satisfies
where is the subdifferential.
Fix . The subgradient optimality condition for
at is
If the source condition holds, choose . Since , the displayed inclusion holds. Because the objective is a convex function, this condition is necessary and sufficient for global minimality.
Conversely, if the stated range condition holds for some , optimality supplies with
Thus , which is precisely the source condition. This proves the equivalence.
For a convex function , its subdifferential is
For ,
The defining inequality shows directly that exactly when is a global minimizer.
Strict convexity means
for distinct and . Two distinct minimizers would make the midpoint have a strictly smaller value, so a minimizer is unique.
A function on an interval is a convex function if
Jensen inequality states that for an integrable random variable , when all terms are defined,
The function is convex on because
Applying Jensen's inequality to the uniform distribution on the positive numbers gives
Writing and multiplying by yields
A vector is a subgradient of the convex function at when
for every . The set of all such vectors is the subdifferential . The proximal operator satisfies
More generally, exactly when .
Each map is an affine function. For , the pointwise maximum of convex functions satisfies
so is convex.
A vector is a subgradient of a convex function at when
for every . Choose any active index . Then
and hence . More generally, every convex combination of the active vectors is a subgradient, and in fact
The first-order subgradient inequality for a differentiable convex function gives
and, after interchanging and ,
Adding and rearranging yields
Thus the gradient is a monotone operator.
The function is strictly convex. Subject to , the Jensen inequality gives
with equality only at
This is the minimum of the relative potential, so it is a stable relative equilibrium.
Let near an equally spaced configuration, and put
The gap perturbations are
whose sum vanishes. A Taylor expansion gives the quadratic energies
Thus the linearized equations are
The matrix in parentheses is the Graph Laplacian of . Its constant eigenvector has eigenvalue zero, while every vector whose components sum to zero has eigenvalue three. One orthonormal set of three normal modes and their angular frequencies is therefore
The zero mode is rigid rotation. The positive and equal squared frequencies on the two-dimensional relative subspace prove stability and agree with the general normal modes of three equal masses with a symmetric gap potential.
Coordinate descent cycles through the intercept and coefficient coordinates, minimizing the convex ridge objective in one coordinate while holding the others fixed. Given current coefficients, update
For coordinate , form the partial residual
and update it by the exact one-dimensional minimizer
Repeated sweeps converge to the unique fitted value because the objective is a convex function; with it is strictly convex in .
The Euclidean projection onto a convex set is nonexpansive, and because the optimum is feasible. Therefore the projected subgradient method satisfies
The subgradient inequality gives , while Lipschitz continuity of the finite convex function gives . Hence
Summing this telescoping inequality for , and then bounding the smallest term by the average, yields
Writing , the right-hand side is minimized by the constant step size
Substitution gives
If , the initial point is already optimal and the result is immediate.
Write for row of , and define affine functions
Then
The pointwise maximum of convex functions is convex, so is a convex function. Moreover,
so has Lipschitz continuity. The simpler bound is also valid.
Let be the active set. The subdifferential is
the convex hull of all active slopes. In particular, choosing any active index gives the subgradient
At a tie, every convex combination of the tied slopes is also valid.
For a proper lower-semicontinuous convex function , its proximal operator is
The squared norm is strongly convex, so the minimizer is unique. The subdifferential sum rule gives the necessary and sufficient condition
More generally,
The subgradient inversion rule for the convex conjugate says exactly when . Hence
which is precisely the proximal optimality condition
Since , this proves the generalized Moreau decomposition
For a convex function , its subdifferential at is
Its elements are the subgradients of at .
Suppose and minimize . Their midpoint is also a minimizer because is a convex function. The group penalty is convex, while the squared Euclidean norm is strictly convex in the fitted value. If , strict convexity would make the midpoint objective strictly smaller than the minimum. Hence , so the fitted values are unique.
Let and put and . For each alphabet symbol , apply the log-sum inequality to , and the corresponding values. Summing over gives
which is joint convexity in .
An element is a -minimizing solution when
It satisfies the source condition in variational regularization when there is a such that
where is the subdifferential of the convex functional .