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.

Kamm K., Pagliarani S., Pascucci A. (2023). Numerical solution of kinetic SPDEs via stochastic Magnus expansion. MATHEMATICS AND COMPUTERS IN SIMULATION, 207, 189-208 [10.1016/j.matcom.2022.12.029].

Numerical solution of kinetic SPDEs via stochastic Magnus expansion

Kamm K.
;
Pagliarani S.;Pascucci A.
2023

Abstract

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.
2023
Kamm K., Pagliarani S., Pascucci A. (2023). Numerical solution of kinetic SPDEs via stochastic Magnus expansion. MATHEMATICS AND COMPUTERS IN SIMULATION, 207, 189-208 [10.1016/j.matcom.2022.12.029].
Kamm K.; Pagliarani S.; Pascucci A.
File in questo prodotto:
File Dimensione Formato  
1-s2.0-S037847542200516X-main.pdf

accesso aperto

Tipo: Versione (PDF) editoriale
Licenza: Licenza per Accesso Aperto. Creative Commons Attribuzione (CCBY)
Dimensione 1.06 MB
Formato Adobe PDF
1.06 MB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11585/916049
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 1
  • ???jsp.display-item.citation.isi??? 1
social impact