TY - JOUR
T1 - Fast CBC construction of randomly shifted lattice rules achieving O (n -1 + δ) convergence for unbounded integrands over ℝs in weighted spaces with POD weights
AU - Nichols, James A.
AU - Kuo, Frances Y.
PY - 2014/8
Y1 - 2014/8
N2 - This paper provides the theoretical foundation for the component-by- component (CBC) construction of randomly shifted lattice rules that are tailored to integrals over ℝs arising from practical applications. For an integral of the form ∫ℝsf(y)∝j=1φ( yj)dy with a univariate probability density φ, our general strategy is to first map the integral into the unit cube [0,1] s using the inverse of the cumulative distribution function of φ, and then apply quasi-Monte Carlo (QMC) methods. However, the transformed integrand in the unit cube rarely falls within the standard QMC setting of Sobolev spaces of functions with mixed first derivatives. Therefore, a non-standard function space setting for integrands over ℝs, previously considered by Kuo, Sloan, Wasilkowski and Waterhouse (2010), is required for the analysis. Motivated by the needs of three applications, the present paper extends the theory of the aforementioned paper in several non-trivial directions, including a new error analysis for the CBC construction of lattice rules with general non-product weights, the introduction of an unanchored variant for the setting, the use of coordinate-dependent weight functions in the norm, and the strategy for fast CBC construction with POD ("product and order dependent") weights.
AB - This paper provides the theoretical foundation for the component-by- component (CBC) construction of randomly shifted lattice rules that are tailored to integrals over ℝs arising from practical applications. For an integral of the form ∫ℝsf(y)∝j=1φ( yj)dy with a univariate probability density φ, our general strategy is to first map the integral into the unit cube [0,1] s using the inverse of the cumulative distribution function of φ, and then apply quasi-Monte Carlo (QMC) methods. However, the transformed integrand in the unit cube rarely falls within the standard QMC setting of Sobolev spaces of functions with mixed first derivatives. Therefore, a non-standard function space setting for integrands over ℝs, previously considered by Kuo, Sloan, Wasilkowski and Waterhouse (2010), is required for the analysis. Motivated by the needs of three applications, the present paper extends the theory of the aforementioned paper in several non-trivial directions, including a new error analysis for the CBC construction of lattice rules with general non-product weights, the introduction of an unanchored variant for the setting, the use of coordinate-dependent weight functions in the norm, and the strategy for fast CBC construction with POD ("product and order dependent") weights.
KW - General weights
KW - POD weights
KW - Quasi-Monte Carlo methods
KW - Randomly shifted lattice rules
KW - Unbounded integrands
UR - http://www.scopus.com/inward/record.url?scp=84901850673&partnerID=8YFLogxK
U2 - 10.1016/j.jco.2014.02.004
DO - 10.1016/j.jco.2014.02.004
M3 - Article
SN - 0885-064X
VL - 30
SP - 444
EP - 468
JO - Journal of Complexity
JF - Journal of Complexity
IS - 4
ER -