Numerical solution of kinetic SPDEs via stochastic Magnus expansion
Fuente:
arXiv
Enregistré dans:
| Auteurs principaux: | , , |
|---|---|
| Format: | Preprint |
| Publié: |
2022
|
| Sujets: | |
| Accès en ligne: | |
| Tags: |
Ajouter un tag
Pas de tags, Soyez le premier à ajouter un tag!
|
| _version_ | 1866911813586649088 |
|---|---|
| author | Kamm, Kevin Pagliarani, Stefano Pascucci, Andrea |
| author_facet | Kamm, Kevin Pagliarani, Stefano Pascucci, Andrea |
| contents | In this paper, we show how the Itô-stochastic Magnus expansion can be used to efficiently solve stochastic partial differential equations (SPDE) with two space variables numerically. To this end, we will first discretize the SPDE in space only by utilizing finite difference methods and vectorize the resulting equation exploiting its sparsity. As a benchmark, we will apply it to the case of the stochastic Langevin equation with constant coefficients, where an explicit solution is available, and compare the Magnus scheme with the Euler-Maruyama scheme. We will see that the Magnus expansion is superior in terms of both accuracy and especially computational time by using a single GPU and verify it in a variable coefficient case. Notably, we will see speed-ups of order ranging form 20 to 200 compared to the Euler-Maruyama scheme, depending on the accuracy target and the spatial resolution. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2207_09776 |
| institution | arXiv |
| publishDate | 2022 |
| record_format | arxiv |
| spellingShingle | Numerical solution of kinetic SPDEs via stochastic Magnus expansion Kamm, Kevin Pagliarani, Stefano Pascucci, Andrea Numerical Analysis Probability In this paper, we show how the Itô-stochastic Magnus expansion can be used to efficiently solve stochastic partial differential equations (SPDE) with two space variables numerically. To this end, we will first discretize the SPDE in space only by utilizing finite difference methods and vectorize the resulting equation exploiting its sparsity. As a benchmark, we will apply it to the case of the stochastic Langevin equation with constant coefficients, where an explicit solution is available, and compare the Magnus scheme with the Euler-Maruyama scheme. We will see that the Magnus expansion is superior in terms of both accuracy and especially computational time by using a single GPU and verify it in a variable coefficient case. Notably, we will see speed-ups of order ranging form 20 to 200 compared to the Euler-Maruyama scheme, depending on the accuracy target and the spatial resolution. |
| title | Numerical solution of kinetic SPDEs via stochastic Magnus expansion |
| topic | Numerical Analysis Probability |
| url | https://arxiv.org/abs/2207.09776 |