Citation:
Mancini, S., Bernal, F., Acebrón, J. A. (2016). An Efficient Algorithm for Accelerating Monte Carlo Approximations of the Solution to Boundary Value Problems. Journal of Scientific Computing, 66(2), 577–597.
Sponsor:
This work was supported by national funds through Fundação para a Ciência e a Tecnologia (FCT) under grants UID/CEC/50021/2013, and PTDC/EIA-CCO/098910/2008. FB also acknowledges FCT funding under grant SFRH/BPD/79986/2011.
Keywords:
Monte Carlo method
,
Romberg extrapolation
,
Bounded difusión
,
Feynman-Kac formula
,
First exit time
,
Parallel computing
The numerical approximation of boundary value problems by means of a probabilistic representations often has the drawback that the Monte Carlo estimate of the solution is substantially biased due to the presence of the domain boundary. We introduce a scheme, wThe numerical approximation of boundary value problems by means of a probabilistic representations often has the drawback that the Monte Carlo estimate of the solution is substantially biased due to the presence of the domain boundary. We introduce a scheme, which we have called the leading-term Monte Carlo regression, which seeks to remove that bias by replacing a ’cloud’ of Monte Carlo estimates—carried out at different discretization levels—for the usual single Monte Carlo estimate. The practical result of our scheme is an acceleration of the Monte Carlo method. Theoretical analysis of the proposed scheme, confirmed by numerical experiments, shows that the achieved speedup can be well over 100.[+][-]