GalSim-developers / GalSim

The modular galaxy image simulation toolkit. Documentation:
http://galsim-developers.github.io/GalSim/
Other
224 stars 105 forks source link

drawFFT image type #1203

Closed jecampagne closed 1 year ago

jecampagne commented 1 year ago

Hello, looking at def drawFFT_finish(self, image, kimage, wrap_size, add_to_image) I was wandering if in this line

real_image = Image(breal, dtype=float)

the dtype would be more the "image dtype". I know that np.float32 is the default but... My 1cp.

rmjarvis commented 1 year ago

No, that's intentional. We do the intermediate calculation in double precision regardless of the dtype to store the result.