Mass Preserving Mappings and Image Registration
MICCAI '01 Proceedings of the 4th International Conference on Medical Image Computing and Computer-Assisted Intervention
Optimal Mass Transport for Registration and Warping
International Journal of Computer Vision
A convergent monotone difference scheme for motion of level sets by mean curvature
Numerische Mathematik
Consistency of Generalized Finite Difference Schemes for the Stochastic HJB Equation
SIAM Journal on Numerical Analysis
SIAM Journal on Numerical Analysis
Journal of Computational Physics
Journal of Scientific Computing
Space deformations, surface deformations and the opportunities in-between
Journal of Computer Science and Technology
SIAM Journal on Numerical Analysis
Moving Mesh Generation Using the Parabolic Monge-Ampère Equation
SIAM Journal on Scientific Computing
The Monge-Ampère equation: Various forms and numerical solution
Journal of Computational Physics
Fast finite difference solvers for singular solutions of the elliptic Monge-Ampère equation
Journal of Computational Physics
Quadratic Finite Element Approximations of the Monge-Ampère Equation
Journal of Scientific Computing
Numerical solution of the Optimal Transportation problem using the Monge-Ampère equation
Journal of Computational Physics
Journal of Computational and Applied Mathematics
Hi-index | 0.01 |
The elliptic Monge-Ampère equation is a fully nonlinear partial differential equation that 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. Novel solution methods are required for stability and convergence to the weak (viscosity) solution. In this article we build a wide stencil finite difference discretization for the Monge-Ampère equation. The scheme is monotone, so the Barles-Souganidis theory allows us to prove that the solution of the scheme converges to the unique viscosity solution of the equation. Solutions of the scheme are found using a damped Newton's method. We prove convergence of Newton's method and provide a systematic method to determine a starting point for the Newton iteration. Computational results are presented in two and three dimensions, which demonstrates the speed and accuracy of the method on a number of exact solutions, which range in regularity from smooth to nondifferentiable.