Primal-dual active-set method for solving the unilateral pricing problem of American better-of options on two assets

In this paper, an efficient numerical algorithm is proposed for the valuation of unilateral American better-of options with two underlying assets. The pricing model can be described as a backward parabolic variational inequality with variable coefficients on a two-dimensional unbounded domain. It can be transformed into a one-dimensional bounded free boundary problem by some conventional transformations and the far-field truncation technique. With appropriate boundary conditions on the free boundary, a bounded linear complementary problem corresponding to the option pricing is established. Furthermore, the full discretization scheme is obtained by applying the backward Euler method and the finite element method in temporal and spatial directions, respectively. Based on the symmetric positive definite property of the discretized matrix, the value of the option and the free boundary are obtained simultaneously by the primal-dual active-set method. The error estimation is established by the variational theory. Numerical experiments are carried out to verify the efficiency of our method at the end.