We present an iterative procedure for computing the optimal Bermudan stopping time, hence the Bermudan Snell envelope. The method produces an increasing sequence of approximations of the Snell envelope from below, which coincide with the Snell envelope after finitely many steps. Then, by duality, the method induces a convergent sequence of upper bounds as well. In a Markovian setting the presented procedure allows to calculate approximative solutions with only a few nestings of conditional expectations and is therefore tailor-made for a plain Monte Carlo implementation. The method may be considered generic for all discrete optimal stopping problems. The power of the procedure is demonstrated for Bermudan swaptions in a full factor LIBOR market model. Copyright Springer-Verlag Berlin/Heidelberg 2006Bermudan options, optimal stopping, Monte Carlo simulation, LIBOR market model,

We present a generic non-nested Monte Carlo procedure for computing true upper bounds for Bermudan products, given an approximation of the Snell envelope. The pleonastic "true" stresses that, by construction, the estimator is biased above the Snell envelope. The key idea is a regression estimator for the Doob martingale part of the approximative Snell envelope, which preserves the martingale property. The so constructed martingale can be employed for computing tight dual upper bounds without nested simulation. In general, this martingale can also be used as a control variate for simulation of conditional expectations. In this context, we develop a variance reduced version of the nested primal-dual estimator (Andersen and Broadie, 2004). Numerical experiments indicate the efficiency of the proposed algorithms.

No abstract

Abstract. We introduce a multiple curve framework that combines tractable dynamics and semi-analytic pricing formulas with positive interest rates and basis spreads. Negatives rates and positive spreads can also be accommodated in this framework. The dynamics of OIS and LIBOR rates are specified following the methodology of the affine LIBOR models and are driven by the wide and flexible class of affine processes. The affine property is preserved under forward measures, which allows us to derive Fourier pricing formulas for caps, swaptions and basis swaptions. A model specification with dependent LIBOR rates is developed, that allows for an efficient and accurate calibration to a system of caplet prices.

Abstract. New classes of stochastic differential equations can now be studied using rough path theory (e.g. Lyons et al. [LCL07] or ). In this paper we investigate, from a numerical analysis point of view, stochastic differential equations driven by Gaussian noise in the aforementioned sense. Our focus lies on numerical implementations, and more specifically on the saving possible via multilevel methods. Our analysis relies on a subtle combination of pathwise estimates, Gaussian concentration, and multilevel ideas. Numerical examples are given which both illustrate and confirm our findings.