Constructing lattice rules based on weighted degree of exactness and worst case error

被引:16
作者
Cools, Ronald [1 ]
Kuo, Frances Y. [2 ]
Nuyens, Dirk [1 ,2 ]
机构
[1] Katholieke Univ Leuven, Dept Comp Sci, Louvain, Belgium
[2] Univ New S Wales, Sch Math & Stat, Sydney, NSW, Australia
基金
澳大利亚研究理事会;
关键词
Multivariate numerical integration; Lattice rules; Degree of exactness; BY-COMPONENT CONSTRUCTION; TRIGONOMETRIC DEGREE; MULTIVARIATE INTEGRATION; CUBATURE FORMULAS; CONVERGENCE; ALGORITHMS; KOROBOV; SPACES;
D O I
10.1007/s00607-009-0076-1
中图分类号
TP301 [理论、方法];
学科分类号
081202 ;
摘要
Recall that an integration rule is said to have a trigonometric degree of exactness m if it integrates exactly all trigonometric polynomials of degree a parts per thousand currency sign m. In this paper we focus on high dimensions, say, d a parts per thousand << 6. We introduce three notions of weighted degree of exactness, where we use weights to characterize the anisotropicness of the integrand with respect to successive coordinate directions. Unlike in the classical unweighted setting, the minimal number of integration points needed to achieve a prescribed weighted degree of exactness no longer grows exponentially with d provided that the weights decay sufficiently fast. We present a component-by-component algorithm for the construction of a rank-1 lattice rule such that (i) it has a prescribed weighted degree of exactness, and (ii) its worst case error achieves the optimal rate of convergence in a weighted Korobov space. Then we introduce a modified, more practical, version of this algorithm which maximizes the weighted degree of exactness in each step of the construction. Both algorithms are illustrated by numerical results.
引用
收藏
页码:63 / 89
页数:27
相关论文
共 19 条
[1]  
Becker G.F., 1893, GEOL SOC AM BULL, V4, P13
[2]   Five- and six-dimensional lattice rules generated by structured matrices [J].
Cools, R ;
Govaert, H .
JOURNAL OF COMPLEXITY, 2003, 19 (06) :715-729
[3]   Smolyak's construction of cubature formulas of arbitrary trigonometric degree [J].
Cools, R ;
Novak, E ;
Ritter, K .
COMPUTING, 1999, 62 (02) :147-162
[4]   Different quality indexes for lattice rules [J].
Cools, R ;
Reztsov, A .
JOURNAL OF COMPLEXITY, 1997, 13 (02) :235-258
[5]   Three- and four-dimensional K-optimal lattice rules of moderate trigonometric degree [J].
Cools, R ;
Lyness, JN .
MATHEMATICS OF COMPUTATION, 2001, 70 (236) :1549-1567
[6]  
Cools R., 1997, Acta Numerica, V6, P1, DOI 10.1017/S0962492900002701
[7]   Minimal cubature formulae of trigonometric degree [J].
Cools, R ;
Sloan, IH .
MATHEMATICS OF COMPUTATION, 1996, 65 (216) :1583-1600
[8]   A Belgian view on lattice rules [J].
Cools, Ronald ;
Nuyens, Dirk .
MONTE CARLO AND QUASI-MONTE CARLO METHODS 2006, 2008, :3-21
[9]   On the convergence rate of the component-by-component construction of good lattice rules [J].
Dick, J .
JOURNAL OF COMPLEXITY, 2004, 20 (04) :493-522
[10]   The construction of good extensible rank-1 lattices [J].
Dick, Josef ;
Pillichshammer, Friedrich ;
Waterhouse, Benjamin J. .
MATHEMATICS OF COMPUTATION, 2008, 77 (264) :2345-2373