We propose a definition of the Dirichlet energy (which is roughly speaking the integral of the square of the gradient) for mappings mu : Omega -> (P(D), W ) defined over a subset Omega of R^p and valued in the space P(D) of probability measures on a compact convex subset D of R^q endowed with the quadratic Wasserstein distance. Our definition relies on a straightforward generalization of the Benamou-Brenier formula (already introduced by Brenier) but is also equivalent to the definition of Koorevaar, Schoen and Jost as limit of approximate Dirichlet energies, and to the definition of Reshetnyak of Sobolev spaces valued in metric spaces. We study harmonic mappings, i.e. minimizers of the Dirichlet energy provided that the values on the boundary d Omega are fixed. The notion of constant-speed geodesics in the Wasserstein space is recovered by taking for Omega a segment of R. As the Wasserstein space (P(D), W ) is positively curved in the sense of Alexandrov we cannot apply the theory of Koorevaar, Schoen and Jost and we use instead arguments based on optimal transport. We manage to get existence of harmonic mappings provided that the boundary values are Lipschitz on d Omega, uniqueness is an open question. If Omega is a segment of R, it is known that a curve valued in the Wasserstein space P(D) can be seen as a superposition of curves valued in D. We show that it is no longer the case in higher dimensions: a generic mapping Omega -> P(D) cannot be represented as the superposition of mappings Omega -> D. We are able to show a Ishihara-type property: the composition F(mu) of a function F : P(D) -> R convex along generalized geodesics and a harmonic mapping mu : Omega -> P(D) is a subharmonic real-valued function. We also study the special case where we restrict ourselves to a given family of elliptically contoured distributions (a finite-dimensional and geodesically convex submanifold of (P(D), W ) which generalizes the case of Gaussian measures) and show that it boils down to harmonic mappings valued in the Riemannian manifold of symmetric matrices endowed with the distance coming from optimal transport.
Finally, in @cite_2 the authors also study mappings valued in the space of probability measures, but are rather interested in the bounded variation norm (the integral of the norm of the gradient) than in the Dirichlet energy. Their provide applications to the denoising of measure-valued images.
@cite_2
2133027741
These lecture notes highlight the mathematical and computational structure relating to the formulation of, and development of algorithms for, the Bayesian approach to inverse problems in differential equations. This approach is fundamental in the quantification of uncertainty within applications in volving the blending of mathematical models with data. The finite dimensional situation is described first, along with some motivational examples. Then the development of probability measures on separable Banach space is undertaken, using a random series over an infinite set of functions to construct draws; these probability measures are used as priors in the Bayesian approach to inverse problems. Regularity of draws from the priors is studied in the natural Sobolev or Besov spaces implied by the choice of functions in the random series construction, and the Kolmogorov continuity theorem is used to extend regularity considerations to the space of Holder continuous functions. Bayes’ theorem is de rived in this prior setting, and here interpreted as finding conditions under which the posterior is absolutely continuous with respect to the prior, and determining a formula for the Radon-Nikodym derivative in terms of the likelihood of the data. Having established the form of the posterior, we then describe various properties common to it in the infinite dimensional setting. These properties include well-posedness, approximation theory, and the existence of maximum a posteriori estimators. We then describe measure-preserving dynamics, again on the infinite dimensional space, including Markov chain-Monte C arlo and sequential Monte Carlo methods, and measure-preserving reversible stochastic differential equations. By formulating the theory and algorithms on the underlying infinite dimensional space, we obtain a framework suitable for rigorous analysis of the accuracy of reconstructions, of computational complexity, as well as naturally constructing algorithms which perform well under mesh refinement, since they are inherently well-defined in infinite dimensions.