A spectral Monte Carlo method for the Poisson equation

Using a sequential Monte Carlo algorithm, we compute a spectral approximation of the solution of the Poisson equation in dimension 1 and 2. The Feyman-Kac computation of the pointwise solution is achieved using either an integral representation or a modified walk on spheres method. The variances decrease geometrically with the number of steps. A global solution is obtained, accurate up to the interpolation error. Surprisingly, the accuracy depends very little on the absorption layer thickness of the walk on spheres.

[1]  C. Bernardi,et al.  Approximations spectrales de problèmes aux limites elliptiques , 2003 .

[2]  J. Halton.  Sequential Monte Carlo , 1962, Mathematical Proceedings of the Cambridge Philosophical Society.

[3]  A. Friedman.  Stochastic Differential Equations and Applications , 1975 .

[4]  Harald Niederreiter,et al.  Random number generation and Quasi-Monte Carlo methods , 1992, CBMS-NSF regional conference series in applied mathematics.

[5]  S. Maire.  Reducing variance using iterated control variates , 2003 .

[6]  M. Freidlin.  Functional Integration And Partial Differential Equations , 1985 .

[7]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[8]  Chi-Ok Hwang,et al.  A Feynman-Kac path-integral implementation for Poisson's equation using an h-conditioned Green's function , 2003, Math. Comput. Simul..

[9]  Karl K. Sabelfeld.  Monte Carlo Methods in Boundary Value Problems. , 1991 .

[10]  T. E. Booth,et al.  Exponential convergence for Monte Carlo particle transport , 1985 .

[11]  D. Cox,et al.  Exponential convergence of adaptive importance sampling for Markov chains , 2000 .

[12]  Robert Dautray.  Méthodes probabilistes pour les équations de la physique , 1989 .

[13]  Paolo Baldi,et al.  Exact asymptotics for the probability of exit from a domain and applications to simulation , 1995 .

[14]  E. Gobet.  Euler schemes and half-space approximation for the simulation of diffusion in a domain , 2001 .

[15]  Åke Björck,et al.  Numerical methods for least square problems , 1996 .

[16]  M. E. Muller.  Some Continuous Monte Carlo Methods for the Dirichlet Problem , 1956 .

[17]  Christoph W. Ueberhuber,et al.  Computational Integration , 2018, An Introduction to Scientific, Symbolic, and Graphical Computation.