We first develop an efficient algorithm to compute Deltas of interest rate derivatives for a number of standard market models. The computational complexity of the algorithms is shown to be proportional to the number of rates times the number of factors per step. We then show how to extend the...

We incorporate a simple and effective control-variate into Fourier inversion formulas for vanilla option prices. The control-variate used in this paper is the Black-Scholes formula whose volatility parameter is determined in a generic non-arbitrary fashion. We analyze contour dependence both in...

We present a new non-nested approach to computing additive upper bounds for callable derivatives using Monte Carlo simulation. It relies on the regression of Greeks computed using adjoint methods. We also show that it is is possible to early terminate paths once points of optimal exercise have...

We introduce a new approach to computing sensitivities of discontinuous integrals. The methodology is generic in that it only requires knowledge of the simulation scheme and the location of the integrand's singularities. The methodology is proven to be optimal in terms of minimizing the variance...

In this paper, we present an efficient approach to compute the first and the second order price sensitivities in the Heston model using the algorithmic differentiation approach. Issues related to the applicability of the pathwise method are discussed in this paper as most existing numerical...

Credit value adjustment (CVA) and related charges have emerged as important risk factors following the Global Financial Crisis. These charges depend on uncertain future values of underlying products, and are usually computed by Monte Carlo simulation. For products that cannot be valued...

We introduce a new class of numerical schemes for discretizing processes driven by Brownian motions. These allow the rapid computation of sensitivities of discontinuous integrals using pathwise methods even when the underlying densities post-discretization are singular. The two new methods...

We develop an efficient algorithm to implement the adjoint method that computes sensitivities of an interest rate derivative (IRD) with respect to different underlying rates in the co-terminal swap-rate market model. The order of computation per step of the new method is shown to be proportional...

We develop new Monte Carlo techniques based on stratifying the stock's hitting-times to the barrier for the pricing and Delta calculations of discretely-monitored barrier options using the Black-Scholes model. We include a new algorithm for sampling an Inverse Gaussian random variable such that...