Path-dependent option pricing with two-dimensional PDE using MPDATA

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Magnuszewski, Paweł, Arabas, Sylwester
Formato: Preprint
Publicado: 2025
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866909629047373824
author Magnuszewski, Paweł
Arabas, Sylwester
author_facet Magnuszewski, Paweł
Arabas, Sylwester
contents In this paper, we discuss a simple yet robust PDE method for evaluating path-dependent Asian-style options using the non-oscillatory forward-in-time second-order MPDATA finite-difference scheme. The valuation methodology involves casting the Black-Merton-Scholes equation as a transport problem by first transforming it into a homogeneous advection-diffusion PDE via variable substitution, and then expressing the diffusion term as an advective flux using the pseudo-velocity technique. As a result, all terms of the Black-Merton-Sholes equation are consistently represented using a single high-order numerical scheme for the advection operator. We detail the additional steps required to solve the two-dimensional valuation problem compared to MPDATA valuations of vanilla instruments documented in a prior study. Using test cases employing fixed-strike instruments, we validate the solutions against Monte Carlo valuations, as well as against an approximate analytical solution in which geometric instead of arithmetic averaging is used. The analysis highlights the critical importance of the MPDATA corrective steps that improve the solution over the underlying first-order "upwind" step. The introduced valuation scheme is robust: conservative, non-oscillatory, and positive-definite; yet lucid: explicit in time, engendering intuitive stability-condition interpretation and inflow/outflow boundary-condition heuristics. MPDATA is particularly well suited for two-dimensional problems as it is not a dimensionally split scheme. The documented valuation workflow also constitutes a useful two-dimensional case for testing advection schemes featuring both Monte Carlo solutions and analytic bounds. An implementation of the introduced valuation workflow, based on the PyMPDATA package and the Numba Just-In-Time compiler for Python, is provided as free and open source software.
format Preprint
id arxiv_https___arxiv_org_abs_2505_24435
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Path-dependent option pricing with two-dimensional PDE using MPDATA
Magnuszewski, Paweł
Arabas, Sylwester
Computational Finance
In this paper, we discuss a simple yet robust PDE method for evaluating path-dependent Asian-style options using the non-oscillatory forward-in-time second-order MPDATA finite-difference scheme. The valuation methodology involves casting the Black-Merton-Scholes equation as a transport problem by first transforming it into a homogeneous advection-diffusion PDE via variable substitution, and then expressing the diffusion term as an advective flux using the pseudo-velocity technique. As a result, all terms of the Black-Merton-Sholes equation are consistently represented using a single high-order numerical scheme for the advection operator. We detail the additional steps required to solve the two-dimensional valuation problem compared to MPDATA valuations of vanilla instruments documented in a prior study. Using test cases employing fixed-strike instruments, we validate the solutions against Monte Carlo valuations, as well as against an approximate analytical solution in which geometric instead of arithmetic averaging is used. The analysis highlights the critical importance of the MPDATA corrective steps that improve the solution over the underlying first-order "upwind" step. The introduced valuation scheme is robust: conservative, non-oscillatory, and positive-definite; yet lucid: explicit in time, engendering intuitive stability-condition interpretation and inflow/outflow boundary-condition heuristics. MPDATA is particularly well suited for two-dimensional problems as it is not a dimensionally split scheme. The documented valuation workflow also constitutes a useful two-dimensional case for testing advection schemes featuring both Monte Carlo solutions and analytic bounds. An implementation of the introduced valuation workflow, based on the PyMPDATA package and the Numba Just-In-Time compiler for Python, is provided as free and open source software.
title Path-dependent option pricing with two-dimensional PDE using MPDATA
topic Computational Finance
url https://arxiv.org/abs/2505.24435