Andy Miller

when sampling fluxes (and holding the shape parameters constant), the spatial extent of stars and galaxies do not need to be recomputed (expected ~ 3x speedup for stars and ~...

Should we track the log-likelihood's computed for each sample? If we let the sampler object own a trace of samples, it can also own a trace of log likelihoods.

This is an ensemble of samplers at different temperatures - each sampler can just be an existing Sampyl sampler. This object would just keep track of the cooling schedule.

we have names of variables and collections of state samples - it'd be nice if we could have a very generic method for summarizing them on the command line (and...

Add a standard test for each sampler to make sure MCMC sample means/variances are approaching the right answer. A mixture of gaussians would be a good fit here.