A Unified Framework for the Error Analysis of Physics-Informed Neural Networks
Marius Zeinhofer, Rami Masri, Kent-André Mardal
TL;DR
This work develops a unified, coercivity based error analysis for physics-informed neural networks (PINNs) solving linear PDEs across elliptic, parabolic, hyperbolic, Stokes, and PDE constrained optimization problems. It introduces an abstract least-squares energy framework and a non-Galerkin Céa type analysis that yields both a priori and a posteriori error estimates by combining energy, duality, and interpolation arguments. A key finding is that $L^2$ penalties for boundary and initial data can degrade the decay rate of the error in natural Sobolev norms, while encoding constraints directly into the neural representation avoids this issue. The authors provide sharp coercivity and continuity results for Poisson, Darcy, elasticity, Stokes, parabolic, hyperbolic, and source recovery problems and corroborate the theory with three- to four-dimensional numerical experiments using recent optimization advances. Together these results offer rigorous guarantees and practical guidance for designing PINNs that solve linear PDEs efficiently with provable accuracy.
Abstract
We prove a priori and a posteriori error estimates for physics-informed neural networks (PINNs) for linear PDEs. We analyze elliptic equations in primal and mixed form, elasticity, parabolic, hyperbolic and Stokes equations; and a PDE constrained optimization problem. For the analysis, we propose an abstract framework in the common language of bilinear forms, and we show that coercivity and continuity lead to error estimates. The obtained estimates are sharp and reveal that the $L^2$ penalty approach for initial and boundary conditions in the PINN formulation weakens the norm of the error decay. Finally, utilizing recent advances in PINN optimization, we present numerical examples that illustrate the ability of the method to achieve accurate solutions.
