amaurea / enlib

5 stars 12 forks source link

Bug? lensing.lens_map seems to not lens correctly if cdelt[0] is negative #30

Closed msyriac closed 6 years ago

msyriac commented 6 years ago

If you specify a geometry like this

box = [[y0,x0],[y1,x1]]

such that x0>x1 (so that cdelt[0] < 0)

then lens_map will lens a map with such a geometry incorrectly. It will introduce the right lensed 2-pt power, but the resulting lensed map will not correlate correctly with the input potential (the cross of the reconstruction with the input is mostly zero). A similar result probably holds along the y direction, but I haven't tested that.

I'll try to collect more details here.

msyriac commented 6 years ago

Fixed with Sigurd's new extent patches.