A multigrid scheme for 3D Monge–Ampère equations*

Jun Liu, Brittany D. Froese, Adam M. Oberman, Mingqing Xiao

Research output: Contribution to journalArticlepeer-review

7 Scopus citations


The elliptic Monge–Ampère equation is a fully nonlinear partial differential equation which has been the focus of increasing attention from the scientific computing community. Fast three-dimensional solvers are needed, for example in medical image registration but are not yet available. We build fast solvers for smooth solutions in three dimensions using a nonlinear full-approximation storage multigrid method. Starting from a second-order accurate centred finite difference approximation, we present a nonlinear Gauss–Seidel iterative method which has a mechanism for selecting the convex solution of the equation. The iterative method is used as an effective smoother, combined with the full-approximation storage multigrid method. Numerical experiments are provided to validate the accuracy of the finite difference scheme and illustrate the computational efficiency of the proposed multigrid solver.

Original languageEnglish (US)
Pages (from-to)1850-1866
Number of pages17
JournalInternational Journal of Computer Mathematics
Issue number9
StatePublished - Sep 2 2017

All Science Journal Classification (ASJC) codes

  • Computer Science Applications
  • Computational Theory and Mathematics
  • Applied Mathematics


  • FAS multigrid method
  • Gauss–Seidel iteration
  • Monge–Ampère equation
  • finite difference method
  • nonlinear partial differential equations


Dive into the research topics of 'A multigrid scheme for 3D Monge–Ampère equations*'. Together they form a unique fingerprint.

Cite this