Picard-Based Acceleration of Newton Continuation for Mean Field Game PDE Systems
Abstract
We develop a method that uses Picard iterations to accelerate Newton continuation for the semi-implicit finite-difference discretization of forward-backward partial differential equation (PDE) systems arising in mean field games (MFGs). We first investigate the computational properties of the Picard and Newton methods, which are widely used separately in the MFG literature but whose comparative cost and robustness across different regimes remain insufficiently explored and documented. The Picard method uses an outer fixed-point iteration that alternates a forward Fokker-Planck solve and a backward Hamilton-Jacobi-Bellman solve. The Newton method instead applies Newton's method directly to the coupled nonlinear space-time system. Across one- and two-dimensional MFG benchmarks, Picard offers substantial computational savings in favorable regimes, but may fail at sufficiently low viscosity or require strong damping under temporal shocks, making Newton continuation preferable. With parameter continuation in the viscosity parameter, the Newton method is more robust in these regimes, at the cost of larger coupled linear systems. We relate these trade-offs to the residuals, Jacobian blocks, and sparsity structures produced by separable, local nonseparable, and nonlocal Hamiltonians. We then demonstrate how to combine inexpensive Picard iterations with Newton continuation in a hybrid method to reduce the total computational cost for a two- dimensional double-well problem.
Related
- Acceleration of Newton's method for nonlinear algebraic systems by preflattening and elimination techniques
- Grad-div stabilized finite element Picard-Newton iterative methods for the Navier-Stokes equations
- On a blended Anderson / NGMRES acceleration for the Navier-Stokes Picard iteration: convergence analysis and adaptive depth selection
- Matching the Lower Bounds: Stochastic Contracting Cubic Newton and Its Optimal Acceleration
- Hybrid accelerated iterative scheme based on the Picard–Ishikawa approach for systems of nonlinear equations
BibTeX
@misc{lauriere2026picardbased,
title = {{Picard-Based Acceleration of Newton Continuation for Mean Field Game PDE Systems}},
author = {Mathieu Lauriere and Andrew Shi},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.11050},
eprint = {2609.11050},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.11050},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Picard-Based Acceleration of Newton Continuation for Mean Field Game PDE Systems
AU - Mathieu Lauriere
AU - Andrew Shi
PY - 2026
DA - 2026/09/15
PB - arXiv
DO - 10.48550/arxiv.2609.11050
UR - https://arxiv.org/abs/2609.11050
AB - We develop a method that uses Picard iterations to accelerate Newton continuation for the semi-implicit finite-difference discretization of forward-backward partial differential equation (PDE) systems arising in mean field games (MFGs). We first investigate the computational properties of the Picard and Newton methods, which are widely used separately in the MFG literature but whose comparative cost and robustness across different regimes remain insufficiently explored and documented. The Picard method uses an outer fixed-point iteration that alternates a forward Fokker-Planck solve and a backward Hamilton-Jacobi-Bellman solve. The Newton method instead applies Newton's method directly to the coupled nonlinear space-time system. Across one- and two-dimensional MFG benchmarks, Picard offers substantial computational savings in favorable regimes, but may fail at sufficiently low viscosity or require strong damping under temporal shocks, making Newton continuation preferable. With parameter continuation in the viscosity parameter, the Newton method is more robust in these regimes, at the cost of larger coupled linear systems. We relate these trade-offs to the residuals, Jacobian blocks, and sparsity structures produced by separable, local nonseparable, and nonlocal Hamiltonians. We then demonstrate how to combine inexpensive Picard iterations with Newton continuation in a hybrid method to reduce the total computational cost for a two- dimensional double-well problem.
ER -