There are currently 2 approaches used to making flux-weighted, psf-convolved velocity maps. One uses galsim to do psf-convolution and the other just generates a pixel matrix and does the convolution in real-space with scipy. Ideally they should agree but the last time I checked there were some problems with the galsim version.
I think we ideally want to keep the functionality of both around for testing, but this can probably be done in a cleaner way.
There are currently 2 approaches used to making flux-weighted, psf-convolved velocity maps. One uses galsim to do psf-convolution and the other just generates a pixel matrix and does the convolution in real-space with scipy. Ideally they should agree but the last time I checked there were some problems with the galsim version.
I think we ideally want to keep the functionality of both around for testing, but this can probably be done in a cleaner way.