ResearchTrend.AI
  • Papers
  • Communities
  • Events
  • Blog
  • Pricing
Papers
Communities
Social Events
Terms and Conditions
Pricing
Parameter LabParameter LabTwitterGitHubLinkedInBlueskyYoutube

© 2025 ResearchTrend.AI, All rights reserved.

  1. Home
  2. Papers
  3. 1007.0842
134
47
v1v2v3v4 (latest)

Higher order scrambled digital nets achieve the optimal rate of the root mean square error for smooth integrands

6 July 2010
J. Dick
ArXiv (abs)PDFHTML
Abstract

We study numerical approximations of integrals ∫[0,1]sf(\bsx) d\bsx\int_{[0,1]^s} f(\bsx) \,\mathrm{d} \bsx∫[0,1]s​f(\bsx)d\bsx by averaging the function at some sampling points. Monte Carlo (MC) sampling yields a convergence of the root mean square error (RMSE) of order N−1/2N^{-1/2}N−1/2 (where NNN is the number of samples). Quasi-Monte Carlo (QMC) sampling on the other hand achieves a convergence of order N−1+εN^{-1+\varepsilon}N−1+ε, for any ε>0\varepsilon >0ε>0. Randomized QMC (RQMC), a combination of MC and QMC, achieves a RMSE of order N−3/2+εN^{-3/2+\varepsilon}N−3/2+ε. A combination of RQMC with local antithetic sampling achieves a convergence of the RMSE of order N−3/2−1/s+εN^{-3/2-1/s+\varepsilon}N−3/2−1/s+ε (where s≥1s \ge 1s≥1 is the dimension). QMC, RQMC and RQMC with local antithetic sampling require that the integrand has some smoothness (for instance, bounded variation). Stronger smoothness assumptions on the integrand do not improve the convergence of the above algorithms further. This paper introduces a new RQMC algorithm, for which we prove that it achieves a convergence of the RMSE of order N−α−1/2+εN^{-\alpha-1/2+\varepsilon}N−α−1/2+ε if the integrand has square integrable partial mixed derivatives up to order α\alphaα in each variable. Known lower bounds show that this rate of convergence cannot be improved. We provide numerical examples for which the RMSE converges approximately with order N−5/2N^{-5/2}N−5/2 and N−7/2N^{-7/2}N−7/2, in accordance with the theoretical upper bound.

View on arXiv
Comments on this paper