π€ AI Summary
To address the low sampling efficiency and unreliable accuracy in Physics-Informed Neural Networks (PINNs), this paper proposes an adaptive sampling method grounded in the Hessian matrix of the PDE residual. The method establishes, for the first time, a provable integral approximation error bound for PINN sampling strategies. By integrating the geometric structure encoded in the residualβs Hessian with a novel weighted numerical quadrature scheme, it overcomes theoretical limitations inherent in uniform or heuristic sampling approaches. Leveraging automatic differentiation and physics-informed modeling, the proposed framework significantly enhances solution accuracy and convergence speed across multiple 1D and 2D PDE benchmarks. Under identical computational cost, it reduces the PDE residual by 30β50% while providing rigorous theoretical guarantees on approximation quality.
π Abstract
Despite considerable scientific advances in numerical simulation, efficiently solving PDEs remains a complex and often expensive problem. Physics-informed Neural Networks (PINN) have emerged as an efficient way to learn surrogate solvers by embedding the PDE in the loss function and minimizing its residuals using automatic differentiation at so-called collocation points. Originally uniformly sampled, the choice of the latter has been the subject of recent advances leading to adaptive sampling refinements for PINNs. In this paper, leveraging a new quadrature method for approximating definite integrals, we introduce a provably accurate sampling method for collocation points based on the Hessian of the PDE residuals. Comparative experiments conducted on a set of 1D and 2D PDEs demonstrate the benefits of our method.