Title: Weighted compound integration rules with higher order convergence for all N
Authors: Hickernell, Fred J. *
Kritzer, Peter *
Kuo, Frances Y. *
Nuyens, Dirk * # ×
Issue Date: 2012
Publisher: Springer New York LLC
Series Title: Numerical Algorithms vol:59 issue:2 pages:161-183
Abstract: Quasi-Monte Carlo integration rules, which are equal-weight sample averages of function values, have been popular for approximating multivariate integrals due to their superior convergence rate of order close to 1/N or better, compared to the order 1/sqrt(N) of simple Monte Carlo algorithms. For practical applications, it is desirable to be able to increase the total number of sampling points N one or several at a time until a desired accuracy is met, while keeping all existing evaluations. We show that although a convergence rate of order close to 1/N can be achieved for all values of N (e.g., by using a good lattice sequence), it is impossible to get better than order 1/N convergence for all values of N by adding equally-weighted sampling points in this manner. We then prove that a convergence of order N^(-α) for α > 1 can be achieved by weighting the sampling points, that is, by using a weighted compound integration rule. We apply our theory to lattice sequences and present some numerical results. The same theory also applies to digital sequences.
ISSN: 1017-1398
Publication status: published
KU Leuven publication type: IT
Appears in Collections:NUMA, Numerical Analysis and Applied Mathematics Section
* (joint) first author
× corresponding author
# (joint) last author

Files in This Item:
File Description Status SizeFormat
stop-anywhere.pdfPreprint Published 740KbAdobe PDFView/Open


All items in Lirias are protected by copyright, with all rights reserved.

© Web of science