Efficiency of Some Predictor-Corrector Methods with Fourth-Order Compact Scheme for a System of Free Boundary Options

The trade-off between numerical accuracy and computational cost is always an important factor to consider when pricing options numerically, due to the inherent irregularity and existence of non-linearity in many models. In this work, we first present fast and accurate (1,2) and (2,2) predictor-corrector methods with a fourth-order compact finite difference scheme for pricing coupled system of the non-linear free boundary option pricing problem consisting of the option value and delta sensitivity. To predict the optimal exercise boundary, we set up a high-order boundary scheme, which is strategically derived using a combination of the fourth-order Robin boundary scheme and the fourth-order compact finite difference scheme near boundary. Furthermore, we implement a three-step high-order correction scheme for computing interior values of the option value and delta sensitivity. The discrete matrix system of this correction scheme has a tri-diagonal structure and strictly diagonal dominance. This nice feature allows for the implementation of the Thomas algorithm, thereby enabling fast computation. The optimal exercise boundary value is also corrected in each of the three correction steps with the derived Robin boundary scheme. Our implementations are fast on both coarse and very refined grids and provide highly accurate numerical approximations. Moreover, we recover a reasonable convergence rate. Further extensions to high-order predictor two-step corrector schemes are elaborated.