Where are the logs?
share
The commonly quoted error rates for QMC integration with an infinite low discrepancy sequence is O(n^-1log(n)^r) with r=d for extensible sequences and r=d-1 otherwise. Such rates hold uniformly over all d dimensional integrands of Hardy-Krause variation one when using n evaluation points. Implicit in those bounds is that for any sequence of QMC points, the integrand can be chosen to depend on n. In this paper we show that rates with any r<(d-1)/2 can hold when f is held fixed as n→∞. This is accomplished following a suggestion of Erich Novak to use some unpublished results of Trojan from the 1980s as given in the information based complexity monograph of Traub, Wasilkowski and Woźniakowski. The proof is made by applying a technique of Roth with the theorem of Trojan. The proof is non constructive and we do not know of any integrand of bounded variation in the sense of Hardy and Krause for which the QMC error exceeds (log n)^1+ϵ/n for infinitely many n when using a digital sequence such as one of Sobol's. An empirical search when d=2 for integrands designed to exploit known weaknesses in certain point sets showed no evidence that r>1 is needed. An example with d=3 and n up to 2^100 might possibly require r>1.
READ FULL TEXT