Radial basis function partition of unity operator splitting method for pricing multi-asset American options

The operator splitting method in combination with finite differences has been shown to be an efficient approach for pricing American options numerically. Here, the operator splitting formulation is extended to the radial basis function partition of unity method. An approach that has previously often been used together with radial basis function methods to deal with the free boundary arising in American option pricing is to solve a penalised version of the Black–Scholes equation. It is shown that the operator splitting technique outperforms the penalty approach when used with the radial basis function partition of unity method. Numerical experiments are performed for one, two and three underlying assets. The advantage of the operator splitting technique grows with the number of dimensions.