Mike Jarvis

Interesting. Definitely worth trying, but I think you're right that it will be more applicable to space PSFs without the atmospheric component. I've been focusing on ground PSFs for now,...

@jmeyers314 @aaronroodman Would either of you like to take a look at this before merging? I'll be doing Convolution next if you're rather wait for that as a more realistic...

I'm game. Thoughts on what to use for the truth? Maybe a GalSim atmospheric screen and then use GP for the interpolation of Kolmogorov for the primary component. Would that...

Hi Josh, I finally got back to this after 8 months of letting your very good suggestion sit unaddressed. I added a test using Atmospheric + Optical screens for the...

Honestly, I think this is a Python bug. Python thinks that something using /= requires a setter, but numpy has overridden that to work in place. You can check that...

Pull request welcome. ;-)

`image += 1` That seems correct, no? Was this surprising behavior from your perspective? Looks to me like that is working as intended.

It looks like they may have escaped documentation, but here is the [code](https://github.com/GalSim-developers/GalSim/blob/releases/2.5/galsim/image.py#L1603) and [tests](https://github.com/GalSim-developers/GalSim/blob/releases/2.5/tests/test_image.py#L1379-L2027) [tests](https://github.com/GalSim-developers/GalSim/blob/releases/2.5/tests/test_image.py#L2755-L3242) of this functionality.

Calling numpy functions on images directly doesn't seem like an appropriate API. But arithmetic seems quite intuitive to me. YMMV.

Sorry, I didn't follow this. What do you think is being computed wrong? The point of the whole rmin, a, b thing is just so that the right number of...