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

被引:0
|
作者
Ronald Cools
Frances Y. Kuo
Dirk Nuyens
机构
[1] K.U. Leuven,Department of Computer Science
[2] University of New South Wales,School of Mathematics and Statistics
来源
Computing | 2010年 / 87卷
关键词
Multivariate numerical integration; Lattice rules; Degree of exactness; 65D30 (Numerical integration); 65D32 (Quadrature and cubature formulas);
D O I
暂无
中图分类号
学科分类号
摘要
Recall that an integration rule is said to have a trigonometric degree of exactness m if it integrates exactly all trigonometric polynomials of degree ≤ m. In this paper we focus on high dimensions, say, d ≫ 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
页数:26
相关论文
共 38 条