A Parallel Algorithm for solving BSDEs - Application to the pricing and hedging of American options

We present a parallel algorithm for solving backward stochastic differential equations (BSDEs in short) which are very useful theoretic tools to deal with many financial problems ranging from option pricing option to risk management. Our algorithm based on Gobet and Labart (2010) exploits the link between BSDEs and non linear partial differential equations (PDEs in short) and hence enables to solve high dimensional non linear PDEs. In this work, we apply it to the pricing and hedging of American options in high dimensional local volatility models, which remains very computationally demanding. We have tested our algorithm up to dimension 10 on a cluster of 512 CPUs and we obtained linear speedups which proves the scalability of our implementationComment: 25 page

Efficient hierarchical approximation of high-dimensional option pricing problems

A major challenge in computational finance is the pricing of options that depend on a large number of risk factors. Prominent examples are basket or index options where dozens or even hundreds of stocks constitute the underlying asset and determine the dimensionality of the corresponding degenerate parabolic equation. The objective of this article is to show how an efficient discretisation can be achieved by hierarchical approximation as well as asymptotic expansions of the underlying continuous problem. The relation to a number of state-of-the-art methods is highlighted

A Parallel Algorithm for solving BSDEs - Application to the pricing and hedging of American options

We present a parallel algorithm for solving backward stochastic differential equations (BSDEs in short) which are very useful theoretic tools to deal with many financial problems ranging from option pricing option to risk management. Our algorithm based on Gobet and Labart (2010) exploits the link between BSDEs and non linear partial differential equations (PDEs in short) and hence enables to solve high dimensional non linear PDEs. In this work, we apply it to the pricing and hedging of American options in high dimensional local volatility models, which remains very computationally demanding. We have tested our algorithm up to dimension 10 on a cluster of 512 CPUs and we obtained linear speedups which proves the scalability of our implementationbackward stochastic differential equations, parallel computing, Monte- Carlo methods, non linear PDE, American options, local volatility model.

Numerical methods for option pricing.

This thesis aims to introduce some fundamental concepts underlying option valuation theory including implementation of computational tools. In many cases analytical solution for option pricing does not exist, thus the following numerical methods are used: binomial trees, Monte Carlo simulations and finite difference methods. First, an algorithm based on Hull and Wilmott is written for every method. Then these algorithms are improved in different ways. For the binomial tree both speed and memory usage is significantly improved by using only one vector instead of a whole price storing matrix. Computational time in Monte Carlo simulations is reduced by implementing a parallel algorithm (in C) which is capable of improving speed by a factor which equals the number of processors used. Furthermore, MatLab code for Monte Carlo was made faster by vectorizing simulation process. Finally, obtained option values are compared to those obtained with popular finite difference methods, and it is discussed which of the algorithms is more appropriate for which purpose

An Irregular Grid Approach for Pricing High-Dimensional American Options

We propose and test a new method for pricing American options in a high-dimensional setting.The method is centred around the approximation of the associated complementarity problem on an irregular grid.We approximate the partial differential operator on this grid by appealing to the SDE representation of the underlying process and computing the root of the transition probability matrix of an approximating Markov chain.Experimental results in five dimensions are presented for four different payoff functions.option pricing;inequality;markov chains

Simulation-based Valuation and Counterparty Exposure Estimation of American Options

Valuing American options is a central problem in option pricing since the early-exercise feature is very common among financial or insurance derivatives products. For high-dimensional American options, Monte Carlo simulation is generally regarded as the only viable approach to price them, and this is the focus of our work. We propose a new regression-based Monte Carlo algorithm for pricing American options. This method typically generates an upper bound of the option value. It is computationally efficient and generates accurate price estimates. To improve the convergence rate, we apply a bias reduction technique to the least-squares Monte Carlo estimators of American option value. It works by subtracting a bias approximation from the original option value estimators at each exercise opportunity. The bias approximation is derived using large sample properties of the least-squares regression estimators. The resulting expression is easy to evaluate, and is applicable to any payoff structures and underlying processes. Numerical results show that this technique can significantly reduce the bias. However, it introduces non-negligible computational costs, thus careful treatment is required when it is adopted in practice. Finally, we extend the least-squares Monte Carlo algorithm to estimate the counterparty exposures of American options. The new algorithm is termed optimized least-squares Monte Carlo (OLSM), which is combined with variance reduction techniques, initial state dispersion and multiple bucketing to enhance its performance. The biggest advantage of OLSM is that it avoids nested simulations, allowing for the computation of risk measures on various time horizons under a reasonable computational budget

Challenges in Developing Great Quasi-Monte Carlo Software

Quasi-Monte Carlo (QMC) methods have developed over several decades. With the explosion in computational science, there is a need for great software that implements QMC algorithms. We summarize the QMC software that has been developed to date, propose some criteria for developing great QMC software, and suggest some steps toward achieving great software. We illustrate these criteria and steps with the Quasi-Monte Carlo Python library (QMCPy), an open-source community software framework, extensible by design with common programming interfaces to an increasing number of existing or emerging QMC libraries developed by the greater community of QMC researchers