make sure high-res map indexing is not nonsense

CMB-S4 · Sep 19, 2023 · 07f5175 · 07f5175
1 parent 1b9087f
commit 07f5175
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/maps/tests/healpix_maps.py b/maps/tests/healpix_maps.py
@@ -221,3 +221,9 @@
 lon, lat = maps.azel.convert_radec_to_gal(alpha, delta)
 pix = xeq.angles_to_pixels(lon, lat)
 assert(np.allclose(np.asarray(pix), np.asarray(xeq)))
+
+# high-res
+hx = maps.HealpixSkyMap(16384)
+pix1 = hx.angle_to_pixel(0, -45 * core.G3Units.deg)
+pix2 = hp.ang2pix(16384, 3 * np.pi / 4, 0)
+assert(pix1 == pix2)