187 research outputs found
The Monge-Ampere equation: various forms and numerical methods
We present three novel forms of the Monge-Ampere equation, which is used,
e.g., in image processing and in reconstruction of mass transportation in the
primordial Universe. The central role in this paper is played by our Fourier
integral form, for which we establish positivity and sharp bound properties of
the kernels. This is the basis for the development of a new method for solving
numerically the space-periodic Monge-Ampere problem in an odd-dimensional
space. Convergence is illustrated for a test problem of cosmological type, in
which a Gaussian distribution of matter is assumed in each localised object,
and the right-hand side of the Monge-Ampere equation is a sum of such
distributions.Comment: 24 pages, 2 tables, 5 figures, 32 references. Submitted to J.
Computational Physics. Times of runs added, multiple improvements of the
manuscript implemented
Fast finite difference solvers for singular solutions of the elliptic Monge-Amp\`ere equation
The elliptic Monge-Ampere equation is a fully nonlinear Partial Differential
Equation which originated in geometric surface theory, and has been applied in
dynamic meteorology, elasticity, geometric optics, image processing and image
registration. Solutions can be singular, in which case standard numerical
approaches fail. In this article we build a finite difference solver for the
Monge-Ampere equation, which converges even for singular solutions. Regularity
results are used to select a priori between a stable, provably convergent
monotone discretization and an accurate finite difference discretization in
different regions of the computational domain. This allows singular solutions
to be computed using a stable method, and regular solutions to be computed more
accurately. The resulting nonlinear equations are then solved by Newton's
method. Computational results in two and three dimensions validate the claims
of accuracy and solution speed. A computational example is presented which
demonstrates the necessity of the use of the monotone scheme near
singularities.Comment: 23 pages, 4 figures, 4 tables; added arxiv links to references, added
coment
Entropic and displacement interpolation: a computational approach using the Hilbert metric
Monge-Kantorovich optimal mass transport (OMT) provides a blueprint for
geometries in the space of positive densities -- it quantifies the cost of
transporting a mass distribution into another. In particular, it provides
natural options for interpolation of distributions (displacement interpolation)
and for modeling flows. As such it has been the cornerstone of recent
developments in physics, probability theory, image processing, time-series
analysis, and several other fields. In spite of extensive work and theoretical
developments, the computation of OMT for large scale problems has remained a
challenging task. An alternative framework for interpolating distributions,
rooted in statistical mechanics and large deviations, is that of Schroedinger
bridges (entropic interpolation). This may be seen as a stochastic
regularization of OMT and can be cast as the stochastic control problem of
steering the probability density of the state-vector of a dynamical system
between two marginals. In this approach, however, the actual computation of
flows had hardly received any attention. In recent work on Schroedinger bridges
for Markov chains and quantum evolutions, we noted that the solution can be
efficiently obtained from the fixed-point of a map which is contractive in the
Hilbert metric. Thus, the purpose of this paper is to show that a similar
approach can be taken in the context of diffusion processes which i) leads to a
new proof of a classical result on Schroedinger bridges and ii) provides an
efficient computational scheme for both, Schroedinger bridges and OMT. We
illustrate this new computational approach by obtaining interpolation of
densities in representative examples such as interpolation of images.Comment: 20 pages, 7 figure
- …