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);
            }