We present a new method for sampling the Lévy area for a two-dimensional Wiener process conditioned on its endpoints. An efficient sampler for the Lévy area is required to implement a strong Milstein numerical scheme to approximate the solution of a stochastic differential equation driven by a...