draft for fast loading galfile

jacobic · May 26, 2021 · 022c98f · 022c98f
1 parent 2c7a88f
commit 022c98f
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/redmapper/galaxy.py b/redmapper/galaxy.py
@@ -37,6 +37,7 @@ def read_cat(path, fn, columns=None):
 
 
 def trim_cat(cat, border, nside, _hpix, nside_cutref=512):
+
     if len(_hpix) == 1 and nside > 0 and border > 0.0:
         # Prepare border cuts if necessary
         boundaries = hp.boundaries(nside, _hpix[0], step=nside_cutref // nside)
@@ -419,7 +420,7 @@ def from_galfile(cls, filename, zredfile=None, nside=0, hpix=[], border=0.0, tru
         # New more memory efficient reading...
         ngals = np.sum(tab.ngals[indices])
 
-        kw_border = dict(border=border, nside=nside, _hpix=hpix, nside_cutref=512)
+        kw_border = dict(border=border, nside=nside, _hpix=_hpix, nside_cutref=512)
 
         # read the files
         if len(_hpix) == 0 or nside == 0: