We present efficient partial differential equation methods for continuous time mean‐variance portfolio allocation problems when the underlying risky asset follows a jump‐diffusion. The standard formulation of mean‐variance optimal portfolio allocation problems, where the total wealth is the underlying stochastic process, gives rise to a one‐dimensional (1D) nonlinear Hamilton–Jacobi–Bellman (HJB) partial integrodifferential equation (PIDE) with the control present in the integrand of the jump term, and thus is difficult to solve efficiently. To preserve the efficient handling of the jump term, we formulate the asset allocation problem as a 2D impulse control problem, 1D for each asset in the portfolio, namely the bond and the stock. We then develop a numerical scheme based on a semi‐Lagrangian timestepping method, which we show to be monotone, consistent, and stable. Hence, assuming a strong comparison property holds, the numerical solution is guaranteed to converge to the unique viscosity solution of the corresponding HJB PIDE. The correctness of the proposed numerical framework is verified by numerical examples. We also discuss the effects on the efficient frontier of realistic financial modeling, such as different borrowing and lending interest rates, transaction costs, and constraints on the portfolio, such as maximum limits on borrowing and solvency. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 664–698, 2014

No abstract

We develop a highly efficient MC method for computing plain vanilla European option prices and hedging parameters under a very general jump-diffusion option pricing model which includes stochastic variance and multi-factor Gaussian interest short rate(s). The focus of our MC approach is variance reduction via dimension reduction. More specifically, the option price is expressed as an expectation of a unique solution to a conditional Partial Integro-Differential Equation (PIDE), which is then solved using a Fourier transform technique. Important features of our approach are (1) the analytical tractability of the conditional PIDE is fully determined by that of the Black-Scholes-Merton model augmented with the same jump component as in our model, and (2) the variances associated with all the interest rate factors are completely removed when evaluating the expectation via iterated conditioning applied to only the Brownian motion associated with the variance factor. For certain cases when numerical methods are either needed or preferred, we propose a discrete fast Fourier transform method to numerically solve the conditional PIDE efficiently. Our method can also effectively compute hedging parameters. Numerical results show that the proposed method is highly efficient.

We develop space-time adaptive and high-order methods for valuing American options using a partial differential equation (PDE) approach. The linear complementarity problem arising due to the free boundary is handled by a penalty method. Both finite difference and finite element methods are considered for the space discretization of the PDE, while classical finite differences, such as Crank-Nicolson, are used for the time discretization. The high-order discretization in space is based on an optimal finite element collocation method, the main computational requirements of which are the solution of one tridiagonal linear system at each time step, while the resulting errors at the gridpoints and midpoints of the space partition are fourth-order. To control the space error, we use adaptive gridpoint distribution based on an error equidistribution principle. A time stepsize selector is used to further increase the efficiency of the methods. Numerical examples show that our methods converge fast and provide highly accurate options prices, Greeks, and early exercise boundaries.

No abstract