Finite-order weights imply tractability of multivariate integration

Multivariate integration of high dimension s occurs in many applications. In many such applications, for example in finance, integrands can be well approximated by sums of functions of just a few variables. In this situation the superposition (or effective) dimension is small, and we can model the problem with finite-order weights, where the weights describe the relative importance of each distinct group of variables up to a given order (where the order is the number of variables in a group), and ignore all groups of variables of higher order.In this paper we consider multivariate integration for the anchored and unanchored (nonperiodic) Sobolev spaces equipped with finite-order weights. Our main interest is tractability and strong tractability of QMC algorithms in the worst-case setting. That is, we want to find how the minimal number of function values needed to reduce the initial error by a factor e depends on s and e-1. If there is no dependence on s, and only polynomial dependence on e-1, we have strong tractability, whereas with polynomial dependence on both s and e-1 we have tractability.We show that for the anchored Sobolev space we have strong tractability for arbitrary finite-order weights, whereas for the unanchored Sobolev space we have tractability for all bounded finite-order weights. In both cases, the dependence on e-1 is quadratic. We can improve the dependence on e-1 at the expense of polynomial dependence on s. For finite-order weights, we may achieve almost linear dependence on e-1 with a polynomial dependence on s whose degree is proportional to the order of the weights.We show that these tractability bounds can be achieved by shifted lattice rules with generators computed by the component-by-component (CBC) algorithm. The computed lattice rules depend on the weights. Similar bounds can also be achieved by well-known low discrepancy sequences such as Halton, Sobol and Niederreiter sequences which do not depend on the weights. We prove that these classical low discrepancy sequences lead to error bounds with almost linear dependence on n-1 and polynomial dependence on d. We present explicit worst-case error bounds for shifted lattice rules and for the Niederreiter sequence. Better tractability and error bounds are possible for finite-order weights, and even for general weights if they satisfy certain conditions. We present conditions on general weights that guarantee tractability and strong tractability of multivariate integration.

[1]  Fred J. Hickernell,et al.  On tractability of weighted integration over bounded and unbounded regions in Reals , 2004, Math. Comput..

[2]  J. Halton On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals , 1960 .

[3]  Frances Y. Kuo,et al.  On the step-by-step construction of quasi-Monte Carlo integration rules that achieve strong tractability error bounds in weighted Sobolev spaces , 2002, Math. Comput..

[4]  I. Sobol On the distribution of points in a cube and the approximate evaluation of integrals , 1967 .

[5]  Henryk Wozniakowski,et al.  Tractability of Multivariate Integration for Weighted Korobov Classes , 2001, J. Complex..

[6]  Fred J. Hickernell,et al.  Integration and approximation in arbitrary dimensions , 2000, Adv. Comput. Math..

[7]  Fred J. Hickernell,et al.  A generalized discrepancy and quadrature error bound , 1998, Math. Comput..

[8]  I. H. SLOAN,et al.  Constructing Randomly Shifted Lattice Rules in Weighted Sobolev Spaces , 2002, SIAM J. Numer. Anal..

[9]  I. Sloan Lattice Methods for Multiple Integration , 1994 .

[10]  Henryk Wozniakowski,et al.  Tractability of Integration in Non-periodic and Periodic Weighted Tensor Product Hilbert Spaces , 2002, J. Complex..

[11]  A. Owen,et al.  Valuation of mortgage-backed securities using Brownian bridges to reduce effective dimension , 1997 .

[12]  Ian H. Sloan,et al.  Why Are High-Dimensional Finance Problems Often of Low Effective Dimension? , 2005, SIAM J. Sci. Comput..

[13]  Fred J. Hickernell,et al.  The error bounds and tractability of quasi-Monte Carlo algorithms in infinite dimension , 2002, Math. Comput..

[14]  Henryk Wozniakowski,et al.  Liberating the weights , 2004, J. Complex..

[15]  E. Novak,et al.  Foundations of Computational Mathematics: When are integration and discrepancy tractable? , 2001 .

[16]  Henryk Wozniakowski,et al.  Intractability Results for Integration and Discrepancy , 2001, J. Complex..

[17]  Joseph F. Traub,et al.  Faster Valuation of Financial Derivatives , 1995 .

[18]  Henryk Wozniakowski,et al.  Weighted Tensor Product Algorithms for Linear Multivariate Problems , 1999, J. Complex..

[19]  I. Sloan,et al.  When Does Monte Carlo Depend Polynomially on the Number of Variables , 2004 .

[20]  Kai-Tai Fang,et al.  The effective dimension and quasi-Monte Carlo integration , 2003, J. Complex..

[21]  Harald Niederreiter,et al.  Random number generation and Quasi-Monte Carlo methods , 1992, CBMS-NSF regional conference series in applied mathematics.

[22]  Henryk Wozniakowski,et al.  When Are Quasi-Monte Carlo Algorithms Efficient for High Dimensional Integrals? , 1998, J. Complex..

[23]  Henryk Wozniakowski,et al.  Good Lattice Rules in Weighted Korobov Spaces with General Weights , 2006, Numerische Mathematik.

[24]  Xiaoqun Wang,et al.  Strong tractability of multivariate integration using quasi-Monte Carlo algorithms , 2003, Math. Comput..