cmbstack 0.1.0__tar.gz → 0.1.1__tar.gz

{cmbstack-0.1.0 → cmbstack-0.1.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cmbstack
-Version: 0.1.0
+Version: 0.1.1
 Summary: Stack patches of the CMB temperature sky around local maxima
 Author-email: Isaac Alexis López Paredes <isaacalexis.lopezparedes@students.uniroma2.eu>, Anushka Sanjay Tilekar <anushka.tilekar.23@alumni.ucl.ac.uk>
 License-Expression: MIT
@@ -11,6 +11,7 @@ Classifier: Operating System :: OS Independent
 Requires-Python: >=3.8
 Description-Content-Type: text/markdown
 Requires-Dist: numpy
+Requires-Dist: pandas
 Requires-Dist: healpy
 Requires-Dist: matplotlib
 

cmbstack-0.1.1/cmbstack/__init__.py

@@ -0,0 +1,5 @@
+__version__ = "0.1.1"
+
+from cmbstack.pipeline import StackingPipeline
+
+__all__ = ["StackingPipeline", "maps", "stacking"]

{cmbstack-0.1.0 → cmbstack-0.1.1}/cmbstack/maps.py

@@ -13,6 +13,7 @@ power spectrum (temperature TT, lensing convergence, y-map, ...).
 
 import numpy as np
 import healpy as hp
+import pandas as pd
 
 
 def dl_to_cl(ell, dl,lmax):
@@ -46,12 +47,13 @@ def dl_to_cl(ell, dl,lmax):
 
     return cl
 
-# Function 1
-def load_cl(path):
+
+def load_cl(path, column="TT"):
     """Load a power-spectrum file and return C_ell for the requested spectrum.
 
-    The expected file columns are ell, Dl_TT, Dl_TE, Dl_EE, Dl_BB, Dl_dd, with
-    D_ell in uK^2. The chosen column is converted from D_ell to C_ell via
+    The expected file columns are ell, Dl_TT, Dl_TE, Dl_EE, Dl_BB, Dl_pp, with
+    D_ell in uK^2 (with names L, TT, TE, EE, BB, PP).
+    The chosen column is converted from D_ell to C_ell via
     :func:`dl_to_cl` before being returned.
 
     Parameters
@@ -59,7 +61,7 @@ def load_cl(path):
     path : str
         Path to the whitespace-delimited spectrum file.
     column : str, optional
-        Which spectrum to return: one of "TT", "TE", "EE", "BB", "dd".
+        Which spectrum to return: one of "TT", "TE", "EE", "BB", "PP".
         Default "TT".
 
     Returns
@@ -68,14 +70,22 @@ def load_cl(path):
         C_ell array indexed from ell=0, suitable for passing to
         :func:`simulate_map`.
     """
-    ell, dl = np.loadtxt(path,usecols=(0,1),unpack=True)
-    ell = ell.astype(int)
-    lmax = ell.max()  
-    cl = dl_to_cl(ell,dl,lmax)
+
+    df = pd.read_csv(
+        path,
+        sep=r"\s+",
+        comment="#",
+        header=None,
+        names=["L", "TT", "TE", "EE", "BB", "PP"],
+    )
+    ell = df['L'].values.astype(int)
+    dl = df[column].values
+    lmax = ell.max()
+    cl = dl_to_cl(ell, dl, lmax)
 
     return cl
 
-# Function 2
+
 def simulate_map(cl, nside=128, seed=None):
     """Simulate a Gaussian random HEALPix map from a power spectrum.
 

{cmbstack-0.1.0 → cmbstack-0.1.1}/cmbstack.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: cmbstack
-Version: 0.1.0
+Version: 0.1.1
 Summary: Stack patches of the CMB temperature sky around local maxima
 Author-email: Isaac Alexis López Paredes <isaacalexis.lopezparedes@students.uniroma2.eu>, Anushka Sanjay Tilekar <anushka.tilekar.23@alumni.ucl.ac.uk>
 License-Expression: MIT
@@ -11,6 +11,7 @@ Classifier: Operating System :: OS Independent
 Requires-Python: >=3.8
 Description-Content-Type: text/markdown
 Requires-Dist: numpy
+Requires-Dist: pandas
 Requires-Dist: healpy
 Requires-Dist: matplotlib
 

{cmbstack-0.1.0 → cmbstack-0.1.1}/cmbstack.egg-info/SOURCES.txt

@@ -10,4 +10,5 @@ cmbstack.egg-info/dependency_links.txt
 cmbstack.egg-info/requires.txt
 cmbstack.egg-info/top_level.txt
 tests/test_find_peaks.py
-tests/test_round.py
+tests/test_round.py
+tests/test_runtime.py

{cmbstack-0.1.0 → cmbstack-0.1.1}/cmbstack.egg-info/requires.txt

@@ -1,3 +1,4 @@
 numpy
+pandas
 healpy
 matplotlib

{cmbstack-0.1.0 → cmbstack-0.1.1}/pyproject.toml

@@ -19,6 +19,7 @@ classifiers = [
 requires-python = ">=3.8"
 dependencies = [
     "numpy",
+    "pandas",
     "healpy",
     "matplotlib",
 ]

cmbstack-0.1.1/tests/test_runtime.py

@@ -0,0 +1,3 @@
+# pStats
+# cProfiler
+

cmbstack-0.1.0/cmbstack/__init__.py

@@ -1,5 +0,0 @@
-__version__ = "0.1.0"
-
-from cmbstack.pipeline import StackingPipeline
-
-__all__ = ["StackingPipeline", "maps", "stacking"]