These results differ in a few cases from the ones I sent you earlier. I reran several of the tests to be sure that they were correct

I also have done the iterations count you asked for in several cases, though I did not try to include it here.

The DM-4368.odt file contains the latest results.

Review complete. I have nothing new to add, but I'll paste some of my statements from our previously off-line conversation here for posterity:

It's interesting that the largest PSF images caused that many more failures (suggesting that it takes more iterations to fit more pixels on average), but looking at the other sizes I don't see a clear trend in terms of PSF image size vs. failure rate.

I'm actually encouraged to see the the CModel speed with Full is only ~8x slower than with SingleGaussian; I expected it to be much worse than that, actually.

Mostly, it's clear that I need to do some work on speeding up ShapeletPsfApprox to make it usable in practice, even for simple models. I'm pretty confident I can do that if I can just get some time to work on it. (I'm much less optimistic about being able to speed up CModel, which is why I've been paying more attention to how it scales with the PSF approximation complexity).

I've added a couple of Full and DoubleShapelet model logs which have printouts of time and number of iterations from the size of the history recorder catalog. All but 2 are outer iterations.

Verbose logs, including number of iterations