I thought it would make sense to have it e.g to fit point sources or test convolution kernels. But at some point
I realized that it is not as easy to implement, because models are evaluated most times on a grid. I wonder if there is a well defined way to handle this? So, I actually want to implement it, but I still need to find out how...
I may ask this on astropy-dev. Thanks for reviewing!