M. Croci (Oxford), M. B. Giles (Oxford), P. E. Farrell (Oxford), M. E. Rognes (Simula) MCQMC July 4, 2018
|
|
- Amberly Powell
- 5 years ago
- Views:
Transcription
1 Efficient wite noise sampling and coupling for multilevel Monte Carlo M. Croci (Oxford), M. B. Giles (Oxford), P. E. Farrell (Oxford), M. E. Rognes (Simula) MCQMC July 4, 2018
2 Overview Introduction Wite noise sampling Numerical results Conclusions and furter work
3 Overview Introduction Wite noise sampling Numerical results Conclusions and furter work
4 Motivation Te motivation of our researc is te sampling of lognormal Gaussian fields. A Matérn Gaussian field (approximately) satisfies a linear elliptic SPDE of te form Lu = W, x D, ω Ω + BCs, were u = u(x, ω) and W is spatial wite noise. Oter approaces can be used (wit pros and cons), but we will not discuss tem ere.
5 Motivation Te motivation of our researc is te sampling of lognormal Gaussian fields. A Matérn Gaussian field (approximately) satisfies a linear elliptic SPDE of te form Lu = W, x D, ω Ω + BCs, were u = u(x, ω) and W is spatial wite noise. Oter approaces can be used (wit pros and cons), but we will not discuss tem ere. Te same tecniques can be used to solve a more general class of SPDEs, e.g. N(u) + Lu = W, x D, ω Ω + BCs. In tis case solving means to compute E[P(u)] for some functional P of te solution.
6 Motivation Te motivation of our researc is te sampling of lognormal Gaussian fields. A Matérn Gaussian field (approximately) satisfies a linear elliptic SPDE of te form Lu = W, x D, ω Ω + BCs, were u = u(x, ω) and W is spatial wite noise. Oter approaces can be used (wit pros and cons), but we will not discuss tem ere. Te same tecniques can be used to solve a more general class of SPDEs, e.g. N(u) + Lu = W, x D, ω Ω + BCs. In tis case solving means to compute E[P(u)] for some functional P of te solution. Common applications: finance, geology, meteorology, biology...
7 Motivation Te motivation of our researc is te sampling of lognormal Gaussian fields. A Matérn Gaussian field (approximately) satisfies a linear elliptic SPDE of te form Lu = W, x D, ω Ω + BCs, were u = u(x, ω) and W is spatial wite noise. Oter approaces can be used (wit pros and cons), but we will not discuss tem ere. Te same tecniques can be used to solve a more general class of SPDEs, e.g. N(u) + Lu = W, x D, ω Ω + BCs. In tis case solving means to compute E[P(u)] for some functional P of te solution. Common applications: finance, geology, meteorology, biology... Main issue: sampling W is ard!
8 Motivation Te motivation of our researc is te sampling of lognormal Gaussian fields. A Matérn Gaussian field (approximately) satisfies a linear elliptic SPDE of te form Lu = W, x D, ω Ω + BCs, were u = u(x, ω) and W is spatial wite noise. Oter approaces can be used (wit pros and cons), but we will not discuss tem ere. Te same tecniques can be used to solve a more general class of SPDEs, e.g. N(u) + Lu = W, x D, ω Ω + BCs. In tis case solving means to compute E[P(u)] for some functional P of te solution. Common applications: finance, geology, meteorology, biology... Main issue: sampling Te efficient sampling of W is ard! W is te focus of tis talk.
9 Wite noise (1D) WARNING! Point evaluation not defined!
10 University of Oxford Matematical Institute Wite noise (2D) WARNING! Point evaluation not defined!
11 University of Oxford Matematical Institute Wite noise (2D) WARNING! Point evaluation not defined! IDEA! Avoid point evaluation by integrating W.
12 University of Oxford Matematical Institute Wite Noise (practical definition) Definition (Spatial Wite Noise W ) R For any φ L2 (D), define W, φi := D W φ dx. For any φi, φj L2 (D), bi = W, φi i, bj = W, φj i are zero-mean Gaussian random variables, wit, Z E[bi bj ] = φi φj dx =: Mij, b N (0, M). (1) D
13 Finite element (FEM) framework Wen solving SPDEs (see 1 st slide) wit FEM, we get (for linear problems) Discrete weak form: find u V s.t. for all v V, Were V = span({φ i } n i=0 ), (e.g. wit Lagrange elements). a(u, v ) = W, v, (2)
14 Finite element (FEM) framework Wen solving SPDEs (see 1 st slide) wit FEM, we get (for linear problems) Discrete weak form: find u V s.t. for all v V, Were V = span({φ i } n i=0 ), (e.g. wit Lagrange elements). FEM linear system: u = i u iφ i, u = [u 0,..., u n ] T, were te entries of b are given by, wit b N (0, M) as before. M is te mass matrix of V. a(u, v ) = W, v, (2) Au = b(ω), (3) W, φ i (ω) = b i (ω), (4)
15 Multilevel Monte Carlo [Giles 2008] + FEM framework For MLMC, we ave two approximation levels l and l 1. For any particular ω Ω, we need to solve: find u l V l, ul 1 V l 1 s.t. for all v l V l, v l 1 V l 1, a(u l 1 a(u l, v l ) = W, v l (ω), (5), v l 1 ) = W, v l 1 (ω). (6)
16 Multilevel Monte Carlo [Giles 2008] + FEM framework For MLMC, we ave two approximation levels l and l 1. For any particular ω Ω, we need to solve: find u l V l, ul 1 V l 1 s.t. for all v l V l, v l 1 V l 1, a(u l 1 Tis yields te linear system [ A l 0 0 A l 1 a(u l, v l ) = W, v l (ω), (5), v l 1 ) = W, v l 1 ] [ u l u l 1 ] = (ω). (6) [ b l b l 1 ] = b,
17 Multilevel Monte Carlo [Giles 2008] + FEM framework For MLMC, we ave two approximation levels l and l 1. For any particular ω Ω, we need to solve: find u l V l, ul 1 V l 1 s.t. for all v l V l, v l 1 V l 1, a(u l 1 Tis yields te linear system [ A l 0 0 A l 1 a(u l, v l ) = W, v l (ω), (5), v l 1 ) = W, v l 1 ] [ u l u l 1 ] = (ω). (6) [ b l b l 1 ] = b, were b N (0, M). Let V l = span({φl i }n l l 1 i=0 ) and V = span({φ l 1 i [ ] M l M l,l 1 M = (M l,l 1 ) T M l 1, M l,l 1 ij = φ l i φ l 1 j dx. } n l 1 i=0 ), ten
18 Multilevel Monte Carlo [Giles 2008] + FEM framework For MLMC, we ave two approximation levels l and l 1. For any particular ω Ω, we need to solve: find u l V l, ul 1 V l 1 s.t. for all v l V l, v l 1 V l 1, a(u l 1 Tis yields te linear system [ A l 0 0 A l 1 a(u l, v l ) = W, v l (ω), (5), v l 1 ) = W, v l 1 ] [ u l u l 1 ] = (ω). (6) [ b l b l 1 ] = b, were b N (0, M). Let V l = span({φl i }n l l 1 i=0 ) and V = span({φ l 1 i [ ] M l M l,l 1 M = (M l,l 1 ) T M l 1, M l,l 1 ij = φ l i φ l 1 j dx. } n l 1 i=0 NOTE: we do not require te FEM approximation subspaces to be nested! ), ten
19 Overview Introduction Wite noise sampling Numerical results Conclusions and furter work
20 Te callenge SAMPLING PROBLEM 1: single level realisations: sample b N (0, M), were M is te mass matrix of V. SAMPLING PROBLEM 2: coupled realisations: sample b N (0, M), were M is te block mass matrix given by V l and V l 1.
21 How to sample b? Sampling b is ard!
22 How to sample b? Sampling b is ard! Naïve approac - Factorise M = HH T (cubic complexity!) and set b = Hz, wit z N (0, I ). E[bb T ] = E[Hz(Hz) T ] = HE[zz T ]H T = HIH T = M.
23 How to sample b? Sampling b is ard! Naïve approac - Factorise M = HH T (cubic complexity!) and set b = Hz, wit z N (0, I ). E[bb T ] = E[Hz(Hz) T ] = HE[zz T ]H T = HIH T = M. - Works well if M diagonal: previous work used eiter mass-lumping [Lindgren, Rue and Lindström 2009], piecewise constant elements [Osborn, Vassilevski and Villa 2017] or a piecewise constant approximation of wite noise [Drzisga, et al. 2017, Du and Zang 2002]. - We do not require M to be diagonal (and we do not approximate wite noise). - We can sample b wit linear complexity.
24 How to sample b? Sampling b is ard! Naïve approac - Factorise M = HH T (cubic complexity!) and set b = Hz, wit z N (0, I ). E[bb T ] = E[Hz(Hz) T ] = HE[zz T ]H T = HIH T = M. - Works well if M diagonal: previous work used eiter mass-lumping [Lindgren, Rue and Lindström 2009], piecewise constant elements [Osborn, Vassilevski and Villa 2017] or a piecewise constant approximation of wite noise [Drzisga, et al. 2017, Du and Zang 2002]. - We do not require M to be diagonal (and we do not approximate wite noise). - We can sample b wit linear complexity. IDEA! H does not need to be square, maybe we can find a more efficient factorisation!
25 Wite noise sampling: single level realisations SAMPLING PROBLEM 1: need to sample b N (0, M). Exploit te FEM assembly M 1 0 M = L T. 0 M.. 2 L = LT diag e (M e )L. (7)
26 Wite noise sampling: single level realisations SAMPLING PROBLEM 1: need to sample b N (0, M). Exploit te FEM assembly N (0, M) b = L T b 1 b 2. = L T vstack e (b e ) (8)
27 Wite noise sampling: single level realisations SAMPLING PROBLEM 1: need to sample b N (0, M). Exploit te FEM assembly - Eac b e can be sampled as b e = H e z e wit z e N (0, I ) and H e H T e = M e. - b = L T vstack e (b e ) is N (0, M) since E[bb T ] = L T E[vstack e (b e )vstack e (b e ) T ]L = L T diag e (H e )diag e (H T e )L = L T diag e (M e )L = M. - If te mapping to te FEM reference element is affine (e.g. Lagrange elements on simplices) we ave tat M e / e = const on eac element and only one local factorisation is needed. Tis approac is trivially parallelisable!
28 Wite noise sampling: coupled realisations SAMPLING PROBLEM 2: need to sample b N (0, M), were M is now te block mixed mass matrix. Definition (Supermes, [Farrell 2009]) Let A and B be two (possibly non-nested) meses. Teir supermes S is one of teir common refinements. A and B are bot nested witin S.
29 Wite noise sampling: coupled realisations SAMPLING PROBLEM 2: need to sample b N (0, M), were M is now te block mixed mass matrix. - Factorise locally, tis time on eac supermes element. - Sample b on S, ten interpolate/project te result onto A and B (tis step can be performed locally). - Since A and B are nested witin S, tis operation is exact. Note tat A and B need not be nested. Previous work on wite noise coupling for MLMC used eiter a nested ierarcy [Drzisga et al. 2017, Osborn et al. 2017] or an algebraically constructed ierarcy of agglomerated meses [Osborn, Vassilevski and Villa 2017].
30 Complexity overview offline cost online cost (per sample) memory storage single level 0 (or O(m 3 N)) O(m 3 N) (or O(m 2 N)) O(m 2 ) (or O(m 2 N)) single l. (affine) O(m 3 ) O(m 2 N) O(m 2 ) coupled 0 (or O(m 3 N S )) O(m 3 N S ) (or O(m 2 N S )) O(m 2 ) (or O(m 2 N S )) coupled (affine) O(m 3 ) O(m 2 N S ) O(m 2 ) Table: Memory and cost complexity of our wite noise sampling strategy. In te non-affine case te cost per sample can be lowered by precomputing and storing te local factorisations (see entries in blue). N S is te number of supermes elements. In our experience wit MLMC, N S c d N l and c d = 2 (1D), c d = 2.5 (2D), c d = 45 (3D).
31 Overview Introduction Wite noise sampling Numerical results Conclusions and furter work
32 Numerical results: convergence of P(u) = u 2 L 2 (D) Consider te linear elliptic SPDE [Lindgren, Rue and Lindström 2009], [Bolin, Kircner and Kovács 2017], ( I κ 2 ) k u(x, ω) = η W, x D R d, ω Ω, ν = 2k d/2 > 0. We compute FEM solutions {u l }l=8 l=1 wit a non-nested ierarcy of subspaces {V l }l=8 l=
33 Numerical results: covariance convergence C(r) = E[u(x)u(y)] = ν 2 ν 1 Γ(ν) (κr)ν K ν (κr), r = x y 2, κ =, x, y D, λ
34 Overview Introduction Wite noise sampling Numerical results Conclusions and furter work
35 Conclusions and furter work Outlook - Wite noise is an extremely non-smoot object and is defined troug its integral. - We can sample single level wite noise realisations efficiently. - We can couple wite noise between different FEM approximation subspaces. A supermes construction is not needed in te nested case. - Te overall order of complexity is linear in te number of elements of te supermes and it can be trivially parallelised. Standard tecniques usually ave cubic complexity. Furter work: extensions to QMC and MLQMC. Paper: ttps://arxiv.org/abs/
36 References - Tank you for listening! [1] M. Croci, M. B. Giles, M. E. Rognes, and P. E. Farrell. Efficient wite noise sampling and coupling for multilevel Monte Carlo wit non-nested meses. Preprint, URL ttps://arxiv.org/abs/ [2] D. Bolin, K. Kircner, and M. Kovács. Numerical solutions of fractional elliptic stocastic PDEs wit spatial wite noise. Preprint, [3] D. Drzisga, B. Gmeiner, U. Ude, R. Sceicl, and B. Wolmut. Sceduling massively parallel multigrid for multilevel Monte Carlo metods. SIAM Journal of Scientific Computing, 39(5): , [4] Q. Du and T. Zang. Numerical approximation of some linear stocastic partial differential equations driven by special additive noises. SIAM Journal on Numerical Analysis, 40(4): , [5] M. B. Giles. Multilevel Monte Carlo pat simulation. Operations Researc, 56(3): , [6] F. Lindgren, H. Rue, and J. Lindström. An explicit link between Gaussian fields and Gaussian Markov random fields: te stocastic partial differential equation approac. Journal of te Royal Statistical Society: Series B (Statistical Metodology), 73(4): , [7] S. Osborn, P. S. Vassilevski, and U. Villa. A multilevel, ierarcical sampling tecnique for spatially correlated random fields. Preprint, [8] S. Osborn, P. Zulian, T. Benson, U. Villa, R. Krause, and P. S. Vassilevski. Scalable ierarcical PDE sampler for generating spatially correlated random fields using non-matcing meses. Preprint, [9] A. J. Waten. Realistic eigenvalue bounds for te Galerkin mass matrix. IMA Journal of Numerical Analysis, 7: , 1987.
WYSE Academic Challenge 2004 Sectional Mathematics Solution Set
WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationAn Adaptive Model Switching and Discretization Algorithm for Gas Flow on Networks
Procedia Computer Science 1 (21) (212) 1 1 1331 134 Procedia Computer Science www.elsevier.com/locate/procedia International Conference on Computational Science, ICCS 21 An Adaptive Model Switcing and
More informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationImproved Monte Carlo methods for computing failure probabilities of porous media flow systems
Improved Monte Carlo metods for computing failure probabilities of porous media flow systems Fritjof Fagerlund 1, Fredrik Hellman 2, Axel Målqvist 3, Auli Niemi 1 August 31, 2015 Abstract We study improvements
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationdu(l) 5 dl = U(0) = 1 and 1.) Substitute for U an unspecified trial function into governing equation, i.e. dx + = =
Consider an ODE of te form: Finite Element Metod du fu g d + wit te following Boundary Conditions: U(0) and du(l) 5 dl.) Substitute for U an unspecified trial function into governing equation, i.e. ^ U
More informationIntegral Calculus, dealing with areas and volumes, and approximate areas under and between curves.
Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationHow to solve the stochastic partial differential equation that gives a Matérn random field using the finite element method
arxiv:1803.03765v1 [stat.co] 10 Mar 2018 How to solve the stochastic partial differential equation that gives a Matérn random field using the finite element method Haakon Bakka bakka@r-inla.org March 13,
More informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationOptimal iterative solvers for linear nonsymmetric systems and nonlinear systems with PDE origins: Balanced black-box stopping tests
Optimal iterative solvers for linear nonsymmetric systems and nonlinear systems wit PDE origins: Balanced black-box stopping tests Pranjal, Prasad and Silvester, David J. 2018 MIMS EPrint: 2018.13 Mancester
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationN igerian Journal of M athematics and Applications V olume 23, (2014), 1 13
N igerian Journal of M atematics and Applications V olume 23, (24), 3 c N ig. J. M at. Appl. ttp : //www.kwsman.com CONSTRUCTION OF POLYNOMIAL BASIS AND ITS APPLICATION TO ORDINARY DIFFERENTIAL EQUATIONS
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationMultigrid Methods for Obstacle Problems
Multigrid Metods for Obstacle Problems by Cunxiao Wu A Researc Paper presented to te University of Waterloo in partial fulfillment of te requirement for te degree of Master of Matematics in Computational
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationConvergence and Descent Properties for a Class of Multilevel Optimization Algorithms
Convergence and Descent Properties for a Class of Multilevel Optimization Algoritms Stepen G. Nas April 28, 2010 Abstract I present a multilevel optimization approac (termed MG/Opt) for te solution of
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationComparison of multimesh hp-fem to interpolation and projection methods for spatial coupling of reactor thermal and neutron diffusion calculations
Comparison of multimes p-fem to interpolation and projection metods for spatial coupling of reactor termal and neutron diffusion calculations Lenka Dubcova a, Pavel Solin a,b,, Glen Hansen c, HyeongKae
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationFinancial Econometrics Prof. Massimo Guidolin
CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationINTRODUCTION TO CALCULUS LIMITS
Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and
More informationTwo new developments in Multilevel Monte Carlo methods
Two new developments in Multilevel Monte Carlo methods Mike Giles Abdul-Lateef Haji-Ali Oliver Sheridan-Methven Mathematical Institute, University of Oxford Rob Scheichl s Farewell Conference University
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationPreface. Here are a couple of warnings to my students who may be here to get a copy of what happened on a day that you missed.
Preface Here are my online notes for my course tat I teac ere at Lamar University. Despite te fact tat tese are my class notes, tey sould be accessible to anyone wanting to learn or needing a refreser
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationKey words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness
PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationTowards Perfectly Matched Layers for time-dependent space fractional PDEs
Towards Perfectly Matced Layers for time-dependent space fractional PDEs Xavier Antoine, Emmanuel Lorin To cite tis version: Xavier Antoine, Emmanuel Lorin. Towards Perfectly Matced Layers for time-dependent
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationLearning based super-resolution land cover mapping
earning based super-resolution land cover mapping Feng ing, Yiang Zang, Giles M. Foody IEEE Fellow, Xiaodong Xiuua Zang, Siming Fang, Wenbo Yun Du is work was supported in part by te National Basic Researc
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationA LOCAL FOURIER ANALYSIS FRAMEWORK FOR FINITE-ELEMENT DISCRETIZATIONS OF SYSTEMS OF PDES
A LOCAL FOURIER ANALYSIS FRAMEWORK FOR FINITE-ELEMENT DISCRETIZATIONS OF SYSTEMS OF PDES SCOTT P. MACLACHLAN AND CORNELIS W. OOSTERLEE Abstract. Since teir popularization in te late 1970s and early 1980s,
More informationarxiv: v1 [math.pr] 28 Dec 2018
Approximating Sepp s constants for te Slepian process Jack Noonan a, Anatoly Zigljavsky a, a Scool of Matematics, Cardiff University, Cardiff, CF4 4AG, UK arxiv:8.0v [mat.pr] 8 Dec 08 Abstract Slepian
More informationA SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types
More informationThird order Approximation on Icosahedral Great Circle Grids on the Sphere. J. Steppeler, P. Ripodas DWD Langen 2006
Tird order Approximation on Icosaedral Great Circle Grids on te Spere J. Steppeler, P. Ripodas DWD Langen 2006 Deasirable features of discretisation metods on te spere Great circle divisions of te spere:
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationLearning to Reject Sequential Importance Steps for Continuous-Time Bayesian Networks
Learning to Reject Sequential Importance Steps for Continuous-Time Bayesian Networks Jeremy C. Weiss University of Wisconsin-Madison Madison, WI, US jcweiss@cs.wisc.edu Sriraam Natarajan Indiana University
More informationChemnitz Scientific Computing Preprints
G. Of G. J. Rodin O. Steinbac M. Taus Coupling Metods for Interior Penalty Discontinuous Galerkin Finite Element Metods and Boundary Element Metods CSC/11-02 Cemnitz Scientific Computing Preprints Impressum:
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA
ROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA Y. EFENDIEV, J. GALVIS, S. KI KANG, AND R.D. LAZAROV Abstract. In tis paper, we study robust iterative solvers for
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationResearch Article Smoothing Analysis of Distributive Red-Black Jacobi Relaxation for Solving 2D Stokes Flow by Multigrid Method
Matematical Problems in Engineering Volume 205, Article ID 57298, 7 pages ttp://dx.doi.org/0.55/205/57298 Researc Article Smooting Analysis of Distributive Red-Black Jacobi Relaxation for Solving 2D Stokes
More informationA Study on Using Hierarchical Basis Error Estimates in Anisotropic Mesh Adaptation for the Finite Element Method
A Study on Using Hierarcical Basis Error Estimates in Anisotropic Mes Adaptation for te Finite Element Metod Lennard Kamenski arxiv:06.603v3 [mat.na] 20 Jan 202 Abstract: A common approac for generating
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationNumerical Analysis of the Double Porosity Consolidation Model
XXI Congreso de Ecuaciones Diferenciales y Aplicaciones XI Congreso de Matemática Aplicada Ciudad Real 21-25 septiembre 2009 (pp. 1 8) Numerical Analysis of te Double Porosity Consolidation Model N. Boal
More information= 0 and states ''hence there is a stationary point'' All aspects of the proof dx must be correct (c)
Paper 1: Pure Matematics 1 Mark Sceme 1(a) (i) (ii) d d y 3 1x 4x x M1 A1 d y dx 1.1b 1.1b 36x 48x A1ft 1.1b Substitutes x = into teir dx (3) 3 1 4 Sows d y 0 and states ''ence tere is a stationary point''
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationAPPLICATION OF A DIRAC DELTA DIS-INTEGRATION TECHNIQUE TO THE STATISTICS OF ORBITING OBJECTS
APPLICATION OF A DIRAC DELTA DIS-INTEGRATION TECHNIQUE TO THE STATISTICS OF ORBITING OBJECTS Dario Izzo Advanced Concepts Team, ESTEC, AG Noordwijk, Te Neterlands ABSTRACT In many problems related to te
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
More informationEffect of the Dependent Paths in Linear Hull
1 Effect of te Dependent Pats in Linear Hull Zenli Dai, Meiqin Wang, Yue Sun Scool of Matematics, Sandong University, Jinan, 250100, Cina Key Laboratory of Cryptologic Tecnology and Information Security,
More informationMore on generalized inverses of partitioned matrices with Banachiewicz-Schur forms
More on generalized inverses of partitioned matrices wit anaciewicz-scur forms Yongge Tian a,, Yosio Takane b a Cina Economics and Management cademy, Central University of Finance and Economics, eijing,
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationImproved Algorithms for Largest Cardinality 2-Interval Pattern Problem
Journal of Combinatorial Optimization manuscript No. (will be inserted by te editor) Improved Algoritms for Largest Cardinality 2-Interval Pattern Problem Erdong Cen, Linji Yang, Hao Yuan Department of
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
More informationFirst-Order Fractionally Integrated Non-Separable Spatial Autoregressive (FINSSAR(1,1)) Model and Some of its Properties
First-Order Fractionally Integrated Non-Separable Spatial Autoregressive (FINSSAR(,)) Model and Some of its Properties Alireza Godsi Department of Matematics, Faculty of Science University Putra, Malaysia.
More information(4.2) -Richardson Extrapolation
(.) -Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Suppose tat lim G 0 and lim F L. Te function F is said to converge to L as
More informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationMath 31A Discussion Notes Week 4 October 20 and October 22, 2015
Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationarxiv: v1 [physics.flu-dyn] 3 Jun 2015
A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows arxiv:156.132v1 [pysics.flu-dyn] 3 Jun 215 S. Dong Center for Computational & Applied Matematics Department
More informationStationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series
Stationary Gaussian Markov Processes As Limits of Stationary Autoregressive Time Series Lawrence D. Brown, Pilip A. Ernst, Larry Sepp, and Robert Wolpert August 27, 2015 Abstract We consider te class,
More informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationHaplotyping. Biostatistics 666
Haplotyping Biostatistics 666 Previously Introduction to te E-M algoritm Approac for likeliood optimization Examples related to gene counting Allele frequency estimation recessive disorder Allele frequency
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationMultigrid Methods for Discretized PDE Problems
Towards Metods for Discretized PDE Problems Institute for Applied Matematics University of Heidelberg Feb 1-5, 2010 Towards Outline A model problem Solution of very large linear systems Iterative Metods
More informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
More informationIntroduction to PDEs and Numerical Methods Tutorial 11. 2D elliptic equations
Platzalter für Bild, Bild auf Titelfolie inter das Logo einsetzen Introduction to PDEs and Numerical Metods Tutorial 11. 2D elliptic equations Dr. Noemi Friedman, 3. 1. 215. Overview Introduction (classification
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationSolving Poisson s equations by the Discrete Least Square meshless method
Boundary Elements and Oter Mes eduction Metods XXV 3 Solving Poisson s equations by te Discrete Least Square mesless metod H. Arzani & M. H. Afsar Said aaee University, Lavizan, eran, ran Department of
More information