In the numerical evaluation of geodetic convolution integrals, whether by quadrature or discrete/fast Fourier transform (D/FFT) techniques, the integration kernel is sometimes computed at the centre of the discretised grid cells. For singular kernels a common case in physical geodesy this approximation produces significant errors near the computation point, where the kernel changes rapidly acro...