Just started pasting convolution code in here. It needs to be generalized, etc.
The existing code for convolution depends upon having a PSF kernel. I have renamed this function accordingly and will next insert some more general code.