Mike Jarvis
Mike Jarvis
Sure. We would welcome a pull request with this feature. :) None of the developers AFAIK have any research interest in varying the inner scale, so you would need to...
This would be an excellent addition Michael. We would definitely like this for some of the LSST SRM tasks I think. To be clear, do you have the math for...
Great. The TelescopePSFLSST should probably go in the lsst module (the galsim/lsst/ directory). And you won't be surprised to hear me say that we'll want more unit tests than you...
I knew about the non-thread-safeness of FFT.cpp. Specifically the fftw calls were known to not be thread safe, so the mutex locks there are probably right. I'm curious though what...
I see Yes, each thread would be using the same SericInfo object, so that setup could be racing. In fact, you probably want to move that mutex lock up a...
Hm. There are probably similar things in SersicInfo that also need locks. e.g. stepK() looks like it wouldn't race, but it would be possible for one thread to overwrite R...
I noticed when writing text for Troxel's paper that the IPC coefficients [here](https://wfirst.gsfc.nasa.gov/science/201907/WFIRST_WFI_Reference_Information_20190626.pptx) don't match what we are using. I couldn't actually find the reference for our coefficients, but at...
I think the number that is most relevant to have documented for people to reference is: At what S/N does photon shooting become faster than fft. (Or vice versa depending...
The other independent variables that could affect the relative timing are the pixel size (relative to some size value of the profile) and the image size. Probably should start with...