In this article, we propose a new numerical approach to high-dimensional partial differential equations (PDEs) arising in the valuation of exotic derivative securities. The proposed method is extended from [22] and uses principal component analysis (PCA) of the underlying process in combination with a Taylor expansion of the value function into solutions to low-dimensional PDEs. The approximation is related to anchored analysis of variance (ANOVA) decompositions and is expected to be accurate whenever the covariance matrix has one or few dominating eigenvalues. A main purpose of the present article is to give a careful analysis of the numerical accuracy and computational complexity compared to state-of-the-art Monte Carlo methods on the example of Bermudan swaptions and Ratchet floors, which are considered difficult benchmark problems. We are able to demonstrate that for problems with medium to high dimensionality and moderate time horizons the presented PDE method delivers results comparable in accuracy to the MC methods considered here in similar or (often significantly) faster runtime.

No abstract

We study an expansion method for high-dimensional parabolic PDEs which constructs accurate approximate solutions by decomposition into solutions to lower-dimensional PDEs, and which is particularly effective if there are a low number of dominant principal components. The focus of the present article is the derivation of sharp error bounds for the constant coefficient case and a first and second order approximation. We give a precise characterisation when these bounds hold for (non-smooth) option pricing applications and provide numerical results demonstrating that the practically observed convergence speed is in agreement with the theoretical predictions.

No abstract

No abstract