Quasi-random integration in high dimensions

  • Authors:
  • George Takhtamyshev;Bart Vandewoestyne;Ronald Cools

  • Affiliations:
  • Boyer & Associates: Member of the Alpha Millennium Group, 1009 S. Stewart Ave., Lombard IL 60148, USA;Department of Computer Science, KULeuven, Belgium;Department of Computer Science, KULeuven, Belgium

  • Venue:
  • Mathematics and Computers in Simulation
  • Year:
  • 2007

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper, we show that the Sobol' and Richtmyer sequences can be effectively used for numerical integration of functions having up to 1000 variables. The results of integration obtained with the two sequences are compared and the parameters C and @a from the convergence model C/N^@a are estimated, where N is the number of points used. For all the tests done, the Sobol' sequence demonstrated somewhat better convergence, but for many practical values of N the relative error is higher than for Richtmyer sequences due to the large value of C. Constructing Sobol' sequences also takes considerably more time than constructing Richtmyer sequences.