A PRIMAL–DUAL ALGORITHM FOR BSDES

We generalize the primal-dual methodology, which is popular in the pricing of early-exercise options, to a backward dynamic programming equation associated with time discretization schemes of (reflected) backward stochastic differential equations (BSDEs). Taking as an input some approximate solution of the backward dynamic program, which was pre-computed, e.g., by least-squares Monte Carlo, our methodology allows to construct a confidence interval for the unknown true solution of the time discretized (reflected) BSDE at time 0. We numerically demonstrate the practical applicability of our method in two five-dimensional nonlinear pricing problems where tight price bounds were previously unavailable.

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

[2]  K. Ritter,et al.  Adaptive Wavelet Methods for Elliptic Stochastic Partial Differential Equations , 2022 .

[3]  Winfried Sickel,et al.  Best m-term approximation and Lizorkin-Triebel spaces , 2011, J. Approx. Theory.

[4]  Robert Denk,et al.  A forward scheme for backward SDEs , 2007 .

[5]  H. Yserentant,et al.  On the Complexity of the Electronic Schrödinger Equation , 2022 .

[6]  Armin Iske,et al.  Curvature analysis of frequency modulated manifolds in dimensionality reduction , 2011 .

[7]  Marco Avellaneda,et al.  Reducing variance in the numerical solution of BSDEs , 2013 .

[8]  M. Hansen,et al.  n-term approximation rates and Besov regularity for elliptic PDEs on polyhedral domains , 2012 .

[9]  Simen Kvaal,et al.  Multiconfigurational time-dependent Hartree method to describe particle loss due to absorbing boundary conditions , 2011, 1102.3899.

[10]  M. Gunzburger,et al.  A SPARSE-GRID METHOD FOR MULTI-DIMENSIONAL BACKWARD STOCHASTIC DIFFERENTIAL EQUATIONS * , 2013 .

[11]  D. Belomestny,et al.  TRUE UPPER BOUNDS FOR BERMUDAN PRODUCTS VIA NON‐NESTED MONTE CARLO , 2009 .

[12]  Michael Gnewuch,et al.  On weighted Hilbert spaces and integration of functions of infinitely many variables , 2014, J. Complex..

[13]  Wolfgang Dahmen,et al.  Adaptive Petrov-Galerkin methods for first order transport equations , 2011 .

[14]  E. Ullmann,et al.  Preconditioning Stochastic Galerkin Saddle Point Problems , 2022 .

[15]  Petru A. Cioica,et al.  On the $L_q(L_p)$-regularity and Besov smoothness of stochastic parabolic equations on bounded Lipschitz domains , 2013, 1301.1180.

[16]  Klaus Ritter,et al.  A Local Refinement Strategy for Constructive Quantization of Scalar SDEs , 2013, Found. Comput. Math..

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

[18]  André Uschmajew,et al.  Local Convergence of the Alternating Least Squares Algorithm for Canonical Tensor Approximation , 2012, SIAM J. Matrix Anal. Appl..

[19]  J. Ballani,et al.  Black box approximation of tensors in hierarchical Tucker format , 2013 .

[20]  Felix Lindner,et al.  Singular behavior of the solution to the stochastic heat equation on a polygonal domain , 2013, Stochastic Partial Differential Equations: Analysis and Computations.

[21]  Steffen Dereich,et al.  A multilevel Monte Carlo algorithm for Lévy-driven stochastic differential equations , 2011 .

[22]  R. Rockafellar.  Convex Analysis: (pms-28) , 1970 .

[23]  Stephan Dahlke,et al.  A note on quarkonial systems and multilevel partition of unity methods , 2013 .

[24]  G. Kutyniok,et al.  Irregular Shearlet Frames: Geometry and Approximation Properties , 2010, 1002.2657.

[25]  Jean-Paul Laurent,et al.  An overview of the valuation of collateralized derivative contracts , 2014 .

[26]  Wolfgang Dahmen,et al.  Classification algorithms using adaptive partitioning , 2014, 1411.0839.

[27]  Jin Ma,et al.  Representations and regularities for solutions to BSDEs with reflections , 2005 .

[28]  Michael Döhler,et al.  Nonequispaced Hyperbolic Cross Fast Fourier Transform , 2010, SIAM J. Numer. Anal..

[29]  Harry Yserentant,et al.  The mixed regularity of electronic wave functions multiplied by explicit correlation factors , 2011 .

[30]  Wang-Q Lim,et al.  Image Separation Using Shearlets , 2011 .

[31]  B. Bouchard,et al.  Discrete time approximation of decoupled Forward-Backward SDE with jumps , 2008 .

[32]  W. Hackbusch,et al.  Black Box Low Tensor-Rank Approximation Using Fiber-Crosses , 2009 .

[33]  Reinhold Schneider,et al.  Error estimates for the Coupled Cluster method , 2013 .

[34]  Howard C. Elman,et al.  Efficient Iterative Solvers for Stochastic Galerkin Discretizations of Log-Transformed Random Diffusion Problems , 2012, SIAM J. Sci. Comput..

[35]  Wolfgang Dahmen,et al.  Fast high-dimensional approximation with sparse occupancy trees , 2011, J. Comput. Appl. Math..

[36]  Pierre Henry-Labordere,et al.  Uncertain Volatility Model: A Monte-Carlo Approach , 2010 .

[37]  Lutz Kämmerer.  Reconstructing hyperbolic cross trigonometric polynomials by sampling along generated sets , 2012 .

[38]  Steffen Dereich,et al.  Multilevel Monte Carlo algorithms for L\'{e}vy-driven SDEs with Gaussian correction , 2011, 1101.1369.

[39]  Erich Novak,et al.  Optimal approximation of elliptic problems by linear and nonlinear mappings IV: Errors in L2 and other norms , 2004, J. Complex..

[40]  B. Bouchard,et al.  Discrete time approximation for continuously and discretely re ected BSDEs , 2007 .

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

[42]  G. Kutyniok,et al.  Construction of Compactly Supported Shearlet Frames , 2010, 1003.5481.

[43]  O. Ernst,et al.  ON THE CONVERGENCE OF GENERALIZED POLYNOMIAL CHAOS EXPANSIONS , 2011 .

[44]  Vivek F. Farias,et al.  Pathwise Optimization for Optimal Stopping Problems , 2012, Manag. Sci..

[45]  Wolfgang Dahmen,et al.  Convergence Rates for Greedy Algorithms in Reduced Basis Methods , 2010, SIAM J. Math. Anal..

[46]  S. Peng,et al.  Backward Stochastic Differential Equations in Finance , 1997 .

[47]  E. Novak,et al.  On the power of function values for the approximation problem in various settings , 2010, 1011.3682.

[48]  Jianfeng Zhang.  A numerical scheme for BSDEs , 2004 .

[49]  Christian Bender,et al.  Dual pricing of multi-exercise options under volume constraints , 2011, Finance Stochastics.

[50]  Harry Yserentant,et al.  The mixed regularity of electronic wave functions in fractional order and weighted Sobolev spaces , 2012, Numerische Mathematik.

[51]  K. Bredies,et al.  Regularization with non-convex separable constraints , 2009 .

[52]  Dirk A. Lorenz,et al.  Minimization of Non-smooth, Non-convex Functionals by Iterative Thresholding , 2015, J. Optim. Theory Appl..

[53]  Armin Iske,et al.  On Groupoid C∗-Algebras, Persistent Homology and Time-Frequency Analysis , 2011 .

[54]  Daniela Rosca,et al.  A New Hybrid Method for Image Approximation Using the Easy Path Wavelet Transform , 2011, IEEE Transactions on Image Processing.

[55]  D. Duffie,et al.  Recursive valuation of defaultable securities and the timing of resolution of uncertainty , 1996 .

[56]  Erich Novak,et al.  The Curse of Dimensionality for Monotone and Convex Functions of Many Variables , 2010, 1011.3680.

[57]  Christian Bender,et al.  Error Criteria for Numerical Solutions of Backward SDEs , 2010 .

[58]  Wang-Q Lim,et al.  Compactly supported shearlets are optimally sparse , 2010, J. Approx. Theory.

[59]  Wolfgang Dahmen,et al.  Adaptive Near-Optimal Rank Tensor Approximation for High-Dimensional Operator Equations , 2013, Foundations of Computational Mathematics.

[60]  Henryk Wozniakowski,et al.  The curse of dimensionality for numerical integration of smooth functions , 2012, Math. Comput..

[61]  Thorsten Rohwedder,et al.  The continuous Coupled Cluster formulation for the electronic Schrödinger equation , 2013 .

[62]  Denis Belomestny,et al.  Pricing Bermudan options using nonparametric regression: optimal rates of convergence for lower estimates , 2009, 0907.5599.

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

[64]  Reinhold Schneider,et al.  On manifolds of tensors of fixed TT-rank , 2012, Numerische Mathematik.

[65]  G. Pagès,et al.  A quantization algorithm for solving multidimensional discrete-time optimal stopping problems , 2003 .

[66]  Dajana Conte,et al.  Mathematical Modelling and Numerical Analysis an Error Analysis of the Multi-configuration Time-dependent Hartree Method of Quantum Dynamics , 2022 .

[67]  Martin B. Haugh,et al.  Pricing American Options: A Duality Approach , 2001, Oper. Res..

[68]  Peng Sun,et al.  Information Relaxations and Duality in Stochastic Dynamic Programs , 2010, Oper. Res..

[69]  Mark Broadie,et al.  A Primal-Dual Simulation Algorithm for Pricing Multi-Dimensional American Options , 2001 .

[70]  Michael Griebel,et al.  An Adaptive Sparse Grid Semi-Lagrangian Scheme for First Order Hamilton-Jacobi Bellman Equations , 2012, Journal of Scientific Computing.

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

[72]  Markus Hansen,et al.  On tensor products of quasi-Banach spaces , 2010 .

[73]  Wang-Q Lim,et al.  Compactly Supported Shearlets , 2010, 1009.4359.

[74]  Karsten Urban,et al.  On Wavelet-Galerkin Methods for Semilinear Parabolic Equations with Additive Noise , 2013 .

[75]  P. Maass,et al.  An analysis of electrical impedance tomography with applications to Tikhonov regularization , 2012 .

[76]  H. Harbrecht,et al.  On the low-rank approximation by the pivoted Cholesky decomposition , 2012 .

[77]  Andreas Neuenkirch,et al.  Multilevel Monte Carlo Quadrature of Discontinuous Payoffs in the Generalized Heston Model Using Malliavin Integration by Parts , 2015, SIAM J. Financial Math..

[78]  John Schoenmakers,et al.  Optimal Dual Martingales, Their Analysis, and Application to New Algorithms for Bermudan Products , 2011, SIAM J. Financial Math..

[79]  Winfried Sickel,et al.  Best m-term aproximation and tensor product of Sobolev and Besov spaces-the case of non-compact embeddings , 2010 .

[80]  D. Brigo,et al.  Funding, Collateral and Hedging: Uncovering the Mechanics and the Subtleties of Funding Valuation Adjustments , 2012, 1210.3811.

[81]  Steffen Dereich,et al.  Foundations of Computational Mathematics, Budapest 2011: On the Complexity of Computing Quadrature Formulas for SDEs , 2012 .

[82]  Karsten Urban,et al.  Adaptive Wavelet Methods on Unbounded Domains , 2012, Journal of Scientific Computing.

[83]  D. Rudolf,et al.  Explicit error bounds for Markov chain Monte Carlo , 2011, 1108.3201.

[84]  Klaus Ritter,et al.  On embeddings of weighted tensor product Hilbert spaces , 2015, J. Complex..

[85]  Lars Grasedyck,et al.  Polynomial Approximation in Hierarchical Tucker Format by Vector – Tensorization , 2010 .

[86]  Fred J. Hickernell,et al.  Multi-level Monte Carlo algorithms for infinite-dimensional integration on RN , 2010, J. Complex..

[87]  Wang-Q Lim,et al.  Shearlets on Bounded Domains , 2010, 1007.3039.

[88]  G. Teschke,et al.  Inversion of the noisy Radon transform on SO(3) by Gabor frames and sparse recovery principles , 2011 .

[89]  Sadegh Jokar,et al.  Sparse recovery and Kronecker products , 2010, 2010 44th Annual Conference on Information Sciences and Systems (CISS).

[90]  M. Joshi,et al.  Effective Sub-Simulation-Free Upper Bounds for the Monte Carlo Pricing of Callable Derivatives and Various Improvements to Existing Methodologies , 2012 .

[91]  Paul Glasserman,et al.  Simulation for American Options: Regression Now or Regression Later? , 2004 .

[92]  Reinhold Schneider,et al.  Optimization problems in contracted tensor networks , 2011, Comput. Vis. Sci..

[93]  Lars Grasedyck,et al.  F ¨ Ur Mathematik in Den Naturwissenschaften Leipzig a Projection Method to Solve Linear Systems in Tensor Format a Projection Method to Solve Linear Systems in Tensor Format , 2022 .

[94]  Samuel N. Cohen,et al.  A general theory of finite state Backward Stochastic Difference Equations , 2008, 0810.4957.

[95]  Henryk Wozniakowski,et al.  Discontinuous information in the worst case and randomized settings , 2011, 1106.2945.

[96]  Ianwei,et al.  Compressive Video Sampling with Approximate Message Passing Decoding , 2011 .

[97]  Reinhold Schneider,et al.  An analysis for the DIIS acceleration method used in quantum chemistry calculations , 2011 .

[98]  Konstantin Grella,et al.  Sparse tensor spherical harmonics approximation in radiative transfer , 2011, J. Comput. Phys..

[99]  Stephan Dahlke,et al.  Spatial Besov regularity for semilinear stochastic partial differential equations on bounded Lipschitz domains , 2012, Int. J. Comput. Math..

[100]  Wolfgang Hackbusch,et al.  Tensorisation of vectors and their efficient convolution , 2011, Numerische Mathematik.

[101]  Wolfgang Dahmen,et al.  Adaptive Petrov-Galerkin Methods for First Order Transport Equations , 2011, SIAM J. Numer. Anal..

[102]  Wolfgang Hackbusch,et al.  An Introduction to Hierarchical (H-) Rank and TT-Rank of Tensors with Examples , 2011, Comput. Methods Appl. Math..

[103]  Harry Yserentant,et al.  A Short Theory of the Rayleigh–Ritz Method , 2013, Comput. Methods Appl. Math..

[104]  Henryk Wozniakowski,et al.  Approximation of infinitely differentiable multivariate functions is intractable , 2009, J. Complex..

[105]  Jianwei Ma,et al.  A Review of Curvelets and Recent Applications , 2009 .

[106]  Armin Iske,et al.  Optimally Sparse Image Representation by the Easy Path Wavelet Transform , 2012, Int. J. Wavelets Multiresolution Inf. Process..

[107]  K. Ritter,et al.  Variable Subspace Sampling and Multi-level Algorithms , 2009 .

[108]  Yaacov Z. Bergman.  Option Pricing with Differential Interest Rates , 1995 .

[109]  Stéphane Crépey,et al.  Counterparty Risk and Funding: The Four Wings of the TVA , 2012, 1210.5046.

[110]  Gerlind Plonka-Hoch,et al.  Curvelet-Wavelet Regularized Split Bregman Iteration for Compressed Sensing , 2011, Int. J. Wavelets Multiresolution Inf. Process..

[111]  G. Teschke,et al.  Compressive sensing principles and iterative sparse recovery for inverse and ill-posed problems , 2010 .

[112]  M. Joshi.  Analyzing the Bias in the Primal-Dual Upper Bound Method for Early Exercisable Derivatives: Bounds, Estimation and Removal , 2015 .

[113]  Klaus Ritter,et al.  Derandomization of the Euler scheme for scalar stochastic differential equations , 2012, J. Complex..

[114]  Hermann G. Matthies,et al.  Efficient low-rank approximation of the stochastic Galerkin matrix in tensor formats , 2014, Comput. Math. Appl..

[115]  Tobias Jahnke,et al.  On Reduced Models for the Chemical Master Equation , 2011, Multiscale Model. Simul..

[116]  D. Crisan,et al.  Robust filtering: Correlated noise and multidimensional observation , 2012, 1201.1858.

[117]  Daniel Rudolf,et al.  Error bounds for computing the expectation by Markov chain Monte Carlo , 2009, Monte Carlo Methods Appl..

[118]  Tobias Jahnke,et al.  Error Bound for Piecewise Deterministic Processes Modeling Stochastic Reaction Systems , 2012, Multiscale Model. Simul..

[119]  Armin Iske,et al.  Adaptive ADER Methods Using Kernel-Based Polyharmonic Spline WENO Reconstruction , 2010, SIAM J. Sci. Comput..

[120]  M. Fornasier,et al.  Multilevel Preconditioning for Adaptive Sparse Optimization , 2009 .

[121]  Klaus Ritter,et al.  Spatial Besov Regularity for Stochastic Partial Differential Equations on Lipschitz Domains , 2010, 1011.1814.

[122]  Lutz Kämmerer,et al.  On the stability of the hyperbolic cross discrete Fourier transform , 2011, Numerische Mathematik.

[123]  Wolfgang Dahmen,et al.  DOUBLE GREEDY ALGORITHMS: REDUCED BASIS METHODS FOR TRANSPORT DOMINATED PROBLEMS ∗ , 2013, 1302.5072.

[124]  K. Ritter,et al.  On the convergence analysis of Rothe ’ s method , 2022 .

[125]  H. Johnson.  Options on the Maximum or the Minimum of Several Assets , 1987, Journal of Financial and Quantitative Analysis.

[126]  L. Rogers.  Monte Carlo valuation of American options , 2002 .

[127]  Erwan Faou,et al.  Computing Semiclassical Quantum Dynamics with Hagedorn Wavepackets , 2009, SIAM J. Sci. Comput..

[128]  Andreas Zeiser,et al.  Wavelet Approximation in Weighted Sobolev Spaces of Mixed Order with Applications to the Electronic Schrödinger Equation , 2012 .

[129]  Denis Belomestny,et al.  Multilevel dual approach for pricing American style derivatives , 2012, Finance Stochastics.

[130]  D. Potts,et al.  Approximation of multivariate functions by trigonometric polynomials based on rank-1 lattice sampling , 2013 .

[131]  Wolfgang Dahmen,et al.  Compressed Sensing and Electron Microscopy , 2010 .

[132]  Lars Grasedyck,et al.  Tree Adaptive Approximation in the Hierarchical Tensor Format , 2014, SIAM J. Sci. Comput..

[133]  Gerlind Plonka-Hoch,et al.  The Easy Path Wavelet Transform: A New Adaptive Wavelet Transform for Sparse Representation of Two-Dimensional Data , 2008, Multiscale Model. Simul..

[134]  Emmanuel Gobet,et al.  L2-time regularity of BSDEs with irregular terminal functions , 2010 .

[135]  Reinhold Schneider,et al.  Numerical analysis of augmented plane wave methods for full-potential electronic structure calculations , 2015 .

[136]  P. Maass,et al.  Function spaces and optimal currents in impedance tomography , 2011 .

[137]  G. Teschke,et al.  Accelerated projected steepest descent method for nonlinear inverse problems with sparsity constraints , 2010 .

[138]  Lutz Kämmerer,et al.  Interpolation lattices for hyperbolic cross trigonometric polynomials , 2012, J. Complex..

[139]  Stephan Dahlke,et al.  An adaptive wavelet method for parameter identification problems in parabolic partial differential equations , 2022 .

[140]  Gabriele Steidl,et al.  Shearlet Coorbit Spaces: Compactly Supported Analyzing Shearlets, Traces and Embeddings , 2011 .

[141]  Reinhold Schneider,et al.  Numerical analysis of Gaussian approximations in quantum chemistry , 2012 .

[142]  Elisabeth Ullmann,et al.  Stochastic Galerkin Matrices , 2010, SIAM J. Matrix Anal. Appl..

[143]  E. Novak,et al.  Optimal Order of Convergence and (In)Tractability of Multivariate Approximation of Smooth Functions , 2009 .

[144]  Wolfgang Dahmen,et al.  Adaptivity and variational stabilization for convection-diffusion equations∗ , 2012 .

[145]  Lars Grasedyck,et al.  Hierarchical Singular Value Decomposition of Tensors , 2010, SIAM J. Matrix Anal. Appl..

[146]  Torsten Görner,et al.  Efficient and accurate computation of spherical mean values at scattered center points , 2012 .

[147]  Christian Bender,et al.  Primal and Dual Pricing of Multiple Exercise Options in Continuous Time , 2011, SIAM J. Financial Math..

[148]  Tobias Jahnke,et al.  Solving chemical master equations by adaptive wavelet compression , 2010, J. Comput. Phys..

[149]  S. Dereich,et al.  Constructive quantization: Approximation by empirical measures , 2011, 1108.5346.

[150]  Stephan Dahlke,et al.  Piecewise tensor product wavelet bases by extensions and approximation rates , 2013, Math. Comput..

[151]  D. Lorenz,et al.  Greedy solution of ill-posed problems: error bounds and exact inversion , 2009, 0904.0154.

[152]  Winfried Sickel,et al.  On Besov regularity of solutions to nonlinear elliptic partial differential equations , 2020, Nonlinear Analysis.

[153]  Fred J. Hickernell,et al.  Deterministic multi-level algorithms for infinite-dimensional integration on RN , 2011, J. Complex..

[154]  Wolfgang Dahmen,et al.  Approximation of High-Dimensional Rank One Tensors , 2013, Constructive Approximation.

[155]  Gitta Kutyniok,et al.  Data Separation by Sparse Representations , 2011, Compressed Sensing.

[156]  R. Schilling,et al.  Weak Order for the Discretization of the Stochastic Heat Equation Driven by Impulsive Noise , 2009, 0911.4681.

[157]  Stephan Dahlke,et al.  Adaptive wavelet methods and sparsity reconstruction for inverse heat conduction problems , 2010, Adv. Comput. Math..

[158]  Emmanuel Gobet,et al.  Error expansion for the discretization of backward stochastic differential equations , 2006, math/0602503.

[159]  Armin Iske,et al.  Optimal representation of piecewise Hölder smooth bivariate functions by the Easy Path Wavelet Transform , 2013, J. Approx. Theory.

[160]  G. Teschke,et al.  A compressive Landweber iteration for solving ill-posed inverse problems , 2008 .

[161]  Winfried Sickel,et al.  Best m-Term Approximation and Sobolev–Besov Spaces of Dominating Mixed Smoothness—the Case of Compact Embeddings , 2012 .

[162]  R. Schneider,et al.  The Alternating Linear Scheme for Tensor Optimisation in the TT Format , 2022 .

[163]  Reinhold Schneider,et al.  Dynamical Approximation by Hierarchical Tucker and Tensor-Train Tensors , 2013, SIAM J. Matrix Anal. Appl..

[164]  Denis Belomestny,et al.  SOLVING OPTIMAL STOPPING PROBLEMS VIA EMPIRICAL DUAL OPTIMIZATION , 2013, 1309.2125.

[165]  Patrick Cheridito,et al.  BS\Delta Es and BSDEs with non-Lipschitz drivers: Comparison, convergence and robustness , 2010, 1002.0175.

[166]  Christian Bender,et al.  A Posteriori Estimates for Backward SDEs , 2013, SIAM/ASA J. Uncertain. Quantification.

[167]  Hermann G. Matthies,et al.  Efficient Analysis of High Dimensional Data in Tensor Formats , 2012 .

[168]  A. Uschmajew,et al.  LOCAL CONVERGENCE OF ALTERNATING SCHEMES FOR OPTIMIZATION OF CONVEX PROBLEMS IN THE TT FORMAT , 2011 .