A game-theoretical interpretation for a doubly nonlinear parabolic equation
Pith reviewed 2026-05-10 15:15 UTC · model grok-4.3
The pith
A new asymptotic mean value formula for the p-Laplacian defines a dynamic programming principle whose solutions match viscosity solutions and arise as values of a two-player stochastic game.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We introduce a game-theoretical framework for the doubly nonlinear parabolic equation |∂_t u|^{p-2} ∂_t u - Δ_p u = 0. A key feature to our approach is a new asymptotic mean value formula for the p-Laplacian that is robust even when the gradient vanishes and is independent of the sign of the p-Laplacian. This new AMVF leads naturally to a dynamic programming principle whose solutions converge to the viscosity solution of the boundary value problem for the differential equation. In addition, solutions to the DPP coincide with value functions for a stochastic, two-players, zero-sum game that we introduce and analyze here.
What carries the argument
The new asymptotic mean value formula for the p-Laplacian, which remains valid at vanishing gradients and independent of sign, and which directly yields both the dynamic programming principle and the stochastic game.
If this is right
- Solutions of the dynamic programming principle converge to the viscosity solution of the boundary value problem.
- Solutions of the dynamic programming principle are exactly the value functions of the stochastic two-player zero-sum game.
- The game supplies a probabilistic representation of solutions to the doubly nonlinear equation.
- The construction applies for all p greater than 2.
Where Pith is reading between the lines
- Game-simulation algorithms might serve as practical numerical schemes for approximating solutions to the PDE.
- Similar robust mean-value formulas, if found for other nonlinear operators, could yield game interpretations for a wider class of parabolic equations.
- Convergence rates or error estimates between the discrete game values and the continuous viscosity solution remain open for quantitative study.
Load-bearing premise
The proposed asymptotic mean value formula continues to hold and produce a well-defined dynamic programming principle at every point where the spatial gradient is zero.
What would settle it
A concrete function or numerical test case in which the dynamic programming principle solutions fail to converge to a known viscosity solution of the boundary value problem, or in which the mean value formula deviates from the p-Laplacian at a zero-gradient point.
read the original abstract
We introduce a game-theoretical framework for the doubly nonlinear parabolic equation \[ |\partial_t u|^{p-2} \partial_t u - \Delta_p u = 0. \] where $\Delta_p u = \nabla \cdot ( |\nabla u |^{p-2} \nabla u)$ with $p>2$ is the standard $p-$Laplacian. A key feature to our approach is a new asymptotic mean value formula (AMVF) for the $p-$Laplacian that is robust even when the gradient vanishes and is independent of the sign of the $p-$Laplacian. This new AMVF leads naturally to a dynamic programming principle (DPP) whose solutions converge to the viscosity solution of the boundary value problem for the differential equation. In addition, solutions to the DPP coincide with value functions for a stochastic, two-players, zero-sum game that we introduce and analyze here.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a game-theoretical framework for the doubly nonlinear parabolic equation |∂_t u|^{p-2} ∂_t u - Δ_p u = 0 (p>2). It proposes a new asymptotic mean value formula (AMVF) for the p-Laplacian that is claimed to be robust at vanishing gradients and independent of the sign of the operator. This AMVF is used to derive a dynamic programming principle (DPP) whose solutions converge to the viscosity solution of the boundary-value problem; the DPP solutions are further shown to coincide with the value functions of a stochastic two-player zero-sum game introduced in the paper.
Significance. If the uniform error control in the new AMVF holds, the work supplies a probabilistic interpretation and a game-theoretic characterization of solutions to a degenerate parabolic equation, extending mean-value approaches to a setting where standard derivations break down. This could open avenues for numerical approximation via game iterations and clarify the connection between viscosity solutions and stochastic processes for doubly nonlinear problems.
major comments (1)
- [Abstract (AMVF and DPP convergence)] The central convergence claim from the DPP to the viscosity solution rests on the new AMVF holding with an error term that is o(r^2) uniformly in |∇u|, including at points where ∇u = 0. Standard p-Laplacian mean-value formulas typically normalize by |∇u|^{p-2} and degenerate there; without an explicit uniform estimate (independent of |∇u|) in the derivation of the AMVF and its insertion into the DPP, the viscosity inequality may fail to be recovered at critical points, which are the most degenerate locations for p>2. This is load-bearing for the abstract's claim that the DPP converges to the viscosity solution.
minor comments (2)
- Clarify the precise definition of the stochastic game (payoff, stopping time, and strategy spaces) to make the coincidence with DPP solutions fully explicit.
- The abstract states that the AMVF is 'independent of the sign of the p-Laplacian'; a brief remark on how this sign-independence is achieved would aid readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading and the detailed comment on the uniformity of the error term in our asymptotic mean value formula. We address this point below.
read point-by-point responses
-
Referee: The central convergence claim from the DPP to the viscosity solution rests on the new AMVF holding with an error term that is o(r^2) uniformly in |∇u|, including at points where ∇u = 0. Standard p-Laplacian mean-value formulas typically normalize by |∇u|^{p-2} and degenerate there; without an explicit uniform estimate (independent of |∇u|) in the derivation of the AMVF and its insertion into the DPP, the viscosity inequality may fail to be recovered at critical points, which are the most degenerate locations for p>2. This is load-bearing for the abstract's claim that the DPP converges to the viscosity solution.
Authors: We agree that uniformity of the o(r^2) error with respect to |∇u| is essential for recovering the viscosity inequalities at critical points. Our AMVF is derived precisely to achieve this robustness. In Theorem 3.1, the formula is obtained via a direct integral representation over the ball that does not involve division by |∇u|^{p-2}; the remainder is controlled by C r^2 where C depends only on the C^2-norm of the test function and on p, but is independent of |∇u(x)|. The proof proceeds by splitting into regimes |∇u| ≥ δ and |∇u| < δ, with uniform bounds in both cases. This estimate is then used verbatim in the DPP (Section 4) and the convergence argument (Section 5), so the viscosity inequalities hold without additional restrictions at points where ∇u = 0. We will add an explicit corollary after Theorem 3.1 stating the |∇u|-independence of the constant to make this feature more prominent. revision: partial
Circularity Check
No circularity: new AMVF introduced independently, then used to derive DPP and game equivalence
full rationale
The derivation begins with an explicitly new asymptotic mean value formula for the p-Laplacian (robust at vanishing gradients), proceeds to a DPP whose solutions are shown to converge to the viscosity solution of the target PDE, and finally equates those DPP solutions to value functions of a newly defined stochastic game. None of these steps reduce by construction to the inputs, rely on fitted parameters renamed as predictions, or depend on load-bearing self-citations whose content is itself unverified. The central claims remain independent of the paper's own fitted quantities or prior self-referential definitions.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Standard existence and uniqueness theory for viscosity solutions of parabolic PDEs
- domain assumption Well-posedness of the dynamic programming principle for the given game
invented entities (1)
-
The stochastic two-player zero-sum game
no independent evidence
Reference graph
Works this paper leans on
- [1]
-
[2]
Akagi, On some doubly-nonlinear parabolic equations posed inRd,Discrete Contin
G. Akagi, On some doubly-nonlinear parabolic equations posed inRd,Discrete Contin. Dyn. Syst. Ser. S, 16(12):3661–3676, 2023
work page 2023
-
[3]
G. Akagi and G. Schimperna, On a class of doubly nonlinear evolution equations in Musielak–Orlicz spaces, Math. Nachr., 297(7):2686–2729, 2024
work page 2024
-
[4]
G. Barles and P. E. Souganidis. Convergence of approximation schemes for fully nonlinear second order equations.Asymptotic Anal, 4(3):271-283, 1991
work page 1991
- [5]
- [6]
-
[7]
P. Blanc and J. D. Rossi. Game Theory and Partial Differential Equations, Berlin, Boston: De Gruyter. 2019
work page 2019
- [8]
-
[9]
A. Bonfiglioli and E. Lanconelli; Subharmonic functions in sub-Riemannian settings,J. Eur. Math. Soc., 15:387–441, 2013
work page 2013
-
[10]
M. G. Crandall, H. Ishii, and P.-L. Lions, User’s guide to viscosity solutions of second order partial differential equations,Bulletin of the American Mathematical Society (N.S.), 27(1):1–67, 1992
work page 1992
-
[11]
F. del Teso and E. Lindgren, Finite difference schemes for the parabolicp-Laplace equation,SeMA Journal, 80(4):527–547, 2022
work page 2022
-
[12]
F. del Teso, M. Medina, and P. Ochoa, Higher-order asymptotic expansions and finite difference schemes for the fractionalp-Laplacian.Math. Ann. 390:157–203, 2024
work page 2024
-
[13]
F. del Teso and J. D. Rossi. Game Theoretical Asymptotic Mean Value Properties for Non-Homogeneous p-Laplace Problems.Calculus of Variations and Partial Differential Equations. 65, art. 13, 2026
work page 2026
-
[14]
F. del Teso, J. D. Rossi and J. Ruiz-Cases. A Viscosity Framework for Dynamic Programming Principles and Applications.https://arxiv.org/abs/2602.09946 [math.AP], 10 Fec 2026
-
[15]
DiBenedetto.Degenerate parabolic equations
E. DiBenedetto.Degenerate parabolic equations. Universitext. Springer-Verlag, New York, 1993
work page 1993
-
[16]
E. Fabes and N. Garofalo; Mean Value Properties of Solutions to Parabolic Equations with Variable Coefficients.J. Math. Anal. Appl., 121:305–316, 1987
work page 1987
-
[17]
Fulks; A mean value theorem for the heat equation,Proc
W. Fulks; A mean value theorem for the heat equation,Proc. Amer. Math. Soc., 17:6–11, 1966
work page 1966
-
[18]
R. Hynd and E. Lindgren. A doubly nonlinear evolution for the optimal Poincaré inequality.Calc. Var. Partial Differential Equations, 55(4): 2016
work page 2016
-
[19]
R. Hynd and E. Lindgren, Approximation of the least Rayleigh quotient for degreephomogeneous func- tionals,J. Funct. Anal., 272(12):4873–4918, 2017
work page 2017
- [20]
-
[21]
T. Kilpeläinen and P. Lindqvist, On the Dirichlet boundary value problem for a degenerate parabolic equation,SIAM J. Math. Anal., 27(3):661–683, 1996. 37
work page 1996
- [22]
-
[23]
Ü. Kuran. On the mean-value property of harmonic functions,Bull. London Math. Soc., 4:311–312, 1972
work page 1972
-
[24]
E. Le Gruyer, and J.C. Archer. Harmonious Extensions,Siam J. Math. Anal.29(1):279–292, 1998
work page 1998
- [25]
-
[26]
A course on tug-of-war games with random noise
M.Lewicka. A course on tug-of-war games with random noise. Mathematical theory. Universitext, Intro- duction and basic constructions, Springer, Cham, 2020
work page 2020
-
[27]
W. Littman, G. Stampacchia, and H.F. Weinberger. Regular points for elliptic equations with discontinuous coefficients,Ann. Scuola. Norm. Sup. Pisa. Cl. Sci., 17(1–2):43–77, 1963
work page 1963
-
[28]
J. J. Manfredi, M. Parviainen and J. D. Rossi. An asymptotic mean value characterization forp-harmonic functions.Proc. Amer. Math. Soc., 138(3): 881–889, 2010
work page 2010
-
[29]
J. J. Manfredi, M. Parviainen, and J. D. Rossi. An asymptotic mean value characterization for a class of nonlinearparabolic equations related to tug-of-war games.SIAM Jour. Math. Anal.42(5):2058–2081, 2010
work page 2058
- [30]
- [31]
-
[32]
Y. Peres, and S. Sheffield; Tug-of-war with noise: a game theoretic view of thep-Laplacian,Duke Math. J.,145(1):91–120, 2008
work page 2008
-
[33]
Privaloff; Sur les fonctions harmoniques,Mat
I. Privaloff; Sur les fonctions harmoniques,Mat. Sb.32(3):464–471, 1925
work page 1925
-
[34]
N. Suzuki and N.A. Watson, A characterization of heat balls by a mean value property for temperatures, Proc. Amer. Math. Soc., 129(9):2709–2714, 2001
work page 2001
-
[35]
J. L. Vazquez and E. Zuazua, Complexity of large time behaviour of evolution equations with bounded data. Dedicated to the memory of Jacques-Louis Lions.Chinese Ann. Math. Ser. B, 23(2), 293–310, 2002
work page 2002
-
[36]
N. A. Watson; A theory of subtemperatures in several variables,Proc. London Math. Soc., 26(3):385–417, 1973
work page 1973
-
[37]
Williams; Probability with martingales, Cambridge University Press, Cambridge, 1991
D. Williams; Probability with martingales, Cambridge University Press, Cambridge, 1991. Félix del Teso and Carlos Fuertes Morán. Departamento de Matemáticas, Universidad Autónoma de Madrid, ICMAT - Instituto de Ciencias Matemáticas, CSIC-UAM-UC3M-UCM, Campus de Cantoblanco, 28049 Madrid, Spain. E-mail:felix.delteso@uam.esWeb-page:https://sites.google.com/...
work page 1991
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.