Using box-muller with low discrepancy points

  • Authors:
  • Tim Pillards;Ronald Cools

  • Affiliations:
  • Dept. of Computer Science, Katholieke Universiteit Leuven, Heverlee, Belgium;Dept. of Computer Science, Katholieke Universiteit Leuven, Heverlee, Belgium

  • Venue:
  • ICCSA'06 Proceedings of the 2006 international conference on Computational Science and Its Applications - Volume Part V
  • Year:
  • 2006

Quantified Score

Hi-index 0.00

Visualization

Abstract

To use quasi-Monte Carlo methods, the integral is usually first (implicitly) transformed to the unit cube. Integrals weighted with the multivariate normal density are usually transformed to the unit cube with the inverse of the multivariate normal cumulative distribution function. However, other transformations are possible, amongst which the transformation by Box and Muller. The danger in using a non-separable transformation is that it might break the low discrepancy structure which makes quasi-Monte Carlo converge faster than regular Monte Carlo. We examine several transformations visually, theoretically and practically and show that it is sometimes preferable to use other transformations than the inverse cumulative distribution function.