PDE formulation of some SABR/LIBOR market models and its numerical solution with a sparse grid combination technique

SABR models have been used to incorporate stochastic volatility to LIBOR market models (LMM) in order to describe interest rate dynamics and price interest rate derivatives. From the numerical point of view, the pricing of derivatives with SABR/LIBOR market models (SABR/LMMs) is mainly carried out with Monte Carlo simulation. However, this approach could involve excessively long computational times. For first time in the literature, in the present paper we propose an alternative pricing based on partial differential equations (PDEs). Thus, we pose original PDE formulations associated to the SABR/LMMs proposed by Hagan and Lesniewsk (2008), Mercurio and Morini (2009) and Rebonato and White (2008). Moreover, as the PDEs associated to these SABR/LMMs are high dimensional in space, traditional full grid methods (like standard finite differences or finite elements) are not able to price derivatives over more than three or four underlying interest rates. In order to overcome this curse of dimensionality, a sparse grid combination technique is proposed. A comparison between Monte Carlo simulation results and the ones obtained with the sparse grid technique illustrates the performance of the method.