20,676 research outputs found
Hot new directions for quasi-Monte Carlo research in step with applications
This article provides an overview of some interfaces between the theory of
quasi-Monte Carlo (QMC) methods and applications. We summarize three QMC
theoretical settings: first order QMC methods in the unit cube and in
, and higher order QMC methods in the unit cube. One important
feature is that their error bounds can be independent of the dimension
under appropriate conditions on the function spaces. Another important feature
is that good parameters for these QMC methods can be obtained by fast efficient
algorithms even when is large. We outline three different applications and
explain how they can tap into the different QMC theory. We also discuss three
cost saving strategies that can be combined with QMC in these applications.
Many of these recent QMC theory and methods are developed not in isolation, but
in close connection with applications
Some Results on the Complexity of Numerical Integration
This is a survey (21 pages, 124 references) written for the MCQMC 2014
conference in Leuven, April 2014. We start with the seminal paper of Bakhvalov
(1959) and end with new results on the curse of dimension and on the complexity
of oscillatory integrals. Some small errors of earlier versions are corrected
Equivalence Theorems in Numerical Analysis : Integration, Differentiation and Interpolation
We show that if a numerical method is posed as a sequence of operators acting
on data and depending on a parameter, typically a measure of the size of
discretization, then consistency, convergence and stability can be related by a
Lax-Richtmyer type equivalence theorem -- a consistent method is convergent if
and only if it is stable. We define consistency as convergence on a dense
subspace and stability as discrete well-posedness. In some applications
convergence is harder to prove than consistency or stability since convergence
requires knowledge of the solution. An equivalence theorem can be useful in
such settings. We give concrete instances of equivalence theorems for
polynomial interpolation, numerical differentiation, numerical integration
using quadrature rules and Monte Carlo integration.Comment: 18 page
Estimation with Numerical Integration on Sparse Grids
For the estimation of many econometric models, integrals without analytical solutions have to be evaluated. Examples include limited dependent variables and nonlinear panel data models. In the case of one-dimensional integrals, Gaussian quadrature is known to work efficiently for a large class of problems. In higher dimensions, similar approaches discussed in the literature are either very specific and hard to implement or suffer from exponentially rising computational costs in the number of dimensions - a problem known as the "curse of dimensionality" of numerical integration. We propose a strategy that shares the advantages of Gaussian quadrature methods, is very general and easily implemented, and does not suffer from the curse of dimensionality. Monte Carlo experiments for the random parameters logit model indicate the superior performance of the proposed method over simulation techniques
Efficient calculation of the worst-case error and (fast) component-by-component construction of higher order polynomial lattice rules
We show how to obtain a fast component-by-component construction algorithm
for higher order polynomial lattice rules. Such rules are useful for
multivariate quadrature of high-dimensional smooth functions over the unit cube
as they achieve the near optimal order of convergence. The main problem
addressed in this paper is to find an efficient way of computing the worst-case
error. A general algorithm is presented and explicit expressions for base~2 are
given. To obtain an efficient component-by-component construction algorithm we
exploit the structure of the underlying cyclic group.
We compare our new higher order multivariate quadrature rules to existing
quadrature rules based on higher order digital nets by computing their
worst-case error. These numerical results show that the higher order polynomial
lattice rules improve upon the known constructions of quasi-Monte Carlo rules
based on higher order digital nets
- …