What We ll Do... Random
|
|
- Jonas Pitts
- 6 years ago
- Views:
Transcription
1 What We ll Do... Random- number generation Random Number Generation Generating random variates Nonstationary Poisson processes Variance reduction Sequential sampling Designing and executing simulation experiments Random-Number Generators (RNGs) Algorithm to generate independent, identically distributed draws from the continuous UNIF (0, ) distribution f(x) These are called Random Number Generators in simulation 0 Basis for generating observations from all other distributions and random processes Transform random numbers in a way that depends on the desired distribution or process It s essential to have a good RNG There are a lot of bad RNGs this is very tricky Methods, coding are both tricky x Properties of RNG Numbers produced should appear to be distributed uniformly on [0,] and should not exhibit any correlation with each other. Generator should be fast. Should be able to reproduce a given stream of random numbers: For debugging, and To use identical random numbers in simulating different systems in order to obtain a more precise comparison Should be a provision in the generator to easily reproduce separate streams of random numbers. The Nature of RNGs Recursive formula (algorithm) Start with a seed (or seed vector) Do something weird to the seed to get the next one Repeat generate the same sequence Will eventually repeat (cycle) want long cycle length Not really random as in unpredictable Want to design RNGs Long cycle length Good statistical properties (uniform, independent) -- tests Fast Streams (subsegments) many and long (for variance reduction later) Linear Congruential Generators (LCGs) The most common of several different methods Generate a sequence of integers Z, Z 2, Z 3, via the recursion Z i = (a Z i + c) (mod m) If c > 0 : mixed LCGs and If c = 0 : multiplicative LCGs a, c, and m are carefully chosen constants Specify a seed Z 0 to start off mod m means take the remainder of dividing by m as the next Z i All the Z i s are between 0 and m Return the ith random number as U i = Z i / m
2 Example of LCG Parameters m = 63, a = 22, c = 4, Z 0 = 9: Z i = (22 Z i + 4) (mod 63), seed with Z 0 = 9 Another Example of LCG m = 6, a = 5, c = 3, Z 0 = 7 LCG Z i = (5 Z i- + 3)(mod 6) with Z 0 = 7 i 22 Z i +4 Z i U i : :: : : :: : Cycling will repeat forever Cycle length is at most m (could be << m depending on parameters) Pick m BIG But that might not be enough for good statistical properties Z = (22*9 + 4) (mod 63) = (48 + 4) mod 63 = 422 mod 63 = 44 U = 44 / 63 = i Z i U i Issues with LCGs Cycle length: < m Typically, m = 2. billion (= 2 3 ) or more Which used to be a lot not anymore Other parameters chosen so that cycle length = m or m Statistical properties Uniformity, independence There are many tests of RNGs Empirical tests Theoretical tests lattice structure (next slide ) Speed, storage both are usually fine Must be carefully, cleverly coded BIG integers Reproducibility streams (long internal subsequences) with fixed seeds Issues with LCGs (cont d.) Regularity of LCGs (and other kinds of RNGs): For the earlier LCG Plot of U i vs. i Plot of U i+ vs. U i Random Numbers Fall Mainly in the Planes George Marsaglia Design RNGs: dense lattice in high dimensions Other kinds of RNGs longer memory in recursion, combination of several RNGs The Current Arena RNG Uses some of the same ideas as LCG Modulo division, recursive on earlier values But is not an LCG Combines two separate component generators Recursion involves more than just the preceding value Combined multiple recursive generator (CMRG) A n = ( A n A n-3 ) mod B n = (52762 B n B n-3 ) mod Z n = (A n B n ) mod Combine the two Z n / if Z n > 0 U n = / if Z n = 0 Seed = a six-vector of first three A n s, B n s The next random number Two simultaneous recursions The Current Arena RNG Properties Extremely good statistical properties Good uniformity Cycle length = To cycle, all six seeds must match up Only slightly slower than old LCG And RNG is usually a minor part of overall computing time
3 Mixed Generators Other RNGs Choose m = 2 b, where b is the bits in the computer system. 32 bit system has b = 3 (left most bit is sign bit) That will give m = 2 3 > 2. billion. Multiplicative Generators Addition of c is not needed Instead of m = 2 b, let m = a large prime number < 2 3 Combined LCGs Combine 2 or more multiplicative LCGs so as to have good statistical properties and longer period. Tests for Random Numbers Need to test uniformity and independence Frequency test Kolmogorov Smirnov (K-S) test or Chi-Square test to compare distribution of the set of numbers generated to a uniform distribution Runs Test Uses Chi-Square test to compare the runs above and below the mean by comparing actual values with expected values Autocorrelation Test Tests the correlation between numbers and compares the sample correlation with expected correlation of zero Gap Test Counts the number of digits that appear between repetitions of a particular digit and then uses K-S test to compare with the expected size of gaps Generating Random Variates Generating from Discrete Distributions How to transform a uniform distribution between 0 and into draws from the input probability distributions you want for your model? Have: Desired input distribution for model (fitted or specified in some way), and RNG (UNIF (0, )) Want: Transform UNIF (0, ) random numbers into draws from the desired input distribution Method: Mathematical transformations of random numbers to deform them to the desired distribution Specific transform depends on desired distribution Do discrete, continuous distributions separately Example: probability mass function Divide [0, ] Into subintervals of length 0., 0.5, 0.4 Generate U ~ UNIF (0, ) See which subinterval it s in Return X = corresponding value Discrete Generation: Another View Example Plot cumulative distribution function; generate U and plot on vertical axis; read across and down Inverting the CDF Equivalent to earlier method The discrete RV X has the following probability function: k P (X = k) Use the uniform random number U = 0.79 to generate an observation for X 0.00 to 0.30 = 0.30 to 0.65 = to 0.85 = 3 U = 0.79 X = 3
4 Continuous Distributions Exponential The exponential distribution probablity density function (PDF) is given by λe λx x 0 f(x) = 0 x < 0 Cumulative Distribution Function (CDF) is given by For any i F(x) = - e λx 0 x 0 x < 0 E(X i ) = / λ = mean interarrival time. Goal is to develop a procedure for generating values X, X 2, etc λ = mean number of occurrences per time unit (arrival rate) EXPO (5) distribution Density (PDF) Distribution (CDF) Example General algorithm Inverse Transform Technique:. Generate a random number U ~ UNIF(0, ) 2. Set U = F(X) and solve for X = F (U) Solving analytically for X may or may not be simple (or possible) Sometimes use numerical approximation to solve Example (cont d.) Uniform Distribution Set U = F(X) = e X/5 e X/5 = U X/5 = ln ( U) X = 5 ln ( U) Picture (inverting the CDF, as in discrete case): PDF: CDF: f(x) = F(x) = a x b b a 0 otherwise 0 x a b a x < a a x b x > b F(x) = (X a) / (b a) = U X = a + (b a) U X is uniformly distributed on the interval [a, b]. PDF: f(x) = CDF: For 0 X For X 2 Triangular Distribution F(x) = U = X 2 / 2 U = (2 X) 2 2 x 0 x 2 x < x 2 0 otherwise 0 x 0 x 2 / 2 0 < x ((2 x) 2 / 2) < x 2 x > 2 Nonstationary Poisson Processes Many systems have externally originating events affecting them e.g., arrivals of customers If process is stationary over time, usually specify a fixed interevent- time distribution But process could vary markedly in its rate Fast-food lunch rush Freeway rush hours Ignoring nonstationarity can lead to serious model and output errors
5 Nonstationary Poisson Processes Definition Usual model: nonstationary Poisson process: Have a rate function l(t) Number of events in [t, t 2 ] ~ Poisson with mean λ(t) Nonstationary Poisson Processes Estimating the Rate Function Estimation of the rate function Probably the most practical method is piecewise constant Decide on a time interval within which rate is fixed Estimate from data the (constant) rate during each interval Be careful to get the units right Issues: How to estimate rate function? Given an estimate, how to generate during simulation? t Other (more complicated) methods exist in the literature Variance Reduction Random input random output (RIRO) In other words, output has variance Higher output variance means less precise results Would like to eliminate or reduce output variance One (bad) way to eliminate: replace all input random variables by constants (like their mean) Will get rid of random output, but will also invalidate model Thus, best hope is to reduce output variance Easy (brute- force) variance reduction: just simulate some more Terminating: additional replications Steady-state: additional replications or a longer run Variance Reduction (cont d.) But sometimes can reduce variance without more runs Key: unlike physical experiments, can control randomness in computer-simulation experiments via manipulating the RNG Re-use the same random numbers either as they were, in some opposite sense, or for a similar but simpler model Induce certain kinds of correlations to exploit to take advantage to reduce the variance of the output. Several different variance-reduction techniques Classified into categories common random numbers, antithetic variates, control variates, indirect estimation, Usually requires thorough understanding of model, code Will look only at common random numbers in detail Common Random Numbers (CRN) Applies when objective is to compare two (or more) alternative configurations or models Interest is in difference(s) of performance measure(s) across alternatives Example: A. Base case (as is) B. 3.5% increase in business (interarrival-time mean falls from 3 to 2.56 minutes) Same run conditions The Usual Comparison Run case A, make the change to get to case B and run it, then Compare Means via Output Analyzer: Difference is not statistically significant Were the runs of A and B statistically independent? Did we use the same random numbers running A and B?
6 CRN Intuition Get sharper comparison if you subject all alternatives to the same conditions Then observed differences are due to model differences rather than random differences in the conditions For both A and B runs, cause: The same parts arrive at the same times Be assigned same attributes (job type) Have the same process times at each step Then observed differences will be attributable to system differences, not random bounce There isn t any random bounce Synchronization of Random Numbers in CRN One approach is to dedicate a stream of random numbers to each place in the model where variates are needed. extra parameter in r.v. calls Fairly simple but might not ensure complete synchronization Another approach is to assign to each entity, immediately upon its arrival, attribute values for all possible processing times, branching decisions, etc. when the entity needs one of these values, you just read it out of the appropriate attribute of the entity instead of generating it on the spot. Takes lot of computer memory Sequential Sampling Always try to quantify imprecision in results If imprecision is small enough, you re done If not, need to do something to increase precision Just saw one way: variance-reduction techniques Obvious way to increase precision: keep simulating one more step at a time, quit when you achieve desired precision Terminating models: step = another replication Cannot extend length of replications that s part of the model Steady-state models: step = another replication if using truncated replications, or step = some extension of the run if using batch means
CPSC 531: System Modeling and Simulation. Carey Williamson Department of Computer Science University of Calgary Fall 2017
CPSC 531: System Modeling and Simulation Carey Williamson Department of Computer Science University of Calgary Fall 2017 Outline Random number generation Properties of random numbers Linear Congruential
More informationTesting Random- Number Generators
Testing Random- Number Generators Raj Jain Washington University Saint Louis, MO 63131 Jain@cse.wustl.edu These slides are available on-line at: http://www.cse.wustl.edu/~jain/cse574-06/ 27-1 Overview
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 13 Random Numbers and Stochastic Simulation Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright
More informationRandom Number Generators
1/17 Random Number Generators Professor Karl Sigman Columbia University Department of IEOR New York City USA 2/17 Introduction Your computer generates" numbers U 1, U 2, U 3,... that are considered independent
More informationChapter 6: Simulation Using Spread-Sheets (Excel)
Chapter 6: Simulation Using Spread-Sheets (Excel) Refer to Reading Assignments 1 Simulation Using Spread-Sheets (Excel) OBJECTIVES To be able to Generate random numbers within a spreadsheet environment.
More informationComputational Methods. Randomness and Monte Carlo Methods
Computational Methods Randomness and Monte Carlo Methods Manfred Huber 2010 1 Randomness and Monte Carlo Methods Introducing randomness in an algorithm can lead to improved efficiencies Random sampling
More informationProbability Models.S4 Simulating Random Variables
Operations Research Models and Methods Paul A. Jensen and Jonathan F. Bard Probability Models.S4 Simulating Random Variables In the fashion of the last several sections, we will often create probability
More informationRandom-Number Generation
Random-Number Generation Overview Desired properties of a good generator Linear-congruential generators Tausworthe generators Survey of random number generators Seed selection Myths about random number
More informationYou ve already read basics of simulation now I will be taking up method of simulation, that is Random Number Generation
Unit 5 SIMULATION THEORY Lesson 39 Learning objective: To learn random number generation. Methods of simulation. Monte Carlo method of simulation You ve already read basics of simulation now I will be
More informationInput Analysis. Input Analysis: Specifying Model Parameters, Distributions. Deterministic vs. Random Inputs
Input Analysis Input Analysis: Specifying Model Parameters, Distributions Structural modeling: what we ve done so far Logical aspects entities, resources, paths, etc. Quantitative modeling Numerical, distributional
More informationAnalysis of Simulation Results
Analysis of Simulation Results Raj Jain Washington University Saint Louis, MO 63130 Jain@cse.wustl.edu Audio/Video recordings of this lecture are available at: http://www.cse.wustl.edu/~jain/cse574-08/
More information>>> SOLUTIONS <<< a) What is a system? Give a short formal definition. (Hint: what you draw a box around is not a formal definition)
Mid-Term Exam for Simulation (CIS 4930) Summer 2009 >>> SOLUTIONS
More informationNumerical Precision. Or, why my numbers aren t numbering right. 1 of 15
Numerical Precision Or, why my numbers aren t numbering right 1 of 15 What s the deal? Maybe you ve seen this #include int main() { float val = 3.6f; printf( %.20f \n, val); Print a float with
More informationCOMPUTER SCIENCE LARGE PRACTICAL.
COMPUTER SCIENCE LARGE PRACTICAL Page 45 of 100 SURVEY RESULTS Approx. 1/5 of class responded; statistically significant? The majority of you have substantial experience in Java, and all have at least
More informationMonte Carlo Integration COS 323
Monte Carlo Integration COS 323 Last time Interpolatory Quadrature Review formulation; error analysis Newton-Cotes Quadrature Midpoint, Trapezoid, Simpson s Rule Error analysis for trapezoid, midpoint
More informationMonte Carlo Integration COS 323
Monte Carlo Integration COS 323 Integration in d Dimensions? One option: nested 1-D integration f(x,y) g(y) y f ( x, y) dx dy ( ) = g y dy x Evaluate the latter numerically, but each sample of g(y) is
More informationMonte Carlo Integration and Random Numbers
Monte Carlo Integration and Random Numbers Higher dimensional integration u Simpson rule with M evaluations in u one dimension the error is order M -4! u d dimensions the error is order M -4/d u In general
More informationDECISION SCIENCES INSTITUTE. Exponentially Derived Antithetic Random Numbers. (Full paper submission)
DECISION SCIENCES INSTITUTE (Full paper submission) Dennis Ridley, Ph.D. SBI, Florida A&M University and Scientific Computing, Florida State University dridley@fsu.edu Pierre Ngnepieba, Ph.D. Department
More informationWhy is Statistics important in Bioinformatics?
Why is Statistics important in Bioinformatics? Random processes are inherent in evolution and in sampling (data collection). Errors are often unavoidable in the data collection process. Statistics helps
More informationSimulation with Arena
Simulation with Arena Sixth Edition W. David Kelton Professor Department of Operations, Business Analytics, and Information Systems University of Cincinnati Randall P. Sadowski Retired Nancy B. Zupick
More informationLecture 2: Introduction to Numerical Simulation
Lecture 2: Introduction to Numerical Simulation Ahmed Kebaier kebaier@math.univ-paris13.fr HEC, Paris Outline of The Talk 1 Simulation of Random variables Outline 1 Simulation of Random variables Random
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: April 11, 2011. Lecture 1: Introduction and Basic Terms Welcome to the course, time table, assessment, etc..
More informationChapter 4: (0,1) Random Number Generation
Chapter 4: (0,1) Random Number Generation Refer to Text Book: Simulation Modeling and ARENA, Manuel Rossetti, Ch. 2 Operations Research: Applications and Algorithms By Wayne L. Winston,Ch. 21 Operations
More informationPackage simed. November 27, 2017
Version 1.0.3 Title Simulation Education Author Barry Lawson, Larry Leemis Package simed November 27, 2017 Maintainer Barry Lawson Imports graphics, grdevices, methods, stats, utils
More informationSimulation Modeling and Analysis
Simulation Modeling and Analysis FOURTH EDITION Averill M. Law President Averill M. Law & Associates, Inc. Tucson, Arizona, USA www. averill-law. com Boston Burr Ridge, IL Dubuque, IA New York San Francisco
More informationMicroscopic Traffic Simulation
Microscopic Traffic Simulation Lecture Notes in Transportation Systems Engineering Prof. Tom V. Mathew Contents Overview 2 Traffic Simulation Models 2 2. Need for simulation.................................
More informationUNIT 9A Randomness in Computation: Random Number Generators
UNIT 9A Randomness in Computation: Random Number Generators 1 Last Unit Computer organization: what s under the hood 3 This Unit Random number generation Using pseudorandom numbers 4 Overview The concept
More informationChapter 16. Microscopic Traffic Simulation Overview Traffic Simulation Models
Chapter 6 Microscopic Traffic Simulation 6. Overview The complexity of traffic stream behaviour and the difficulties in performing experiments with real world traffic make computer simulation an important
More informationCS 112: Computer System Modeling Fundamentals. Prof. Jenn Wortman Vaughan April 21, 2011 Lecture 8
CS 112: Computer System Modeling Fundamentals Prof. Jenn Wortman Vaughan April 21, 2011 Lecture 8 Quiz #2 Reminders & Announcements Homework 2 is due in class on Tuesday be sure to check the posted homework
More informationFathom Dynamic Data TM Version 2 Specifications
Data Sources Fathom Dynamic Data TM Version 2 Specifications Use data from one of the many sample documents that come with Fathom. Enter your own data by typing into a case table. Paste data from other
More informationLecture: Simulation. of Manufacturing Systems. Sivakumar AI. Simulation. SMA6304 M2 ---Factory Planning and scheduling. Simulation - A Predictive Tool
SMA6304 M2 ---Factory Planning and scheduling Lecture Discrete Event of Manufacturing Systems Simulation Sivakumar AI Lecture: 12 copyright 2002 Sivakumar 1 Simulation Simulation - A Predictive Tool Next
More informationRandom Number Generation and Monte Carlo Methods
James E. Gentle Random Number Generation and Monte Carlo Methods With 30 Illustrations Springer Contents Preface vii 1 Simulating Random Numbers from a Uniform Distribution 1 1.1 Linear Congruential Generators
More informationLinear Congruential Number Generators. A useful, if not important, ability of modern computers is random number
Jagannath Pisharath Newer Math Fall 2003 Linear Congruential Number Generators A useful, if not important, ability of modern computers is random number generation. Without this ability, if you wanted to,
More informationIntroduction to the course
Introduction to the course Lecturer: Dmitri A. Moltchanov E-mail: moltchan@cs.tut.fi http://www.cs.tut.fi/ moltchan/modsim/ http://www.cs.tut.fi/kurssit/tlt-2706/ 1. What is the teletraffic theory? Multidisciplinary
More informationGetting started with simulating data in R: some helpful functions and how to use them Ariel Muldoon August 28, 2018
Getting started with simulating data in R: some helpful functions and how to use them Ariel Muldoon August 28, 2018 Contents Overview 2 Generating random numbers 2 rnorm() to generate random numbers from
More informationTutorial: Random Number Generation
Tutorial: Random Number Generation John Lau johnlau@umail.ucsb.edu Henry Yu henryyu@umail.ucsb.edu June 2018 1 Abstract This tutorial will cover the basics of Random Number Generation. This includes properties
More informationn! = 1 * 2 * 3 * 4 * * (n-1) * n
The Beauty and Joy of Computing 1 Lab Exercise 9: Problem self-similarity and recursion Objectives By completing this lab exercise, you should learn to Recognize simple self-similar problems which are
More information(Refer Slide Time: 00:26)
Programming, Data Structures and Algorithms Prof. Shankar Balachandran Department of Computer Science and Engineering Indian Institute Technology, Madras Module 07 Lecture 07 Contents Repetitive statements
More informationCDA6530: Performance Models of Computers and Networks. Chapter 8: Statistical Simulation --- Discrete-Time Simulation
CDA6530: Performance Models of Computers and Networks Chapter 8: Statistical Simulation --- Discrete-Time Simulation Simulation Studies Models with analytical formulas Calculate the numerical solutions
More informationDiscrete-Event Simulation:
Discrete-Event Simulation: A First Course Section 32: Multi-Stream Lehmer RNGs Section 32: Multi-Stream Lehmer RNGs Discrete-Event Simulation c 2006 Pearson Ed, Inc 0-13-142917-5 Section 32: Multi-Stream
More informationSection 3.2: Multi Stream Lehmer RNGs
Section 32: Multi Stream Lehmer RNGs Discrete-Event Simulation: A First Course c 2006 Pearson Ed, Inc 0-13-142917-5 Discrete-Event Simulation: A First Course Section 32: Multi Stream Lehmer RNGs 1/ 19
More informationNotes for Lecture 10
COS 533: Advanced Cryptography Lecture 10 (October 16, 2017) Lecturer: Mark Zhandry Princeton University Scribe: Dylan Altschuler Notes for Lecture 10 1 Motivation for Elliptic Curves Diffie-Hellman For
More informationDiscrete-Event Simulation:
Discrete-Event Simulation: A First Course Section 2.2: Lehmer Random Number Generators: Implementation Section 2.2: Discrete-Event Simulation c 2006 Pearson Ed., Inc. 0-13-142917-5 Section 2.2: Lehmer
More information% Close all figure windows % Start figure window
CS1112 Fall 2016 Project 3 Part A Due Monday 10/3 at 11pm You must work either on your own or with one partner. If you work with a partner, you must first register as a group in CMS and then submit your
More informationCPSC 467b: Cryptography and Computer Security
CPSC 467b: Cryptography and Computer Security Michael J. Fischer Lecture 7 January 30, 2012 CPSC 467b, Lecture 7 1/44 Public-key cryptography RSA Factoring Assumption Computing with Big Numbers Fast Exponentiation
More informationA Quick Introduction to R
Math 4501 Fall 2012 A Quick Introduction to R The point of these few pages is to give you a quick introduction to the possible uses of the free software R in statistical analysis. I will only expect you
More informationSIMULATION OF A SINGLE-SERVER QUEUEING SYSTEM
SIMULATION OF A SINGLE-SERVER QUEUEING SYSTEM Will show how to simulate a specific version of the single-server queuing system Though simple, it contains many features found in all simulation models 1-
More informationSimulation. Outline. Common Mistakes in Simulation (3 of 4) Common Mistakes in Simulation (2 of 4) Performance Modeling Lecture #8
Introduction (1 of 3) The best advice to those about to embark on a very large simulation is often the same as Punch s famous advice to those about to marry: Don t! Bratley, Fox and Schrage (1986) Simulation
More informationA Rant on Queues. Van Jacobson. July 26, MIT Lincoln Labs Lexington, MA
A Rant on Queues Van Jacobson July 26, 2006 MIT Lincoln Labs Lexington, MA Unlike the phone system, the Internet supports communication over paths with diverse, time varying, bandwidth. This means we often
More information! Addition! Multiplication! Bigger Example - RSA cryptography
! Addition! Multiplication! Bigger Example - RSA cryptography Modular Arithmetic Modular Exponentiation Primality Testing (Fermat s little theorem) Probabilistic algorithm Euclid s Algorithm for gcd (greatest
More informationChapter 1. Introduction
Chapter 1 Introduction A Monte Carlo method is a compuational method that uses random numbers to compute (estimate) some quantity of interest. Very often the quantity we want to compute is the mean of
More informationThe Plan: Basic statistics: Random and pseudorandom numbers and their generation: Chapter 16.
Scientific Computing with Case Studies SIAM Press, 29 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit IV Monte Carlo Computations Dianne P. O Leary c 28 What is a Monte-Carlo method?
More informationis bad when combined with a LCG with a small multiplier. In section 2 this observation is examined. Section 3 gives portable implementations for LCGs
Version: 25.11.92 (To appear in ACM TOMS) A Portable Uniform Random Number Generator Well Suited for the Rejection Method W. Hormann and G. Deringer University of Economics and Business Administration
More informationA simple OMNeT++ queuing experiment using different random number generators
A simple OMNeT++ queuing experiment using different random number generators Bernhard Hechenleitner and Karl Entacher December 5, 2002 Abstract We apply a simple queuing-experiment using parallel streams
More informationChapter 3. Bootstrap. 3.1 Introduction. 3.2 The general idea
Chapter 3 Bootstrap 3.1 Introduction The estimation of parameters in probability distributions is a basic problem in statistics that one tends to encounter already during the very first course on the subject.
More informationData and File Structures Chapter 11. Hashing
Data and File Structures Chapter 11 Hashing 1 Motivation Sequential Searching can be done in O(N) access time, meaning that the number of seeks grows in proportion to the size of the file. B-Trees improve
More informationKINEMATICS FOR ANIMATION. Rémi Ronfard, Animation, M2R MOSIG
KINEMATICS FOR ANIMATION Rémi Ronfard, Animation, M2R MOSIG Direct and inverse kinematics in animation Fast numerical methods for inverse kinematics by Bill Baxter My Adventures with Inverse Kinematics
More informationCongruence Arithmetic
Module 4 Congruence Arithmetic Popper 4 Introduction to what is like Modulus choices Partitions by modulus Mod 5 Mod 7 Mod 30 Modular Arithmetic Addition Subtraction Multiplication INTEGERS! Mod 12 Cayley
More informationWHOLE NUMBER AND DECIMAL OPERATIONS
WHOLE NUMBER AND DECIMAL OPERATIONS Whole Number Place Value : 5,854,902 = Ten thousands thousands millions Hundred thousands Ten thousands Adding & Subtracting Decimals : Line up the decimals vertically.
More informationReproducibility in Stochastic Simulation
Reproducibility in Stochastic Simulation Prof. Michael Mascagni Department of Computer Science Department of Mathematics Department of Scientific Computing Graduate Program in Molecular Biophysics Florida
More informationTABLES AND HASHING. Chapter 13
Data Structures Dr Ahmed Rafat Abas Computer Science Dept, Faculty of Computer and Information, Zagazig University arabas@zu.edu.eg http://www.arsaliem.faculty.zu.edu.eg/ TABLES AND HASHING Chapter 13
More information(Refer Slide Time: 00:02:00)
Computer Graphics Prof. Sukhendu Das Dept. of Computer Science and Engineering Indian Institute of Technology, Madras Lecture - 18 Polyfill - Scan Conversion of a Polygon Today we will discuss the concepts
More informationIntegration. Volume Estimation
Monte Carlo Integration Lab Objective: Many important integrals cannot be evaluated symbolically because the integrand has no antiderivative. Traditional numerical integration techniques like Newton-Cotes
More informationPackage rstream. R topics documented: June 21, Version Date Title Streams of Random Numbers
Version 1.3.5 Date 2017-06-21 Title Streams of Random Numbers Package rstream June 21, 2017 Author Maintainer Depends R (>= 2.0.0) Imports methods, stats Unified object oriented interface for multiple
More informationIntroduction to Programming in C Department of Computer Science and Engineering\ Lecture No. #02 Introduction: GCD
Introduction to Programming in C Department of Computer Science and Engineering\ Lecture No. #02 Introduction: GCD In this session, we will write another algorithm to solve a mathematical problem. If you
More informationWhy Use Graphs? Test Grade. Time Sleeping (Hrs) Time Sleeping (Hrs) Test Grade
Analyzing Graphs Why Use Graphs? It has once been said that a picture is worth a thousand words. This is very true in science. In science we deal with numbers, some times a great many numbers. These numbers,
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit IV Monte Carlo
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit IV Monte Carlo Computations Dianne P. O Leary c 2008 1 What is a Monte-Carlo
More informationMean Value Analysis and Related Techniques
Mean Value Analysis and Related Techniques 34-1 Overview 1. Analysis of Open Queueing Networks 2. Mean-Value Analysis 3. Approximate MVA 4. Balanced Job Bounds 34-2 Analysis of Open Queueing Networks Used
More informationModeling and Simulation (An Introduction)
Modeling and Simulation (An Introduction) 1 The Nature of Simulation Conceptions Application areas Impediments 2 Conceptions Simulation course is about techniques for using computers to imitate or simulate
More informationExcel 2010 with XLSTAT
Excel 2010 with XLSTAT J E N N I F E R LE W I S PR I E S T L E Y, PH.D. Introduction to Excel 2010 with XLSTAT The layout for Excel 2010 is slightly different from the layout for Excel 2007. However, with
More informationLecture 2: Analyzing Algorithms: The 2-d Maxima Problem
Lecture 2: Analyzing Algorithms: The 2-d Maxima Problem (Thursday, Jan 29, 1998) Read: Chapter 1 in CLR. Analyzing Algorithms: In order to design good algorithms, we must first agree the criteria for measuring
More information36 Modular Arithmetic
36 Modular Arithmetic Tom Lewis Fall Term 2010 Tom Lewis () 36 Modular Arithmetic Fall Term 2010 1 / 10 Outline 1 The set Z n 2 Addition and multiplication 3 Modular additive inverse 4 Modular multiplicative
More information(Refer Slide Time: 02:59)
Numerical Methods and Programming P. B. Sunil Kumar Department of Physics Indian Institute of Technology, Madras Lecture - 7 Error propagation and stability Last class we discussed about the representation
More informationImplementing Coroutines with call/cc. Producer/Consumer using Coroutines
Implementing Coroutines with call/cc Producer/Consumer using Coroutines Coroutines are a very handy generalization of subroutines. A coroutine may suspend its execution and later resume from the point
More informationCPSC 340: Machine Learning and Data Mining
CPSC 340: Machine Learning and Data Mining Fundamentals of learning (continued) and the k-nearest neighbours classifier Original version of these slides by Mark Schmidt, with modifications by Mike Gelbart.
More informationChapter 1. Math review. 1.1 Some sets
Chapter 1 Math review This book assumes that you understood precalculus when you took it. So you used to know how to do things like factoring polynomials, solving high school geometry problems, using trigonometric
More informationPractice Problems (Con t) The ALU performs operation x and puts the result in the RR The ALU operand Register B is loaded with the contents of Rx
Microprogram Control Practice Problems (Con t) The following microinstructions are supported by each CW in the CS: RR ALU opx RA Rx RB Rx RB IR(adr) Rx RR Rx MDR MDR RR MDR Rx MAR IR(adr) MAR Rx PC IR(adr)
More informationApplied Cryptography and Network Security
Applied Cryptography and Network Security William Garrison bill@cs.pitt.edu 6311 Sennott Square Lecture #8: RSA Didn t we learn about RSA last time? During the last lecture, we saw what RSA does and learned
More informationCS302 Topics: * More on Simulation and Priority Queues * Random Number Generation
CS302 Topics: * More on Simulation and Priority Queues * Random Number Generation Thursday, Oct. 5, 2006 Greetings from Salvador, Brazil!! (i.e., from my recent trip to present a plenary lecture on multi-robot
More informationCMSC 341 Lecture 16/17 Hashing, Parts 1 & 2
CMSC 341 Lecture 16/17 Hashing, Parts 1 & 2 Prof. John Park Based on slides from previous iterations of this course Today s Topics Overview Uses and motivations of hash tables Major concerns with hash
More informationSimulation Input Data Modeling
Introduction to Modeling and Simulation Simulation Input Data Modeling OSMAN BALCI Professor Department of Computer Science Virginia Polytechnic Institute and State University (Virginia Tech) Blacksburg,
More information(Refer Slide Time: 00:51)
Programming, Data Structures and Algorithms Prof. Shankar Balachandran Department of Computer Science and Engineering Indian Institute Technology, Madras Module 10 E Lecture 24 Content Example: factorial
More informationA-SSE.1.1, A-SSE.1.2-
Putnam County Schools Curriculum Map Algebra 1 2016-2017 Module: 4 Quadratic and Exponential Functions Instructional Window: January 9-February 17 Assessment Window: February 20 March 3 MAFS Standards
More information6. Control Statements II
Visibility Declaration in a block is not visible outside of the block. 6. Control Statements II Visibility, Local Variables, While Statement, Do Statement, Jump Statements main block int main () int i
More informationMidterm Exam. CS381-Cryptography. October 30, 2014
Midterm Exam CS381-Cryptography October 30, 2014 Useful Items denotes exclusive-or, applied either to individual bits or to sequences of bits. The same operation in Python is denoted ˆ. 2 10 10 3 = 1000,
More informationOverview of the Simulation Process. CS1538: Introduction to Simulations
Overview of the Simulation Process CS1538: Introduction to Simulations Simulation Fundamentals A computer simulation is a computer program that models the behavior of a physical system over time. Program
More informationIntroduction to Information Retrieval (Manning, Raghavan, Schutze)
Introduction to Information Retrieval (Manning, Raghavan, Schutze) Chapter 3 Dictionaries and Tolerant retrieval Chapter 4 Index construction Chapter 5 Index compression Content Dictionary data structures
More informationTo complete the computer assignments, you ll use the EViews software installed on the lab PCs in WMC 2502 and WMC 2506.
An Introduction to EViews The purpose of the computer assignments in BUEC 333 is to give you some experience using econometric software to analyse real-world data. Along the way, you ll become acquainted
More informationMontana City School GRADE 5
Montana City School GRADE 5 Montana Standard 1: Students engage in the mathematical processes of problem solving and reasoning, estimation, communication, connections and applications, and using appropriate
More informationRANDOM NUMBERS GENERATION
Chapter 4 RANDOM NUMBERS GENERATION M. Ragheb 10/2/2015 4.1. INTRODUCTION The use and generation of random numbers uniformly distributed over the unit interval: [0, 1] is a unique feature of the Monte
More informationRANDOM NUMBERS GENERATION
Chapter 4 RANDOM NUMBERS GENERATION M. Ragheb 9//2013 4.1. INTRODUCTION The use and generation of random numbers uniformly distributed over the unit interval: [0, 1] is a unique feature of the Monte Carlo
More informationFast Automated Estimation of Variance in Discrete Quantitative Stochastic Simulation
Fast Automated Estimation of Variance in Discrete Quantitative Stochastic Simulation November 2010 Nelson Shaw njd50@uclive.ac.nz Department of Computer Science and Software Engineering University of Canterbury,
More informationPlanar Graphs and Surfaces. Graphs 2 1/58
Planar Graphs and Surfaces Graphs 2 1/58 Last time we discussed the Four Color Theorem, which says that any map can be colored with at most 4 colors and not have two regions that share a border having
More informationTesting parallel random number generators
Parallel Computing 29 (2003) 69 94 www.elsevier.com/locate/parco Testing parallel random number generators Ashok Srinivasan a, Michael Mascagni b, *, David Ceperley c a Department of Computer Science,
More informationMAT 243 Test 2 SOLUTIONS, FORM A
MAT 243 Test 2 SOLUTIONS, FORM A 1. [15 points] Calculate the following quantities: a. 17 mod 4 Solution: 17 17 4 = 17 4 4 = 1. 4 b. 17 div 4 17 Solution: = 4. 4 c. (( 1) mod 12) mod (27 div 5) Solution:
More informationHi everyone. Starting this week I'm going to make a couple tweaks to how section is run. The first thing is that I'm going to go over all the slides
Hi everyone. Starting this week I'm going to make a couple tweaks to how section is run. The first thing is that I'm going to go over all the slides for both problems first, and let you guys code them
More informationReading Assignment. Lazy Evaluation
Reading Assignment Lazy Evaluation MULTILISP: a language for concurrent symbolic computation, by Robert H. Halstead (linked from class web page Lazy evaluation is sometimes called call by need. We do an
More informationVARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung
POLYTECHNIC UNIVERSITY Department of Computer and Information Science VARIANCE REDUCTION TECHNIQUES IN MONTE CARLO SIMULATIONS K. Ming Leung Abstract: Techniques for reducing the variance in Monte Carlo
More informationEECS 121: Coding for Digital Communication & Beyond Fall Lecture 22 December 3. Introduction
EECS 121: Coding for Digital Communication & Beyond Fall 2013 Lecture 22 December 3 Lecturer: K. V. Rashmi Scribe: Ajay Shanker Tripathi 22.1 Context Introduction Distributed storage is a deeply relevant
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More information