Hochdimensionale Integration
|
|
- George Powell
- 6 years ago
- Views:
Transcription
1 Oliver Ernst Institut für Numerische Mathematik und Optimierung Hochdimensionale Integration 4-tägige Vorlesung im Wintersemester 2/ im Rahmen des Moduls Ausgewählte Kapitel der Numerik
2 Contents. Introduction. An Example.2 A Selection of Strategies 2. Monte Carlo Integration 2. Convergence and Accuracy 2.2 Sampling Methods 2.3 Variance Reduction Methods 3. Sparse Grids 4. Quasi-Monte Carlo Integration 5. Extensions 5. ANOVA Decomposition 5.2 Concentration of Measure Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/
3 Inhalt. Introduction. An Example.2 A Selection of Strategies 2. Monte Carlo Integration 2. Convergence and Accuracy 2.2 Sampling Methods 2.3 Variance Reduction Methods 3. Sparse Grids 4. Quasi-Monte Carlo Integration 5. Extensions 5. ANOVA Decomposition 5.2 Concentration of Measure Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 2
4 Object of Study In these lectures we will study numerical methods for approximating definite integrals [,] d f(x ) dx. [, ] d d-dimensional unit cube, d large f : [, ] d R continuous. Today: Why is this useful? What methods are out there? Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 3
5 Inhalt. Introduction. An Example.2 A Selection of Strategies 2. Monte Carlo Integration 2. Convergence and Accuracy 2.2 Sampling Methods 2.3 Variance Reduction Methods 3. Sparse Grids 4. Quasi-Monte Carlo Integration 5. Extensions 5. ANOVA Decomposition 5.2 Concentration of Measure Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 4
6 Example : Pricing Derivatives [Traub & Passkov (995)]: Goldman Sachs CMO Fannie Mae REMIC Trust CMO : Collateralized Mortgage Obligation REMIC : Real Estate Mortgage Investment Conduit Pool of mortgages with 3-year maturity and monthly payment, i.e., 36 cash flows Cash flows distributed across tranches A, B, C, D, E, G, H, J, R, Z. Notation C : monthly payment on underlying pool of mortgages, a 36 k+ : remaining annuity after month k, i k : interest rate in month k, random w k : percentage prepaying in month k, random a k = + v + + v k, v = + i, k =,..., 36, i : current monthly interest rate. Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 5
7 Example : Pricing Derivatives Interest rate model: i k = K e iξ k i k = K k i e ξ+ +ξ k, {ξ k } 36 k= i.i.d. random variables, ξ k N(, σ 2 ), σ = 4 4. Prepayment model: w k = w k (ξ,..., ξ k ) = K + K 2 arctan(k 3 i k + K 4 ), K j constants Cash flow in month k: M k = M k (ξ,..., ξ k ) = C ( w (ξ ) )( w 2(ξ, ξ 2) ) ( w k (ξ,..., ξ k ) ) + w k (ξ,..., ξ k )a 36 k+ Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 6
8 Example : Pricing Derivatives Cash flow distributed to tranches according to CMO rules. Portion of M k directed to tranche T in month k described by function G k,t = G k,t (ξ,..., ξ k ) (continuous, composition of smooth functions with min, complicated). Present value of tranche T in month k: multiply G k,t with discount factor u k = u k (ξ,..., ξ k ) := v v (ξ )... v k (ξ,..., ξ k ), v j (ξ,..., ξ j ) = + i j (ξ,..., ξ j ) = + K j j,..., 359. eξ+ +ξj Sum up values over all months: present value PV of tranche T given by 36 P V T (ξ,..., ξ 36 ) = G k,t (ξ,..., ξ k ) u k (ξ,..., ξ k ). k= Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 7
9 Example : Pricing Derivatives Expected value E(P V T ) = E ( P V T (ξ,..., ξ 36 ) ) after change of variables y i := y i (x i ), x i = 2πσ yi e t2 2σ given by E(P V T ) = P V T (y (x ),..., y 36 (x 36 )) dx dx 36 [,] 36 dt, Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 8
10 Inhalt. Introduction. An Example.2 A Selection of Strategies 2. Monte Carlo Integration 2. Convergence and Accuracy 2.2 Sampling Methods 2.3 Variance Reduction Methods 3. Sparse Grids 4. Quasi-Monte Carlo Integration 5. Extensions 5. ANOVA Decomposition 5.2 Concentration of Measure Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 9
11 Numerical Quadrature Numerical quadrature (integration) of univariate functions is essentially a solved problem. Denoting the desired integral of f by I(f) := f(x) dx quadrature formulas for approximating I(f) have the form I(f) Q N (f) := with nodes {ξ j } N j= and weights {w j} N j=. N w j f(ξ j ) j= Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/
12 Numerical Quadrature Some Simple Quadrature Rules Trapezoidal Rule Composite Trapezoidal Rule Q T N (f) = h f() N + 2 j=2 N-Point Gauß-Legendre Rule Q T (f) = Q GL N (f) = f(jh) + f() 2 N j= f() + f(). 2 w j 2 f, h = N, N 2. ( ) ξj + with {ξ j } N j= [, ] the zeros of the Legendre polynomial of degree N and w j chosen such that the quadrature formula is interpolatory. 2 Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/
13 Numerical Quadrature Some Quadrature Rules Example: f(x) = cos x. 2 Composite Trapezoidal Gauss Legendre 4 (If Q N f)/if N Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 2
14 Multivariate Quadrature, d = 2 Simplest idea: product rule of n-point D quadrature rule: I(f) = f(x ) dx = f(x, y) dy dx [,] 2 x= y= n w j f(x, η j ) dx x= k= j= j= n n w k w j f(ξ k, η j ) Result: N = n 2 -point quadrature rule in 2 dimensions. Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 3
15 Numerical Quadrature Benchmark Problems A widely used set of benchmark problems for multivariate integration schemes [Genz (984)], [Genz (987)] uses the functions ( ) d f (x ) = cos 2πu + a i x i Oscillatory f 2 (x ) = f 3 (x ) = i= i= d a 2 Product Peak i + (x i u i ) 2 ( f 4 (x ) = exp + ( d i= a i x i ) d+ ) d a 2 i (x i u i ) 2 i= Corner Peak Gaussian Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 4
16 Numerical Quadrature Benchmark Problems (cont d) f 5 (x ) = exp ( ) d a i x i u i i= { x > u x d > u d, f 6 (x ) = ( d ) exp i= a ix i, otherwise. C Function Discontinuous u, u = [u,..., u d ] chosen randomly from [, ], inneffective parameters integration problem more difficult with increasing a Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 5
17 f Numerical Quadrature Benchmark Problems Oscillatory: a = 2*rand(d,); u= x y Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 6
18 f Numerical Quadrature Benchmark Problems Product Peak: a = 2*ones(d,); u =.5*ones(d,); y.4.2 Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ x.6.8
19 Numerical Quadrature Benchmark Problems Corner Peak: a = 2*ones(d,); u =.5*ones(d,); f y x Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 8
20 f Numerical Quadrature Benchmark Problems Gaussian: a = 5*ones(d,); u =.5*ones(d,); y.4.2 Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ x.6.8
21 f Numerical Quadrature Benchmark Problems C Function: a = 3*ones(d,); u =.5*ones(d,); y.4.2 Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ x.6.8
22 f Numerical Quadrature Benchmark Problems Discontinuous: a = 5*ones(d,); u =.5*ones(d,); y.8 Oliver Ernst (TU Freiberg) Hochdimensionale Integration x Wintersemester 2/ 2
23 Numerical Quadrature 2D Benchmark Examples With Product Gauss-Legendre (If Q N f)/if (If Q N f)/if 8 (If Q N f)/if N Corner Peak N Gaussian N C Function Summary: For small d numerical integration is tractable. Issues not mentioned: Geometry Adaptivity Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 22
24 Numerical Quadrature in Higher Dimensions For the CMO valuation example d = 36. Assuming integrand smooth (it is not) such that 2 nodes in each direction suffice, a product rule requires function evaluations. N = The supercomputer currently leading the Top 5 list runs at a theoretical peak performance of 2.3 Petaflops, i.e., operations per second. This machine would require 93 seconds, or 3 85 years. The estimated numer of atoms in the observable universe: 8. This fundamental limitation is known as the curse of dimensionality [Bellman (957)]. Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 23
25 Numerical Quadrature in Higher Dimensions Monte Carlo Integration [Neumann, Ulam & Metropolis (946)] Idea: Choose N quadrature nodes {ξ j } N j= [, ]d sampled from a uniform distribution and approximate I(f) = f(x ) dx Q N (f) := N f(ξ [,] N j ). d j= Terminology: I(f) : mean of f, σ 2 (f) := I(f 2 ) [I(f)] 2 : variance of f 64 random nodes. Known: error estimate Figure random points Carlo method is that it does not suffer from th curse of dimensionality: in particular the O(N /2 convergence rate, while slow and erratic, does no depend on the dimension d so long as f is squar I(f) Q N (f) σ(f) (independent of d!), σ(f) Q N (f 2 ) (Q N f) 2. N Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 24
26 Numerical Quadrature in Higher Dimensions Monte Carlo Integration D Examples 2 MC error MC error estimate 2 MC error MC error estimate If Q N f /2 If Q N f 3 4 / N f(x) = cos x N f(x) = x 2 Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 25
27 Numerical Quadrature in Higher Dimensions Monte Carlo Integration 2D Examples 2 MC error MC error estimate 2 MC error MC error estimate If Q N f /2 If Q N f 3 4 / N Corner Peak N C Function Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 26
28 ntroduction concludes with a tion Numerical of the 36-dimensional Quadrature in Higher Dimensions above, since many features of mmon Sparse to problems Gridsfrom math- Figure. Product rule with 64 points. oncerned with the valuation of ge-backed securities held by a [Smolyak (963)] mers of the bank borrow money rs. Each month every customer ay the loan, and of course re- sparse-grid methods, which are generalizations of Idea: Organize points in a product rule in a hierarchical way and only use a construction first devised by Smolyak. Figure 2 certain levels of nodes. shows an example of a regular sparse grid with 49 rly will reduce its value to the odel, the proportion of those ay will depend on the interest e higher the interest rate, the e to repay the loan. The interto follow a (geometric) Brownth-by-month changes in the indom variables, so the present of mortgages is a (suitably dissional expected value, because ble repayment occasions. This -dimensional Euclidean space into an integral over the 36- be by an appropriate variable y other high-dimensional probluding options of all varieties) sional expected values, with the ng either from discretization in cause there are multiple assets acteristics, or both. One feasible strategy is to organize the points of a product rule in a hierarchical way and use only a few levels of points. This is the principle behind points. Figure 2. Sparse grid with 49 points. Smolyak sparse grid with 49 points Modern sparse-grid methods are dimensionadaptive: they find the important dimensions automatically and use more integration points in those dimensions. For details on sparse-grid methods, we refer readers to the recent survey article re Possible? e can approximate an integral by Bungartz and Griebel [2]. Another possible strategy is the Monte Carlo Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 27
29 ntroduction concludes with a tion Numerical of the 36-dimensional Quadrature in Higher Dimensions above, since many features of mmon Sparse to problems Gridsfrom math- Figure. Product rule with 64 points. oncerned with the valuation of ge-backed securities held by a [Smolyak (963)] mers of the bank borrow money rs. Each month every customer ay the loan, and of course re- sparse-grid methods, which are generalizations of Idea: Organize points in a product rule in a hierarchical way and only use a construction first devised by Smolyak. Figure 2 certain levels of nodes. shows an example of a regular sparse grid with 49 rly will reduce its value to the odel, the proportion of those ay will depend on the interest e higher the interest rate, the e to repay the loan. The interto follow a (geometric) Brownth-by-month changes in the indom variables, so the present of mortgages is a (suitably dissional expected value, because ble repayment occasions. This -dimensional Euclidean space into an integral over the 36- be by an appropriate variable y other high-dimensional probluding options of all varieties) sional expected values, with the ng either from discretization in cause there are multiple assets acteristics, or both. One feasible strategy is to organize the points of a product rule in a hierarchical way and use only a few levels of points. This is the principle behind points. Figure 2. Sparse grid with 49 points. Smolyak sparse grid with 49 points Modern sparse-grid methods are dimensionadaptive: they find the important dimensions automatically and use more integration points in those dimensions. For details on sparse-grid methods, we refer readers to the recent survey article re Possible? e can approximate an integral by Bungartz and Griebel [2]. Another possible strategy is the Monte Carlo Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 28
30 Quasi-Monte Carlo methods are equal weight rules, just like the Monte Carlo method, except Numerical Quadrature in Higher Dimensions Figure 5. Fibonacci lattice rule with 5 In this article we focus on lattice rule that the points x,...,x N are now designed in a Quasi-Monte continuing we should acknowledge a disa clever way Carlo to be Methods more uniformly distributed size, each rectangle than will include exactly one point. random of all deterministic methods, such as sp points so that a convergence Details rate on close both totheory and construction of nets and O(N ) is possible. or quasi-monte Carlo, when compared to t (Note, however, that sequences the implied can be found in the book of Niederreiter on but d.) [3]. Figure design 4 the nodes to be more evenly Carlo method namely, that they come Idea: constant Use equal can depend weights exponentially like in MC shows the first 64 points of a 2-dimensional Lattice Sobol rules (deterministically) distributed than random. are a any different practical kind information of quasi-monte about the erro sequence, priori estimates involving, for example, h the first example of the Carlo now method. widelythe points x,...,x N of a lattice rivatives of f are essentially never useful renowned concept of (t,m,s)-nets rule and are so (t,s) regular - that they form a group under tical error estimation.) This has led to a Can achieve sequences O(/N) established rate, by Niederreiter. but constant the operation generally of addition dependent modulo the integers. d. Figure 5 shows a lattice rule with 55 points. interest in hybrid methods that are essen terministic but that also have some el randomness, thereby seeking to capture efits of both approaches. We shall see an Figure random points. of a hybrid method later in the article. More on Lattice Rules There are many kinds of lattice rules (ind the product of left-rectangle rules is a lat but for our purposes it is enough to con the oldest and simplest kind, technical now as a rank- lattice rule, which takes (2) Q N f = N ({ f k z }) Figure 4. First Sobol 64 sequence points of 2D with Sobol 64 Figure sequence. 5. Fibonacci lattice rule with with points.. N N k= Informally, the basic points are now designed idea in ais to have In this the article right we focus points on lattice rules. Before Here z Z number of points in various subcubes. continuing For example, if in Figure 4 we divide the unit of all square deterministic into methods, such as sparse-grid we should acknowledge d is the a disadvantage generating vector, and t indicate that each component is to be rep its fractional part in [, ). In this case the strips of size by /64, then there or is quasi-monte exactly onecarlo, when compared to the Monte group formed by the points is the cyc point in each of the 64 strips, with Carlo any method namely, point on that they come without generated by {z/n}. Without loss of gene method is that it does not suffer from the of dimensionality: in particular the O(N /2 ) rgence rate, while slow and erratic, does not d on the dimension d so long as f is square able. Furthermore, it is cheap and easy to pron effective error estimate, since the first term the variance can be estimated by making use same function values as already used for apmating If. asi-monte Carlo methods are equal weight just like the Monte Carlo method, except he points x,...,x N way to be more uniformly distributed than m points so that a convergence rate close to ) is possible. (Note, however, that the implied ant can depend exponentially on d.) Figure 4 Oliver Ernst (TU Freiberg) Hochdimensionale any practical Integration information about the Wintersemester error. (The 2/ a 29
31 Numerical Quadrature in Higher Dimensions Effective Dimension Consider function f(x ) = f(x, x 2, x 3, x 4 ) = x cos x 2 + x 3. What is the dimensionality of f? (2, 3, 4?) Generally: can write any d-dimensional function as the 2 d -term sum f(x ) = f α (x α ), α {,...,d} where f α depends on only α of the original d variables. Uniquess from requirement f α (x α ) dx j = j α (ANOVA decomposition). Idea: Approximate f by neglecting unimportant terms for efficient quadrature.. Oliver Ernst (TU Freiberg) Hochdimensionale Integration Wintersemester 2/ 3
Hochdimensionale Integration
Oliver Ernst Institut für Numerische Mathematik und Optimierung Hochdimensionale Integration 14-tägige Vorlesung im Wintersemester 2010/11 im Rahmen des Moduls Ausgewählte Kapitel der Numerik Contents
More informationLifting the Curse of Dimensionality
Lifting the Curse of Dimensionality Frances Y. Kuo and Ian H. Sloan Introduction Richard Bellman [1] coined the phrase the curse of dimensionality to describe the extraordinarily rapid growth in the difficulty
More informationProgress 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 informationQuasi-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 informationGlobal Focus on Knowledge Lecture Series. Finance and Math. Practical Problems of Finance and Integral Computations of Higher Dimensions
Global Focus on Knowledge Lecture Series Finance and Math The 9 th Lecture Practical Problems of Finance and Integral Computations of Higher Dimensions How to Compute Integration Math Application New Mathematiques
More informationSlow Growth for Gauss Legendre Sparse Grids
Slow Growth for Gauss Legendre Sparse Grids John Burkardt, Clayton Webster April 4, 2014 Abstract A sparse grid for multidimensional quadrature can be constructed from products of 1D rules. For multidimensional
More informationDifferentiation and Integration
Differentiation and Integration (Lectures on Numerical Analysis for Economists II) Jesús Fernández-Villaverde 1 and Pablo Guerrón 2 February 12, 2018 1 University of Pennsylvania 2 Boston College Motivation
More informationSlow Exponential Growth for Clenshaw Curtis Sparse Grids jburkardt/presentations/... slow growth paper.
Slow Exponential Growth for Clenshaw Curtis Sparse Grids http://people.sc.fsu.edu/ jburkardt/presentations/...... slow growth paper.pdf John Burkardt, Clayton Webster April 30, 2014 Abstract A widely used
More informationAccuracy, Precision and Efficiency in Sparse Grids
John, Information Technology Department, Virginia Tech.... http://people.sc.fsu.edu/ jburkardt/presentations/ sandia 2009.pdf... Computer Science Research Institute, Sandia National Laboratory, 23 July
More informationOrthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals
1/33 Orthogonal Polynomials, Quadratures & Sparse-Grid Methods for Probability Integrals Dr. Abebe Geletu May, 2010 Technische Universität Ilmenau, Institut für Automatisierungs- und Systemtechnik Fachgebiet
More informationPascal s Triangle on a Budget. Accuracy, Precision and Efficiency in Sparse Grids
Covering : Accuracy, Precision and Efficiency in Sparse Grids https://people.sc.fsu.edu/ jburkardt/presentations/auburn 2009.pdf... John Interdisciplinary Center for Applied Mathematics & Information Technology
More informationQuasi- 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 informationOn ANOVA expansions and strategies for choosing the anchor point
On ANOVA expansions and strategies for choosing the anchor point Zhen Gao and Jan S. Hesthaven August 31, 2010 Abstract The classic Lebesgue ANOVA expansion offers an elegant way to represent functions
More information4. Numerical Quadrature. Where analytical abilities end... continued
4. Numerical Quadrature Where analytical abilities end... continued Where analytical abilities end... continued, November 30, 22 1 4.3. Extrapolation Increasing the Order Using Linear Combinations Once
More informationDimension Adaptive Tensor Product Quadrature
Dimension Adaptive Tensor Product Quadrature Thomas Gerstner and Michael Griebel Department for Applied Mathematics University of Bonn {gerstner,griebel}@iam.uni-bonn.de Abstract We consider the numerical
More informationSection 6.6 Gaussian Quadrature
Section 6.6 Gaussian Quadrature Key Terms: Method of undetermined coefficients Nonlinear systems Gaussian quadrature Error Legendre polynomials Inner product Adapted from http://pathfinder.scar.utoronto.ca/~dyer/csca57/book_p/node44.html
More informationNumerical Probabilistic Analysis under Aleatory and Epistemic Uncertainty
Numerical Probabilistic Analysis under Aleatory and Epistemic Uncertainty Boris S. Dobronets Institute of Space and Information Technologies, Siberian Federal University, Krasnoyarsk, Russia BDobronets@sfu-kras.ru
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationPart 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 informationTutorial 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 informationPHYS-4007/5007: Computational Physics Course Lecture Notes Appendix G
PHYS-4007/5007: Computational Physics Course Lecture Notes Appendix G Dr. Donald G. Luttermoser East Tennessee State University Version 7.0 Abstract These class notes are designed for use of the instructor
More informationANALYSIS OF THE ANOMALY OF ran10 GENERATOR IN MONTE CARLO PRICING OF FINANCIAL DERIVATIVES
Journal of the Operations Research Society of Japan Vol. 41, No. 3, September 1998 1998 The Operations Research Society of Japan ANALYSIS OF THE ANOMALY OF ran10 GENERATOR IN MONTE CARLO PRICING OF FINANCIAL
More informationScenario generation for stochastic optimization problems via the sparse grid method
Scenario generation for stochastic optimization problems via the sparse grid method Michael Chen Sanjay Mehrotra Dávid Papp March 7, 0 Abstract We study the use of sparse grids methods for the scenario
More informationFast evaluation of mixed derivatives and calculation of optimal weights for integration. Hernan Leovey
Fast evaluation of mixed derivatives and calculation of optimal weights for integration Humboldt Universität zu Berlin 02.14.2012 MCQMC2012 Tenth International Conference on Monte Carlo and Quasi Monte
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Lecture 9: Sensitivity Analysis ST 2018 Tobias Neckel Scientific Computing in Computer Science TUM Repetition of Previous Lecture Sparse grids in Uncertainty Quantification
More informationA note on multivariate Gauss-Hermite quadrature
A note on multivariate Gauss-Hermite quadrature Peter Jäckel 1th May 5 1 Introduction Gaussian quadratures are an ingenious way to approximate the integral of an unknown function f(x) over a specified
More informationApproximation of High-Dimensional Rank One Tensors
Approximation of High-Dimensional Rank One Tensors Markus Bachmayr, Wolfgang Dahmen, Ronald DeVore, and Lars Grasedyck March 14, 2013 Abstract Many real world problems are high-dimensional in that their
More informationAdvanced Monte Carlo Methods - Computational Challenges
Advanced Monte Carlo Methods - Computational Challenges Ivan Dimov IICT, Bulgarian Academy of Sciences Sofia, February, 2012 Aims To cover important topics from probability theory and statistics. To present
More informationFrom Random Numbers to Monte Carlo. Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that
From Random Numbers to Monte Carlo Random Numbers, Random Walks, Diffusion, Monte Carlo integration, and all that Random Walk Through Life Random Walk Through Life If you flip the coin 5 times you will
More information3. Numerical Quadrature. Where analytical abilities end...
3. Numerical Quadrature Where analytical abilities end... Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 32 3.1. Preliminary Remarks The Integration Problem Numerical quadrature denotes numerical
More informationCS 450 Numerical Analysis. Chapter 8: Numerical Integration and Differentiation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationIntroduction to Discrepancy Theory
Introduction to Discrepancy Theory Aleksandar (Sasho) Nikolov University of Toronto Sasho Nikolov (U of T) Discrepancy 1 / 23 Outline 1 Monte Carlo 2 Discrepancy and Quasi MC 3 Combinatorial Discrepancy
More informationSparse Grid Collocation for Uncertainty Quantification
Sparse Grid Collocation for Uncertainty Quantification John Burkardt Department of Scientific Computing Florida State University Ajou University 14 August 2013 http://people.sc.fsu.edu/ jburkardt/presentations/......
More informationLinear Models for Regression. Sargur Srihari
Linear Models for Regression Sargur srihari@cedar.buffalo.edu 1 Topics in Linear Regression What is regression? Polynomial Curve Fitting with Scalar input Linear Basis Function Models Maximum Likelihood
More informationCh. 03 Numerical Quadrature. Andrea Mignone Physics Department, University of Torino AA
Ch. 03 Numerical Quadrature Andrea Mignone Physics Department, University of Torino AA 2017-2018 Numerical Quadrature In numerical analysis quadrature refers to the computation of definite integrals. y
More informationApplication of QMC methods to PDEs with random coefficients
Application of QMC methods to PDEs with random coefficients a survey of analysis and implementation Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with Ivan Graham
More informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More information19 : Slice Sampling and HMC
10-708: Probabilistic Graphical Models 10-708, Spring 2018 19 : Slice Sampling and HMC Lecturer: Kayhan Batmanghelich Scribes: Boxiang Lyu 1 MCMC (Auxiliary Variables Methods) In inference, we are often
More informationAPPLIED 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 informationUncertainty Propagation
Setting: Uncertainty Propagation We assume that we have determined distributions for parameters e.g., Bayesian inference, prior experiments, expert opinion Ṫ 1 = 1 - d 1 T 1 - (1 - ")k 1 VT 1 Ṫ 2 = 2 -
More informationQuasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient
Quasi-optimal and adaptive sparse grids with control variates for PDEs with random diffusion coefficient F. Nobile, L. Tamellini, R. Tempone, F. Tesei CSQI - MATHICSE, EPFL, Switzerland Dipartimento di
More informationScientific Computing: Monte Carlo
Scientific Computing: Monte Carlo Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 April 5th and 12th, 2012 A. Donev (Courant Institute)
More informationQuadrature for Uncertainty Analysis Stochastic Collocation. What does quadrature have to do with uncertainty?
Quadrature for Uncertainty Analysis Stochastic Collocation What does quadrature have to do with uncertainty? Quadrature for Uncertainty Analysis Stochastic Collocation What does quadrature have to do with
More informationScientific Computing: Numerical Integration
Scientific Computing: Numerical Integration Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Fall 2015 Nov 5th, 2015 A. Donev (Courant Institute) Lecture
More informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationChapter 4: Monte-Carlo Methods
Chapter 4: Monte-Carlo Methods A Monte-Carlo method is a technique for the numerical realization of a stochastic process by means of normally distributed random variables. In financial mathematics, it
More informationDimension-adaptive sparse grid for industrial applications using Sobol variances
Master of Science Thesis Dimension-adaptive sparse grid for industrial applications using Sobol variances Heavy gas flow over a barrier March 11, 2015 Ad Dimension-adaptive sparse grid for industrial
More informationSTRONGLY NESTED 1D INTERPOLATORY QUADRATURE AND EXTENSION TO ND VIA SMOLYAK METHOD
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 2018, Glasgow, UK STRONGLY NESTED 1D INTERPOLATORY QUADRATURE AND
More informationMonte Carlo Method for Numerical Integration based on Sobol s Sequences
Monte Carlo Method for Numerical Integration based on Sobol s Sequences Ivan Dimov and Rayna Georgieva Institute of Information and Communication Technologies, Bulgarian Academy of Sciences, 1113 Sofia,
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationUncertainty Quantification in Computational Science
DTU 2010 - Lecture I Uncertainty Quantification in Computational Science Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu Objective of lectures The main objective of these lectures are To offer
More informationInterpolation and function representation with grids in stochastic control
with in classical Sparse for Sparse for with in FIME EDF R&D & FiME, Laboratoire de Finance des Marchés de l Energie (www.fime-lab.org) September 2015 Schedule with in classical Sparse for Sparse for 1
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC55 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 0-) Time Allowed: Hours INSTRUCTIONS TO CANDIDATES. This examination paper contains FIVE questions and comprises
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
More informationHot new directions for QMC research in step with applications
Hot new directions for QMC research in step with applications Frances Kuo f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia Frances Kuo @ UNSW Australia p.1 Outline Motivating example
More informationNew 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 informationSample of Ph.D. Advisory Exam For MathFinance
Sample of Ph.D. Advisory Exam For MathFinance Students who wish to enter the Ph.D. program of Mathematics of Finance are required to take the advisory exam. This exam consists of three major parts. The
More informationNumerical integration and importance sampling
2 Numerical integration and importance sampling 2.1 Quadrature Consider the numerical evaluation of the integral I(a, b) = b a dx f(x) Rectangle rule: on small interval, construct interpolating function
More informationChapter 4: Monte Carlo Methods. Paisan Nakmahachalasint
Chapter 4: Monte Carlo Methods Paisan Nakmahachalasint Introduction Monte Carlo Methods are a class of computational algorithms that rely on repeated random sampling to compute their results. Monte Carlo
More informationA 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 informationReduction of Variance. Importance Sampling
Reduction of Variance As we discussed earlier, the statistical error goes as: error = sqrt(variance/computer time). DEFINE: Efficiency = = 1/vT v = error of mean and T = total CPU time How can you make
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationTRB 2005: For Presentation and Publication Paper # Final Submission: April 1, 2005 Word Count: 8,191
Simulation Estimation of Mixed Discrete Choice Models Using Randomized Quasi-Monte Carlo Sequences: A Comparison of Alternative Sequences, Scrambling Methods, and Uniform-to-Normal Variate Transformation
More informationImportance 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 informationIn manycomputationaleconomicapplications, one must compute thede nite n
Chapter 6 Numerical Integration In manycomputationaleconomicapplications, one must compute thede nite n integral of a real-valued function f de ned on some interval I of
More informationMathematics. Pre-Leaving Certificate Examination, Paper 1 Higher Level Time: 2 hours, 30 minutes. 300 marks L.17 NAME SCHOOL TEACHER
L.7 NAME SCHOOL TEACHER Pre-Leaving Certificate Examination, 205 Name/v Printed Checke To: Update Name/v Comple Paper Higher Level Time: 2 hours, 30 minutes 300 marks For examiner Question 2 School stamp
More informationAlgebra Readiness. Curriculum (445 topics additional topics)
Algebra Readiness This course covers the topics shown below; new topics have been highlighted. Students navigate learning paths based on their level of readiness. Institutional users may customize the
More informationMachine Learning. Nonparametric Methods. Space of ML Problems. Todo. Histograms. Instance-Based Learning (aka non-parametric methods)
Machine Learning InstanceBased Learning (aka nonparametric methods) Supervised Learning Unsupervised Learning Reinforcement Learning Parametric Non parametric CSE 446 Machine Learning Daniel Weld March
More informationMulti-Element Probabilistic Collocation Method in High Dimensions
Multi-Element Probabilistic Collocation Method in High Dimensions Jasmine Foo and George Em Karniadakis Division of Applied Mathematics, Brown University, Providence, RI 02912 USA Abstract We combine multi-element
More informationYou are allowed two hours to answer this question paper. All questions are compulsory.
Examination Question and Answer Book Write here your full examination number Centre Code: Hall Code: Desk Number: Foundation Level 3c Business Mathematics FBSM 0 May 00 Day 1 late afternoon INSTRUCTIONS
More informationMath 101 Final Exam Review Solutions. Eric Schmutz
Math 101 Final Exam Review Solutions Eric Schmutz Problem 1. Write an equation of the line passing through (,7) and (-1,1). Let (x 1, y 1 ) = (, 7) and (x, y ) = ( 1, 1). The slope is m = y y 1 x x 1 =
More informationCalifornia Algebra 1
California Algebra 1 This course covers the topics shown below. Students navigate learning paths based on their level of readiness. Institutional users may customize the scope and sequence to meet curricular
More informationMiddle School Math Course 2
Middle School Math Course 2 This course covers the topics shown below. Students navigate learning paths based on their level of readiness. Institutional users may customize the scope and sequence to meet
More informationSample Institution Memphis, TN
Morgan Keegan FIXED INCOME RESEARCH Morgan Keegan & Company, Inc. Members New York Stock Exchange Memphis, TN Asset & Liability Analysis MORGAN KEEGAN & COMPANY, INC. Memphis, Tennessee June 2, 2008 The
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationGENERALIZED ANNUITIES AND ASSURANCES, AND INTER-RELATIONSHIPS. BY LEIGH ROBERTS, M.Sc., ABSTRACT
GENERALIZED ANNUITIES AND ASSURANCES, AND THEIR INTER-RELATIONSHIPS BY LEIGH ROBERTS, M.Sc., A.I.A ABSTRACT By the definition of generalized assurances and annuities, the relation is shown to be the simplest
More informationSparse Grid Collocation for Uncertainty Quantification
Sparse Grid Collocation for Uncertainty Quantification John Burkardt Department of Scientific Computing Florida State University 13 February 2014 http://people.sc.fsu.edu/ jburkardt/presentations/......
More informationSampling Distributions Allen & Tildesley pgs and Numerical Recipes on random numbers
Sampling Distributions Allen & Tildesley pgs. 345-351 and Numerical Recipes on random numbers Today I explain how to generate a non-uniform probability distributions. These are used in importance sampling
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
More informationIntegration. Topic: Trapezoidal Rule. Major: General Engineering. Author: Autar Kaw, Charlie Barker.
Integration Topic: Trapezoidal Rule Major: General Engineering Author: Autar Kaw, Charlie Barker 1 What is Integration Integration: The process of measuring the area under a function plotted on a graph.
More informationLow 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 informationAdvanced Algorithms for Multidimensional Sensitivity Studies of Large-scale Air Pollution Models Based on Sobol Sequences
Advanced Algorithms for Multidimensional Sensitivity Studies of Large-scale Air Pollution Models Based on Sobol Sequences I. Dimov a,, R. Georgieva a, Tz. Ostromsky a, Z. Zlatev b a Department of Parallel
More informationCSE446: non-parametric methods Spring 2017
CSE446: non-parametric methods Spring 2017 Ali Farhadi Slides adapted from Carlos Guestrin and Luke Zettlemoyer Linear Regression: What can go wrong? What do we do if the bias is too strong? Might want
More informationA Polynomial Chaos Approach to Robust Multiobjective Optimization
A Polynomial Chaos Approach to Robust Multiobjective Optimization Silvia Poles 1, Alberto Lovison 2 1 EnginSoft S.p.A., Optimization Consulting Via Giambellino, 7 35129 Padova, Italy s.poles@enginsoft.it
More informationIntroduction to Uncertainty Quantification in Computational Science Handout #3
Introduction to Uncertainty Quantification in Computational Science Handout #3 Gianluca Iaccarino Department of Mechanical Engineering Stanford University June 29 - July 1, 2009 Scuola di Dottorato di
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationSTRONG 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 informationThis chapter illustrates some of the basic features of Mathematica useful for finite mathematics and business calculus. Example
The Mathematics Companion for Finite Mathematics and Business Calculus is a dictionary-like reference guide for learning and applying mathematical ideas, techniques, and formulas with the help of Mathematica,
More informationSlow Growth for Sparse Grids
Slow Growth for Sparse Grids John Burkardt, Clayton Webster, Guannan Zhang... http://people.sc.fsu.edu/ jburkardt/presentations/... slow growth 2014 savannah.pdf... SIAM UQ Conference 31 March - 03 April
More informationAPPM/MATH Problem Set 6 Solutions
APPM/MATH 460 Problem Set 6 Solutions This assignment is due by 4pm on Wednesday, November 6th You may either turn it in to me in class or in the box outside my office door (ECOT ) Minimal credit will
More informationModern Methods of Data Analysis - WS 07/08
Modern Methods of Data Analysis Lecture III (29.10.07) Contents: Overview & Test of random number generators Random number distributions Monte Carlo The terminology Monte Carlo-methods originated around
More information6 Markov Chain Monte Carlo (MCMC)
6 Markov Chain Monte Carlo (MCMC) The underlying idea in MCMC is to replace the iid samples of basic MC methods, with dependent samples from an ergodic Markov chain, whose limiting (stationary) distribution
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationMonte Carlo Methods for Statistical Inference: Variance Reduction Techniques
Monte Carlo Methods for Statistical Inference: Variance Reduction Techniques Hung Chen hchen@math.ntu.edu.tw Department of Mathematics National Taiwan University 3rd March 2004 Meet at NS 104 On Wednesday
More informationScenario generation Contents
Scenario generation W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Page 1 of 39 Tutorial, ICSP12, Halifax, August 14, 2010 (1) Aproximation issues (2) Generation
More informationMFM Practitioner Module: Quantitative Risk Management. John Dodson. September 23, 2015
MFM Practitioner Module: Quantitative Risk Management September 23, 2015 Mixtures Mixtures Mixtures Definitions For our purposes, A random variable is a quantity whose value is not known to us right now
More informationKrylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms
Krylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms 2. First Results and Algorithms Michael Eiermann Institut für Numerische Mathematik und Optimierung Technische
More information$690,440,549 Freddie Mac. Multiclass Certificates, Series 3450
Offering Circular Supplement (To Offering Circular Dated December 31, 2007) $690,440,549 Freddie Mac Multiclass Certificates, Series 3450 Offered Classes: REMIC Classes shown below and MACR Classes shown
More informationFast integration using quasi-random numbers. J.Bossert, M.Feindt, U.Kerzel University of Karlsruhe ACAT 05
Fast integration using quasi-random numbers J.Bossert, M.Feindt, U.Kerzel University of Karlsruhe ACAT 05 Outline Numerical integration Discrepancy The Koksma-Hlawka inequality Pseudo- and quasi-random
More informationNumerical integration - I. M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V
Numerical integration - I M. Peressi - UniTS - Laurea Magistrale in Physics Laboratory of Computational Physics - Unit V deterministic methods in 1D equispaced points (trapezoidal, Simpson...), others...
More information