Numerical computation of Theta in a jump-diffusion model by integration by parts

Using the Malliavin calculus in time inhomogeneous jump-diffusion models, we obtain an expression for the sensitivity Theta of an option price (with respect to maturity) as the expectation of the option payoff multiplied by a stochastic weight. This expression is used to design efficient numerical algorithms that are compared with traditional finite-difference methods for the computation of Theta. Our proof can be viewed as a generalization of Dupire's integration by parts to arbitrary and possibly non-smooth payoff functions. In the time homogeneous case, Theta admits an expression from the Black-Scholes PDE in terms of Delta and Gamma but the representation formula obtained in this way is different from ours. Numerical simulations are presented in order to compare the efficiency of the finite-difference and Malliavin methods.