pith. sign in

arxiv: 1411.3823 · v2 · pith:77XN7DDOnew · submitted 2014-11-14 · 🧮 math.NA

Open type quasi-Monte Carlo integration based on Halton sequences in weighted Sobolev spaces

classification 🧮 math.NA
keywords errorboldsymboladichaltonintegrationsequencesalgorithmsbounds
0
0 comments X
read the original abstract

In this paper, we study quasi-Monte Carlo (QMC) integration in weighted Sobolev spaces. In contrast to many previous results the QMC algorithms considered here are of open type, i.e., they are extensible in the number of sample points without having to discard the samples already used. As the underlying integration nodes we consider randomized Halton sequences in prime bases $\boldsymbol{p}=(p_1,...,p_s)$ for which we study the root mean square (RMS) worst-case error. The randomization method is a $\boldsymbol{p}$-adic shift which is based on $\boldsymbol{p}$-adic arithmetic. The obtained error bounds are optimal in the order of magnitude of the number of sample nodes. Furthermore we obtain conditions on the coordinate weights under which the error bounds are independent of the dimension $s$. In terms of the field of Information-Based Complexity this means that the corresponding QMC rule achieves a strong polynomial tractability error bound. Our findings on the RMS worst-case error of randomized Halton sequences can be carried over to the RMS $L_2$-discrepancy. Except for the $\boldsymbol{p}$-adic shift our results are fully constructive and no search algorithms (such as the component-by-component algorithm) are required.

This paper has not been read by Pith yet.

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.