Generalize to multi-component shapelet models

[griffinfoster]

The shapelet decomposition works well with Gaussian like sources, but sources with double lobes, or many components tend to no do very well. Instead, the fitter/solver could decompose an image into multiple shapelet models.