Non-rigid image registration is a key technology in medical image computing. It is about establishing spatial correspondence between two or more images. The registration problem is ill-posed. To ensure well-posedness a regularizer is introduced. In the present work diffusion regularisation is considered. The fundamental contribution is the introduction of a novel numerical scheme that is fast, efficient, accurate and has a low implementation complexity.