Skip to content

fix_lenspyx_result allows 1e-5 tolerance on lenspyx and pixell pixels #318

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
msyriac merged 3 commits into from
Jan 16, 2026
+2 −2
Merged

fix_lenspyx_result allows 1e-5 tolerance on lenspyx and pixell pixels #318

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
078a738
fix_lenspyx_result allows 1e-5 tolerance on lenspyx and pixell pixels
zatkins2 Nov 11, 2025
4efff9d
underscore before _fix_lenspyx_result
zatkins2 Nov 11, 2025
64a0958
Merge branch 'master' into quick_newlens_fix
zatkins2 Dec 25, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions pixell/lensing.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -189,7 +189,7 @@ def _fix_lenspyx_result(lenspyx_result, lenspyx_geom_info, shape, wcs):
# in the full-sky geometry
_, _phi0_ind = enmap.sky2pix(fs_shape, fs_wcs, [0, phi0[0]])
phi0_ind = np.round(_phi0_ind).astype(int)
assert phi0_ind == _phi0_ind, \
assert np.allclose(phi0_ind, _phi0_ind, rtol=0, atol=1e-5), \
('we cannot handle the case of a non-integer pixel with cut and paste '
'but could roll the whole array by a fractional pixel using ffts, this '
'needs to be implemented')
Expand Down Expand Up @@ -332,7 +332,7 @@ def _lens_map_curved_lenspyx(shape, wcs, phi_alm, cmb_alm, phi_ainfo=None,
cmb_obs = lenspyx.alm2lenmap(cmb_alm, d_alm, geom_info,
epsilon=epsilon, verbose=verbose,
nthreads=nthreads, pol=pol)
cmb_obs = fix_lenspyx_result(cmb_obs, geom_info, oshape, wcs)
cmb_obs = _fix_lenspyx_result(cmb_obs, geom_info, oshape, wcs)
cmb_obs = cmb_obs.astype(dtype=dtype, copy=False)

# possibly get extra outputs
Expand Down
Loading