bbf 0.7.1__tar.gz → 0.7.2__tar.gz

{bbf-0.7.1 → bbf-0.7.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: bbf
-Version: 0.7.1
+Version: 0.7.2
 Summary: Fast computation of broadband fluxes and magnitudes
 Keywords: astronomy,astrophysics
 Author-Email: Nicolas Regnault <nicolas.regnault@lpnhe.in2p3.fr>
@@ -28,10 +28,17 @@ Requires-Dist: scikit-sparse<0.5
 Requires-Dist: scipy>=1.14
 Requires-Dist: sncosmo>=2.12
 Requires-Dist: gaiahealpixcache>=0.4.0
+Provides-Extra: hdf5
+Requires-Dist: tables; extra == "hdf5"
 Provides-Extra: test
 Requires-Dist: pytest>=6.0; extra == "test"
 Requires-Dist: pytest-cov; extra == "test"
 Requires-Dist: lemaitre-bandpasses>=0.4; extra == "test"
+Requires-Dist: tables; extra == "test"
+Provides-Extra: docs
+Requires-Dist: sphinx>=7.0; extra == "docs"
+Requires-Dist: myst-parser>=2.0; extra == "docs"
+Requires-Dist: furo; extra == "docs"
 Description-Content-Type: text/markdown
 
 # Broadband fluxes (`bbf`)
@@ -116,6 +123,14 @@ To work on the sources and, e.g. run the tests locally:
   VERBOSE=1 pip install --no-build-isolation -Ceditable.rebuild=true -ve .
   ```
 
+To build the documentation locally (API reference generated from the docstrings, plus
+this README and the `docs/audit` pages):
+
+  ```bash
+  pixi run -e docs docs
+  # open docs/_build/html/index.html
+  ```
+
 
 ## Installing the Lemaitre bandpasses
 

{bbf-0.7.1 → bbf-0.7.2}/README.md

@@ -80,6 +80,14 @@ To work on the sources and, e.g. run the tests locally:
   VERBOSE=1 pip install --no-build-isolation -Ceditable.rebuild=true -ve .
   ```
 
+To build the documentation locally (API reference generated from the docstrings, plus
+this README and the `docs/audit` pages):
+
+  ```bash
+  pixi run -e docs docs
+  # open docs/_build/html/index.html
+  ```
+
 
 ## Installing the Lemaitre bandpasses
 

{bbf-0.7.1 → bbf-0.7.2}/bbf/atmosphere.py

@@ -151,22 +151,36 @@ def get_airmass_extinction(source="palomar"):
         )
         return spec
 
-    if source.lower() == "palomareac4":
+    if source.lower() == "palomareac4median":
         lbda, data = pd.read_csv(
             importlib.resources.files(__package__)
             / "data"
-            / "palomar_extinction_eac4.dat"
+            / "palomar_extinction_eac4_median.dat"
         ).values.T
-        spec = ExtinctionSpectrum(lbda, data, "EAC4 average 2017-2025")
+        spec = ExtinctionSpectrum(lbda, data, "EAC4 median 2017-2025")
         return spec
 
-    if source.lower() == "palomarmerra2":
+    if source.lower() == "palomarmerra2median":
         lbda, data = pd.read_csv(
             importlib.resources.files(__package__)
             / "data"
-            / "palomar_extinction_merra2.dat"
+            / "palomar_extinction_merra2_median.dat"
         ).values.T
-        spec = ExtinctionSpectrum(lbda, data, "MERRA2 average 2017-2025")
+        spec = ExtinctionSpectrum(lbda, data, "MERRA2 median 2017-2025")
+        return spec
+
+    if source.lower() == "palomareac4":
+        lbda, data = pd.read_csv(
+            importlib.resources.files(__package__) / 'data' / 'palomar_extinction_eac4.dat').values.T
+        spec = ExtinctionSpectrum(lbda, data,
+                                  "EAC4 average DR2.5")
+        return spec
+
+    if source.lower() == "palomarmerra2":
+        lbda, data = pd.read_csv(
+            importlib.resources.files(__package__) / 'data' / 'palomar_extinction_merra2.dat').values.T
+        spec = ExtinctionSpectrum(lbda, data,
+                                  "MERRA2 average DR2.5")
         return spec
 
     if source.lower() == "buton":

{bbf-0.7.1 → bbf-0.7.2}/bbf/bspline/projector.py

@@ -30,9 +30,16 @@ class Projector:
 
         self._factor, self._J = self._get_projector(self._grid)
 
-    def _get_projector(self, x):
+    def _get_projector(self, x, beta=1.0e-6):
+        # beta=1e-6: matches every other cholesky_AAt call site in the codebase
+        # (docs/bug_fix_plan.md B-9). Without it, a dataset that does not fully
+        # cover the basis (any basis function untouched by any sample point) makes
+        # the Gram matrix exactly rank-deficient and this raised
+        # sksparse.cholmod.CholmodNotPositiveDefiniteError outright, unlike the
+        # one_by_one/BSpline.linear_fit projection path (also beta=1e-6), which
+        # degrades gracefully instead.
         J = self.basis.eval(x).tocsr()
-        factor = sksparse.cholmod.cholesky_AAt(J.T)
+        factor = sksparse.cholmod.cholesky_AAt(J.T, beta=beta)
         return factor, J
 
     def _compatible_arrays(self, y, x):