Superdiffusion and Lévy Flights. A Particle Transport Monte Carlo Simulation Code
|
|
- Evan Bennett
- 5 years ago
- Views:
Transcription
1 Superdiffusion and Lévy Flights A Particle Transport Monte Carlo Simulation Code Eduardo J. Nunes-Pereira Centro de Física Escola de Ciências Universidade do Minho Page 1 of 49
2 ANOMALOUS TRANSPORT Definitions and Examples SUPERDIFFUSION AND LÉVY FLIGHTS Serial code and opportunities for parallelization DYNAMICS OF SUPERDIFFUSION Numerical algorithms Page 2 of 49
3 Page 3 of 49
4 BROWNIAN MOTION ONE OF PARADIGMS OF XX CENTURY Page 4 of 49
5 BROWNIAN MOTION DIFFUSION EQUATION Page 5 of 49
6 ANOMALOUS TRANSPORT WHY ANOMALOUS? NOT VALID DIFFUSION EQ. Page 6 of 49
7 Brownian motion Characteristic scale Mean Free Path Diffusion Equation Diffusion Coefficient Anomalous Transport No characteristic scale (fractal self-similar) No second moment (sometimes no mean free path) Diffusion Equation not valid No Diffusion Coefficient Page 7 of 49
8 Anomalous Transport No characteristic scale Jump size distribution Is a power law (all scales; self similar fractal) Page 8 of 49
9 ANOMALOUS TRANSPORT (SUPERDIFFUSION) PATHLENGTH DISTRIBUTION IS A POWER LAW LÉVY FLIGHT FLUORESCENT LAMPS Page 9 of 49
10 ANOMALOUS TRANSPORT (SUPERDIFFUSION) SUPERDIFUSION IS FASTER (SUBDIFFUSION IS SLOWER; NOT SHOWN) Page 10 of 49
11 POWER LAWS (ARE EVERYWHERE) Page 11 of 49
12 POWER LAWS (ARE EVERYWHERE) Page 12 of 49
13 ANOMALOUS TRANSPORT EXAMPLE 1 FORAGING BEHAVIOR IN SCARCE RESOURCES ENVIRONMENTS Page 13 of 49
14 ANOMALOUS TRANSPORT EXAMPLE 1 FORAGING BEHAVIOR IN SCARCE RESOURCES ENVIRONMENTS Page 14 of 49
15 ANOMALOUS TRANSPORT EXAMPLE 1 FORAGING BEHAVIOR IN SCARCE RESOURCES ENVIRONMENTS Page 15 of 49
16 ANOMALOUS TRANSPORT EXAMPLE 1 ANIMAL TRACKING Page 16 of 49
17 ANOMALOUS TRANSPORT EXAMPLE 2 POPULATION DYNAMICS AND EPIDEMICS Page 17 of 49
18 ANOMALOUS TRANSPORT EXAMPLE 2 POPULATION DYNAMICS AND EPIDEMICS Page 18 of 49
19 ANOMALOUS TRANSPORT EXAMPLE 2 POPULATION DYNAMICS AND EPIDEMICS Page 19 of 49
20 ANOMALOUS TRANSPORT EXAMPLE 3 EXTREME STATISTICS AND NATURAL HAZARDS Page 20 of 49
21 ANOMALOUS TRANSPORT EXAMPLE 4 ANOMALOUS TRASPORT FOR LIGHT SUPERDIFFUSIVE RADIATIVE TRANSPORT Page 21 of 49
22 ANOMALOUS TRANSPORT EXAMPLE 4 ANOMALOUS TRASPORT FOR LIGHT SUPERDIFFUSIVE RADIATIVE TRANSPORT Page 22 of 49
23 ANOMALOUS TRANSPORT Definitions and Examples SUPERDIFFUSION AND LÉVY FLIGHTS Serial code and opportunities for parallelization DYNAMICS OF SUPERDIFFUSION Numerical algorithms Page 23 of 49
24 THEORY VERIFIED EXPERIMENTALLY Page 24 of 49 SUPERDIFFUSIVE RADIATIVE TRANSPORT
25 LE CAMEMBERT EXPERIMENT (with R.KAISER, NICE) Transmission profile under superdiffusion SUPERDIFFUSIVE RADIATIVE TRANSPORT WITH POWER LAWS (LÉVY FLIGHTS) PARTICLE (LIGHT) TRANSPORT MONTE CARLO SIMULATION Page 25 of 49
26 SUPERDIFFUSIVE RADIATIVE TRANSPORT WITH POWER LAWS (LÉVY FLIGHTS) PARTICLE (LIGHT) TRANSPORT MONTE CARLO SIMULATION PARTICLE (LIGHT) TRANSPORT MONTE CARLO SIMULATION GENERATE TRAJECTORIES AND KEEP TRACK AND ACCUMULATE MONTE CARLO MASSIVE USE OF RANDOM NUMBERS SLOW CONVERGENCE STATISTICS (MONTE CARLO) AND POWER LAWS (EVENTS IN THE TAIL OF THE DISTRIBUTION DOMINANT) HIGH NUMBER OF TRAJECTORIES Page 26 of 49
27 SERIAL RANDOM NUMBER GENERATORS PORTABLE, PSEUDO-RANDOM NUMBERS GENERATORS STEP 1: UNIFORM DEVIATES IN [0,1[ x is distributed uniformly in [0,1[ STEP 2 (EXAMPLE): EXPONENTIAL DEVIATES Page 27 of 49
28 SERIAL RANDOM NUMBER GENERATORS PORTABLE, PSEUDO-RANDOM NUMBERS GENERATORS STEP 3: GENERAL TRANSFORMATION METHOD (FOR ANY STATISTICAL DISTRIBUTION) Solve with Indefinite integral of desired distribution SOLVE ANALYTICALLY OR NUMERICALLY STEP 2 (EXAMPLE): EXPONENTIAL DEVIATES ANALYTICAL Page 28 of 49
29 SERIAL RANDOM NUMBER GENERATORS PORTABLE, PSEUDO-RANDOM NUMBERS GENERATORS STEP 3: GENERAL TRANSFORMATION METHOD (FOR ANY STATISTICAL DISTRIBUTION) SOLVE ANALYTICALLY OR NUMERICALLY Page 29 of 49
30 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT STEP 1: UNIFORM DEVIATES IN [0,1[ x is distributed uniformly in [0,1[ BASIC ALGORITHM LINEAR CONGRENTIAL GENERATOR INTEGER ARITHMETIC / IS THE DESIRED REAL RANDOM NUMBER SEED + SET OF MAGIC NUMBERS, FOR MAXIMUM REPETITION PERIOD Page 30 of 49
31 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT Page 31 of 49 ran2 NUMERICAL RECIPES 2 GENERATORS PLUS ADDITIONAL SHUFFLING PROCEDURE (BREAK SEQUENTIAL CORRELATION)
32 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT ANALYTICAL GENERAL TRANSFORMATION METHOD IN MONTE CARLO TRAJECTORIES SPHERICAL COORDINATES THE DIRECTION OF THE NEXT STEP IS ANALYTICAL HOWEVER, THE STEP LENGTH IS NOT ANALYTICAL Page 32 of 49
33 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT NUMERICAL GENERAL TRANSFORMATION METHOD IN MONTE CARLO TRAJECTORIES THE SINGLE STEP F(y) IS COMPUTED NUMERICALLY ONCE IN A SET OF DISCRETE DATA POINTS Page 33 of 49
34 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT CASE I DOPPLER RANDOM DEVIATES THEN A NATURAL CUBIC SPLINE INTERPOLATION IMPLEMENTS THE TRANSFORMATION METHOD Page 34 of 49
35 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT CASE II VOIGT RANDOM DEVIATES THEN A NATURAL CUBIC SPLINE INTERPOLATION IMPLEMENTS THE TRANSFORMATION METHOD Page 35 of 49
36 RANDOM NUMBER GENERATORS IN SERIAL CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT THEN A NATURAL CUBIC SPLINE INTERPOLATION IMPLEMENTS THE TRANSFORMATION METHOD THE CUBIC SPLINE INTERPOLATION IS A 2 STEP ALGORITHM: STEP 1: SOLVE LINEAR TRIDIAGONAL SYSTEM (ONCE; NO TIME COST) STEP 2: INTERPOLATE BY TABLE LOOKUP AND BISSECTION (REPEAT ) THE MONTE CARLO CODE IS EXPECTED TO SPENT MOST OF THE TIME: GENERATING UNIFORM RANDOM NUMBERS IN [0,1[ IMPLEMENTING NUMERICALLY THE TRANSFORMATION METHOD (BISSECTION) Page 36 of 49
37 OPPORTUNITIES FOR PARALLELIZATION (PROVIDE FEEDBACK) MONTE CARLO PARTICLE TRANSPORT CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT: SERIAL CODE IN FORTRAN 90 BOTTLENECK RANDOM NUMBERS (UNIFORM DEVIATES & TRANSFORMATION METHOD) EACH TIME A RANDOM NUMBER IS NEEDED, IS GENERATED FIGURE OF MERIT TO STOP SIMULATION TOTAL NUMBER TRAJECTORIES MONTE CARLO PARTICLE TRANSPORT CODE FOR SUPERDIFFUSIVE RADIATIVE TRANSPORT ALGOTITHM: GENERATE TRAJECTORIES ACUMULATE POSITIONS REPEAT UNTIL TOTAL NUMBER OF TRAJECTORIES Page 37 of 49
38 OPPORTUNITIES FOR PARALLELIZATION (PROVIDE FEEDBACK) OPPORTUNITY #1 MODEL MPI MASTER/SLAVE MODEL WITH DYNAMIC LOAD BALANCE IN HETEROGENEOUS SYSTEMS TIME IS LINEAR IN # TRAJECTORIES => LAUNCH REPLICAS OF SMALL MONTE CARLO SIMULATIONS IN EACH COMPUTATIONAL NODE AND, BASED ON FEEDBACK, PARTITION TOTAL NUMBER OF TRAJECTORIES TO OBTAIN OVERALL LOAD BALANCE DISADVANTAGE: EACH NODE IS INDEPENDENT AND DOES NOT SHARES RANDOM NUMBER SEQUENCES TESTED IN PROTOTYPE CODE Page 38 of 49
39 OPPORTUNITIES FOR PARALLELIZATION (PROVIDE FEEDBACK) OPPORTUNITY #2 PARALLEL RANDOM NUMBER GENERATORS DATA SHARING THE UNIFORM RANDOM NUMBER SEQUENCE IS GLOBALLY ACCESSIBLE FROM ALL COMPUTING NODES POSSIBLY COMBINE WITH OPPORTUNITY #1; MASTER/SLAVE MODEL WITH DYNAMIC LOAD BALANCE IN HETEROGENEOUS SYSTEMS NO EXPERIENCE WITH PARALLEL RANDOM NUMBER GENERATORS Page 39 of 49
40 OPPORTUNITIES FOR PARALLELIZATION (PROVIDE FEEDBACK) OPPORTUNITY #3 HYBRID CPU/GPU MODELS LEVEL #1: LINEAR CONGRUENTIAL GENERATOR USES INTEGER ARITHMETIC (IS THIS STILL RELEVANT TO EXPLORE, BASED IN LATEST HARDWARE?) LEVEL #2: RANDOM NUMBER GENERATOR BY NUMERICAL INTERPOLATION IS A LOOKUP ALGORITHM IN AN ORDERED TABLE; COULD BE MORE EFFICIENTLY RECODED IN GPU? Page 40 of 49
41 ANOMALOUS TRANSPORT Definitions and Examples SUPERDIFFUSION AND LÉVY FLIGHTS Serial code and opportunities for parallelization DYNAMICS OF SUPERDIFFUSION Numerical algorithms Page 41 of 49
42 GOAL ALGORITHM/NUMERICAL CODE DYNAMICS OF ANOMALOUS TRANSPORT Page 42 of 49
43 APPLICABLE FOR OTHER POWER LAW CASES Page 43 of 49 SUPERDIFFUSION (POWER LAWS)
44 MASTER EQUATION INTEGRO-DIFFERENTIAL EQUATION DESCRIBES DYNAMICS 1. LINEAR REGIME AND EXPANSION IN # OF JUMPS (MONTE CARLO CODE) 2. LINEAR REGIME AND EXPANSION IN EIGENVALUES/EIGENVECTORS 3. LINEAR/NON LINEAR REGIME AND INTEGRATION MASTER EQUATION IN TIME AND SPACE Page 44 of 49
45 1. LINEAR EXPANSION IN # OF JUMPS + MONTE CARLO SIMULATION NUMERICAL ALGORITHMS SPATIAL DISTRIBUTIONS IN # JUMPS (MONTE CARLO) TEMPORAL DISTRIBUTIONS ARE ANALYTICAL (DECOUPLED FROM SPACE) MONTE CARLO SIMULATION OF TRAJECTORIES (ONLY SPACE) BOTTLENECK RANDOM NUMBER GENERATION/MANIPULATION Page 45 of 49
46 2. LINEAR REGIME AND EIGENVALUE/EIGENVECTOR DECOMPOSITION + MONTE CARLO SIMULATION (FOR RANDOM CONFIGURATIONS) TRANSITION MATRIX FOR EACH CONFIGURATION (EIGENVALUES/EIGENVECTORS) + CONFIGURATIONAL ENSEMBLES (MONTE CARLO REPLICAS OF DIFFERENT CONFIGURATIONS) CORE COMPUTACIONAL LINEAR ALGEBRA (EIGENVALUES/EIGENVECTORS) + RANDOM NUMBERS CPU or CPU/GPU HYBRID LIBRARIES Page 46 of 49
47 3. LINEAR/NON LINEAR REGIME USING FRACTIONAL DERIVATIVES MODELS FOR SPACE AND TIME DISCRETIZATION THE MASTER EQUATION CAN BE WRITTE AS A GENERALIZATION OF THE DIFFUSION EQUATION, BUT NOW USING (FRACTIONAL) DERIVATIVE ORDER FRACTIONAL CALCULUS / FRACTIONAL DIFFUSION / FRACTIONAL KINETICS Page 47 of 49
48 3. LINEAR/NON LINEAR REGIME USING FRACTIONAL DERIVATIVES MODELS FOR SPACE AND TIME TIME SPACE (SUPERDIFFUSION) CORE COMPUTATIONAL LINEAR ALGEBRA HOT TOPIC AND matlab CODE Page 48 of 49
49 3. LINEAR/NON LINEAR REGIME USING FRACTIONAL DERIVATIVES MODELS FOR SPACE AND TIME DISCRETIZATION Page 49 of 49
Monte Carlo Method for Solving Inverse Problems of Radiation Transfer
INVERSE AND ILL-POSED PROBLEMS SERIES Monte Carlo Method for Solving Inverse Problems of Radiation Transfer V.S.Antyufeev. ///VSP/// UTRECHT BOSTON KÖLN TOKYO 2000 Contents Chapter 1. Monte Carlo modifications
More informationStochastic Simulation: Algorithms and Analysis
Soren Asmussen Peter W. Glynn Stochastic Simulation: Algorithms and Analysis et Springer Contents Preface Notation v xii I What This Book Is About 1 1 An Illustrative Example: The Single-Server Queue 1
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 informationJoe Wingbermuehle, (A paper written under the guidance of Prof. Raj Jain)
1 of 11 5/4/2011 4:49 PM Joe Wingbermuehle, wingbej@wustl.edu (A paper written under the guidance of Prof. Raj Jain) Download The Auto-Pipe system allows one to evaluate various resource mappings and topologies
More informationFast Multipole Method on the GPU
Fast Multipole Method on the GPU with application to the Adaptive Vortex Method University of Bristol, Bristol, United Kingdom. 1 Introduction Particle methods Highly parallel Computational intensive Numerical
More informationContents Computing with Formulas
Contents 1 Computing with Formulas... 1 1.1 The First Programming Encounter: a Formula... 1 1.1.1 Using a Program as a Calculator... 2 1.1.2 About Programs and Programming... 2 1.1.3 Tools for Writing
More informationSequential Monte Carlo Adaptation in Low-Anisotropy Participating Media. Vincent Pegoraro Ingo Wald Steven G. Parker
Sequential Monte Carlo Adaptation in Low-Anisotropy Participating Media Vincent Pegoraro Ingo Wald Steven G. Parker Outline Introduction Related Work Monte Carlo Integration Radiative Energy Transfer SMC
More informationMonte Carlo Method on Parallel Computing. Jongsoon Kim
Monte Carlo Method on Parallel Computing Jongsoon Kim Introduction Monte Carlo methods Utilize random numbers to perform a statistical simulation of a physical problem Extremely time-consuming Inherently
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 informationIntroduction to Computational Mathematics
Introduction to Computational Mathematics Introduction Computational Mathematics: Concerned with the design, analysis, and implementation of algorithms for the numerical solution of problems that have
More informationDEPARTMENT - Mathematics. Coding: N Number. A Algebra. G&M Geometry and Measure. S Statistics. P - Probability. R&P Ratio and Proportion
DEPARTMENT - Mathematics Coding: N Number A Algebra G&M Geometry and Measure S Statistics P - Probability R&P Ratio and Proportion YEAR 7 YEAR 8 N1 Integers A 1 Simplifying G&M1 2D Shapes N2 Decimals S1
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 informationRandom Numbers Random Walk
Random Numbers Random Walk Computational Physics Random Numbers Random Walk Outline Random Systems Random Numbers Monte Carlo Integration Example Random Walk Exercise 7 Introduction Random Systems Deterministic
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationWhat is Monte Carlo Modeling*?
What is Monte Carlo Modeling*? Monte Carlo Modeling is a statisitcal method used here to simulate radiative transfer by simulating photon (or more exactly light rays/beams) interaction with a medium. MC
More informationMonte Carlo Methods; Combinatorial Search
Monte Carlo Methods; Combinatorial Search Parallel and Distributed Computing Department of Computer Science and Engineering (DEI) Instituto Superior Técnico November 22, 2012 CPD (DEI / IST) Parallel and
More informationCITS2401 Computer Analysis & Visualisation
FACULTY OF ENGINEERING, COMPUTING AND MATHEMATICS CITS2401 Computer Analysis & Visualisation SCHOOL OF COMPUTER SCIENCE AND SOFTWARE ENGINEERING Topic 13 Revision Notes CAV review Topics Covered Sample
More informationCombinatorial Search; Monte Carlo Methods
Combinatorial Search; Monte Carlo Methods Parallel and Distributed Computing Department of Computer Science and Engineering (DEI) Instituto Superior Técnico May 02, 2016 CPD (DEI / IST) Parallel and Distributed
More information1. Mathematical Modelling
1. describe a given problem with some mathematical formalism in order to get a formal and precise description see fundamental properties due to the abstraction allow a systematic treatment and, thus, solution
More informationAdvanced Numerical Methods for Numerical Weather Prediction
Advanced Numerical Methods for Numerical Weather Prediction Francis X. Giraldo Naval Research Laboratory Monterey, CA 93943-5502 phone: (831) 656-4882 fax: (831) 656-4769 e-mail: giraldo@nrlmry.navy.mil
More informationMonte Carlo Simulations
Monte Carlo Simulations DESCRIPTION AND APPLICATION Outline Introduction Description of Method Cost Estimating Example Other Considerations Introduction Most interesting things are probabilistic (opinion)
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 informationTutorial: Modeling Liquid Reactions in CIJR Using the Eulerian PDF transport (DQMOM-IEM) Model
Tutorial: Modeling Liquid Reactions in CIJR Using the Eulerian PDF transport (DQMOM-IEM) Model Introduction The purpose of this tutorial is to demonstrate setup and solution procedure of liquid chemical
More informationGauss for Econometrics: Simulation
Gauss for Econometrics: Simulation R.G. Pierse 1. Introduction Simulation is a very useful tool in econometric modelling. It allows the economist to examine the properties of models and estimators when
More informationA methodology for the rigorous verification of plasma simulation codes
A methodology for the rigorous verification of plasma simulation codes Fabio Riva P. Ricci, C. Beadle, F.D. Halpern, S. Jolliet, J. Loizu, J. Morales, A. Mosetto, P. Paruta, C. Wersal École Polytechnique
More informationHow to perform HPL on CPU&GPU clusters. Dr.sc. Draško Tomić
How to perform HPL on CPU&GPU clusters Dr.sc. Draško Tomić email: drasko.tomic@hp.com Forecasting is not so easy, HPL benchmarking could be even more difficult Agenda TOP500 GPU trends Some basics about
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 informationRandom Walks. Dieter W. Heermann. Monte Carlo Methods. Dieter W. Heermann (Monte Carlo Methods) Random Walks / 22
Random Walks Dieter W. Heermann Monte Carlo Methods 2009 Dieter W. Heermann (Monte Carlo Methods) Random Walks 2009 1 / 22 Outline 1 Introduction 2 Random Walk 3 Self-Avoiding Random Walks Simple Sampling
More informationThe Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer
The Spherical Harmonics Discrete Ordinate Method for Atmospheric Radiative Transfer K. Franklin Evans Program in Atmospheric and Oceanic Sciences University of Colorado, Boulder Computational Methods in
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Dick K.P. Yue Center for Ocean Engineering Massachusetts Institute of Technology Room 5-321, 77 Massachusetts Ave, Cambridge, MA 02139 phone:
More informationIntel Performance Libraries
Intel Performance Libraries Powerful Mathematical Library Intel Math Kernel Library (Intel MKL) Energy Science & Research Engineering Design Financial Analytics Signal Processing Digital Content Creation
More informationModule 4. Computer-Aided Design (CAD) systems
Module 4. Computer-Aided Design (CAD) systems Nowadays the design of complex systems is unconceivable without computers. The fast computers, the sophisticated developing environments and the well elaborated
More informationAccelerated tokamak transport simulations
Accelerated tokamak transport simulations via Neural-Network based regression of TGLF turbulent energy, particle and momentum fluxes by Teobaldo Luda 1 O. Meneghini 2, S. Smith 2, G. Staebler 2 J. Candy
More informationHigh Performance Computing. Introduction to Parallel Computing
High Performance Computing Introduction to Parallel Computing Acknowledgements Content of the following presentation is borrowed from The Lawrence Livermore National Laboratory https://hpc.llnl.gov/training/tutorials
More informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationNumerical Methods in Physics Lecture 2 Interpolation
Numerical Methods in Physics Pat Scott Department of Physics, Imperial College November 8, 2016 Slides available from http://astro.ic.ac.uk/pscott/ course-webpage-numerical-methods-201617 Outline The problem
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 informationSTEPHEN WOLFRAM MATHEMATICADO. Fourth Edition WOLFRAM MEDIA CAMBRIDGE UNIVERSITY PRESS
STEPHEN WOLFRAM MATHEMATICADO OO Fourth Edition WOLFRAM MEDIA CAMBRIDGE UNIVERSITY PRESS Table of Contents XXI a section new for Version 3 a section new for Version 4 a section substantially modified for
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 informationInterpolation and Basis Fns
CS148: Introduction to Computer Graphics and Imaging Interpolation and Basis Fns Topics Today Interpolation Linear and bilinear interpolation Barycentric interpolation Basis functions Square, triangle,,
More informationcomputational Fluid Dynamics - Prof. V. Esfahanian
Three boards categories: Experimental Theoretical Computational Crucial to know all three: Each has their advantages and disadvantages. Require validation and verification. School of Mechanical Engineering
More informationarxiv: v2 [math.na] 11 Jun 2014
Parallel stochastic methods for PDE based grid generation Alexander Bihlo and Ronald D. Haynes Department of Mathematics and Statistics, Memorial University of Newfoundland, St. John s (NL), AC 5S7, Canada
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture XV (04.02.08) Contents: Function Minimization (see E. Lohrmann & V. Blobel) Optimization Problem Set of n independent variables Sometimes in addition some constraints
More informationTime Series Analysis by State Space Methods
Time Series Analysis by State Space Methods Second Edition J. Durbin London School of Economics and Political Science and University College London S. J. Koopman Vrije Universiteit Amsterdam OXFORD UNIVERSITY
More informationParallel Computing. Slides credit: M. Quinn book (chapter 3 slides), A Grama book (chapter 3 slides)
Parallel Computing 2012 Slides credit: M. Quinn book (chapter 3 slides), A Grama book (chapter 3 slides) Parallel Algorithm Design Outline Computational Model Design Methodology Partitioning Communication
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 informationToday. Participating media. Participating media. Rendering Algorithms: Participating Media and. Subsurface scattering
Today Rendering Algorithms: Participating Media and Subsurface Scattering Introduction Rendering participating media Rendering subsurface scattering Spring 2009 Matthias Zwicker Participating media Participating
More informationDesigning Parallel Programs. This review was developed from Introduction to Parallel Computing
Designing Parallel Programs This review was developed from Introduction to Parallel Computing Author: Blaise Barney, Lawrence Livermore National Laboratory references: https://computing.llnl.gov/tutorials/parallel_comp/#whatis
More informationDiscrete Event Simulation & VHDL. Prof. K. J. Hintz Dept. of Electrical and Computer Engineering George Mason University
Discrete Event Simulation & VHDL Prof. K. J. Hintz Dept. of Electrical and Computer Engineering George Mason University Discrete Event Simulation Material from VHDL Programming with Advanced Topics by
More informationVideo 11.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar
Video 11.1 Vijay Kumar 1 Smooth three dimensional trajectories START INT. POSITION INT. POSITION GOAL Applications Trajectory generation in robotics Planning trajectories for quad rotors 2 Motion Planning
More informationDirect Numerical Simulation of Turbulent Boundary Layers at High Reynolds Numbers.
Direct Numerical Simulation of Turbulent Boundary Layers at High Reynolds Numbers. G. Borrell, J.A. Sillero and J. Jiménez, Corresponding author: guillem@torroja.dmt.upm.es School of Aeronautics, Universidad
More informationIntroduction. Chapter 1
Chapter 1 Introduction In the physical sciences we often encounter problems of evaluating various properties of a given function f (x). Typicaloperationsaredifferentiation,integrationandfinding the roots
More informationComputational Study of Protein Diffusion in a Membrane. By Kate Schneider Advisor Dr. Ken Ritchie
Computational Study of Protein Diffusion in a Membrane By Kate Schneider Advisor Dr. Ken Ritchie 1 Cell RBC Membrane Red blood cell membrane. Embedded Proteins Membrane protects the cell. Lipids have two
More informationMassively Parallel Computation for Three-Dimensional Monte Carlo Semiconductor Device Simulation
L SIMULATION OF SEMICONDUCTOR DEVICES AND PROCESSES Vol. 4 Edited by W. Fichtner, D. Aemmer - Zurich (Switzerland) September 12-14,1991 - Hartung-Gorre Massively Parallel Computation for Three-Dimensional
More informationMotion Synthesis and Editing. Yisheng Chen
Motion Synthesis and Editing Yisheng Chen Overview Data driven motion synthesis automatically generate motion from a motion capture database, offline or interactive User inputs Large, high-dimensional
More informationProf. Fanny Ficuciello Robotics for Bioengineering Trajectory planning
Trajectory planning to generate the reference inputs to the motion control system which ensures that the manipulator executes the planned trajectories path and trajectory joint space trajectories operational
More informationInterpolation and Basis Fns
CS148: Introduction to Computer Graphics and Imaging Interpolation and Basis Fns Topics Today Interpolation Linear and bilinear interpolation Barycentric interpolation Basis functions Square, triangle,,
More informationImproved Convergence Rates in Implicit Monte Carlo Simulations Through Stratified Sampling
Improved Convergence Rates in Implicit Monte Carlo Simulations Through Stratified Sampling ANS Winter Conference 2013 Alex Long and Ryan McClarren Texas A&M University Alex Long (Texas A&M) ANS Winter
More informationParallel Algorithm Design. Parallel Algorithm Design p. 1
Parallel Algorithm Design Parallel Algorithm Design p. 1 Overview Chapter 3 from Michael J. Quinn, Parallel Programming in C with MPI and OpenMP Another resource: http://www.mcs.anl.gov/ itf/dbpp/text/node14.html
More informationComputer Simulations for a Fractional Calculus Derived Internet Traffic Model
Paper 21 Civil-Comp Press, 26. Proceedings of the Fifth International Conference on Engineering Computational Technology, B.H.V. Topping, G. Montero and R. Montenegro, (Editors), Civil-Comp Press, Stirlingshire,
More informationCENTRALITIES. Carlo PICCARDI. DEIB - Department of Electronics, Information and Bioengineering Politecnico di Milano, Italy
CENTRALITIES Carlo PICCARDI DEIB - Department of Electronics, Information and Bioengineering Politecnico di Milano, Italy email carlo.piccardi@polimi.it http://home.deib.polimi.it/piccardi Carlo Piccardi
More informationNovember , Universität Leipzig, CompPhys12
Parallel streams of pseudorandom numbers for Monte Carlo simulations: Using most reliable algorithms and applying parallelism of modern CPUs and GPUs L.Yu. Barash Landau Institute for Theoretical Physics
More informationNIH Public Access Author Manuscript Opt Lett. Author manuscript; available in PMC 2012 January 27.
NIH Public Access Author Manuscript Published in final edited form as: Opt Lett. 2011 June 1; 36(11): 2095 2097. Fast perturbation Monte Carlo method for photon migration in heterogeneous turbid media
More informationIntroduction to Parallel Computing
Introduction to Parallel Computing Bootcamp for SahasraT 7th September 2018 Aditya Krishna Swamy adityaks@iisc.ac.in SERC, IISc Acknowledgments Akhila, SERC S. Ethier, PPPL P. Messina, ECP LLNL HPC tutorials
More informationThe Course Structure for the MCA Programme
The Course Structure for the MCA Programme SEMESTER - I MCA 1001 Problem Solving and Program Design with C 3 (3-0-0) MCA 1003 Numerical & Statistical Methods 4 (3-1-0) MCA 1007 Discrete Mathematics 3 (3-0-0)
More informationIntroduction to parallel Computing
Introduction to parallel Computing VI-SEEM Training Paschalis Paschalis Korosoglou Korosoglou (pkoro@.gr) (pkoro@.gr) Outline Serial vs Parallel programming Hardware trends Why HPC matters HPC Concepts
More informationOpen Compute Stack (OpenCS) Overview. D.D. Nikolić Updated: 20 August 2018 DAE Tools Project,
Open Compute Stack (OpenCS) Overview D.D. Nikolić Updated: 20 August 2018 DAE Tools Project, http://www.daetools.com/opencs What is OpenCS? A framework for: Platform-independent model specification 1.
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationIntroduction to High Performance Computing
Introduction to High Performance Computing Gregory G. Howes Department of Physics and Astronomy University of Iowa Iowa High Performance Computing Summer School University of Iowa Iowa City, Iowa 25-26
More informationParallel Systems. Project topics
Parallel Systems Project topics 2016-2017 1. Scheduling Scheduling is a common problem which however is NP-complete, so that we are never sure about the optimality of the solution. Parallelisation is a
More informationArrow Tech Course Catalog
Arrow Tech Course Catalog Advanced Training for the Weapons Designer COURSE PAGE DURATION (DAYS) Basic PRODAS 3 3 Basic PRODAS (Extended) 5 5 Advanced User PRODAS 7 3 BALANS (Projectile Balloting Analysis)
More informationAssessing the Quality of the Natural Cubic Spline Approximation
Assessing the Quality of the Natural Cubic Spline Approximation AHMET SEZER ANADOLU UNIVERSITY Department of Statisticss Yunus Emre Kampusu Eskisehir TURKEY ahsst12@yahoo.com Abstract: In large samples,
More informationDESIGNING OF STREAM CIPHER ARCHITECTURE USING THE CELLULAR AUTOMATA
DESIGNING OF STREAM CIPHER ARCHITECTURE USING THE CELLULAR AUTOMATA 1 Brundha K A MTech Email: 1 brundha1905@gmail.com Abstract Pseudo-random number generators (PRNGs) are a key component of stream ciphers
More informationECS 289 / MAE 298, Lecture 15 Mar 2, Diffusion, Cascades and Influence, Part II
ECS 289 / MAE 298, Lecture 15 Mar 2, 2011 Diffusion, Cascades and Influence, Part II Diffusion and cascades in networks (Nodes in one of two states) Viruses (human and computer) contact processes epidemic
More information"Fast, High-Fidelity, Multi-Spacecraft Trajectory Simulation for Space Catalogue Applications"
"Fast, High-Fidelity, Multi-Spacecraft Trajectory Simulation for Space Catalogue Applications" Ryan P. Russell Assistant Professor Nitin Arora Ph.D. Candidate School of Aerospace Engineering Georgia Institute
More informationGeometric and Solid Modeling. Problems
Geometric and Solid Modeling Problems Define a Solid Define Representation Schemes Devise Data Structures Construct Solids Page 1 Mathematical Models Points Curves Surfaces Solids A shape is a set of Points
More information6.375 Group 2: Adaptive PIV Final Report
6.375 Group 2: Adaptive PIV Final Report Robin Deits rdeits@csail.mit.edu May 15, 2013 Contents 1 Overview 2 2 Background: PIV 2 3 Adaptive PIV 2 4 Target Application 3 5 Implementation 3 5.1 Module: PIV.............................
More informationSupersedes: S-S-01. (a) acceptance sampling of discrete units presented for inspection in lots;
Specifications Category: STATISTICAL METHOD Specification: S-S-01 (rev.1) Page: 1 of 13 Specifications for Random Sampling and Randomization 1.0 Scope 1.1 These specifications define algorithms for random
More informationPLB-HeC: A Profile-based Load-Balancing Algorithm for Heterogeneous CPU-GPU Clusters
PLB-HeC: A Profile-based Load-Balancing Algorithm for Heterogeneous CPU-GPU Clusters IEEE CLUSTER 2015 Chicago, IL, USA Luis Sant Ana 1, Daniel Cordeiro 2, Raphael Camargo 1 1 Federal University of ABC,
More informationVerification and Validation of X-Sim: A Trace-Based Simulator
http://www.cse.wustl.edu/~jain/cse567-06/ftp/xsim/index.html 1 of 11 Verification and Validation of X-Sim: A Trace-Based Simulator Saurabh Gayen, sg3@wustl.edu Abstract X-Sim is a trace-based simulator
More informationParallelization of DQMC Simulations for Strongly Correlated Electron Systems
Parallelization of DQMC Simulations for Strongly Correlated Electron Systems Che-Rung Lee Dept. of Computer Science National Tsing-Hua University Taiwan joint work with I-Hsin Chung (IBM Research), Zhaojun
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 informationModeling with Uncertainty Interval Computations Using Fuzzy Sets
Modeling with Uncertainty Interval Computations Using Fuzzy Sets J. Honda, R. Tankelevich Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, CO, U.S.A. Abstract A new method
More information1\C 1 I)J1mptMIll 'betl~flbejlli
1\C 1 )J1mptMll 'betl~flbejll l~~1hlbd ada]t6gh -or~«ejf. '~illrlf J~..6 '~~~ll!~ 4iJ~ "Mf:i',nl.Nqr2l' ~':l:mj}.i~:tv t.~l '\h Dr. N.Homsup, Abstract n this paper, two high-order FDTD schemes are developed
More informationLectures & Excercises
TLTE.3120 Computer Simulation in Communication and Systems (5 ECTS) http://www.uva.fi/~timan/tlte3120 Lecture 1 9.9.2015 Timo Mantere Professor, Embedded systems University of Vaasa http://www.uva.fi/~timan
More informationHYPERDRIVE IMPLEMENTATION AND ANALYSIS OF A PARALLEL, CONJUGATE GRADIENT LINEAR SOLVER PROF. BRYANT PROF. KAYVON 15618: PARALLEL COMPUTER ARCHITECTURE
HYPERDRIVE IMPLEMENTATION AND ANALYSIS OF A PARALLEL, CONJUGATE GRADIENT LINEAR SOLVER AVISHA DHISLE PRERIT RODNEY ADHISLE PRODNEY 15618: PARALLEL COMPUTER ARCHITECTURE PROF. BRYANT PROF. KAYVON LET S
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationGAMES Webinar: Rendering Tutorial 2. Monte Carlo Methods. Shuang Zhao
GAMES Webinar: Rendering Tutorial 2 Monte Carlo Methods Shuang Zhao Assistant Professor Computer Science Department University of California, Irvine GAMES Webinar Shuang Zhao 1 Outline 1. Monte Carlo integration
More informationResearch on time optimal trajectory planning of 7-DOF manipulator based on genetic algorithm
Acta Technica 61, No. 4A/2016, 189 200 c 2017 Institute of Thermomechanics CAS, v.v.i. Research on time optimal trajectory planning of 7-DOF manipulator based on genetic algorithm Jianrong Bu 1, Junyan
More informationVector Visualization
Vector Visualization Vector Visulization Divergence and Vorticity Vector Glyphs Vector Color Coding Displacement Plots Stream Objects Texture-Based Vector Visualization Simplified Representation of Vector
More informationWRITING, STRUCTURING AND DEVELOPING MONTE CARLO RADIATIVE TRANSFER CODES ANTONIA BEVAN, UCL ST ANDREWS MONTE CARLO SUMMER SCHOOL 2017
WRITING, STRUCTURING AND DEVELOPING MONTE CARLO RADIATIVE TRANSFER CODES ANTONIA BEVAN, UCL ST ANDREWS MONTE CARLO SUMMER SCHOOL 2017 Me in Oct 2012 just after I started my PhD Herschel, Planck and Spitzer
More informationMaster-Worker pattern
COSC 6397 Big Data Analytics Master Worker Programming Pattern Edgar Gabriel Fall 2018 Master-Worker pattern General idea: distribute the work among a number of processes Two logically different entities:
More informationOptimized Scientific Computing:
Optimized Scientific Computing: Coding Efficiently for Real Computing Architectures Noah Kurinsky SASS Talk, November 11 2015 Introduction Components of a CPU Architecture Design Choices Why Is This Relevant
More informationScalable Video Coding
Introduction to Multimedia Computing Scalable Video Coding 1 Topics Video On Demand Requirements Video Transcoding Scalable Video Coding Spatial Scalability Temporal Scalability Signal to Noise Scalability
More informationValidation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons
Validation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons Stores Separation Introduction Flight Test Expensive, high-risk, sometimes catastrophic
More informationMaster-Worker pattern
COSC 6397 Big Data Analytics Master Worker Programming Pattern Edgar Gabriel Spring 2017 Master-Worker pattern General idea: distribute the work among a number of processes Two logically different entities:
More informationHigh performance 2D Discrete Fourier Transform on Heterogeneous Platforms. Shrenik Lad, IIIT Hyderabad Advisor : Dr. Kishore Kothapalli
High performance 2D Discrete Fourier Transform on Heterogeneous Platforms Shrenik Lad, IIIT Hyderabad Advisor : Dr. Kishore Kothapalli Motivation Fourier Transform widely used in Physics, Astronomy, Engineering
More informationCombos: Advection-Diffusion and operator splitting
96 Chapter 7 Combos: Advection-Diffusion and operator splitting Up to this point we have treated advection and diffusion as two separate processes and considered appropriate numerical techniques for solving
More informationEECS 556 Image Processing W 09. Interpolation. Interpolation techniques B splines
EECS 556 Image Processing W 09 Interpolation Interpolation techniques B splines What is image processing? Image processing is the application of 2D signal processing methods to images Image representation
More informationCS 470 Spring Parallel Algorithm Development. (Foster's Methodology) Mike Lam, Professor
CS 470 Spring 2018 Mike Lam, Professor Parallel Algorithm Development (Foster's Methodology) Graphics and content taken from IPP section 2.7 and the following: http://www.mcs.anl.gov/~itf/dbpp/text/book.html
More information