Random and Quasi-Random Point Sets

Size: px
Start display at page:

Download "Random and Quasi-Random Point Sets"

Transcription

1 Peter Hellekalek Gerhard Larcher (Editors) Random and Quasi-Random Point Sets With contributions by JözsefBeck Peter Hellekalek Fred J. Hickemell Gerhard Larcher Pierre L'Ecuyer Harald Niederreiter Shu Tezuka Chaoping Xing Springer

2 Contents Preface v From Probabilistic Diophantine Approximation to Quadratic Fields 1 Jozsef Beck 1 Part I: Super Irregularity 1 2 Part II: Probabilistic Diophantine Approximation Local Case: Inhomogeneous Pell Inequalities - Hyperbolas Beyond Quadratic Irrationals Global Case: Lattice Points in Tilted Rectangles Simultaneous Case 13 3 Part III: Quadratic Fields and Continued Fractions Cesaro Mean of ^{na 1 / 2 } and Quadratic Fields Hardy-Littlewood Lemma Part IV: Class Number One Problems An Attempt to Reduce the Yokoi's Conjecture to a Finite Amount of Computation 27 5 Part V: Cesaro Mean of ^J{na} - 1/2) 32 6 References 46 On the Assessment of Random and Quasi-Random Point Sets 49 Peter Hellekalek 1 Introduction 49 2 Chapter for the Practitioner Assessing RNGs Correlation Analysis for RNGs I Correlation Analysis for RNGs II Theory vs. Practice I: Leap-Frog Streams Theory vs. Practice II: Parallel Monte Carlo Integration Assessing LDPs Good Lattice Points GLPs vs. (t,m,s)-nets Conclusion 69 3 Mathematical Preliminaries Haar and Walsh Series Integration Lattices 74 4 Uniform Distribution Modulo One The Definition of Uniformly Distributed Sequences 76

3 x Contents 4.2 Weyl Sums and Weyl's Criterion Remarks 78 5 The Spectral Test Definition Properties Examples Geometrie Interpretation Remarks 88 6 The Weighted Spectral Test Definition Examples and Properties Remarks 93 7 Discrepancy Definition The Inequality of Erdös-Turän-Koksma Remarks 98 8 Summary 99 9 Acknowledgements References 101 Lattice Rules: How Well Do They Measure Up? 109 Fred J. Hickernell 1 Introduction Some Basic Properties of Lattice Rules A General Approach to Worst-Case and Average-Case Error Analysis Worst-Case Quadrature Error for Reproducing Kernel Hubert Spaces A More General Worst-Case Quadrature Error Analysis Average-Case Quadrature Error Analysis Examples of Other Discrepancies The ANOVA Decomposition A Generalization of P a (L) with Weights The Periodic Bernoulli Discrepancy - Another Generalization of P a (L) The Non-Periodic Bernoulli Discrepancy The Star Discrepancy The Unanchored Discrepancy The Wrap-Around Discrepancy The Symmetrie Discrepancy Shift-Invariant Kernels and Discrepancies Discrepancy Bounds Upper Bounds for P a (L) A Lower Bound on üjr, a,i(p) 145

4 Contents xi 6.3 Quadrature Rules with Different Weights Copy Rules Discrepancies of Integration Lattices and Nets The Expected Discrepancy of Randomized (0, m, s)-nets Infinite Sequences of Embedded Lattices Tractability of High Dimensional Quadrature Quadrature in Arbitrarily High Dimensions The Effective Dimension of an Integrand Discussion and Conclusion References 163 Digital Point Sets: Analysis and Application 167 Gerhard Larcher 1 Introduction The Concept and Basic Properties of Digital Point Sets Discrepancy Bounds for Digital Point Sets Special Classes of Digital Point Sets and Quality Bounds Digital Sequences Based on Formal Laurent Series and Non-Archimedean Diophantine Approximation Analysis of Pseudo-Random-Number Generators by Digital Nets The Digital Lattice Rule Outlook and Open Research Topics References 217 Random Number Generators: Selection Criteria and Testing 223 Pierre L 'Ecuyer and Peter Hellekalek 1 Introduction Design Principles and Figures of Merit A Roulette Wheel Sampling from $ ( The Lattice Structure of MRG's Equidistribution for Regulär Partitions in Cubic Boxes Other Measures of Divergence Empirical Statistical Tests What are the Good Tests? Two-Level Tests Collections of Empirical Tests Examples of Empirical Tests Serial Tests of Equidistribution Tests Based on Close Points in Space Collections of Small RNGs Small Linear Congruential Generators 237

5 xii Contents 5.2 Explicit Inversive Congruential Generators Compound Cubic Congruential Generators Systematic Testing for Small RNGs Serial Tests of Equidistribution for LCGs Serial Tests of Equidistribution for Nonlinear Generators A Summary of the Serial Tests Results Close-Pairs Tests for LCGs Close-Pairs Tests for Nonlinear Generators A Summary of the Close-Pairs Tests Results How Do Real-Life Generators Fare in These Tests? Acknowledgements References 259 Nets, (t, s)-sequences, and Algebraic Geometry 267 Harald Niederreiter and Chaoping Xing 1 Introduction Basic Concepts The Digital Method Background on Algebraic Curves over Finite Fields Construction of (t, s)-sequences New Constructions of (t, m, s)-nets New Algebraic Curves with Many Rational Points References 299 Financial Applications of Monte Carlo and Quasi-Monte Carlo Methods 303 Shu Tezuka 1 Introduction Monte Carlo Methods for Finance Applications Preliminaries for Derivative Pricing Variance Reduction Techniques Caveats for Computer Implementation Speeding Up by Quasi-Monte Carlo Methods What are Quasi-Monte Carlo Methods? Generalized Faure Sequences Numerical Experiments Discussions Future Topics Monte Carlo Simulations for American Options Research Issues Related to Quasi-Monte Carlo Methods References 329

Tutorial on quasi-monte Carlo methods

Tutorial on quasi-monte Carlo methods Tutorial on quasi-monte Carlo methods Josef Dick School of Mathematics and Statistics, UNSW, Sydney, Australia josef.dick@unsw.edu.au Comparison: MCMC, MC, QMC Roughly speaking: Markov chain Monte Carlo

More information

QUASI-MONTE CARLO VIA LINEAR SHIFT-REGISTER SEQUENCES. Pierre L Ecuyer Christiane Lemieux

QUASI-MONTE CARLO VIA LINEAR SHIFT-REGISTER SEQUENCES. Pierre L Ecuyer Christiane Lemieux Proceedings of the 1999 Winter Simulation Conference P. A. Farrington, H. B. Nembhard, D. T. Sturrock, and G. W. Evans, eds. QUASI-MONTE CARLO VIA LINEAR SHIFT-REGISTER SEQUENCES Pierre L Ecuyer Christiane

More information

f(x n ) [0,1[ s f(x) dx.

f(x n ) [0,1[ s f(x) dx. ACTA ARITHMETICA LXXX.2 (1997) Dyadic diaphony by Peter Hellekalek and Hannes Leeb (Salzburg) 1. Introduction. Diaphony (see Zinterhof [13] and Kuipers and Niederreiter [6, Exercise 5.27, p. 162]) is a

More information

Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed?

Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed? Low Discrepancy Sequences in High Dimensions: How Well Are Their Projections Distributed? Xiaoqun Wang 1,2 and Ian H. Sloan 2 1 Department of Mathematical Sciences, Tsinghua University, Beijing 100084,

More information

QUASI-RANDOM NUMBER GENERATION IN THE VIEW OF PARALLEL IMPLEMENTATION. 1. Introduction

QUASI-RANDOM NUMBER GENERATION IN THE VIEW OF PARALLEL IMPLEMENTATION. 1. Introduction Trends in Mathematics Information Center for Mathematical Sciences Volume 6, Number 2, December, 2003, Pages 155 162 QUASI-RANDOM NUMBER GENERATION IN THE VIEW OF PARALLEL IMPLEMENTATION CHI-OK HWANG Abstract.

More information

Pseudo-Random Numbers Generators. Anne GILLE-GENEST. March 1, Premia Introduction Definitions Good generators...

Pseudo-Random Numbers Generators. Anne GILLE-GENEST. March 1, Premia Introduction Definitions Good generators... 14 pages 1 Pseudo-Random Numbers Generators Anne GILLE-GENEST March 1, 2012 Contents Premia 14 1 Introduction 2 1.1 Definitions............................. 2 1.2 Good generators..........................

More information

Random Numbers. Pierre L Ecuyer

Random Numbers. Pierre L Ecuyer 1 Random Numbers Pierre L Ecuyer Université de Montréal, Montréal, Québec, Canada Random numbers generators (RNGs) are available from many computer software libraries. Their purpose is to produce sequences

More information

Improved Discrepancy Bounds for Hybrid Sequences. Harald Niederreiter. RICAM Linz and University of Salzburg

Improved Discrepancy Bounds for Hybrid Sequences. Harald Niederreiter. RICAM Linz and University of Salzburg Improved Discrepancy Bounds for Hybrid Sequences Harald Niederreiter RICAM Linz and University of Salzburg MC vs. QMC methods Hybrid sequences The basic sequences Deterministic discrepancy bounds The proof

More information

Random Number Generation: RNG: A Practitioner s Overview

Random Number Generation: RNG: A Practitioner s Overview Random Number Generation: A Practitioner s Overview Department of Computer Science Department of Mathematics Department of Scientific Computing Florida State University, Tallahassee, FL 32306 USA E-mail:

More information

On the existence of higher order polynomial. lattices with large figure of merit ϱ α.

On the existence of higher order polynomial. lattices with large figure of merit ϱ α. On the existence of higher order polynomial lattices with large figure of merit Josef Dick a, Peter Kritzer b, Friedrich Pillichshammer c, Wolfgang Ch. Schmid b, a School of Mathematics, University of

More information

typical random point set. There are several ways of dening a discrepancy, many of them leading to an error bound of the form je n j V (f)d(p n ) for a

typical random point set. There are several ways of dening a discrepancy, many of them leading to an error bound of the form je n j V (f)d(p n ) for a QUASI-MONTE CARLO VIA LINEAR SHIFT-REGISTER SEQUENCES Pierre L'Ecuyer and Christiane Lemieux Departement d'informatique et de Recherche Operationnelle Universite de Montreal, C.P. 6128, Succ. Centre-Ville

More information

Progress in high-dimensional numerical integration and its application to stochastic optimization

Progress in high-dimensional numerical integration and its application to stochastic optimization Progress in high-dimensional numerical integration and its application to stochastic optimization W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Short course,

More information

On the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules

On the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules On the Behavior of the Weighted Star Discrepancy Bounds for Shifted Lattice Rules Vasile Sinescu and Pierre L Ecuyer Abstract We examine the question of constructing shifted lattice rules of rank one with

More information

Construction of Equidistributed Generators Based on Linear Recurrences Modulo 2

Construction of Equidistributed Generators Based on Linear Recurrences Modulo 2 Construction of Equidistributed Generators Based on Linear Recurrences Modulo 2 Pierre L Ecuyer and François Panneton Département d informatique et de recherche opérationnelle Université de Montréal C.P.

More information

Digital lattice rules

Digital lattice rules Digital lattice rules Multivariate integration and discrepancy estimates A thesis presented to The University of New South Wales in fulfilment of the thesis requirement for the degree of Doctor of Philosophy

More information

Variance and discrepancy with alternative scramblings

Variance and discrepancy with alternative scramblings Variance and discrepancy with alternative scramblings ART B. OWEN Stanford University This paper analyzes some schemes for reducing the computational burden of digital scrambling. Some such schemes have

More information

Uniform Random Number Generators

Uniform Random Number Generators JHU 553.633/433: Monte Carlo Methods J. C. Spall 25 September 2017 CHAPTER 2 RANDOM NUMBER GENERATION Motivation and criteria for generators Linear generators (e.g., linear congruential generators) Multiple

More information

A Coding Theoretic Approach to Building Nets with Well-Equidistributed Projections

A Coding Theoretic Approach to Building Nets with Well-Equidistributed Projections A Coding Theoretic Approach to Building Nets with Well-Equidistributed Projections Yves Edel 1 and Pierre L Ecuyer 2 1 Mathematisches Institut der Universität Im Neuenheimer Feld 288, 69120 Heidelberg,

More information

New Multilevel Algorithms Based on Quasi-Monte Carlo Point Sets

New Multilevel Algorithms Based on Quasi-Monte Carlo Point Sets New Multilevel Based on Quasi-Monte Carlo Point Sets Michael Gnewuch Institut für Informatik Christian-Albrechts-Universität Kiel 1 February 13, 2012 Based on Joint Work with Jan Baldeaux (UTS Sydney)

More information

QUASI-MONTE CARLO METHODS FOR SIMULATION. Pierre L Ecuyer

QUASI-MONTE CARLO METHODS FOR SIMULATION. Pierre L Ecuyer Proceedings of the 2003 Winter Simulation Conference S Chick, P J Sánchez, D Ferrin, and D J Morrice, eds QUASI-MONTE CARLO METHODS FOR SIMULATION Pierre L Ecuyer Département d Informatique et de Recherche

More information

Resolution-Stationary Random Number Generators

Resolution-Stationary Random Number Generators Resolution-Stationary Random Number Generators Francois Panneton Caisse Centrale Desjardins, 1 Complexe Desjardins, bureau 2822 Montral (Québec), H5B 1B3, Canada Pierre L Ecuyer Département d Informatique

More information

Quasi-Monte Carlo methods with applications in finance

Quasi-Monte Carlo methods with applications in finance Finance Stoch (2009) 13: 307 349 DOI 10.1007/s00780-009-0095-y Quasi-Monte Carlo methods with applications in finance Pierre L Ecuyer Received: 30 May 2008 / Accepted: 20 August 2008 / Published online:

More information

Quasi- Monte Carlo Multiple Integration

Quasi- Monte Carlo Multiple Integration Chapter 6 Quasi- Monte Carlo Multiple Integration Introduction In some sense, this chapter fits within Chapter 4 on variance reduction; in some sense it is stratification run wild. Quasi-Monte Carlo methods

More information

An Invitation to Modern Number Theory. Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD

An Invitation to Modern Number Theory. Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD An Invitation to Modern Number Theory Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Foreword Preface Notation xi xiii xix PART 1. BASIC NUMBER THEORY

More information

APPLIED MATHEMATICS REPORT AMR04/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR. J. Dick and F. Pillichshammer

APPLIED MATHEMATICS REPORT AMR04/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR. J. Dick and F. Pillichshammer APPLIED MATHEMATICS REPORT AMR4/2 DIAPHONY, DISCREPANCY, SPECTRAL TEST AND WORST-CASE ERROR J. Dick and F. Pillichshammer January, 24 Diaphony, discrepancy, spectral test and worst-case error Josef Dick

More information

Discrepancy Bounds for Deterministic Acceptance-Rejection Samplers

Discrepancy Bounds for Deterministic Acceptance-Rejection Samplers Discrepancy Bounds for Deterministic Acceptance-Rejection Samplers Houying Zhu Athesisinfulfilmentoftherequirementsforthedegreeof Doctor of Philosophy School of Mathematics and Statistics Faculty of Science

More information

Dyadic diaphony of digital sequences

Dyadic diaphony of digital sequences Dyadic diaphony of digital sequences Friedrich Pillichshammer Abstract The dyadic diaphony is a quantitative measure for the irregularity of distribution of a sequence in the unit cube In this paper we

More information

Kernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.

Kernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt. SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University

More information

Updated Tables of Parameters of (T, M, S)-Nets

Updated Tables of Parameters of (T, M, S)-Nets Updated Tables of Parameters of (T, M, S)-Nets Andrew T. Clayman, 1 K. Mark Lawrence, 2 Gary L. Mullen, 3 Harald Niederreiter, 4 N. J. A. Sloane 5 1 Raytheon Systems Company, State College Operations,

More information

Uniform Distribution and Quasi-Monte Carlo Methods

Uniform Distribution and Quasi-Monte Carlo Methods Workshop on Uniform Distribution and Quasi-Monte Carlo Methods October 14-18, 2013 as part of the Radon Special Semester 2013 on Applications of Algebra and Number Theory Metric number theory, lacunary

More information

A LeVeque-type lower bound for discrepancy

A LeVeque-type lower bound for discrepancy reprinted from Monte Carlo and Quasi-Monte Carlo Methods 998, H. Niederreiter and J. Spanier, eds., Springer-Verlag, 000, pp. 448-458. A LeVeque-type lower bound for discrepancy Francis Edward Su Department

More information

STRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE CARLO ALGORITHMS

STRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE CARLO ALGORITHMS MATHEMATICS OF COMPUTATION Volume 72, Number 242, Pages 823 838 S 0025-5718(02)01440-0 Article electronically published on June 25, 2002 STRONG TRACTABILITY OF MULTIVARIATE INTEGRATION USING QUASI MONTE

More information

2 FRED J. HICKERNELL the sample mean of the y (i) : (2) ^ 1 N The mean square error of this estimate may be written as a sum of two parts, a bias term

2 FRED J. HICKERNELL the sample mean of the y (i) : (2) ^ 1 N The mean square error of this estimate may be written as a sum of two parts, a bias term GOODNESS OF FIT STATISTICS, DISCREPANCIES AND ROBUST DESIGNS FRED J. HICKERNELL Abstract. The Cramer{Von Mises goodness-of-t statistic, also known as the L 2 -star discrepancy, is the optimality criterion

More information

EXPONENTIAL SUMS OVER THE SEQUENCES OF PRN S PRODUCED BY INVERSIVE GENERATORS

EXPONENTIAL SUMS OVER THE SEQUENCES OF PRN S PRODUCED BY INVERSIVE GENERATORS Annales Univ. Sci. Budapest. Sect. Comp. 48 018 5 3 EXPONENTIAL SUMS OVER THE SEQUENCES OF PRN S PRODUCED BY INVERSIVE GENERATORS Sergey Varbanets Odessa Ukraine Communicated by Imre Kátai Received February

More information

Uniform Random Number Generation

Uniform Random Number Generation Uniform Random Number Generation Pierre L Ecuyer Département d Informatique et de Recherche Opérationnelle, Université de Montréal, C.P. 6128, Succ. Centre-Ville, Montréal (Québec), H9S 5B8, Canada. http://www.iro.umontreal.ca/

More information

Incomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators

Incomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators ACTA ARITHMETICA XCIII.4 (2000 Incomplete exponential sums over finite fields and their applications to new inversive pseudorandom number generators by Harald Niederreiter and Arne Winterhof (Wien 1. Introduction.

More information

Constructions of digital nets using global function fields

Constructions of digital nets using global function fields ACTA ARITHMETICA 105.3 (2002) Constructions of digital nets using global function fields by Harald Niederreiter (Singapore) and Ferruh Özbudak (Ankara) 1. Introduction. The theory of (t, m, s)-nets and

More information

Part III. Quasi Monte Carlo methods 146/349

Part III. Quasi Monte Carlo methods 146/349 Part III Quasi Monte Carlo methods 46/349 Outline Quasi Monte Carlo methods 47/349 Quasi Monte Carlo methods Let y = (y,...,y N ) be a vector of independent uniform random variables in Γ = [0,] N, u(y)

More information

uniform distribution theory

uniform distribution theory Uniform Distribution Theory 6 (2011), no.1, 185 200 uniform distribution theory CONSTRUCTIONS OF UNIFORMLY DISTRIBUTED SEQUENCES USING THE b-adic METHOD Peter Hellekalek Harald Niederreiter ABSTRACT. For

More information

Introduction to (randomized) quasi-monte Carlo

Introduction to (randomized) quasi-monte Carlo 1 aft Introduction to (randomized) quasi-monte Carlo Dr Pierre L Ecuyer MCQMC Conference, Stanford University, August 2016 Program Monte Carlo, Quasi-Monte Carlo, Randomized quasi-monte Carlo QMC point

More information

Randomized Quasi-Monte Carlo: An Introduction for Practitioners

Randomized Quasi-Monte Carlo: An Introduction for Practitioners Randomized Quasi-Monte Carlo: An Introduction for Practitioners Pierre L Ecuyer Abstract We survey basic ideas and results on randomized quasi-monte Carlo (RQMC) methods, discuss their practical aspects,

More information

Walsh Series and Transforms

Walsh Series and Transforms Walsh Series and Transforms Theory and Applications by B. Golubov Moscow Institute of Engineering, A. Efimov Moscow Institute of Engineering, and V. Skvortsov Moscow State University, W KLUWER ACADEMIC

More information

Randomized Quasi-Monte Carlo Simulation of Markov Chains with an Ordered State Space

Randomized Quasi-Monte Carlo Simulation of Markov Chains with an Ordered State Space Randomized Quasi-Monte Carlo Simulation of Markov Chains with an Ordered State Space Pierre L Ecuyer 1, Christian Lécot 2, and Bruno Tuffin 3 1 Département d informatique et de recherche opérationnelle,

More information

Polynomial Integration Lattices

Polynomial Integration Lattices Polynomial Integration Lattices Pierre L Ecuyer Département d informatique et de recherche opérationnelle Université de Montréal, C.P. 6128, Succ. Centre-Ville Montréal (Québec), H3C 3J7, Canada lecuyer@iro.umontreal.ca

More information

Generalized Halton Sequences in 2008: A Comparative Study

Generalized Halton Sequences in 2008: A Comparative Study Generalized alton Sequences in 2008: A Comparative Study ENRI FAURE Institut Mathématiques de Luminy and Université Paul Cézanne and CRISTIANE LEMIEUX University of Waterloo alton sequences have always

More information

Optimal Randomized Algorithms for Integration on Function Spaces with underlying ANOVA decomposition

Optimal Randomized Algorithms for Integration on Function Spaces with underlying ANOVA decomposition Optimal Randomized on Function Spaces with underlying ANOVA decomposition Michael Gnewuch 1 University of Kaiserslautern, Germany October 16, 2013 Based on Joint Work with Jan Baldeaux (UTS Sydney) & Josef

More information

NEW SEQUENCES FROM OLD: A SEARCH FOR LOW DISCREPANCY SEQUENCES. Jarvis Haupt

NEW SEQUENCES FROM OLD: A SEARCH FOR LOW DISCREPANCY SEQUENCES. Jarvis Haupt EW SEQUECES FROM OLD: A SEARCH FOR LOW DISCREPACY SEQUECES Jarvis Haupt Electrical and Computer Engineering Department University of Wisconsin, Madison, WI, USA jdhaupt@wisc.edu ABSTRACT Aside from their

More information

Analytic Number Theory

Analytic Number Theory American Mathematical Society Colloquium Publications Volume 53 Analytic Number Theory Henryk Iwaniec Emmanuel Kowalski American Mathematical Society Providence, Rhode Island Contents Preface xi Introduction

More information

Ian Sloan and Lattice Rules

Ian Sloan and Lattice Rules www.oeaw.ac.at Ian Sloan and Lattice Rules P. Kritzer, H. Niederreiter, F. Pillichshammer RICAM-Report 2016-34 www.ricam.oeaw.ac.at Ian Sloan and Lattice Rules Peter Kritzer, Harald Niederreiter, Friedrich

More information

Random Number Generation and Quasi-Monte Carlo

Random Number Generation and Quasi-Monte Carlo Random Number Generation and Quasi-Monte Carlo Pierre L Ecuyer Université de Montréal, Canada, and Inria Rennes, France November 2014 Keywords: random number generator, pseudorandom numbers, linear generator,

More information

Uniform Random Binary Floating Point Number Generation

Uniform Random Binary Floating Point Number Generation Uniform Random Binary Floating Point Number Generation Prof. Dr. Thomas Morgenstern, Phone: ++49.3943-659-337, Fax: ++49.3943-659-399, tmorgenstern@hs-harz.de, Hochschule Harz, Friedrichstr. 57-59, 38855

More information

A Randomized Quasi-Monte Carlo Simulation Method for Markov Chains

A Randomized Quasi-Monte Carlo Simulation Method for Markov Chains A Randomized Quasi-Monte Carlo Simulation Method for Markov Chains Pierre L Ecuyer GERAD and Département d Informatique et de Recherche Opérationnelle Université de Montréal, C.P. 6128, Succ. Centre-Ville,

More information

Numerical Methods in Economics MIT Press, Chapter 9 Notes Quasi-Monte Carlo Methods. Kenneth L. Judd Hoover Institution.

Numerical Methods in Economics MIT Press, Chapter 9 Notes Quasi-Monte Carlo Methods. Kenneth L. Judd Hoover Institution. 1 Numerical Methods in Economics MIT Press, 1998 Chapter 9 Notes Quasi-Monte Carlo Methods Kenneth L. Judd Hoover Institution October 26, 2002 2 Quasi-Monte Carlo Methods Observation: MC uses random sequences

More information

Exploring Monte Carlo Methods

Exploring Monte Carlo Methods Exploring Monte Carlo Methods William L Dunn J. Kenneth Shultis AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO ELSEVIER Academic Press Is an imprint

More information

QMC methods in quantitative finance. and perspectives

QMC methods in quantitative finance. and perspectives , tradition and perspectives Johannes Kepler University Linz (JKU) WU Research Seminar What is the content of the talk Valuation of financial derivatives in stochastic market models using (QMC-)simulation

More information

2 P. L'Ecuyer and R. Simard otherwise perform well in the spectral test, fail this independence test in a decisive way. LCGs with multipliers that hav

2 P. L'Ecuyer and R. Simard otherwise perform well in the spectral test, fail this independence test in a decisive way. LCGs with multipliers that hav Beware of Linear Congruential Generators with Multipliers of the form a = 2 q 2 r Pierre L'Ecuyer and Richard Simard Linear congruential random number generators with Mersenne prime modulus and multipliers

More information

APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H.

APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION. I.H. Sloan, X. Wang and H. APPLIED MATHEMATICS REPORT AMR04/16 FINITE-ORDER WEIGHTS IMPLY TRACTABILITY OF MULTIVARIATE INTEGRATION I.H. Sloan, X. Wang and H. Wozniakowski Published in Journal of Complexity, Volume 20, Number 1,

More information

Numerical Analysis for Statisticians

Numerical Analysis for Statisticians Kenneth Lange Numerical Analysis for Statisticians Springer Contents Preface v 1 Recurrence Relations 1 1.1 Introduction 1 1.2 Binomial CoefRcients 1 1.3 Number of Partitions of a Set 2 1.4 Horner's Method

More information

A Randomized Quasi-Monte Carlo Simulation Method for Markov Chains

A Randomized Quasi-Monte Carlo Simulation Method for Markov Chains A Randomized Quasi-Monte Carlo Simulation Method for Markov Chains Pierre L Ecuyer GERAD and Département d Informatique et de Recherche Opérationnelle Université de Montréal, C.P. 6128, Succ. Centre-Ville,

More information

Tae-Soo Kim and Young-Kyun Yang

Tae-Soo Kim and Young-Kyun Yang Kangweon-Kyungki Math. Jour. 14 (2006), No. 1, pp. 85 93 ON THE INITIAL SEED OF THE RANDOM NUMBER GENERATORS Tae-Soo Kim and Young-Kyun Yang Abstract. A good arithmetic random number generator should possess

More information

Workshop on Heterogeneous Computing, 16-20, July No Monte Carlo is safe Monte Carlo - more so parallel Monte Carlo

Workshop on Heterogeneous Computing, 16-20, July No Monte Carlo is safe Monte Carlo - more so parallel Monte Carlo Workshop on Heterogeneous Computing, 16-20, July 2012 No Monte Carlo is safe Monte Carlo - more so parallel Monte Carlo K. P. N. Murthy School of Physics, University of Hyderabad July 19, 2012 K P N Murthy

More information

PART I INTRODUCTION The meaning of probability Basic definitions for frequentist statistics and Bayesian inference Bayesian inference Combinatorics

PART I INTRODUCTION The meaning of probability Basic definitions for frequentist statistics and Bayesian inference Bayesian inference Combinatorics Table of Preface page xi PART I INTRODUCTION 1 1 The meaning of probability 3 1.1 Classical definition of probability 3 1.2 Statistical definition of probability 9 1.3 Bayesian understanding of probability

More information

arxiv: v2 [math.nt] 27 Mar 2018

arxiv: v2 [math.nt] 27 Mar 2018 arxiv:1710.09313v [math.nt] 7 Mar 018 The Champernowne constant is not Poissonian Ísabel Pirsic, Wolfgang Stockinger Abstract We say that a sequence (x n n N in [0,1 has Poissonian pair correlations if

More information

SOFTWARE FOR UNIFORM RANDOM NUMBER GENERATION: DISTINGUISHING THE GOOD AND THE BAD. Pierre L Ecuyer

SOFTWARE FOR UNIFORM RANDOM NUMBER GENERATION: DISTINGUISHING THE GOOD AND THE BAD. Pierre L Ecuyer SOFTWARE FOR UNIFORM RANDOM NUMBER GENERATION: DISTINGUISHING THE GOOD AND THE BAD Pierre L Ecuyer Département d Informatique et de Recherche Opérationnelle Université demontréal, C.P. 6128, Succ. Centre-Ville

More information

2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized

2 K. ENTACHER 2 Generalized Haar function systems In the following we x an arbitrary integer base b 2. For the notations and denitions of generalized BIT 38 :2 (998), 283{292. QUASI-MONTE CARLO METHODS FOR NUMERICAL INTEGRATION OF MULTIVARIATE HAAR SERIES II KARL ENTACHER y Deartment of Mathematics, University of Salzburg, Hellbrunnerstr. 34 A-52 Salzburg,

More information

A MONTE CARLO ALGORITHM FOR WEIGHTED INTEGRATION OVER R d

A MONTE CARLO ALGORITHM FOR WEIGHTED INTEGRATION OVER R d MATHEMATICS OF COMPUTATION Volume 73, Number 246, Pages 813 825 S 25-5718(3)1564-3 Article electronically published on August 19, 23 A MONTE CARLO ALGORITHM FOR WEIGHTED INTEGRATION OVER R d PIOTR GAJDA,

More information

Low-discrepancy sequences obtained from algebraic function fields over finite fields

Low-discrepancy sequences obtained from algebraic function fields over finite fields ACTA ARITHMETICA LXXII.3 (1995) Low-discrepancy sequences obtained from algebraic function fields over finite fields by Harald Niederreiter (Wien) and Chaoping Xing (Hefei) 1. Introduction. We present

More information

F 2 -Linear Random Number Generators

F 2 -Linear Random Number Generators F 2 -Linear Random Number Generators Pierre L Ecuyer and François Panneton Abstract Random number generators based on linear recurrences modulo 2 are among the fastest long-period generators currently

More information

Joshua N. Cooper Department of Mathematics, Courant Institute of Mathematics, New York University, New York, NY 10012, USA.

Joshua N. Cooper Department of Mathematics, Courant Institute of Mathematics, New York University, New York, NY 10012, USA. CONTINUED FRACTIONS WITH PARTIAL QUOTIENTS BOUNDED IN AVERAGE Joshua N. Cooper Department of Mathematics, Courant Institute of Mathematics, New York University, New York, NY 10012, USA cooper@cims.nyu.edu

More information

High Dimensional Integration: New Weapons Fighting the Curse of Dimensionality

High Dimensional Integration: New Weapons Fighting the Curse of Dimensionality High Dimensional Integration: ew Weapons Fighting the Curse of Dimensionality Peter Zinterhof Technical Report 26-4 October 26 Department of Computer Sciences Jakob-Haringer-Straße 2 52 Salzburg Austria

More information

THEORY OF GROUP REPRESENTATIONS AND APPLICATIONS

THEORY OF GROUP REPRESENTATIONS AND APPLICATIONS THEORY OF GROUP REPRESENTATIONS AND APPLICATIONS ASIM 0. BARUT Institute for Theoretical Physics, University of Colorado, Boulder, Colo., U.S.A. RYSZARD RATJZKA Institute for Nuclear Research, Warszawa,

More information

2008 Winton. Review of Statistical Terminology

2008 Winton. Review of Statistical Terminology 1 Review of Statistical Terminology 2 Formal Terminology An experiment is a process whose outcome is not known with certainty The experiment s sample space S is the set of all possible outcomes. A random

More information

Simulation Estimation of Mixed Discrete Choice Models Using Randomized and Scrambled Halton Sequences

Simulation Estimation of Mixed Discrete Choice Models Using Randomized and Scrambled Halton Sequences Simulation Estimation of Mixed Discrete Choice Models Using Randomized and Scrambled Halton Sequences Chandra R. Bhat Dept of Civil Engineering, ECJ 6.8, Univ. of Texas at Austin, Austin, TX, 78712 Phone:

More information

Monte Carlo Methods. Handbook of. University ofqueensland. Thomas Taimre. Zdravko I. Botev. Dirk P. Kroese. Universite de Montreal

Monte Carlo Methods. Handbook of. University ofqueensland. Thomas Taimre. Zdravko I. Botev. Dirk P. Kroese. Universite de Montreal Handbook of Monte Carlo Methods Dirk P. Kroese University ofqueensland Thomas Taimre University ofqueensland Zdravko I. Botev Universite de Montreal A JOHN WILEY & SONS, INC., PUBLICATION Preface Acknowledgments

More information

A NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1)-SEQUENCES. Peter Kritzer 1. Abstract

A NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1)-SEQUENCES. Peter Kritzer 1. Abstract INTEGERS: ELECTRONIC JOURNAL OF COMBINATORIAL NUMBER THEORY 5(3 (2005, #A11 A NEW UPPER BOUND ON THE STAR DISCREPANCY OF (0,1-SEQUENCES Peter Kritzer 1 Department of Mathematics, University of Salzburg,

More information

Manfred Einsiedler Thomas Ward. Ergodic Theory. with a view towards Number Theory. ^ Springer

Manfred Einsiedler Thomas Ward. Ergodic Theory. with a view towards Number Theory. ^ Springer Manfred Einsiedler Thomas Ward Ergodic Theory with a view towards Number Theory ^ Springer 1 Motivation 1 1.1 Examples of Ergodic Behavior 1 1.2 Equidistribution for Polynomials 3 1.3 Szemeredi's Theorem

More information

FAST RANDOM NUMBER GENERATORS BASED ON LINEAR RECURRENCES MODULO 2: OVERVIEW AND COMPARISON. Pierre L Ecuyer François Panneton

FAST RANDOM NUMBER GENERATORS BASED ON LINEAR RECURRENCES MODULO 2: OVERVIEW AND COMPARISON. Pierre L Ecuyer François Panneton Proceedings of the 2005 Winter Simulation Conference M. E. Kuhl, N. M. Steiger, F. B. Armstrong, and J. A. Joines, eds. FAST RANDOM NUMBER GENERATORS BASED ON LINEAR RECURRENCES MODULO 2: OVERVIEW AND

More information

Effective dimension of some weighted pre-sobolev spaces with dominating mixed partial derivatives

Effective dimension of some weighted pre-sobolev spaces with dominating mixed partial derivatives Effective dimension of some weighted pre-sobolev spaces with dominating mixed partial derivatives Art B. Owen Stanford University August 2018 Abstract This paper considers two notions of effective dimension

More information

Ordinary Differential Equations and Smooth Dynamical Systems

Ordinary Differential Equations and Smooth Dynamical Systems D.V. Anosov S.Kh. Aranson V.l. Arnold I.U. Bronshtein V.Z. Grines Yu.S. Il'yashenko Ordinary Differential Equations and Smooth Dynamical Systems With 25 Figures Springer I. Ordinary Differential Equations

More information

Low-discrepancy constructions in the triangle

Low-discrepancy constructions in the triangle Low-discrepancy constructions in the triangle Department of Statistics Stanford University Joint work with Art B. Owen. ICERM 29th October 2014 Overview Introduction 1 Introduction The Problem Motivation

More information

The Discrepancy Function and the Small Ball Inequality in Higher Dimensions

The Discrepancy Function and the Small Ball Inequality in Higher Dimensions The Discrepancy Function and the Small Ball Inequality in Higher Dimensions Dmitriy Georgia Institute of Technology (joint work with M. Lacey and A. Vagharshakyan) 2007 Fall AMS Western Section Meeting

More information

CONTENTS. Preface Preliminaries 1

CONTENTS. Preface Preliminaries 1 Preface xi Preliminaries 1 1 TOOLS FOR ANALYSIS 5 1.1 The Completeness Axiom and Some of Its Consequences 5 1.2 The Distribution of the Integers and the Rational Numbers 12 1.3 Inequalities and Identities

More information

Quasi-Monte Carlo integration over the Euclidean space and applications

Quasi-Monte Carlo integration over the Euclidean space and applications Quasi-Monte Carlo integration over the Euclidean space and applications f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with James Nichols (UNSW) Journal of Complexity 30

More information

IMPROVEMENT ON THE DISCREPANCY. 1. Introduction

IMPROVEMENT ON THE DISCREPANCY. 1. Introduction t m Mathematical Publications DOI: 10.2478/tmmp-2014-0016 Tatra Mt. Math. Publ. 59 (2014 27 38 IMPROVEMENT ON THE DISCREPANCY OF (t es-sequences Shu Tezuka ABSTRACT. Recently a notion of (t es-sequences

More information

A Belgian view on lattice rules

A Belgian view on lattice rules A Belgian view on lattice rules Ronald Cools Dept. of Computer Science, KU Leuven Linz, Austria, October 14 18, 2013 Atomium Brussels built in 1958 height 103m figure = 2e coin 5 10 6 in circulation Body

More information

Malvin H. Kalos, Paula A. Whitlock. Monte Carlo Methods. Second Revised and Enlarged Edition WILEY- BLACKWELL. WILEY-VCH Verlag GmbH & Co.

Malvin H. Kalos, Paula A. Whitlock. Monte Carlo Methods. Second Revised and Enlarged Edition WILEY- BLACKWELL. WILEY-VCH Verlag GmbH & Co. Malvin H. Kalos, Paula A. Whitlock Monte Carlo Methods Second Revised and Enlarged Edition WILEY- BLACKWELL WILEY-VCH Verlag GmbH & Co. KGaA v I Contents Preface to the Second Edition IX Preface to the

More information

B. Maddah ENMG 622 Simulation 11/11/08

B. Maddah ENMG 622 Simulation 11/11/08 B. Maddah ENMG 622 Simulation 11/11/08 Random-Number Generators (Chapter 7, Law) Overview All stochastic simulations need to generate IID uniformly distributed on (0,1), U(0,1), random numbers. 1 f X (

More information

Workshop on the Occasion of Harald Niederreiter's 70th Birthday: Applications of Algebra and Number Theory

Workshop on the Occasion of Harald Niederreiter's 70th Birthday: Applications of Algebra and Number Theory Workshop on the Occasion of Harald Niederreiter's 70th Birthday: Applications of Algebra and Number Theory June 23-27, 2014 Harald Niederreiter in Marseille Monday Tuesday Wednesday Thursday Friday 08:15

More information

Sum-discrepancy test on pseudorandom number generators

Sum-discrepancy test on pseudorandom number generators Sum-discrepancy test on pseudorandom number generators Makoto Matsumoto a,, Takuji Nishimura b a Faculty of Science, Hiroshima University, Hiroshima 739-8526, JAPAN b Faculty of Science, Yamagata University,

More information

IN QUASI-MONTE CARLO INTEGRATION

IN QUASI-MONTE CARLO INTEGRATION MATHEMATICS OF COMPUTATION VOLUME, NUMBER 0 JULY, PAGES - GOOD PARAMETERS FOR A CLASS OF NODE SETS IN QUASI-MONTE CARLO INTEGRATION TOM HANSEN, GARY L. MULLEN, AND HARALD NIEDERREITER Dedicated to the

More information

Importance Sampling to Accelerate the Convergence of Quasi-Monte Carlo

Importance Sampling to Accelerate the Convergence of Quasi-Monte Carlo Importance Sampling to Accelerate the Convergence of Quasi-Monte Carlo Wolfgang Hörmann, Josef Leydold Department of Statistics and Mathematics Wirtschaftsuniversität Wien Research Report Series Report

More information

Necessity of low effective dimension

Necessity of low effective dimension Necessity of low effective diension Art B. Owen Stanford University October 2002, Orig: July 2002 Abstract Practitioners have long noticed that quasi-monte Carlo ethods work very well on functions that

More information

The performance of a maximum-period multiple recursive generator (MRG) depends on the choices of the

The performance of a maximum-period multiple recursive generator (MRG) depends on the choices of the Published online ahead of print October 17, 2011 INFORMS Journal on Computing Articles in Advance, pp. 1 12 issn 1091-9856 eissn 1526-5528 http://dx.doi.org/10.1287/ijoc.1110.0477 2011 INFORMS Large-Order

More information

A Classical Introduction to Modern Number Theory

A Classical Introduction to Modern Number Theory Kenneth Ireland Michael Rosen A Classical Introduction to Modern Number Theory Second Edition Springer Contents Preface to the Second Edition Preface v vii CHAPTER 1 Unique Factorization 1 1 Unique Factorization

More information

Distribution of Lattice Points

Distribution of Lattice Points Computing 78, 173 193 (2006) Digital Object Identifier (DOI) 10.1007/s00607-006-0173-3 Distribution of Lattice Points F. Sezgin, Ankara Received January 11, 2006; revised June 21, 2006 Published online:

More information

Uniform and Exponential Random Floating Point Number Generation

Uniform and Exponential Random Floating Point Number Generation Uniform and Exponential Random Floating Point Number Generation Thomas Morgenstern Hochschule Harz, Friedrichstr. 57-59, D-38855 Wernigerode tmorgenstern@hs-harz.de Summary. Pseudo random number generators

More information

Nonuniform Random Variate Generation

Nonuniform Random Variate Generation Nonuniform Random Variate Generation 1 Suppose we have a generator of i.i.d. U(0, 1) r.v. s. We want to generate r.v. s from other distributions, such as normal, Weibull, Poisson, etc. Or other random

More information

Neverending Fractions

Neverending Fractions Neverending Fractions An Introduction to Continued Fractions ( c March 13, 2014: do not circulate) Jonathan Borwein, Alf van der Poorten, Jeffrey Shallit, and Wadim Zudilin Contents Preface page ix

More information

The dichotomy between structure and randomness. International Congress of Mathematicians, Aug Terence Tao (UCLA)

The dichotomy between structure and randomness. International Congress of Mathematicians, Aug Terence Tao (UCLA) The dichotomy between structure and randomness International Congress of Mathematicians, Aug 23 2006 Terence Tao (UCLA) 1 A basic problem that occurs in many areas of analysis, combinatorics, PDE, and

More information

THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES

THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES THE DISCREPANCY OF GENERALIZED VAN-DER-CORPUT-HALTON SEQUENCES MICHAEL DRMOTA Abstract. The purpose of this paper is to provide upper bounds for the discrepancy of generalized Van-der-Corput-Halton sequences

More information