Multistep schemes for solving backward stochastic differential equations on GPU

The goal of this work is to parallelize the multistep scheme for the numerical approximation of the backward stochastic differential equations (BSDEs) in order to achieve both, a high accuracy and a reduction of the computation time as well. In the multistep scheme the computations at each grid point are independent and this fact motivates us to select massively parallel GPU computing using CUDA. In our investigations we identify performance bottlenecks and apply appropriate optimization techniques for reducing the computation time, using a uniform domain. Finally, some examples with financial applications are provided to demonstrate the achieved acceleration on GPUs.

[1]  Weidong Zhao,et al.  A New Kind of Accurate Numerical Method for Backward Stochastic Differential Equations , 2006, SIAM J. Sci. Comput..

[2]  Christian Bender,et al.  Least-Squares Monte Carlo for Backward SDEs , 2012 .

[3]  Hui Liu,et al.  Option Pricing on the GPU with Backward Stochastic Differential Equation , 2011, 2011 Fourth International Symposium on Parallel Architectures, Algorithms and Programming.

[4]  Dan Crisan,et al.  Solving Backward Stochastic Differential Equations Using the Cubature Method: Application to Nonlinear Pricing , 2010, SIAM J. Financial Math..

[5]  E. Gobet,et al.  Rate of convergence of an empirical regression method for solving generalized backward stochastic differential equations , 2006 .

[6]  Shige Peng,et al.  Probabilistic interpretation for systems of quasilinear parabolic partial differential equations , 1991 .

[7]  S. Peng,et al.  Adapted solution of a backward stochastic differential equation , 1990 .

[8]  Irene A. Stegun,et al.  Handbook of Mathematical Functions. , 1966 .

[9]  Céline Labart,et al.  A Parallel Algorithm for solving BSDEs - Application to the pricing and hedging of American options , 2011, ArXiv.

[10]  Carlos Vázquez,et al.  Stratified Regression Monte-Carlo Scheme for Semilinear PDEs and BSDEs with Large Scale Parallelization on GPUs , 2016, SIAM J. Sci. Comput..

[11]  Ying Peng,et al.  Parallel Option Pricing with BSDE Method on GPU , 2010, 2010 Ninth International Conference on Grid and Cloud Computing.

[12]  Christophette Blanchet-Scalliet,et al.  CREDIT RISK PREMIA AND QUADRATIC BSDEs WITH A SINGLE JUMP , 2009, 0907.1221.

[13]  B. Bouchard,et al.  Discrete-time approximation and Monte-Carlo simulation of backward stochastic differential equations , 2004 .

[14]  Michael Gunther,et al.  A multi-step scheme based on cubic spline for solving backward stochastic differential equations , 2018, 1809.00324.

[15]  Nizar Touzi,et al.  A Probabilistic Numerical Method for Fully Nonlinear Parabolic PDEs , 2009, 0905.1863.

[16]  E. Gobet,et al.  A regression-based Monte Carlo method to solve backward stochastic differential equations , 2005, math/0508491.

[17]  Hai-ping Shi.  Backward stochastic differential equations in finance , 2010 .

[18]  Weidong Zhao,et al.  A Stable Multistep Scheme for Solving Backward Stochastic Differential Equations , 2010, SIAM J. Numer. Anal..

[19]  Anne Eyraud-Loisel.  Backward stochastic differential equations with enlarged filtration: Option hedging of an insider trader in a financial market with jumps , 2005 .

[20]  Olof B. Widlund,et al.  Smoothing of initial data and rates of convergence for parabolic difference equations , 1970 .