"On the Solution of Complementarity Problems Arising in American Options Pricing"
L. Feng, V. Linetsky, J.L. Morales, J. Nocedal
<Technical Report, Optimization Center, Northwestern University (2009).

In the Black-Scholes-Merton model, as well as in more general stochastic models in finance, the price of an American option solves a system of partial differential variational inequalities. When these inequalities are discretized, one obtains a linear complementar- ity problem that must be solved at each time step. This paper presents an algorithm for the solution of these types of linear complementarity problems that is signicantly faster than the methods currently used in practice. The new algorithm is a two-phase method that combines the active-set identification properties of the projected Gauss-Seidel (or SOR) iteration with the second-order acceleration of a (recursive) reduced-space phase. We show how to design the algorithm so that it exploits the structure of the linear com- plementarity problems arising in these financial models and present numerical results that show the effectiveness of our approach.
Download (pdf)