u/jbosch/DM-641

This branch replaces the never-fully-developed code for fitting multi-shapelet models to PSF model realizations (which is necessary to convolve galaxies in shapelet space).