This paper is devoted to studying the numerical performance of a power penalty method for a linear parabolic complementarity problem arising from American option valuation. The penalized problem is a nonlinear parabolic partial differential equation (PDE). A fitted finite volume method and an implicit time-stepping scheme are used for respectively the spatial and time discretizations of the PDE. The rate of convergence of the penalty methods with respect to the penalty parameters is investigated both theoretically and numerically. The numerical robustness and computational effectiveness of the penalty method with respect to the market parameters are also studied and compared with those from an existing popular method, PSOR.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.