On efficient weighted integration via a change of variables

• Published in: Numerische Mathematik 2020-11, Vol.146 (3), p.545-570
• Main Authors: Kritzer, P., Pillichshammer, F., Plaskota, L., Wasilkowski, G. W.
• Format: Article
• Language: English
• Subjects: Distribution functions; Mathematical and Computational Engineering; Mathematical and Computational Physics; Mathematical Methods in Physics; Mathematics; Mathematics and Statistics; Monte Carlo simulation; Numerical Analysis; Numerical and Computational Physics; Simulation; Smoothness; Sobolev space; Stress concentration; Theoretical
• ISSN: 0029-599X; 0945-3245
• DOI: 10.1007/s00211-020-01147-7

In this paper, we study the approximation of d -dimensional ϱ -weighted integrals over unbounded domains R + d or R d using a special change of variables, so that quasi-Monte Carlo (QMC) or sparse grid rules can be applied to the transformed integrands over the unit cube. We consider a class of integrands with bounded L p norm of mixed partial derivatives of first order, where p ∈ [ 1 , + ∞ ] . The main results give sufficient conditions on the change of variables ν which guarantee that the transformed integrand belongs to the standard Sobolev space of functions over the unit cube with mixed smoothness of order one. These conditions depend on ϱ and p . The proposed change of variables is in general different than the standard change based on the inverse of the cumulative distribution function. We stress that the standard change of variables leads to integrands over a cube; however, those integrands have singularities which make the application of QMC and sparse grids ineffective. Our conclusions are supported by numerical experiments.