The multi-index Monte Carlo method for semilinear stochastic partial differential equations

Fuente: arXiv
Enregistré dans:
Détails bibliographiques
Auteurs principaux: Haji-Ali, Abdul-Lateef, Hoel, Håkon, Petersson, Andreas
Format: Preprint
Publié: 2025
Sujets:
Accès en ligne:
Tags: Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
_version_ 1866917908462960640
author Haji-Ali, Abdul-Lateef
Hoel, Håkon
Petersson, Andreas
author_facet Haji-Ali, Abdul-Lateef
Hoel, Håkon
Petersson, Andreas
contents Stochastic partial differential equations (SPDEs) are often difficult to solve numerically due to their low regularity and high dimensionality. These challenges limit the practical use of computer-aided studies and pose significant barriers to statistical analysis of SPDEs. In this work, we introduce a highly efficient multi-index Monte Carlo method (MIMC) designed to approximate statistics of mild solutions to semilinear parabolic SPDEs. Key to our approach is the proof of a multiplicative convergence property for coupled solutions generated by an exponential integrator numerical solver, which we incorporate with MIMC. We further describe theoretically how the asymptotic computational cost of MIMC can be bounded in terms of the input accuracy tolerance, as the tolerance goes to zero. Notably, our methodology illustrates that for an SPDE with low regularity, MIMC offers substantial performance improvements over other viable methods. Numerical experiments comparing the performance of MIMC with the multilevel Monte Carlo method on relevant test problems validate our theoretical findings. These results also demonstrate that MIMC significantly outperforms state-of-the-art multilevel Monte Carlo, thereby underscoring its potential as a robust and tractable tool for solving semilinear parabolic SPDEs.
format Preprint
id arxiv_https___arxiv_org_abs_2502_00393
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle The multi-index Monte Carlo method for semilinear stochastic partial differential equations
Haji-Ali, Abdul-Lateef
Hoel, Håkon
Petersson, Andreas
Numerical Analysis
Probability
60H15, 65C05, 60H35, 65Y20, 35K58
Stochastic partial differential equations (SPDEs) are often difficult to solve numerically due to their low regularity and high dimensionality. These challenges limit the practical use of computer-aided studies and pose significant barriers to statistical analysis of SPDEs. In this work, we introduce a highly efficient multi-index Monte Carlo method (MIMC) designed to approximate statistics of mild solutions to semilinear parabolic SPDEs. Key to our approach is the proof of a multiplicative convergence property for coupled solutions generated by an exponential integrator numerical solver, which we incorporate with MIMC. We further describe theoretically how the asymptotic computational cost of MIMC can be bounded in terms of the input accuracy tolerance, as the tolerance goes to zero. Notably, our methodology illustrates that for an SPDE with low regularity, MIMC offers substantial performance improvements over other viable methods. Numerical experiments comparing the performance of MIMC with the multilevel Monte Carlo method on relevant test problems validate our theoretical findings. These results also demonstrate that MIMC significantly outperforms state-of-the-art multilevel Monte Carlo, thereby underscoring its potential as a robust and tractable tool for solving semilinear parabolic SPDEs.
title The multi-index Monte Carlo method for semilinear stochastic partial differential equations
topic Numerical Analysis
Probability
60H15, 65C05, 60H35, 65Y20, 35K58
url https://arxiv.org/abs/2502.00393