Fast and exact projected convolution for non-equidistant grids

Usually, the fast evaluation of a convolution integral $$\int_{{\mathbb{R}}}f(y)g(x-y)\mathrm{d}y$$ requires that the functions f, g are discretised on an equidistant grid in order to apply the fast Fourier transform. Here we discuss the efficient performance of the convolution in locally refined grids. More precisely, the convolution result is projected into some given locally refined grid (Galerkin approximation). Under certain conditions, the overall costs are still $${\mathcal{O}}(N\log N),$$ where N is the sum of the dimensions of the subspaces containing f, g and the resulting function.