A superconvergent monte carlo method for multiple integrals on the grid

  • Authors:
  • Sofiya Ivanovska;Emanouil Atanassov;Aneta Karaivanova

  • Affiliations:
  • Institute for Parallel Processing, Bulgarian Academy of Sciences, Sofia, Bulgaria;Institute for Parallel Processing, Bulgarian Academy of Sciences, Sofia, Bulgaria;Institute for Parallel Processing, Bulgarian Academy of Sciences, Sofia, Bulgaria

  • Venue:
  • ICCS'05 Proceedings of the 5th international conference on Computational Science - Volume Part III
  • Year:
  • 2005

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper we present error and performance analysis of a Monte Carlo variance reduction method for solving multidimensional integrals and integral equations. This method combines the idea of separation of the domain into small subdomains with the approach of importance sampling. The importance separation method is originally described in our previous works [7,9]. Here we present a new variant of this method adding polynomial interpolation in subdomains. We also discuss the performance of the algorithms in comparison with crude Monte Carlo. We propose efficient parallel implementation of the importance separation method for a grid environment and we demonstrate numerical experiments on a heterogeneous grid. Two versions of the algorithm are compared – a Monte Carlo version using pseudorandom numbers and a quasi-Monte Carlo version using the Sobol and Halton low-discrepancy sequences [13,8].