In addition to a CModel plugin for galaxy photometry (DM-240), we should create a plugin to do sampling-based galaxy fitting, using the optimizer as a starting point and the existing AdaptiveImportanceSampler class to do most of the work.

A major blocker for this is the fact that the SourceTable/SourceRecord don't currently provide any way to save multiple samples per object.

This issue may get worked on before it falls into an official sprint, as it's something I want to on my 20% time.