Low-res PSF gaussian fits

[bd-j]

It might be useful to generate multi-gaussian fits to versions of the PSF smoothed by different gaussians. These fits might require fewer components, and can be used when convolving the galaxy gaussian of the same radius, thus cutting down on the total number of gaussians required.

The relevant code is in forcepho.mixtures.fit_f090w though additional code for smoothing the PSF images is necessary. The current set of gaussian radii being used for the galaxies is

minrad, maxrad, dlnr = 0.20, 28.0, np.log(2)
lnradii = np.arange(np.log(minrad), np.log(maxrad), dlnr)

[bd-j]

It would also be good to have PSF fits that allow for negative gaussian components. The current EM-based PSF fitting may not be the best for this.

[bd-j]

This will be useful in more general schemes of paring down the number of gaussians by not convolving every galaxy gaussian with every PSF gaussian.