Currently creates PSF stack, then computes more FFTs to do loss in fourier domain. Should change propagation code to keep spectra (as an option).