Support python 3.14, update NumPy to 2 and Ducc to 0.41

Author: ziotom78

litebird_sim/grasp2alm.py

@@ -57,7 +57,10 @@ def __init__(self, healpix_map):
     def _check_reason(
         self, reason: int, iter_count: int, residual_norm: float, quality: float
     ) -> None:
-        assert reason in (1, 2), (
+        assert reason in (
+            1,
+            2,
+        ), (
             f"pseudo_analysis failed, reason: {REASON_DESCRIPTION[reason]}, {iter_count=}, {residual_norm=}, {quality=}"
         )
 
@@ -95,7 +98,7 @@ def to_alm(
             (3, SphericalHarmonics.num_of_alm_from_lmax(lmax, mmax)),
             dtype=np.complex128,
         )
-        (_, reason, iter_count, residual_norm, quality) = ducc0.sht.pseudo_analysis(
+        _, reason, iter_count, residual_norm, quality = ducc0.sht.pseudo_analysis(
             map=self.map[0].reshape(1, -1),  # Make this a 2D matrix
             alm=alm[0, :].reshape(1, -1),
             lmax=lmax,
@@ -190,7 +193,7 @@ def beam_mapmaker(
         if hit_map[i] > 0:
             output_map[i] /= hit_map[i]
         else:
-            output_map[i] = np.NaN
+            output_map[i] = np.nan
 
 
 @dataclass

pyproject.toml

@@ -30,23 +30,22 @@ authors = [
 ]
 readme = "README.md"
 license = {text = "GPL-3.0"}
-requires-python = ">=3.10,<3.14"
+requires-python = ">=3.10,<=3.14"
 dependencies = [
-    "numba>=0.63.0",
-    "numpy>=1.26.4",
+    "numba>=0.65.0",
+    "numpy>=2.2.6",
     "astropy>=6.1.7",
     "tomlkit>=0.12.1",
     "jinja2>=3.1",
-    "matplotlib>=3.10.3",
+    "matplotlib>=3.10.8",
     "healpy>=1.18.1",
     "pyyaml>=6.0",
     "jplephem>=2.22",
     "requests>=2.32.0",
     "rich>=14.0.0",
-    "ducc0>=0.40.0",
-    "pysm3>=3.4.2",
+    "ducc0>=0.41.0",
+    "pysm3>=3.4.4",
     "h5py>=3.9",
-    "deprecation>=2.1.0",
     "deprecated>=1.2.18",
     "scipy>=1.14.0",
     "libinsdb>=0.9.0",

test/test_binner.py

@@ -54,7 +54,11 @@ def test_accumulate_map_and_info():
     assert np.allclose(res_info, info)
 
     rhs = mapping._extract_map_and_fill_info(info)
-    assert np.allclose(np.linalg.solve(info, rhs), res_map)
+    # Starting from NumPy 2, we must reshape `rhs`
+    # so that it matches the shape of the `info` tensor.
+    # Plus, we add `.squeeze(1)` to remove the (new)
+    # dimension at the end of the result
+    assert np.allclose(np.linalg.solve(info, rhs[..., np.newaxis]).squeeze(-1), res_map)
 
 
 def test_make_binned_map_api_simulation(tmp_path):