Results 1  10
of
82
QuasiRandom Sequences and Their Discrepancies
 SIAM J. Sci. Comput
, 1994
"... Quasirandom (also called low discrepancy) sequences are a deterministic alternative to random sequences for use in Monte Carlo methods, such as integration and particle simulations of transport processes. The error in uniformity for such a sequence of N points in the sdimensional unit cube is meas ..."
Abstract

Cited by 73 (6 self)
 Add to MetaCart
Quasirandom (also called low discrepancy) sequences are a deterministic alternative to random sequences for use in Monte Carlo methods, such as integration and particle simulations of transport processes. The error in uniformity for such a sequence of N points in the sdimensional unit cube is measured by its discrepancy, which is of size (log N) s N \Gamma1 for large N , as opposed to discrepancy of size (log log N) 1=2 N \Gamma1=2 for a random sequence (i.e. for almost any randomlychosen sequence). Several types of discrepancy, one of which is new, are defined and analyzed. A critical discussion of the theoretical bounds on these discrepancies is presented. Computations of discrepancy are presented for a wide choice of dimension s, number of points N and different quasirandom sequences. In particular for moderate or large s, there is an intermediate regime in which the discrepancy of a quasirandom sequence is almost exactly the same as that of a randomly chosen sequence...
Latin Supercube Sampling for Very High Dimensional Simulations
, 1997
"... This paper introduces Latin supercube sampling (LSS) for very high dimensional simulations, such as arise in particle transport, finance and queuing. LSS is developed as a combination of two widely used methods: Latin hypercube sampling (LHS), and QuasiMonte Carlo (QMC). In LSS, the input variables ..."
Abstract

Cited by 69 (7 self)
 Add to MetaCart
This paper introduces Latin supercube sampling (LSS) for very high dimensional simulations, such as arise in particle transport, finance and queuing. LSS is developed as a combination of two widely used methods: Latin hypercube sampling (LHS), and QuasiMonte Carlo (QMC). In LSS, the input variables are grouped into subsets, and a lower dimensional QMC method is used within each subset. The QMC points are presented in random order within subsets. QMC methods have been observed to lose effectiveness in high dimensional problems. This paper shows that LSS can extend the benefits of QMC to much higher dimensions, when one can make a good grouping of input variables. Some suggestions for grouping variables are given for the motivating examples. Even a poor grouping can still be expected to do as well as LHS. The paper also extends LHS and LSS to infinite dimensional problems. The paper includes a survey of QMC methods, randomized versions of them (RQMC) and previous methods for extending Q...
Computer Experiments
, 1996
"... Introduction Deterministic computer simulations of physical phenomena are becoming widely used in science and engineering. Computers are used to describe the flow of air over an airplane wing, combustion of gasses in a flame, behavior of a metal structure under stress, safety of a nuclear reactor, a ..."
Abstract

Cited by 67 (5 self)
 Add to MetaCart
Introduction Deterministic computer simulations of physical phenomena are becoming widely used in science and engineering. Computers are used to describe the flow of air over an airplane wing, combustion of gasses in a flame, behavior of a metal structure under stress, safety of a nuclear reactor, and so on. Some of the most widely used computer models, and the ones that lead us to work in this area, arise in the design of the semiconductors used in the computers themselves. A process simulator starts with a data structure representing an unprocessed piece of silicon and simulates the steps such as oxidation, etching and ion injection that produce a semiconductor device such as a transistor. A device simulator takes a description of such a device and simulates the flow of current through it under varying conditions to determine properties of the device such as its switching speed and the critical voltage at which it switches. A circuit simulator takes a list of devices and the
Recent Advances In Randomized QuasiMonte Carlo Methods
"... We survey some of the recent developments on quasiMonte Carlo (QMC) methods, which, in their basic form, are a deterministic counterpart to the Monte Carlo (MC) method. Our main focus is the applicability of these methods to practical problems that involve the estimation of a highdimensional inte ..."
Abstract

Cited by 59 (12 self)
 Add to MetaCart
We survey some of the recent developments on quasiMonte Carlo (QMC) methods, which, in their basic form, are a deterministic counterpart to the Monte Carlo (MC) method. Our main focus is the applicability of these methods to practical problems that involve the estimation of a highdimensional integral. We review several QMC constructions and dierent randomizations that have been proposed to provide unbiased estimators and for error estimation. Randomizing QMC methods allows us to view them as variance reduction techniques. New and old results on this topic are used to explain how these methods can improve over the MC method in practice. We also discuss how this methodology can be coupled with clever transformations of the integrand in order to reduce the variance further. Additional topics included in this survey are the description of gures of merit used to measure the quality of the constructions underlying these methods, and other related techniques for multidimensional integration. 1 2 1.
QuasiMonte Carlo Integration
 JOURNAL OF COMPUTATIONAL PHYSICS
, 1995
"... The standard Monte Carlo approach to evaluating multidimensional integrals using (pseudo)random integration nodes is frequently used when quadrature methods are too difficult or expensive to implement. As an alternative to the random methods, it has been suggested that lower error and improved con ..."
Abstract

Cited by 42 (6 self)
 Add to MetaCart
The standard Monte Carlo approach to evaluating multidimensional integrals using (pseudo)random integration nodes is frequently used when quadrature methods are too difficult or expensive to implement. As an alternative to the random methods, it has been suggested that lower error and improved convergence may be obtained by replacing the pseudorandom sequences with more uniformly distributed sequences known as quasirandom. In this paper the Halton, Sobol' and Faure quasirandom sequences are compared in computational experiments designed to determine the effects on convergence of certain properties of the integrand, including variance, variation, smoothness and dimension. The results show that variation, which plays an important role in the theoretical upper bound given by the KoksmaHlawka inequality, does not affect convergence; while variance, the determining factor in random Monte Carlo, is shown to provide a rough upper bound, but does not accurately predict performance. In ge...
Monte Carlo and QuasiMonte Carlo methods
 Acta Numerica
, 1998
"... Monte Carlo is one of the most versatile and widely used numerical methods. Its convergence rate, O(N ~ 1 ^ 2), is independent of dimension, which shows Monte Carlo to be very robust but also slow. This article presents an introduction to Monte Carlo methods for integration problems, including conve ..."
Abstract

Cited by 35 (1 self)
 Add to MetaCart
Monte Carlo is one of the most versatile and widely used numerical methods. Its convergence rate, O(N ~ 1 ^ 2), is independent of dimension, which shows Monte Carlo to be very robust but also slow. This article presents an introduction to Monte Carlo methods for integration problems, including convergence theory, sampling methods and variance reduction techniques. Accelerated convergence for Monte Carlo quadrature is attained using quasirandom (also called lowdiscrepancy) sequences, which are a deterministic alternative to random or pseudorandom sequences. The points in a quasirandom sequence are correlated to provide greater uniformity. The resulting quadrature method, called quasiMonte Carlo, has a convergence rate of approximately O((log N^N ' 1). For quasiMonte Carlo, both theoretical error estimates and practical limitations are presented. Although the emphasis in this article is on integration, Monte Carlo simulation of rarefied gas dynamics is also discussed. In the limit of small mean free path (that is, the fluid dynamic limit), Monte Carlo loses its effectiveness because the collisional distance is much less
Extensible Lattice Sequences For QuasiMonte Carlo Quadrature
 SIAM Journal on Scientific Computing
, 1999
"... Integration lattices are one of the main types of low discrepancy sets used in quasiMonte Carlo methods. However, they have the disadvantage of being of fixed size. This article describes the construction of an infinite sequence of points, the first b m of which form a lattice for any nonnegative ..."
Abstract

Cited by 29 (5 self)
 Add to MetaCart
Integration lattices are one of the main types of low discrepancy sets used in quasiMonte Carlo methods. However, they have the disadvantage of being of fixed size. This article describes the construction of an infinite sequence of points, the first b m of which form a lattice for any nonnegative integer m. Thus, if the quadrature error using an initial lattice is too large, the lattice can be extended without discarding the original points. Generating vectors for extensible lattices are found by minimizing a loss function based on some measure of discrepancy or nonuniformity of the lattice. The spectral test used for finding pseudorandom number generators is one important example of such a discrepancy. The performance of the extensible lattices proposed here is compared to that of other methods for some practical quadrature problems.
Monte Carlo Variance of Scrambled Net Quadrature
 SIAM J. Numer. Anal
, 1997
"... . Hybrids of equidistribution and Monte Carlo methods of integration can achieve the superior accuracy of the former while allowing the simple error estimation methods of the latter. This paper studies the variance of one such hybrid, scrambled nets, by applying a multidimensional multiresolution (w ..."
Abstract

Cited by 28 (1 self)
 Add to MetaCart
. Hybrids of equidistribution and Monte Carlo methods of integration can achieve the superior accuracy of the former while allowing the simple error estimation methods of the latter. This paper studies the variance of one such hybrid, scrambled nets, by applying a multidimensional multiresolution (wavelet) analysis to the integrand. The integrand is assumed to be measurable and square integrable but not necessarily of bounded variation. In simple Monte Carlo, every nonconstant term of the multiresolution contributes to the variance of the estimated integral. For scrambled nets, certain lowdimensional and coarse terms do not contribute to the variance. For any integrand in L 2 , the sampling variance tends to zero faster under scrambled net quadrature than under Monte Carlo sampling, as the number of function evaluations n tends to infinity. Some finite n results bound the variance under scrambled net quadrature by a small constant multiple of the Monte Carlo variance, uniformly ove...
On the use of low discrepancy sequences in Monte Carlo methods
 MONTE CARLO METHODS AND APPLICATIONS
, 1996
"... Quasirandom (or low discrepancy) sequences are sequences for which the convergence to the uniform distribution on [0; 1) s occurs rapidly. Such sequences are used in quasiMonte Carlo methods for which the convergence speed, with respect to the N first terms of the sequence, is in O(N \Gamma1 ..."
Abstract

Cited by 28 (1 self)
 Add to MetaCart
Quasirandom (or low discrepancy) sequences are sequences for which the convergence to the uniform distribution on [0; 1) s occurs rapidly. Such sequences are used in quasiMonte Carlo methods for which the convergence speed, with respect to the N first terms of the sequence, is in O(N \Gamma1 (ln N) s ), where s is the mathematical dimension of the problem considered. The disadvantage of these methods is that error bounds, even if they exist theoretically, are inefficient in practice. Nevertheless, to take advantage of these methods for what concerns their convergence speed, we use them as a variance reduction technique, which lead to great improvements with respect to standard Monte Carlo methods. We consider in this paper two different approaches which combine Monte Carlo and quasiMonte Carlo methods. The first one can use every low discrepancy sequence and the second one, called Owen's method, uses only Niederreiter sequences. We prove that the first approach has the same...
Randomized Halton Sequences
 Mathematical and Computer Modelling
, 2000
"... The Halton sequence is a wellknown multidimensional low discrepancy sequence. In this paper, we propose a new method for randomizing the Halton sequence: we randomize the start point of each component of the sequence. This method combines the potential accuracy advantage of Halton sequence in mult ..."
Abstract

Cited by 26 (1 self)
 Add to MetaCart
The Halton sequence is a wellknown multidimensional low discrepancy sequence. In this paper, we propose a new method for randomizing the Halton sequence: we randomize the start point of each component of the sequence. This method combines the potential accuracy advantage of Halton sequence in multidimensional integration with the practical error estimation advantage of Monte Carlo methods. Theoretically, using multiple randomized Halton sequences as a variance reduction technique we can obtain an efficiency improvement over standard Monte Carlo under rather general conditions. Numerical results show that randomized Halton sequences have better performance than not only Monte Carlo, but also randomly shifted Halton sequences and (single long) purely deterministic skipped Halton sequence. Key Words: QuasiMonte Carlo methods, low discrepancy sequences, Monte Carlo methods, numerical integration, variance reduction. AMS 1991 Subject Classification: 65C05, 65D30. This work was suppor...