Efficient Kernel Convolution for Smooth Surfaces without Edge Effects
  One of the most efficient ways to produce unconditional simulations is with the kernel convolution using fast Fourier transform (FFT) [1]. However, when data is located on a surface, this approach is not efficient because data needs to be processed in a three-dimensional enclosing box. This paper describes a novel approach based on integer transformation to reduce the volume of the enclosing box. 
Alexander Gribov
stat.CO