use spatially-varying PSFs

[lspitler]

Modify convolution so it does a direct convolution as the Fourier convolution won't work too well. Have .a hybrid convolution: at a certain radius the Fourier does the PSF wings, while the core is done via direct convolution.