Numerical methods to solve PDE models for pricing business companies in different regimes and implementation in GPUs

In this paper we propose appropriate numerical methods for companies valuation models proposed in [2]. Moreover, small modifications of these models allow to price the debt of the company and obtain the credit spread. The models are formulated in terms of final-boundary value problems associated to Kolmogorov type equations, which in some cases include an additional unilateral constraint on the solution. We also analyze the required boundary conditions so that the final-boundary value problem is well posed. This allows us to remove one of the unnecessary boundary conditions proposed in [2]. The numerical methods are mainly based on the use of characteristics (also known as semilagrangian) schemes in the direction without diffusion combined with implicit second order finite differences schemes in the direction where diffusion is present in the equation. This particular choice of numerical methods allows to develop an original parallelization strategy, which results to be specially highly efficient when using GPUs technologies.

[1]  J. Hugger.  Wellposedness of the boundary value formulation of a fixed strike Asian option , 2006 .

[2]  Eduardo S. Schwartz,et al.  Rational Pricing of Internet Companies Revisited , 2000 .

[3]  A. Pascucci,et al.  MATHEMATICAL ANALYSIS AND NUMERICAL METHODS FOR A PARTIAL DIFFERENTIAL EQUATIONS MODEL GOVERNING A RATCHET CAP PRICING IN THE LIBOR MARKET MODEL , 2011 .

[4]  Cornelis Vuik,et al.  GPU implementation of a Helmholtz Krylov solver preconditioned by a shifted Laplace multigrid method , 2011, J. Comput. Appl. Math..

[5]  N. Clarke,et al.  The parallel solution of early-exercise Asian options with stochastic volatility , 1999 .

[6]  P. Wilmott,et al.  A model for the value of a business, some optimization problems in its operating procedures and the valuation of its debt , 1997 .

[7]  O. Pironneau.  On the transport-diffusion algorithm and its applications to the Navier-Stokes equations , 1982 .

[8]  Paul Wilmott,et al.  The valuation of a firm advertising optimally , 1998 .

[9]  Carlos Vázquez,et al.  A characteristics-finite differences method for the Hobson-Rogers uncertain volatility model , 2010, Math. Comput. Model..

[10]  Eduardo S. Schwartz,et al.  Investment Under Uncertainty. , 1994 .

[11]  P. Wilmott,et al.  Option pricing: Mathematical models and computation , 1994 .

[12]  C. Cryer.  The Solution of a Quadratic Programming Problem Using Systematic Overrelaxation , 1971 .

[13]  Carlos Vázquez,et al.  Numerical Analysis of Convection-Diffusion-Reaction Problems with Higher Order Characteristics/Finite Elements. Part I: Time Discretization , 2006, SIAM J. Numer. Anal..

[14]  George Labahn,et al.  A Semi-Lagrangian Approach for American Asian Options under Jump Diffusion , 2005, SIAM J. Sci. Comput..

[15]  R. C. Merton,et al.  Continuous-Time Finance , 1990 .

[16]  Steven Li.  A valuation model for firms with stochastic earnings , 2003 .

[17]  Emmett Kilgariff,et al.  Fermi GF100 GPU Architecture , 2011, IEEE Micro.

[18]  Jason Sanders,et al.  CUDA by example: an introduction to general purpose GPU programming , 2010 .

[19]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[20]  Obstacle problem for Arithmetic Asian options , 2009, 0910.4257.

[21]  The Price of Earnings for Firms with Volatile Earnings Growth , 2008 .

[22]  Carlos Vázquez,et al.  An upwind numerical approach for an American and European option pricing model , 1998, Appl. Math. Comput..

[23]  Maria R. Nogueiras,et al.  Numerical solution of variational inequalities for pricing Asian options by higher order Lagrange-Galerkin methods , 2006 .

[24]  Emilio Barucci,et al.  Some Results on Partial Differential Equations and Asian Options , 2001 .