foscat 2025.7.1__tar.gz → 2025.7.2__tar.gz

{foscat-2025.7.1/src/foscat.egg-info → foscat-2025.7.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: foscat
-Version: 2025.7.1
+Version: 2025.7.2
 Summary: Generate synthetic Healpix or 2D data using Cross Scattering Transform
 Author-email: Jean-Marc DELOUIS <jean.marc.delouis@ifremer.fr>
 Maintainer-email: Theo Foulquier <theo.foulquier@ifremer.fr>

{foscat-2025.7.1 → foscat-2025.7.2}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "foscat"
-version = "2025.07.1"
+version = "2025.07.2"
 description = "Generate synthetic Healpix or 2D data using Cross Scattering Transform"
 readme = "README.md"
 license = { text = "BSD-3-Clause" }

{foscat-2025.7.1 → foscat-2025.7.2}/src/foscat/FoCUS.py

@@ -35,7 +35,7 @@ class FoCUS:
             mpi_rank=0
     ):
 
-        self.__version__ = "2025.07.1"
+        self.__version__ = "2025.07.2"
         # P00 coeff for normalization for scat_cov
         self.TMPFILE_VERSION = TMPFILE_VERSION
         self.P1_dic = None
@@ -1691,9 +1691,10 @@ class FoCUS:
                 except:
                     lcell_ids=self.to_numpy(cell_ids)
                     idx_map[lcell_ids]=np.arange(lcell_ids.shape[0],dtype='int32')
-                    
+
                 lidx=np.where(idx_map[tmp[:,1]%(12*nside**2)]!=-1)[0]
                 orientation=tmp[lidx,1]//(12*nside**2)
+                orientation2=tmp[lidx,0]//(12*nside**2)
                 tmp=tmp[lidx]
                 wr=wr[lidx]
                 wi=wi[lidx]
@@ -1703,16 +1704,21 @@ class FoCUS:
                 wi[lidx]=0.0
                 tmp[lidx,0]=0
                 tmp[:,1]+=orientation*lcell_ids.shape[0]
+                tmp[:,0]+=orientation2*lcell_ids.shape[0]
                 
                 idx_map=-np.ones([12*nside**2],dtype='int32')
                 idx_map[lcell_ids]=np.arange(cell_ids.shape[0],dtype='int32')
-                lidx=np.where(idx_map[tmp2[:,1]]!=-1)[0]
+                lidx=np.where(idx_map[tmp2[:,1]%(12*nside**2)]!=-1)[0]
+                i_id=tmp2[lidx,1]//(12*nside**2)
+                i_id2=tmp2[lidx,0]//(12*nside**2)
                 tmp2=tmp2[lidx]
                 ws=ws[lidx]
-                tmp2=idx_map[tmp2]
+                tmp2=idx_map[tmp2%(12*nside**2)]
                 lidx=np.where(tmp2[:,0]==-1)[0]
                 ws[lidx]=0.0
                 tmp2[lidx,0]=0
+                tmp2[:,1]+=i_id*lcell_ids.shape[0]
+                tmp2[:,0]+=i_id2*lcell_ids.shape[0]
                 
         else:
             tmp = indice

{foscat-2025.7.1 → foscat-2025.7.2}/src/foscat/HealSpline.py

@@ -6,11 +6,14 @@ import healpy as hp
 
 class heal_spline:
     def __init__(
-        self,
-        level):
+            self,
+            level,
+            gamma=1,
+    ):
         nside=2**level
         self.nside_store=2**(level//2)
         self.spline_tree={}
+        self.gamma=gamma
         
         self.nside=nside
         #compute colatitude
@@ -79,6 +82,8 @@ class heal_spline:
             
         hit=np.bincount(all_idx.flatten(),weights=www.flatten())
         www[hit[all_idx]<threshold]=0.0
+        if self.gamma!=1:
+            www=www**self.gamma
         www=www/np.sum(www,0)[None,:]
         return www,all_idx,heal_idx
         

{foscat-2025.7.1 → foscat-2025.7.2/src/foscat.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: foscat
-Version: 2025.7.1
+Version: 2025.7.2
 Summary: Generate synthetic Healpix or 2D data using Cross Scattering Transform
 Author-email: Jean-Marc DELOUIS <jean.marc.delouis@ifremer.fr>
 Maintainer-email: Theo Foulquier <theo.foulquier@ifremer.fr>