Modulus-based Successive Overrelaxation Iteration Method for Pricing American Options with the Two-asset Black–Scholes and Heston's Models Based on Finite Volume Discretization

In this paper we introduce a new numerical method for the linear complementarity problems (LCPs) arising from two-asset Black–Scholes and Heston's stochastic volatility American options pricing. Based on barycenter dual mesh, a class of finite volume method (FVM) is proposed for the spatial discretization, coupled with the backward Euler and Crank–Nicolson schemes are employed for time stepping of the partial differential equations (PDEs). Then, for the resulting time-dependent LCPs are solved by using an efficient modulus-based successive overrelaxation (MSOR) iteration method. Numerical experiments are carried out to verify the efficiency and usefulness of the proposed method.