An \({hp}\)-Version Time Stepping Spectral Monte Carlo Method for Semi-linear Parabolic Equations

Abstract

SIAM Journal on Numerical Analysis, Volume 64, Issue 5, Page 2051-2078, October 2026.
Abstract. In this paper, we present an [math]-version time-stepping spectral Monte Carlo method for solving semilinear parabolic equations. The key innovation lies in constructing an exponentially accurate stochastic algorithm that integrates a residual iteration scheme on Gauss-type nodes in both temporal and spatial directions with a reconstruction strategy rooted in spectral methods. To address the long-time simulations and initial singularities that are often challenging for traditional stochastic algorithms (e.g., walk-on-spheres method), we further develop an [math]-version time-stepping framework that employs multiple time steps and, respectively, geometric time partitions with linearly increasing polynomial degrees to handle these difficulties. Notably, the proposed algorithm bypasses the need to solve linear systems required by traditional spectral methods and remarkably supports parallel computation at both temporal and spatial grid points. We rigorously establish exponential convergence rates for the multistep method within a finite number of iterations. Extensive numerical experiments are conducted to demonstrate the spectral accuracy and computational efficiency of the proposed method in long-time simulations, problems with initial singularities, and a five-dimensional problem, thereby validating the theoretical results.