public static Complex[,] CalcPaddedFourierConvolution(double[,] psf, int yPadding, int xPadding)
{
var psfPadded = Residuals.Pad(psf, yPadding, xPadding);
FFT.Shift(psfPadded);
var PSFPadded = FFT.Forward(psfPadded, 1.0);
return(PSFPadded);
}