On the convergence rates of generalized conditional gradient method for fully discretized Mean Field Games
Haruka Nakamura, Norikazu Saito
TL;DR
This work studies the convergence of a fully discrete generalized conditional gradient (GCG) method for Mean Field Games (MFG) by reformulating the problem via the Cole–Hopf transform into CH–GCG and applying finite-difference discretization. It proves a discrete maximum principle and derives explicit joint error bounds that couple discretization errors $( au, au_x)$ with iteration errors (in $k$), showing non-uniform iteration dependence but convergence as both mesh finiteness and iteration progress. Higher regularity of CH–GCG components yields improved discretization rates (up to first- or second-order in space-time under suitable conditions), while the iteration term decays polynomially in $k$, with the constants potentially growing with $k$. Numerical experiments in 1D and 2D corroborate the theoretical rates and demonstrate the practical impact of step-size rules on convergence, guiding parameter choices for efficient computation of fully discretized MFG solutions.
Abstract
We study convergence rates of the generalized conditional gradient (GCG) method applied to fully discretized Mean Field Games (MFG) systems. While explicit convergence rates of the GCG method have been established at the continuous PDE level, a rigorous analysis that simultaneously accounts for time-space discretization and iteration errors has been missing. In this work, we discretize the MFG system using finite difference method and analyze the resulting fully discrete GCG scheme. Under suitable structural assumptions on the Hamiltonian and coupling terms, we establish discrete maximum principles and derive explicit error estimates that quantify both discretization errors and iteration errors within a unified framework. Our estimates show how the convergence rates depend on the mesh sizes and the iteration number, and they reveal a non-uniform behavior with respect to the iteration. Moreover, we prove that higher convergence rates can be achieved under additional regularity assumptions on the solution. Numerical experiments are presented to illustrate the theoretical results and to confirm the predicted convergence behavior.
