Exponentially accurate spectral Monte Carlo method for linear PDEs and their error estimates

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Feng, Jiaying, Sheng, Changtao, Xu, Chenglong
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866929723886534656
author Feng, Jiaying
Sheng, Changtao
Xu, Chenglong
author_facet Feng, Jiaying
Sheng, Changtao
Xu, Chenglong
contents This paper introduces a spectral Monte Carlo iterative method (SMC) for solving linear Poisson and parabolic equations driven by $α$-stable Lévy process with $α\in (0,2)$, which was initially proposed and developed by Gobet and Maire in their pioneering works (Monte Carlo Methods Appl 10(3-4), 275--285, 2004, and SIAM J Numer Anal 43(3), 1256--1275, 2005) for the case $α=2$. The novel method effectively integrates multiple computational techniques, including the interpolation based on generalized Jacobi functions (GJFs), space-time spectral methods, control variates techniques, and a novel walk-on-sphere method (WOS). The exponential convergence of the error bounds is rigorously established through finite iterations for both Poisson and parabolic equations involving the integral fractional Laplacian operator. Remarkably, the proposed space-time spectral Monte Carlo method (ST-SMC) for the parabolic equation is unified for both $α\in(0,2)$ and $α=2$. Extensive numerical results are provided to demonstrate the spectral accuracy and efficiency of the proposed method, thereby validating the theoretical findings.
format Preprint
id arxiv_https___arxiv_org_abs_2502_15123
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Exponentially accurate spectral Monte Carlo method for linear PDEs and their error estimates
Feng, Jiaying
Sheng, Changtao
Xu, Chenglong
Numerical Analysis
65N35, 65C05, 60G52, 65M15, 33C45
G.1
This paper introduces a spectral Monte Carlo iterative method (SMC) for solving linear Poisson and parabolic equations driven by $α$-stable Lévy process with $α\in (0,2)$, which was initially proposed and developed by Gobet and Maire in their pioneering works (Monte Carlo Methods Appl 10(3-4), 275--285, 2004, and SIAM J Numer Anal 43(3), 1256--1275, 2005) for the case $α=2$. The novel method effectively integrates multiple computational techniques, including the interpolation based on generalized Jacobi functions (GJFs), space-time spectral methods, control variates techniques, and a novel walk-on-sphere method (WOS). The exponential convergence of the error bounds is rigorously established through finite iterations for both Poisson and parabolic equations involving the integral fractional Laplacian operator. Remarkably, the proposed space-time spectral Monte Carlo method (ST-SMC) for the parabolic equation is unified for both $α\in(0,2)$ and $α=2$. Extensive numerical results are provided to demonstrate the spectral accuracy and efficiency of the proposed method, thereby validating the theoretical findings.
title Exponentially accurate spectral Monte Carlo method for linear PDEs and their error estimates
topic Numerical Analysis
65N35, 65C05, 60G52, 65M15, 33C45
G.1
url https://arxiv.org/abs/2502.15123