Fast finite difference solvers for singular solutions of the elliptic Monge-Ampère equation
arXiv:1006.5748 · doi:10.1016/j.jcp.2010.10.020
Abstract
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.
23 pages, 4 figures, 4 tables; added arxiv links to references, added coments
Cited by in corpus (9)
- Solving the Monge-Ampère Equations for the Inverse Reflector Problem
- Pointwise rates of convergence for the Oliker-Prussner method for the Monge-Ampère equation
- Monotone Mixed Finite Difference Scheme for Monge-Ampère Equation
- A viscosity solution approach to the Monge-Ampere formulation of the Optimal Transportation Problem
- A multigrid scheme for 3D Monge-Ampère equations
- Blended numerical schemes for the advection equation and conservation laws
- Numerical Methods for the 2-Hessian Elliptic Partial Differential Equation
- Confinement-induced drift in Marangoni-driven transport of surfactant: a Lagrangian perspective
- Convergence of a hybrid scheme for the elliptic Monge-Ampere equation