Parallel Multilevel Monte Carlo Simulation
|
|
- Randall Wood
- 5 years ago
- Views:
Transcription
1 Parallel Simulation Mathematisches Institut Goethe-Universität Frankfurt am Main Advances in Financial Mathematics Paris January 7-10, 2014 Simulation
2 Outline 1 Monte Carlo Algorithm Numerical Results Simulation
3 Outline 1 Monte Carlo Algorithm Numerical Results Simulation
4 Option pricing Model: Black-Scholes ds(t) = µs(t)dt + σs(t)dw (t) Euler-Maruyama discretization: Ŝ(t j+1 ) = Ŝ(t j) + µŝ(t j)h + σŝ(t j)z j with z j N(0, h) Martingale approach: V (S, 0) = e rt E [V (S, T )] Monte Carlo simulation: ˆV (S, 0) = e rt 1 N N ˆV ({Ŝ (i) (t 1 ),..., Ŝ (i) (t d )}, T ) i=1 where ˆV is the discretized payoff, e.g. for a lookback option ˆV ({Ŝ(t 1 ),... Ŝ(t d )}, T ) = Ŝ(t d ) min Ŝ(t j ) 1 j d Simulation
5 Discretization Error Monte Carlo Shortcut notation: Mean square error: Y = E[f (S)] and Ŷ = 1 N N f (Ŝ (i) ) i=1 MSE = E[ E[f (S)] Ê[f (Ŝ)] 2 ] = (E[Ŷ Y ]) 2 + Var[Ŷ ] = O(h 2α + N 1 ) if the approximation f (Ŝ) converges with weak order α max E[f (Ŝ(t j))] E[f (S(t j ))] ch α. 1 j d Simulation
6 Error balancing: To achieve an RMSE of O(ɛ) it is necessary to select The corresponding cost C is then N = O(ɛ 2 ) and h = O(ɛ 1/α ) C = N d = O(ɛ 2 1/α ). Solving for ɛ yields the optimal rate of convergence Optimal refinement rule: RMSE = O(ɛ) = O(C α 2α+1 ). d 2d N 4αN Simulation
7 10 1 Asian Option RMSE Bias Std. Deviation C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, d 1 = 1, N 1 = 5, 100 repetitions all rates 1/3 Simulation
8 10 1 RMSE of various options European option Asian option Barrier option RMSE C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, d 1 = 1, N 1 = 5, 100 repetitions Simulation
9 Convergence rates European option: α = 1 Asian option: α = 1 Barrier option: α = 1/2 α RMSE rate α/(2α+1) 1 1/3 1/2 1/4 Improvement: discrete minimum correction for Barrier options: Ŝ min = min Ŝ(t j ) k σ T /d j=1,...,d with k = recovers α = 1 (Kou, 2003) Simulation
10 10 1 RMSE of various options European option Asian option Barrier option Barrier option transf. min. RMSE C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, d 1 = 1, N 1 = 5, 100 repetitions Simulation
11 Outline 1 Monte Carlo Algorithm Numerical Results Simulation
12 Multilevel Monte Carlo (Giles, 2008): simulate asset prices for different mesh widths h 1,..., h L in time small mesh width = low discretization error, but large costs large mesh width = high discretization error, but small costs Rewrite the payoff on the finest level L as a telescope sum E[ˆP L ] = E[ˆP 0 ] + L E[ˆP l ˆP l 1 ]. l=1 where ˆP l is the approximation for mesh width h l = M l T. Simulation
13 Multilevel Monte Carlo Computation: From the estimate of the expectation E[ˆP l ˆP l 1 ]: Ŷ l = 1 N l N l i=1 (ˆP (i) l and the estimate of E[ˆP 0 ]: ˆP (i) l 1 ), for l = 1,..., L, Ŷ 0 = 1 N 0 N 0 i=1 ˆP (i) 0 results the multilevel estimate of E[ˆP L ]: L Ŷ = Ŷ l, l=0 with N l being the number of simulations for level l. Simulation
14 Multilevel Conditions weak convergence strong convergence ] E [ˆP l P c 1 hl α ] V [Ŷl c 2 N 1 l h β l which corresponds to the strong order of the SDE discretization ( max E f (Ŝ(t j)) f (S(t j )) p) 1 p c 3 h β/2 1 j d Multilevel Complexity (Giles, 2008): { O(C 1/2 ) for β 1, RMSE O(C α 2α+1 β ) for β < 1. Simulation
15 Monte Carlo 10 1 Alpha of various options weak error Asian option Lookback option Barrier option Digital option d Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
16 Monte Carlo 10 0 Beta of various options Variance Asian option Lookback option Barrier option Digital option d Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
17 Monte Carlo 10 1 RMSE of various options Asian option Lookback option Barrier option Digital option RMSE C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
18 Summary Option expected computed α β RMSE α β RMSE Asian 1 1 1/ Lookback 1/2 1 1/ Barrier 1/2 1/2 1/ Digital 1 1/2 2/ RMSE = α/(2α + 1 β) if β < 1 and RMSE=1/2 if β 1 Simulation
19 Outline 1 Monte Carlo Algorithm Numerical Results Simulation
20 Dimension- (Dimension-) Algorithm: 1 set N 0, N 1 = determine V l for l = 0,..., L such that V := L l=0 V l estimates the variance and B := (ŶL/(M α 1)) 2 estimates the bias 2 3 if V + B < ɛ 2 stop 4 else if V > B determine V l, l = 0,..., L, which has the largest variance/work and double N l if B > V then L L + 1 and set N L+1 = goto step 2 Simulation
21 Dimension- (Dimension-) Algorithm: 1 set N 0, N 1 = determine V l for l = 0,..., L such that V := L l=0 V l estimates the variance and B := (ŶL/(M α 1)) 2 estimates the bias 2 3 if V + B < ɛ 2 stop 4 else if V > B determine V l, l = 0,..., L, which has the largest variance/work and double N l if B > V set N L+1 = goto step Dimension adaptive MLMC 1600 Dimension adaptive MLMC N(L) N(L) Level L Level L Simulation
22 Dimension- (Dimension-) Algorithm: 1 set N 0, N 1 = determine V l for l = 0,..., L such that V := L l=0 V l estimates the variance and B := (ŶL/(M α 1)) 2 estimates the bias 2 3 if V + B < ɛ 2 stop 4 else if V > B determine V l, l = 0,..., L, which has the largest variance/work and double N l if B > V set N L+1 = goto step Dimension adaptive MLMC 1600 Dimension adaptive MLMC N(L) N(L) Level L Level L Simulation
23 Dimension- (Dimension-) Algorithm: 1 set N 0, N 1 = determine V l for l = 0,..., L such that V := L l=0 V l estimates the variance and B := (ŶL/(M α 1)) 2 estimates the bias 2 3 if V + B < ɛ 2 stop 4 else if V > B determine V l, l = 0,..., L, which has the largest variance/work and double N l if B > V set N L+1 = goto step Dimension adaptive MLMC 6400 Dimension adaptive MLMC N(L) 200 N(L) Level L Level L Simulation
24 Algorithm 1 set N 0, N 1 = determine V l for l = 0,..., L such that V := L l=0 V l estimates the variance and B := (ŶL/(M α 1)) 2 estimates the bias 2 3 if V + B < ɛ 2 stop 4 else if V > B determine V l, l = 0,..., L, which has the largest variance/work and double N l if B > V set N L+1 = goto step 2 Simulation
25 Monte Carlo 10 2 European Option RMSE 10 3 standard method standard error bound adaptive method adaptive error bound C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
26 Monte Carlo 10 1 Asian Option RMSE Bias Std. Deviation C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions all rates 1/2 Simulation
27 Monte Carlo RMSE RMSE of various options Asian option Lookback option Barrier option Digital option RMSE Option MLMC adapt. MLMC Asian Lookback Barrier Digital C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
28 MLMC not optimal for barrier options one idea: rewrite payoff as a product of probabilities with N 1 P = e r (S(t N ) K) + i=0 ( 2(Sn B) + (S n+1 B) + ) p i = 1 exp σ 2 SnT/N 2 also possible for Double Barrier Options but getting more complicated with more conditions in the payoff p i Simulation
29 more general idea: use adaptive path discretization close to the barrier adaptively refine time intervals if the barrier-crossing probability is large Ψ(S i 1, S i ) := P ( S i 1/2 < B ) > w no additional complexity for more complicated options Simulation
30 How to construct the mid points after refinement? version 1: Brownian bridge using the mean of forward and backward Euler-Maruyama estimates S i 1/2 = 1 ( Si 1 + S i 1 r(t i 1/2 t i 1 ) + S i 1 σ(w i 1/2 W i 1 ) ) ( ) S i r(t i t i 1/2 ) + σ(w i W i 1/2 ) version 2: Brownian bridge using an arithmetic Brownian motion S i 1/2 = S i 1 + r(t i 1/2 t i 1 ) + σ(w i 1/2 W i 1 ) Simulation
31 Example Monte Carlo 1.4 Down and Out Call t S(t) W(t) Barrier Simulation
32 Example Monte Carlo 1.4 Down and Out Call S(t) W(t) Barrier t Simulation
33 Example Monte Carlo 1.4 Down and Out Call S(t) W(t) Barrier t Simulation
34 Example Monte Carlo paths for a Barrier option Refinement steps with w= 1.28 and d= paths for a Barrier option Refinement steps with w= 1.28 and d= Number log(number) Steps Steps Simulation
35 Time-adaptive Again rewrite the payoff on the finest level L as a telescope sum E[ˆP w L L ] = E[ˆP w 0 0 ] + L l=1 E[ˆP w l l ˆP w l 1 l 1 ]. where ˆP w l l is the approximation for mesh width h l = M l T and adaptive path discretization parameter w l. Calculate the expectations on the right side with the corresponding Monte Carlo estimator. Simulation
36 Example Monte Carlo 10 1 Alpha of various methods for Barrier Options 10 2 weak error transformed Barrier adapt. refining Barrier Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions d Simulation
37 Example Monte Carlo 10 1 Beta of various methods for Barrier Options transformed Barrier adaptive refining Barrier 10 2 Variance d Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
38 Example Monte Carlo 10 1 RMSE of various methods for Barrier Options transformed Barrier adaptive refining Barrier RMSE 10 2 Method RMSE Barrier 0.13 transf. Barrier 0.29 adapt. refining C Parameters: S 0 = 1, r = 0.05, σ = 0.2, K = 1, T = 1, N 0 = 100, 100 repetitions Simulation
39 Algorithm Numerical Results Outline 1 Monte Carlo Algorithm Numerical Results Simulation
40 Algorithm Numerical Results Parallel programming Properties Computer with more than one processor/kernel needed Linear speed-up in the number of processors possible Efficient for loops that do the same recurring tasks Very useful for Monte Carlo simulation Example: Parallel Monte Carlo with M processors N/M ˆV j (S, 0) = e rt 1 ˆV ({Ŝ (i,j) (t 1 ),..., Ŝ (i,j) (t d )}, T ) N/M i=1 and the aggretated estimator for the option price is ˆV (S, 0) = M j=1 ˆV j (S, 0) Simulation
41 Algorithm Numerical Results Idea Parallelize each sum of the MLMC estimator Ŷ = L l=0 Ŷ l with Ŷ l = 1 N l N l i=1 (ˆP (i) l ˆP (i) l 1 ). Do not overcharge the memory - do not save the whole path of S but only the necessary part for the option value - do not save all option values but calculate the expectation and variance recursively Simulation
42 Algorithm Monte Carlo Algorithm Numerical Results 1 set L := 0, N M0 = 1000/M 2 determine the variances V k,l for k ( = 1,..., M and l = 0,..., L M ) on each processor such that V l = k=1 V k,l /M and V := L l=0 V l 3 define optimal N l, l = 0,..., L as in the standard algorithm and if N l has increased calculate N Ml = ( Nl new Nl old ) /M extra samples on each processor 4 stop if RMSE< ɛ and L 2 5 else set L := L + 1, N ML = 1000/M and go to step 2. Simulation
43 European Option Monte Carlo Algorithm Numerical Results threads 10 threads 1 thread RMSE Time Figure: Convergence rates in time for a european option using 1,10 and 50 threads. Simulation
44 Algorithm Numerical Results European Option ɛ kernels kernels Table: Time factor improvements for a european option with 10 and 50 kernerls compared to 1 kernel. linear speed-up in the number of kernels if the program is running for more than 1 second code written in C++ using MPI Simulation
45 Algorithm Numerical Results Conclusions Summary: Monte Carlo: RMSE rate 1/3 : RMSE rate 1/2 in best case : RMSE rate 1/2 also for barrier options : linear speed-up Extensions: MLMC for Milstein and higher order schemes (Giles, 2007) Multilevel Quasi-Monte Carlo (Giles, Waterhouse, 2009; G., Noll 2012) Multilevel (adaptive) sparse grid integration (G., Heinz, 2012) Simulation
Multilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
More informationMultilevel Monte Carlo for Basket Options
MLMC for basket options p. 1/26 Multilevel Monte Carlo for Basket Options Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance WSC09,
More informationMonte Carlo Methods. Prof. Mike Giles. Oxford University Mathematical Institute. Lecture 1 p. 1.
Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Lecture 1 p. 1 Geometric Brownian Motion In the case of Geometric Brownian Motion ds t = rs t dt+σs
More informationModule 4: Monte Carlo path simulation
Module 4: Monte Carlo path simulation Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Module 4: Monte Carlo p. 1 SDE Path Simulation In Module 2, looked at the case
More informationMultilevel Monte Carlo Simulation
Multilevel Monte Carlo p. 1/48 Multilevel Monte Carlo Simulation Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance Workshop on Computational
More informationComputing Greeks with Multilevel Monte Carlo Methods using Importance Sampling
Computing Greeks with Multilevel Monte Carlo Methods using Importance Sampling Supervisor - Dr Lukas Szpruch Candidate Number - 605148 Dissertation for MSc Mathematical & Computational Finance Trinity
More informationMultilevel path simulation for jump-diffusion SDEs
Multilevel path simulation for jump-diffusion SDEs Yuan Xia, Michael B. Giles Abstract We investigate the extension of the multilevel Monte Carlo path simulation method to jump-diffusion SDEs. We consider
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationMultilevel Monte Carlo for VaR
Multilevel Monte Carlo for VaR Mike Giles, Wenhui Gou, Abdul-Lateef Haji-Ali Mathematical Institute, University of Oxford (BNP Paribas, Hong Kong) (also discussions with Ralf Korn, Klaus Ritter) Advances
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More informationMultilevel Monte Carlo Methods for American Options
Multilevel Monte Carlo Methods for American Options Simon Gemmrich, PhD Kellog College University of Oxford A thesis submitted in partial fulfillment of the MSc in Mathematical Finance November 19, 2012
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationComputational Finance
Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk
More informationMultilevel Change of Measure for Complex Digital Options
Multilevel Change of Measure for Complex Digital Options Jiaxing Wang Somerville College University of Oxford A thesis submitted in partial fulfillment of the MSc in Mathematical Finance Trinity 2014 This
More informationAnalysing multi-level Monte Carlo for options with non-globally Lipschitz payoff
Finance Stoch 2009 13: 403 413 DOI 10.1007/s00780-009-0092-1 Analysing multi-level Monte Carlo for options with non-globally Lipschitz payoff Michael B. Giles Desmond J. Higham Xuerong Mao Received: 1
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationFinancial Mathematics and Supercomputing
GPU acceleration in early-exercise option valuation Álvaro Leitao and Cornelis W. Oosterlee Financial Mathematics and Supercomputing A Coruña - September 26, 2018 Á. Leitao & Kees Oosterlee SGBM on GPU
More informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationVariance Reduction Through Multilevel Monte Carlo Path Calculations
Variance Reduction Through Mutieve Monte Caro Path Cacuations Mike Gies gies@comab.ox.ac.uk Oxford University Computing Laboratory Mutieve Monte Caro p. 1/30 Mutigrid A powerfu technique for soving PDE
More informationComputational Finance. Computational Finance p. 1
Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy
More informationMath 623 (IOE 623), Winter 2008: Final exam
Math 623 (IOE 623), Winter 2008: Final exam Name: Student ID: This is a closed book exam. You may bring up to ten one sided A4 pages of notes to the exam. You may also use a calculator but not its memory
More informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
More informationMonte Carlo Path Simulation and the Multilevel Monte Carlo Method. Krister Janzon
Monte Carlo Path Simulation and the Multilevel Monte Carlo Method Krister Janzon Master s Thesis, 30 ECTS Master of Science Programme in Engineering Physics, 300 ECTS Spring 2018 Monte Carlo Path Simulation
More informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More information"Vibrato" Monte Carlo evaluation of Greeks
"Vibrato" Monte Carlo evaluation of Greeks (Smoking Adjoints: part 3) Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance MCQMC 2008,
More informationStochastic Differential Equations in Finance and Monte Carlo Simulations
Stochastic Differential Equations in Finance and Department of Statistics and Modelling Science University of Strathclyde Glasgow, G1 1XH China 2009 Outline Stochastic Modelling in Asset Prices 1 Stochastic
More informationMonte-Carlo Methods in Financial Engineering
Monte-Carlo Methods in Financial Engineering Universität zu Köln May 12, 2017 Outline Table of Contents 1 Introduction 2 Repetition Definitions Least-Squares Method 3 Derivation Mathematical Derivation
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationStochastic Grid Bundling Method
Stochastic Grid Bundling Method GPU Acceleration Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee London - December 17, 2015 A. Leitao &
More informationStochastic Modelling in Finance
in Finance Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH April 2010 Outline and Probability 1 and Probability 2 Linear modelling Nonlinear modelling 3 The Black Scholes
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationPricing Early-exercise options
Pricing Early-exercise options GPU Acceleration of SGBM method Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee Lausanne - December 4, 2016
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationValuation of performance-dependent options in a Black- Scholes framework
Valuation of performance-dependent options in a Black- Scholes framework Thomas Gerstner, Markus Holtz Institut für Numerische Simulation, Universität Bonn, Germany Ralf Korn Fachbereich Mathematik, TU
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationFinancial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA
Financial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA Rajesh Bordawekar and Daniel Beece IBM T. J. Watson Research Center 3/17/2015 2014 IBM Corporation
More informationValuation of Asian Option. Qi An Jingjing Guo
Valuation of Asian Option Qi An Jingjing Guo CONTENT Asian option Pricing Monte Carlo simulation Conclusion ASIAN OPTION Definition of Asian option always emphasizes the gist that the payoff depends on
More informationMath Computational Finance Option pricing using Brownian bridge and Stratified samlping
. Math 623 - Computational Finance Option pricing using Brownian bridge and Stratified samlping Pratik Mehta pbmehta@eden.rutgers.edu Masters of Science in Mathematical Finance Department of Mathematics,
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationMulti-scale methods for stochastic differential equations
Multi-scale methods for stochastic differential equations by Niklas Zettervall Department of Physics Umeå University February 2012 Abstract Standard Monte Carlo methods are used extensively to solve stochastic
More informationGamma. The finite-difference formula for gamma is
Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas
More informationBarrier Option. 2 of 33 3/13/2014
FPGA-based Reconfigurable Computing for Pricing Multi-Asset Barrier Options RAHUL SRIDHARAN, GEORGE COOKE, KENNETH HILL, HERMAN LAM, ALAN GEORGE, SAAHPC '12, PROCEEDINGS OF THE 2012 SYMPOSIUM ON APPLICATION
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Spring 2010 Computer Exercise 2 Simulation This lab deals with
More informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
More informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationEstimating the Greeks
IEOR E4703: Monte-Carlo Simulation Columbia University Estimating the Greeks c 207 by Martin Haugh In these lecture notes we discuss the use of Monte-Carlo simulation for the estimation of sensitivities
More informationFrom CFD to computational finance (and back again?)
From CFD to computational finance (and back again?) Mike Giles University of Oxford Mathematical Institute MIT Center for Computational Engineering Seminar March 14th, 2013 Mike Giles (Oxford) CFD to finance
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationAsian Option Pricing: Monte Carlo Control Variate. A discrete arithmetic Asian call option has the payoff. S T i N N + 1
Asian Option Pricing: Monte Carlo Control Variate A discrete arithmetic Asian call option has the payoff ( 1 N N + 1 i=0 S T i N K ) + A discrete geometric Asian call option has the payoff [ N i=0 S T
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
More informationModern Methods of Option Pricing
Modern Methods of Option Pricing Denis Belomestny Weierstraß Institute Berlin Motzen, 14 June 2007 Denis Belomestny (WIAS) Modern Methods of Option Pricing Motzen, 14 June 2007 1 / 30 Overview 1 Introduction
More informationEvaluating the Longstaff-Schwartz method for pricing of American options
U.U.D.M. Project Report 2015:13 Evaluating the Longstaff-Schwartz method for pricing of American options William Gustafsson Examensarbete i matematik, 15 hp Handledare: Josef Höök, Institutionen för informationsteknologi
More informationMONTE CARLO EXTENSIONS
MONTE CARLO EXTENSIONS School of Mathematics 2013 OUTLINE 1 REVIEW OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO 3 SUMMARY MONTE CARLO SO FAR... Simple to program
More informationMonte Carlo Based Numerical Pricing of Multiple Strike-Reset Options
Monte Carlo Based Numerical Pricing of Multiple Strike-Reset Options Stavros Christodoulou Linacre College University of Oxford MSc Thesis Trinity 2011 Contents List of figures ii Introduction 2 1 Strike
More informationSparse Grid Quadrature Methods. for. Computational Finance
Sparse Grid Quadrature Methods for Computational Finance Habilitationsschrift an der Mathematisch Naturwissenschaftlichen Fakultät der Rheinischen Friedrich Wilhelms Universität Bonn eingereicht von Thomas
More informationA new PDE approach for pricing arithmetic average Asian options
A new PDE approach for pricing arithmetic average Asian options Jan Večeř Department of Mathematical Sciences, Carnegie Mellon University, Pittsburgh, PA 15213. Email: vecer@andrew.cmu.edu. May 15, 21
More informationMulti-Asset Options. A Numerical Study VILHELM NIKLASSON FRIDA TIVEDAL. Master s thesis in Engineering Mathematics and Computational Science
Multi-Asset Options A Numerical Study Master s thesis in Engineering Mathematics and Computational Science VILHELM NIKLASSON FRIDA TIVEDAL Department of Mathematical Sciences Chalmers University of Technology
More informationToward a coherent Monte Carlo simulation of CVA
Toward a coherent Monte Carlo simulation of CVA Lokman Abbas-Turki (Joint work with A. I. Bouselmi & M. A. Mikou) TU Berlin January 9, 2013 Lokman (TU Berlin) Advances in Mathematical Finance 1 / 16 Plan
More informationAD in Monte Carlo for finance
AD in Monte Carlo for finance Mike Giles giles@comlab.ox.ac.uk Oxford University Computing Laboratory AD & Monte Carlo p. 1/30 Overview overview of computational finance stochastic o.d.e. s Monte Carlo
More informationAccelerated Option Pricing Multiple Scenarios
Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationConditional sampling for barrier option pricing under the Heston model
Conditional sampling for barrier option pricing under the Heston model Nico Achtsis, Ronald Cools, and Dir Nuyens Abstract We propose a quasi-monte Carlo algorithm for pricing noc-out and noc-in barrier
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
More informationFast Convergence of Regress-later Series Estimators
Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser
More informationMAFS Computational Methods for Pricing Structured Products
MAFS550 - Computational Methods for Pricing Structured Products Solution to Homework Two Course instructor: Prof YK Kwok 1 Expand f(x 0 ) and f(x 0 x) at x 0 into Taylor series, where f(x 0 ) = f(x 0 )
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationS (17) DOI: Reference: CAM Journal of Computational and Applied Mathematics
Accepted Manuscript A multi-level dimension reduction Monte-Carlo method for jump-diffusion models Duy-Minh Dang PII: S377-47173169-3 DOI: http://dx.doi.org/1.116/j.cam.17.4.14 Reference: CAM 1193 To appear
More informationMultilevel Monte Carlo Path Simulation
Mutieve Monte Caro Path Simuation Mike Gies gies@comab.ox.ac.uk Oxford University Computing Laboratory First IMA Conference on Computationa Finance Mutieve Monte Caro p. 1/34 Generic Probem Stochastic
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationF1 Acceleration for Montecarlo: financial algorithms on FPGA
F1 Acceleration for Montecarlo: financial algorithms on FPGA Presented By Liang Ma, Luciano Lavagno Dec 10 th 2018 Contents Financial problems and mathematical models High level synthesis Optimization
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationCh 5. Several Numerical Methods
Ch 5 Several Numerical Methods I Monte Carlo Simulation for Multiple Variables II Confidence Interval and Variance Reduction III Solving Systems of Linear Equations IV Finite Difference Method ( 有限差分法
More informationExact Sampling of Jump-Diffusion Processes
1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance
More informationMONTE CARLO METHODS FOR AMERICAN OPTIONS. Russel E. Caflisch Suneal Chaudhary
Proceedings of the 2004 Winter Simulation Conference R. G. Ingalls, M. D. Rossetti, J. S. Smith, and B. A. Peters, eds. MONTE CARLO METHODS FOR AMERICAN OPTIONS Russel E. Caflisch Suneal Chaudhary Mathematics
More informationarxiv: v1 [cs.dc] 14 Jan 2013
A parallel implementation of a derivative pricing model incorporating SABR calibration and probability lookup tables Qasim Nasar-Ullah 1 University College London, Gower Street, London, United Kingdom
More informationImproved Lower and Upper Bound Algorithms for Pricing American Options by Simulation
Improved Lower and Upper Bound Algorithms for Pricing American Options by Simulation Mark Broadie and Menghui Cao December 2007 Abstract This paper introduces new variance reduction techniques and computational
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationGRAPHICAL ASIAN OPTIONS
GRAPHICAL ASIAN OPTIONS MARK S. JOSHI Abstract. We discuss the problem of pricing Asian options in Black Scholes model using CUDA on a graphics processing unit. We survey some of the issues with GPU programming
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
More informationMultilevel Monte Carlo methods for finance
Multilevel Monte Carlo methods for finance Mike Giles Mathematical Institute, University of Oxford Oxford-Man Institute of Quantitative Finance HPCFinance Final Conference March 14, 2016 Mike Giles (Oxford)
More informationThe Binomial Lattice Model for Stocks: Introduction to Option Pricing
1/33 The Binomial Lattice Model for Stocks: Introduction to Option Pricing Professor Karl Sigman Columbia University Dept. IEOR New York City USA 2/33 Outline The Binomial Lattice Model (BLM) as a Model
More informationOptimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics
More informationAN IMPROVED BINOMIAL METHOD FOR PRICING ASIAN OPTIONS
Commun. Korean Math. Soc. 28 (2013), No. 2, pp. 397 406 http://dx.doi.org/10.4134/ckms.2013.28.2.397 AN IMPROVED BINOMIAL METHOD FOR PRICING ASIAN OPTIONS Kyoung-Sook Moon and Hongjoong Kim Abstract. We
More informationModule 2: Monte Carlo Methods
Module 2: Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute MC Lecture 2 p. 1 Greeks In Monte Carlo applications we don t just want to know the expected
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationBinomial model: numerical algorithm
Binomial model: numerical algorithm S / 0 C \ 0 S0 u / C \ 1,1 S0 d / S u 0 /, S u 3 0 / 3,3 C \ S0 u d /,1 S u 5 0 4 0 / C 5 5,5 max X S0 u,0 S u C \ 4 4,4 C \ 3 S u d / 0 3, C \ S u d 0 S u d 0 / C 4
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationMath 239 Homework 1 solutions
Math 239 Homework 1 solutions Question 1. Delta hedging simulation. (a) Means, standard deviations and histograms are found using HW1Q1a.m with 100,000 paths. In the case of weekly rebalancing: mean =
More informationContinuous Time Mean Variance Asset Allocation: A Time-consistent Strategy
Continuous Time Mean Variance Asset Allocation: A Time-consistent Strategy J. Wang, P.A. Forsyth October 24, 2009 Abstract We develop a numerical scheme for determining the optimal asset allocation strategy
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More informationComputational Finance Finite Difference Methods
Explicit finite difference method Computational Finance Finite Difference Methods School of Mathematics 2018 Today s Lecture We now introduce the final numerical scheme which is related to the PDE solution.
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationFINITE DIFFERENCE METHODS
FINITE DIFFERENCE METHODS School of Mathematics 2013 OUTLINE Review 1 REVIEW Last time Today s Lecture OUTLINE Review 1 REVIEW Last time Today s Lecture 2 DISCRETISING THE PROBLEM Finite-difference approximations
More information