Abstract:We propose an iterative method for pricing American options under jumpdiffusion models. A finite difference discretization is performed on the partial integro-differential equation, and the American option pricing problem is formulated as a linear complementarity problem (LCP). Jump-diffusion models include an integral term, which causes the resulting system to be dense. We propose an iteration to solve the LCPs efficiently and prove its convergence. Numerical examples with Kou's and Merton's jump-diffusion mo… Show more
“…The integrals can be discretized using a second-order accurate quadrature formula. Here the linear interpolation is used for w between grid points and exact integration; see [11], for details. Under the Merton model the discretization of the integral leads to a full matrix while under the Bates model it leads to full diagonal blocks.…”
Section: Full Order Modelsmentioning
confidence: 99%
“…A penalty method together with FFT based fast method for evaluating the jump integral was used in [9]. An iterative method was proposed for LCPs with full matrices in [11]. An implicit-explicit (IMEX) method was proposed in [10] to treat the integral term explicitly and the same approach was studied in [13].…”
European options can be priced by solving parabolic partial(-integro) differential equations under stochastic volatility and jump-diffusion models like Heston, Merton, and Bates models. American option prices can be obtained by solving linear complementary problems (LCPs) with the same operators. A finite difference discretization leads to a so-called full order model (FOM). Reduced order models (ROMs) are derived employing proper orthogonal decomposition (POD). The early exercise constraint of American options is enforced by a penalty on subset of grid points. The presented numerical experiments demonstrate that pricing with ROMs can be orders of magnitude faster within a given model parameter variation range.
“…The integrals can be discretized using a second-order accurate quadrature formula. Here the linear interpolation is used for w between grid points and exact integration; see [11], for details. Under the Merton model the discretization of the integral leads to a full matrix while under the Bates model it leads to full diagonal blocks.…”
Section: Full Order Modelsmentioning
confidence: 99%
“…A penalty method together with FFT based fast method for evaluating the jump integral was used in [9]. An iterative method was proposed for LCPs with full matrices in [11]. An implicit-explicit (IMEX) method was proposed in [10] to treat the integral term explicitly and the same approach was studied in [13].…”
European options can be priced by solving parabolic partial(-integro) differential equations under stochastic volatility and jump-diffusion models like Heston, Merton, and Bates models. American option prices can be obtained by solving linear complementary problems (LCPs) with the same operators. A finite difference discretization leads to a so-called full order model (FOM). Reduced order models (ROMs) are derived employing proper orthogonal decomposition (POD). The early exercise constraint of American options is enforced by a penalty on subset of grid points. The presented numerical experiments demonstrate that pricing with ROMs can be orders of magnitude faster within a given model parameter variation range.
“…The authors in 20 derive an analytical formula for the price of European options for any model including local volatility and Poisson jump process by using Malliavin calculus techniques. Various authors apart of 14 used finite difference schemes for PIDEs in [21][22][23][24][25][26][27] . Discretization of the integral term leads to full matrices due to its nonlocal character.…”
A new discretization strategy is introduced for the numerical solution of partial integrodifferential equations appearing in option pricing jump diffusion models. In order to consider the unknown behaviour of the solution in the unbounded part of the spatial domain, a double discretization is proposed. Stability, consistency, and positivity of the resulting explicit scheme are analyzed. Advantages of the method are illustrated with several examples.
“…Many authors used the finite difference (FD) schemes for solving these PIDE problems [2,4,5,17,25,54,74,75,82,85,87]. Dealing with FD methods for such PIDEs, the following challenges should be addressed.…”
Section: Consistency For Integral Equationmentioning
confidence: 99%
“…Tavella and Randall in [82] use an implicit time discretization and propose a stationary rapid convergent iterative method to solve the full matrix problem quoted above, but with poor numerical analysis. A generalization of their iterative method to price American options is proposed in [75].…”
Section: Consistency For Integral Equationmentioning
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.