Trivariate Spline Collocation Methods for Numerical Solution to 3D Monge-Ampère Equation

We use trivariate spline functions for the numerical solution of the Dirichlet problem of the 3D elliptic Monge-Amp\'ere equation. Mainly we use the spline collocation method introduced in \cite{LL21} to numerically solve iterative Poisson equations and use an averaged algorithm to ensure the convergence of the iterations. We shall also establish the rate of convergence under a sufficient condition and provide some numerical evidence to show the numerical rates. Then we present many computational results to demonstrate that this approach works very well. In particular, we tested many known convex solutions as well as nonconvex solutions over convex and nonconvex domains and compared them with several existing numerical methods to show the efficiency and effectiveness of our approach.