A least-squares method for optimal transport using the Monge-Ampère equation

In this article we introduce a novel numerical method to solve the problem of optimal mass transport and the related elliptic Monge-Ampère equation. It is one of the few numerical algorithms capable of solving this problem efficiently with the proper boundary conditions. It scales well with the grid size and has the additional advantage that the target domain may be non-convex. We present the method and several numerical experiments.