A damped Newton algorithm for generated Jacobian equations

Abstract

Generated Jacobian Equations have been introduced by Trudinger (Discrete Contin Dyn Syst A 34(4):1663–1681, 2014) as a generalization of Monge–Ampère equations arising in optimal transport. In this paper, we introduce and study a damped Newton algorithm for solving these equations in the semi-discrete setting, meaning that one of the two measures involved in the problem is finitely supported and the other one is absolutely continuous. We also present a numerical application of this algorithm to the near-field parallel reflector problem arising in non-imaging problems.