WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid
|
|
- Kory Morrison
- 5 years ago
- Views:
Transcription
1 WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1
2 I. OVERVIEW 2
3 The importance of splitting integrators for (d/dt)x = f a (x) + f b (x) + keeps increasing: evolutionary PDEs, geometric integration,... Format: maps that advances one step given by composition ψ h = φ a λh φb µh, where h is the stepsize and φ a h, φb h,... are exact solution flows of (d/dt)x = f a (x), (d/dt)x = f b (x),... 3
4 Standard analysis of ψ h starts by seeing φ a h, φb h,... as (noncommuting) exponentials. BCH formula then used to write ψ h as a single exponential. The exponent provides the modified differential equation of the integrator, ie the (h-dependent) equation whose h flow formally coincides with ψ h. Difference between modified and true (d/dt)x = f a (x) + f b (x) + equations yields information on the properties of the integrator. 4
5 This is an indirect approach based on comparing differential equations rather than maps. It stands apart from usual techniques in numerical ODEs based on expansion of integrator map ψ h. For Runge-Kutta and related integrators the expansion of ψ h is best carried out by means of B-series (Hairer and Wanner 1974). B-series are parameterized by rooted trees. 5
6 Murua and SS (1999) developed a B-series technique to analyze splitting integrators (Hairer, Lubich, Wanner, GI Book, Chapter III.3). Recently Murua and SS suggested word series as an alternative to the B-series. Word series are parameterized by words from an alphabet. Use of word series based on simple algebraic systematic manipulations. 6
7 In numerical analysis, word series may be used to study the order of consistency of splitting integrators, find modified equations, etc. Outside numerical mathematics, word series may be used (Chartier, Murua, SS ) to perform averaging, to find normal forms, to compute integrals of motion, etc. In this talk we present an example of the use of word series to analyze splitting methods for SDEs. (Also: Ito SDEs, modified equations,... ) 7
8 II. DEFINING WORD SERIES 8
9 A MOTIVATING (DETERMINISTIC) EXAMPLE: dx = a A λ a (t) dtf a (x). A finite or infinite index set (alphabet). f a (x) vector field in D dimensions. λ a (t) scalar-valued function. 9
10 Solution with x(0) = x 0 given by (cf Chen series) x(t) = x 0 + J a1 a n (t)f a1 a n (x 0 ), n=1 a 1,...,a n A with J a1 a n (t) given by t 0 λ an (t n ) dt n tn 0 λ an 1 (t n 1 ) dt n 1 t2 0 λ a1 (t 1 ) dt 1, and, recursively, f a1 a n (x) = x f a2 a n (x) f a1 (x). Note separation: λ a, J a1 a n / f a, f a1 a n. 10
11 WORDS AND THEIR BASIS FUNCTIONS W denotes set of all words from the alphabet A (including the empty word ); W n is the set of words with n letters. Associate with w = a 1 a n W its word basis function f w (x) = f a1 a n (x) (with f (x) x). (Similar to elementary differentials.) 11
12 DEFINITION If δ C W (ie δ maps words into scalars), the word series with coefficients δ is the formal series: W δ (x) = δ w f w (x). w W Note that the notion of word series is relative to the f a. For each t, the solution value x(t) corresponds to the coefficients J w (t) built above from the λ a (t). 12
13 SOME PARTICULAR CASES: EXAMPLE I: A consists of a single letter a. For each n there is a single word w = a a with n letters. If furthermore λ a (t) 1, J w (t), w W n is found to be t n /n!. The word series W α(t) (x 0 ) is the Taylor expansion in powers of t of the solution x(t) of (d/dt)x = f a (x), x(0) = x 0. 13
14 EXAMPLE II: A consists of a two letters a, b. For each n there are 2 n words with n letters. If furthermore λ a (t) 1, λ b (t) 1, α w (t) = t n /n! for w W n. The word series W J(t) (x 0 ) is the Taylor expansion of x(t), (d/dt)x = f a (x) + f b (x), written in terms of f a, f b. 14
15 STOCHASTIC SDEs: dx = f a (x)dt + f A db(t) (Stratonovich), fits in this framework with λ a (t) 1 and λ A (t)dt = db(t). Solution formally given by W J(t) (x 0 ) where now J w (t) is a stochastic process. (Stochastic Taylor expansion see Kloeden and Platen, Chapter 5.) 15
16 III. OPERATING WITH WORD SERIES 16
17 THE CONVOLUTION PRODUCT If δ, δ C W, their convolution product δ δ C W is (δ δ ) a1 a n = δ δ a 1 a n + Not commutative. n 1 j=1 δ a1 a j δ a j+1 a n + δ a1 a n δ Associative. Unit: 1 C W with 1 = 1 and 1 w = 0 for w. 17
18 THE GROUP G (Group of characters of shuffle Hopf algebra) Let denote shuffle product of words. (Eg ab c = abc + acb + cab.) The set G of those γ C W that satisfy the so-called shuffle relations: γ = 1 and, for each w, w W, γ w γ w = N γ wj if w w = N w j. j=1 j=1 is a noncommutative group for. 18
19 For fixed t (and ω Ω), the coefficients J w (t), w W satisfy the shuffle relations, ie they give an element of G. For γ G, the series W γ (x) has special properties: (1) W γ (x) acts on the vector space of all word series by composition: W δ ( Wγ (x) ) = W γ δ (x), δ C W. 19
20 (2) If χ : R D R is any smooth observable: χ ( W γ (x) ) = γ w D w (χ)(x), w W where, for a A, D a = f a is the Lie operator associated with f a and for w = a 1... a n, D w is the differential operator D a1... D an. (3) Equivariance with respect to arbitrary changes of variables x = x(x). 20
21 IIII. SPLITTING INTEGRATORS FOR SDES 21
22 Consider integrators for Stratonovich equation dx = f a1 (x)dt + + f an (x)dt +f A1 (x)db A1 + + f AN (x)db AN such that the one-step map ψ h is composition of flows of individual pieces dx = f ai (x)dt or dx = f Aj (x)db Aj (or blocks, eg dx = f a1 (x)dt + f A2 (x)db Aj ). Each flow entering the composition is (for fixed t and ω Ω) a word series with known coefficients in G. 22
23 By property (1), ψ h = W γ(t) (x), where the coefficients γ w (t) are readily found by using the convolution formula. γ w (t) compared with coefficients J w (t) of true solution. For strong local error to be O(h p ), p = 3/2, 2, 5/2,... γ w (t) = J w (t) as for all words of weight p 1/2 (weight = number of lower case letters + number of upper case letters/2). (This assumes all relevant basis functions are 0.) 23
24 By property (3) E (χ ( W γ (t)(x) )) = w W E ( γ w (t) ) D w (χ)(x), For weak local error to be O(h p ), p = 2, 3,..., E ( γ w (t) ) = E ( J w (t) ), for all words of weight p 1. (This assumes all relevant basis functions are 0.) 24
25 V. LANGEVIN DYNAMICS 25
26 Consider dq = M 1 p dt dp = F (q) dt γp dt + σm 1/2 db, M diagonal with diagonal elements m i > 0, i = 1,..., d, γ > 0, B d-dimensional Wiener process. 26
27 Leimkuhler and Matthews (2013) consider several integrators based on the split systems: a: dq = M 1 p dt, dp = 0. b: dq = 0, dp = F (q) dt. o: dq = 0, dp = γp dt + σm 1/2 db. 27
28 They use acronym aboba for the Strang-like splitting ψ h = φ a h/2 φb h/2 φo h φb h/2 φa h/2 ; baoab defined similarly. Both aboba and baoab possess strong local error O(h 3/2 ); weak O(h 3 ). In spite of similarity, baoab turns out to be clearly superior to aboba. 28
29 We use word series based on the vector fields: a: (M 1 p, 0). b: (0, F (q)). c: (0, γp). A j : (0, σ m j e j ), j = 1,..., d. 29
30 Sparsity pattern results in many zero word basis functions. For weight < 3, f w 0 only for A j, A j a, A j c, A j ab, A j ca, A j cc and some purely deterministic words (that cause no trouble). The coefficients of the methods in the word series expansion are easily found to be: 30
31 w weight aboba baoab A j 1/2 A j a 3/2 (h/2)j Aj (h/2)j Aj A j c 3/2 A j ab 5/2 0 (h 2 /4)J Aj A j ca 5/2 (h/2)j Aj c (h/2)j Aj c A j cc 5/2 Only difference in A j ab, where baoab provides a better approximation: correlation truth-approximation is 5/3 74% for baoab, 0 for aboba. 31
32 Symmetrically, for word A j ba, baoab has γ Aj ba = 0, while aboba has γ Aj ba = (h 2 /4)J Aj. But then the basis function is zero and aboba does not benefit from it. The same phenomena appear for longer words. Word analysis reveals that the source of discrepancy is the following algorithmic flaw in aboba: in any given time step, it uses the same value of the force in both kicks of a given step. 32
Formal series and numerical integrators: some history and some new techniques
Formal series and numerical integrators: some history and some new techniques Jesús María Sanz-Serna and Ander Murua Abstract. This paper provides a brief history of B-series and the associated Butcher
More informationComputing ergodic limits for SDEs
Computing ergodic limits for SDEs M.V. Tretyakov School of Mathematical Sciences, University of Nottingham, UK Talk at the workshop Stochastic numerical algorithms, multiscale modelling and high-dimensional
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More informationHigher order weak approximations of stochastic differential equations with and without jumps
Higher order weak approximations of stochastic differential equations with and without jumps Hideyuki TANAKA Graduate School of Science and Engineering, Ritsumeikan University Rough Path Analysis and Related
More informationHigh order Cubature on Wiener space applied to derivative pricing and sensitivity estimations
1/25 High order applied to derivative pricing and sensitivity estimations University of Oxford 28 June 2010 2/25 Outline 1 2 3 4 3/25 Assume that a set of underlying processes satisfies the Stratonovich
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationBernardo D Auria Stochastic Processes /10. Notes. Abril 13 th, 2010
1 Stochastic Calculus Notes Abril 13 th, 1 As we have seen in previous lessons, the stochastic integral with respect to the Brownian motion shows a behavior different from the classical Riemann-Stieltjes
More informationThe Hopf algebraic structure of stochastic expansions and efficient simulation
The Hopf algebraic structure of stochastic expansions and efficient simulation Kurusch Ebrahimi Fard, Alexander Lundervold, Simon J.A. Malham, Hans Munthe Kaas and Anke Wiese York: 15th October 2012 KEF,
More informationThe Koopmanoperator action on a state space function a(x) is to replace it by
Appendix A9 Implementing evolution A9. Koopmania The Koopmanoperator action on a state space function a(x) is to replace it by its downstream value time t later, a(x) a(x(t)), evaluated at the trajectory
More informationAlgebraic Structure of Stochastic Expansions and Efficient Simulation
Algebraic Structure of Stochastic Expansions and Efficient Simulation Anke Wiese Heriot Watt University, Edinburgh, UK SSSC2012: ICMAT, Madrid November 5-9, 2012 Supported by EMS, IMA, LMS Work in collaboration
More informationStochastic Integration (Simple Version)
Stochastic Integration (Simple Version) Tuesday, March 17, 2015 2:03 PM Reading: Gardiner Secs. 4.1-4.3, 4.4.4 But there are boundary issues when (if ) so we can't apply the standard delta function integration
More informationIntroduction - Motivation. Many phenomena (physical, chemical, biological, etc.) are model by differential equations. f f(x + h) f(x) (x) = lim
Introduction - Motivation Many phenomena (physical, chemical, biological, etc.) are model by differential equations. Recall the definition of the derivative of f(x) f f(x + h) f(x) (x) = lim. h 0 h Its
More informationNumerical Integration of SDEs: A Short Tutorial
Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx
More informationTHEWAY in which time evolution acts on observables may be rephrased in
Appendix F Implementing evolution F.1 Koopmania THEWAY in which time evolution acts on observables may be rephrased in the language of functional analysis, by introducing the Koopman operator, whose action
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 3. Calculaus in Deterministic and Stochastic Environments Steve Yang Stevens Institute of Technology 01/31/2012 Outline 1 Modeling Random Behavior
More informationNumerical stroboscopic averaging for ODEs and DAEs
Numerical stroboscopic averaging for ODEs and DAEs M. P. Calvo Universidad de Valladolid, Spain Joint work with Ph. Chartier, A. Murua, J. M. Sanz-Serna 1 I. HIGHLY OSCILLATORY PROBLEMS 2 Consider the
More informationAccurate approximation of stochastic differential equations
Accurate approximation of stochastic differential equations Simon J.A. Malham and Anke Wiese (Heriot Watt University, Edinburgh) Birmingham: 6th February 29 Stochastic differential equations dy t = V (y
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationQuasi-stroboscopic averaging: The non-autonomous case
Outline Quasi-stroboscopic averaging: The non-autonomous case Philippe Chartier INRIA and ENS Cachan Bruz Joint work with Ander Murua and Jesus-Maria Sanz-Serna Geometric Numerical Integration, Oberwolfach
More informationStatistics 992 Continuous-time Markov Chains Spring 2004
Summary Continuous-time finite-state-space Markov chains are stochastic processes that are widely used to model the process of nucleotide substitution. This chapter aims to present much of the mathematics
More informationHairer /Gubinelli-Imkeller-Perkowski
Hairer /Gubinelli-Imkeller-Perkowski Φ 4 3 I The 3D dynamic Φ 4 -model driven by space-time white noise Let us study the following real-valued stochastic PDE on (0, ) T 3, where ξ is the space-time white
More informationFormal solution Chen-Fliess series
Formal solution Chen-Fliess series Ṡ = S(t) u a (t) a, S() = 1 Z = m a Z on algebra Â(Z ) of formal power series in aset Z of noncommuting indeterminates (letters) has the unique solution CF(T, u) = w
More informationDownloaded 12/15/15 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 40, No. 4, pp. 56 537 c 2002 Society for Industrial and Applied Mathematics PREDICTOR-CORRECTOR METHODS OF RUNGE KUTTA TYPE FOR STOCHASTIC DIFFERENTIAL EQUATIONS KEVIN BURRAGE
More informationStochastic Differential Equations
Chapter 5 Stochastic Differential Equations We would like to introduce stochastic ODE s without going first through the machinery of stochastic integrals. 5.1 Itô Integrals and Itô Differential Equations
More informationc 2002 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 4 No. pp. 507 5 c 00 Society for Industrial and Applied Mathematics WEAK SECOND ORDER CONDITIONS FOR STOCHASTIC RUNGE KUTTA METHODS A. TOCINO AND J. VIGO-AGUIAR Abstract. A general
More informationOn Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA
On Stochastic Adaptive Control & its Applications Bozenna Pasik-Duncan University of Kansas, USA ASEAS Workshop, AFOSR, 23-24 March, 2009 1. Motivation: Work in the 1970's 2. Adaptive Control of Continuous
More informationOn rough PDEs. Samy Tindel. University of Nancy. Rough Paths Analysis and Related Topics - Nagoya 2012
On rough PDEs Samy Tindel University of Nancy Rough Paths Analysis and Related Topics - Nagoya 2012 Joint work with: Aurélien Deya and Massimiliano Gubinelli Samy T. (Nancy) On rough PDEs Nagoya 2012 1
More informationQualitative behaviour of numerical methods for SDEs and application to homogenization
Qualitative behaviour of numerical methods for SDEs and application to homogenization K. C. Zygalakis Oxford Centre For Collaborative Applied Mathematics, University of Oxford. Center for Nonlinear Analysis,
More informationAlgebraic structures of B-series
Algebraic structures of B-series Philippe Chartier, Ernst Hairer, and Gilles Vilmart November 10, 2009 Abstract B-series are a fundamental tool in practical and theoretical aspects of numerical integrators
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationLie Symmetry of Ito Stochastic Differential Equation Driven by Poisson Process
American Review of Mathematics Statistics June 016, Vol. 4, No. 1, pp. 17-30 ISSN: 374-348 (Print), 374-356 (Online) Copyright The Author(s).All Rights Reserved. Published by American Research Institute
More informationThe Wiener Itô Chaos Expansion
1 The Wiener Itô Chaos Expansion The celebrated Wiener Itô chaos expansion is fundamental in stochastic analysis. In particular, it plays a crucial role in the Malliavin calculus as it is presented in
More informationCalculus II - Basic Matrix Operations
Calculus II - Basic Matrix Operations Ryan C Daileda Terminology A matrix is a rectangular array of numbers, for example 7,, 7 7 9, or / / /4 / / /4 / / /4 / /6 The numbers in any matrix are called its
More informationSecond order weak Runge Kutta type methods for Itô equations
Mathematics and Computers in Simulation 57 001 9 34 Second order weak Runge Kutta type methods for Itô equations Vigirdas Mackevičius, Jurgis Navikas Department of Mathematics and Informatics, Vilnius
More informationON LIE GROUP CLASSIFICATION OF A SCALAR STOCHASTIC DIFFERENTIAL EQUATION
Journal of Nonlinear Mathematical Physics ISSN: 1402-9251 (Print) 1776-0852 (Online) Journal homepage: http://www.tandfonline.com/loi/tnmp20 ON LIE GROUP CLASSIFICATION OF A SCALAR STOCHASTIC DIFFERENTIAL
More informationExponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationBernardo D Auria Stochastic Processes /12. Notes. March 29 th, 2012
1 Stochastic Calculus Notes March 9 th, 1 In 19, Bachelier proposed for the Paris stock exchange a model for the fluctuations affecting the price X(t) of an asset that was given by the Brownian motion.
More informationAn Overview of the Martingale Representation Theorem
An Overview of the Martingale Representation Theorem Nuno Azevedo CEMAPRE - ISEG - UTL nazevedo@iseg.utl.pt September 3, 21 Nuno Azevedo (CEMAPRE - ISEG - UTL) LXDS Seminar September 3, 21 1 / 25 Background
More information10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections )
c Dr. Igor Zelenko, Fall 2017 1 10. Linear Systems of ODEs, Matrix multiplication, superposition principle (parts of sections 7.2-7.4) 1. When each of the functions F 1, F 2,..., F n in right-hand side
More informationStochastic Integration and Stochastic Differential Equations: a gentle introduction
Stochastic Integration and Stochastic Differential Equations: a gentle introduction Oleg Makhnin New Mexico Tech Dept. of Mathematics October 26, 27 Intro: why Stochastic? Brownian Motion/ Wiener process
More informationSome Topics in Stochastic Partial Differential Equations
Some Topics in Stochastic Partial Differential Equations November 26, 2015 L Héritage de Kiyosi Itô en perspective Franco-Japonaise, Ambassade de France au Japon Plan of talk 1 Itô s SPDE 2 TDGL equation
More informationSection 7.4 Runge-Kutta Methods
Section 7.4 Runge-Kutta Methods Key terms: Taylor methods Taylor series Runge-Kutta; methods linear combinations of function values at intermediate points Alternatives to second order Taylor methods Fourth
More informationTheoretical Tutorial Session 2
1 / 36 Theoretical Tutorial Session 2 Xiaoming Song Department of Mathematics Drexel University July 27, 216 Outline 2 / 36 Itô s formula Martingale representation theorem Stochastic differential equations
More informationMonoidal Categories, Bialgebras, and Automata
Monoidal Categories, Bialgebras, and Automata James Worthington Mathematics Department Cornell University Binghamton University Geometry/Topology Seminar October 29, 2009 Background: Automata Finite automata
More informationarxiv: v1 [math.na] 5 Apr 2016
THE BCH-FORMULA AND ORDER CONDITIONS FOR SPLITTING METHODS arxiv:1604.01190v1 [math.na] 5 Apr 2016 WINFRIED AUZINGER, WOLFGANG HERFORT, OTHMAR KOCH, AND MECHTHILD THALHAMMER Abstract. As an application
More informationON BIALGEBRAS AND HOPF ALGEBRAS OF ORIENTED GRAPHS
Confluentes Mathematici, Vol. 4, No. 1 (2012) 1240003 (10 pages) c World Scientific Publishing Company DOI: 10.1142/S1793744212400038 ON BIALGEBRAS AND HOPF ALGEBRAS OF ORIENTED GRAPHS DOMINIQUE MANCHON
More informationIntroduction to the Numerical Solution of SDEs Selected topics
0 / 23 Introduction to the Numerical Solution of SDEs Selected topics Andreas Rößler Summer School on Numerical Methods for Stochastic Differential Equations at Vienna University of Technology, Austria
More informationManipulating exponential products Instead of working with complicated concatenations of flows
Manipulating exponential products Instead of working with complicated concatenations of flows z(t) = e t9 t2 t8 (f +f 1 +f 2 ) dt... e t1 (f +f 1 f 2 ) dt t1 e (f +f 1 +f 2 ) dt (p) it is desirable to
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationStrong Predictor-Corrector Euler Methods for Stochastic Differential Equations
Strong Predictor-Corrector Euler Methods for Stochastic Differential Equations Nicola Bruti-Liberati 1 and Eckhard Platen July 14, 8 Dedicated to the 7th Birthday of Ludwig Arnold. Abstract. This paper
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationJournal of Computational and Applied Mathematics. Higher-order semi-implicit Taylor schemes for Itô stochastic differential equations
Journal of Computational and Applied Mathematics 6 (0) 009 0 Contents lists available at SciVerse ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationNew Physical Principle for Monte-Carlo simulations
EJTP 6, No. 21 (2009) 9 20 Electronic Journal of Theoretical Physics New Physical Principle for Monte-Carlo simulations Michail Zak Jet Propulsion Laboratory California Institute of Technology, Advance
More informationImproved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3
Improved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3 J.H. Verner August 3, 2004 Abstract. In [5], Jackiewicz and Verner derived formulas for, and tested the implementation of two-step
More informationSDE Coefficients. March 4, 2008
SDE Coefficients March 4, 2008 The following is a summary of GARD sections 3.3 and 6., mainly as an overview of the two main approaches to creating a SDE model. Stochastic Differential Equations (SDE)
More informationAlexander Lundervold
Lie Butcher series and geometric numerical integration on manifolds PhD Thesis Alexander Lundervold Department of Mathematics University of Bergen June, 2011 Acknowledgements This dissertation is submitted
More informationVectors. January 13, 2013
Vectors January 13, 2013 The simplest tensors are scalars, which are the measurable quantities of a theory, left invariant by symmetry transformations. By far the most common non-scalars are the vectors,
More informationValidated Explicit and Implicit Runge-Kutta Methods
Validated Explicit and Implicit Runge-Kutta Methods Alexandre Chapoutot joint work with Julien Alexandre dit Sandretto and Olivier Mullier U2IS, ENSTA ParisTech 8th Small Workshop on Interval Methods,
More information1 Matrices and Systems of Linear Equations. a 1n a 2n
March 31, 2013 16-1 16. Systems of Linear Equations 1 Matrices and Systems of Linear Equations An m n matrix is an array A = (a ij ) of the form a 11 a 21 a m1 a 1n a 2n... a mn where each a ij is a real
More information1 Differentiable manifolds and smooth maps
1 Differentiable manifolds and smooth maps Last updated: April 14, 2011. 1.1 Examples and definitions Roughly, manifolds are sets where one can introduce coordinates. An n-dimensional manifold is a set
More informationBessel-like SPDEs. Lorenzo Zambotti, Sorbonne Université (joint work with Henri Elad-Altman) 15th May 2018, Luminy
Bessel-like SPDEs, Sorbonne Université (joint work with Henri Elad-Altman) Squared Bessel processes Let δ, y, and (B t ) t a BM. By Yamada-Watanabe s Theorem, there exists a unique (strong) solution (Y
More informationAxioms for Set Theory
Axioms for Set Theory The following is a subset of the Zermelo-Fraenkel axioms for set theory. In this setting, all objects are sets which are denoted by letters, e.g. x, y, X, Y. Equality is logical identity:
More informationStochastic Processes and Advanced Mathematical Finance
Steven R. Dunbar Department of Mathematics 23 Avery Hall University of Nebraska-Lincoln Lincoln, NE 68588-13 http://www.math.unl.edu Voice: 42-472-3731 Fax: 42-472-8466 Stochastic Processes and Advanced
More informationIntroduction to numerical simulations for Stochastic ODEs
Introduction to numerical simulations for Stochastic ODEs Xingye Kan Illinois Institute of Technology Department of Applied Mathematics Chicago, IL 60616 August 9, 2010 Outline 1 Preliminaries 2 Numerical
More informationB-Series and their Cousins - A. Contemporary Geometric View
B-Series and their Cousins - A Contemporary Geometric View Hans Munthe-Kaas 1 Department of Mathematics University of Bergen Norway url: hans.munthe-kaas.no mail: hans.munthe-kaas@uib.no SciCADE Bath,
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
More informationSimulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients
Journal of mathematics and computer Science 8 (2014) 28-32 Simulation Method for Solving Stochastic Differential Equations with Constant Diffusion Coefficients Behrouz Fathi Vajargah Department of statistics,
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationMatthew Zyskowski 1 Quanyan Zhu 2
Matthew Zyskowski 1 Quanyan Zhu 2 1 Decision Science, Credit Risk Office Barclaycard US 2 Department of Electrical Engineering Princeton University Outline Outline I Modern control theory with game-theoretic
More informationAn algebraic approach to invariant preserving integators: The case of quadratic and Hamiltonian invariants
Numerische Mathematik manuscript No. (will be inserted by the editor) An algebraic approach to invariant preserving integators: The case of quadratic and Hamiltonian invariants Philippe Chartier 1, Erwan
More informationStochastic Differential Equations.
Chapter 3 Stochastic Differential Equations. 3.1 Existence and Uniqueness. One of the ways of constructing a Diffusion process is to solve the stochastic differential equation dx(t) = σ(t, x(t)) dβ(t)
More informationEquivalence groupoids of classes of linear ordinary differential equations and their group classification
Journal of Physics: Conference Series PAPER OPEN ACCESS Equivalence groupoids of classes of linear ordinary differential equations and their group classification To cite this article: Vyacheslav M Boyko
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationRough Burgers-like equations with multiplicative noise
Rough Burgers-like equations with multiplicative noise Martin Hairer Hendrik Weber Mathematics Institute University of Warwick Bielefeld, 3.11.21 Burgers-like equation Aim: Existence/Uniqueness for du
More informationLocal Parametrization and Puiseux Expansion
Chapter 9 Local Parametrization and Puiseux Expansion Let us first give an example of what we want to do in this section. Example 9.0.1. Consider the plane algebraic curve C A (C) defined by the equation
More informationTHE GEOMETRY OF ITERATED STRATONOVICH INTEGRALS
THE GEOMETRY OF ITERATED STRATONOVICH INTEGRALS CHRISTIAN BAYER Abstract. We give a summary on the geometry of iterated Stratonovich integrals. For this exposition, we always have the connection to stochastic
More informationNumerical methods for solving stochastic differential equations
Mathematical Communications 4(1999), 251-256 251 Numerical methods for solving stochastic differential equations Rózsa Horváth Bokor Abstract. This paper provides an introduction to stochastic calculus
More informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
More informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationA Unified Formulation of Gaussian Versus Sparse Stochastic Processes
A Unified Formulation of Gaussian Versus Sparse Stochastic Processes Michael Unser, Pouya Tafti and Qiyu Sun EPFL, UCF Appears in IEEE Trans. on Information Theory, 2014 Presented by Liming Wang M. Unser,
More information1 Matrices and matrix algebra
1 Matrices and matrix algebra 1.1 Examples of matrices A matrix is a rectangular array of numbers and/or variables. For instance 4 2 0 3 1 A = 5 1.2 0.7 x 3 π 3 4 6 27 is a matrix with 3 rows and 5 columns
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationCALCULUS ON MANIFOLDS
CALCULUS ON MANIFOLDS 1. Introduction. Crash course on the Multivariate Calculus 1.1. Linear algebra. Field of real numbers R and its completeness. Real line R 1 and real spaces R n. Linear functionals,
More informationIntroduction to Tree Algebras in Quantum Field Theory
Introduction to Tree Algebras in Quantum Field Theory M. D. Sheppeard Abstract Renormalization procedures for quantum field theories are today reinterpreted using Hopf algebras and related function algebras.
More informationReminder on basic differential geometry
Reminder on basic differential geometry for the mastermath course of 2013 Charts Manifolds will be denoted by M, N etc. One should think of a manifold as made out of points (while the elements of a vector
More informationA stroboscopic averaging technique for highly-oscillatory Schrödinger equation
A stroboscopic averaging technique for highly-oscillatory Schrödinger equation F. Castella IRMAR and INRIA-Rennes Joint work with P. Chartier, F. Méhats and A. Murua Saint-Malo Workshop, January 26-28,
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationIntroduction to Computational Stochastic Differential Equations
Introduction to Computational Stochastic Differential Equations Gabriel J. Lord Catherine E. Powell Tony Shardlow Preface Techniques for solving many of the differential equations traditionally used by
More informationNorthwestern University Department of Electrical Engineering and Computer Science
Northwestern University Department of Electrical Engineering and Computer Science EECS 454: Modeling and Analysis of Communication Networks Spring 2008 Probability Review As discussed in Lecture 1, probability
More information