Full error analysis of the random deep splitting method for nonlinear parabolic PDEs and PIDEs
Ariel Neufeld, Philipp Schmocker, Sizhou Wu
TL;DR
The paper develops a randomized deep splitting framework that uses random neural networks to solve high-dimensional nonlinear PDEs and PIDEs with jumps, including infinite-activity cases. It provides a full error analysis proving convergence to the unique viscosity solution and derives explicit bounds for both deterministic and random-network variants, with the random version offering controlled generalization error. Numerical experiments in pricing derivatives under default risk demonstrate the method scales to dimensions up to $d=10^4$, delivering solutions in seconds and outperforming deterministic DS and MLP approaches, especially for PIDEs with nonlocal terms. The results establish a scalable, provably convergent approach for high-dimensional nonlinear PIDE problems, combining rigorous theory with practical efficiency for financial engineering applications.
Abstract
In this paper, we present a randomized extension of the deep splitting algorithm introduced in [Beck, Becker, Cheridito, Jentzen, and Neufeld (2021)] using random neural networks suitable to approximately solve both high-dimensional nonlinear parabolic PDEs and PIDEs with jumps having (possibly) infinite activity. We provide a full error analysis of our so-called random deep splitting method. In particular, we prove that our random deep splitting method converges to the (unique viscosity) solution of the nonlinear PDE or PIDE under consideration. Moreover, we empirically analyze our random deep splitting method by considering several numerical examples including both nonlinear PDEs and nonlinear PIDEs relevant in the context of pricing of financial derivatives under default risk. In particular, we empirically demonstrate in all examples that our random deep splitting method can approximately solve nonlinear PDEs and PIDEs in 10'000 dimensions within seconds.
