Table of Contents
Fetching ...

A Globally Convergent Flow for Time-Dependent Mean Field Games and a Solver-Agnostic Framework for Inverse Problems

Hanwei Yan, Xianjin Yang, Jingguo Zhang

TL;DR

A solver-agnostic framework is proposed, in which parameter updates are computed by implicitly differentiating the discrete MFG equations satisfied by the converged MFG solution, rather than by differentiating through a particular forward solver.

Abstract

Mean field games (MFGs) model the limit of large populations of strategically interacting agents, yet both forward and inverse problems remain challenging. For the forward problem, a difficulty is to design numerical methods with global convergence guarantees whose convergence does not depend on careful initialization. For the inverse problem, a difficulty is to decouple parameter optimization from the forward solver, so that parameter updates do not depend on implementation details, and the inverse method does not need to be reformulated when the forward solver is changed. We address both issues as follows. For the forward problem, we propose a monotone Hessian-Riemannian flow for time-dependent MFGs on the feasible manifold of densities. The flow preserves the positivity of densities and is proved to be globally convergent. For the inverse problem, we cast parameter estimation as an outer optimization problem over the unknown coefficients, with the MFG system solved in an inner step for each parameter value. For solving this problem, we consider an adjoint-based gradient method together with a Gauss-Newton acceleration. This leads to a solver-agnostic framework, in which parameter updates are computed by implicitly differentiating the discrete MFG equations satisfied by the converged MFG solution, rather than by differentiating through a particular forward solver. We demonstrate the approach on several stationary and time-dependent MFG examples, where the Gauss-Newton method consistently requires fewer outer iterations than gradient descent.

A Globally Convergent Flow for Time-Dependent Mean Field Games and a Solver-Agnostic Framework for Inverse Problems

TL;DR

A solver-agnostic framework is proposed, in which parameter updates are computed by implicitly differentiating the discrete MFG equations satisfied by the converged MFG solution, rather than by differentiating through a particular forward solver.

Abstract

Mean field games (MFGs) model the limit of large populations of strategically interacting agents, yet both forward and inverse problems remain challenging. For the forward problem, a difficulty is to design numerical methods with global convergence guarantees whose convergence does not depend on careful initialization. For the inverse problem, a difficulty is to decouple parameter optimization from the forward solver, so that parameter updates do not depend on implementation details, and the inverse method does not need to be reformulated when the forward solver is changed. We address both issues as follows. For the forward problem, we propose a monotone Hessian-Riemannian flow for time-dependent MFGs on the feasible manifold of densities. The flow preserves the positivity of densities and is proved to be globally convergent. For the inverse problem, we cast parameter estimation as an outer optimization problem over the unknown coefficients, with the MFG system solved in an inner step for each parameter value. For solving this problem, we consider an adjoint-based gradient method together with a Gauss-Newton acceleration. This leads to a solver-agnostic framework, in which parameter updates are computed by implicitly differentiating the discrete MFG equations satisfied by the converged MFG solution, rather than by differentiating through a particular forward solver. We demonstrate the approach on several stationary and time-dependent MFG examples, where the Gauss-Newton method consistently requires fewer outer iterations than gradient descent.
Paper Structure (23 sections, 7 theorems, 168 equations, 11 figures, 4 tables)

This paper contains 23 sections, 7 theorems, 168 equations, 11 figures, 4 tables.

Key Result

Proposition 1

Assume that the numerical Hamiltonian $g$ satisfies ass:disc_g and that the coupling $f$ satisfies ass:disc_f. Let $Y$ and $\widetilde{Y}$ be two feasible states with the same fixed end-point samples $(M_0,U_{N_T})$. Then $F_h$ in eq:Fh_full_def is strictly monotone with respect to the space-time pa

Figures (11)

  • Figure 1: Numerical results for the inverse problem of the one-dimensional stationary MFG in \ref{['eq:stationaryMFG_oneDim_inverse_effHam']}. (a), (b), (c) are references for $m,u,V$; (d) Log-log plot comparing the GD and GN losses across iterations; (e), (f), (g) recovered $m,u,V$ via GD; (h), (m), (o) errors of $m,u,V$ via GD; (i), (j), (k) recovered $m,u,V$ via GN; (l), (n), (p) errors of $m,u,V$ via GN.
  • Figure 2: The inverse problem of the stationary MFG \ref{['eq:first_order_MFGS_congestion']}: samples for $m$, $u$ and $V$ & corresponding observation points.
  • Figure 3: Numerical results for solving the inverse problem of the MFG system in \ref{['eq:first_order_MFGS_congestion']} using the adjoint-based GD method and the GN method: (a)-(c) show the $L^2$ errors for $m$, $u$, and $V$, respectively, versus their exact counterparts as the number of observation points increases.
  • Figure 4: Numerical results for the inverse problem of the stationary MFG in \ref{['eq:first_order_MFGS_congestion']}. (a), (b), (c) are references for $m,u,V$; (d) log-log plot of the loss comparison for GD and GN across iterations; (e), (f), (g) recovered $m,u,V$ via GD; (h), (m), (o) errors of $m,u,V$ via GD; (i), (j), (k) recovered $m,u,V$ via GN; (l), (n), (p) errors of $m,u,V$ via GN.
  • Figure 5: The inverse problem of the stationary MFG \ref{['eq:2DstationaryMFGinvSmall_nu']}: samples for $m$, $u$ and $V$ & corresponding observation points.
  • ...and 6 more figures

Theorems & Definitions (15)

  • Proposition 1: Strict monotonicity of $F_h$
  • Proposition 2: Explicit Discrete HRF
  • Proposition 3: Global Existence, Positivity, and Mass Preservation
  • Theorem 2.1: Global convergence of the HRF
  • Remark 3.1: RKHS surrogates for continuous observations
  • Proposition 4
  • Proposition 5: Gauss--Newton step
  • Lemma 1: Discrete adjoint structure of the transport operator
  • proof
  • proof : Proof of \ref{['prop:Fh_monotone_MU']}
  • ...and 5 more