PyCosmoMMF 0.0.7__tar.gz → 0.0.8__tar.gz

{PyCosmoMMF-0.0.7 → PyCosmoMMF-0.0.8}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PyCosmoMMF
-Version: 0.0.7
+Version: 0.0.8
 Summary: A package for detecting structures in the Cosmic Web.
 Home-page: https://github.com/James11222/PyCosmoMMF/
 Author: James Sunseri

{PyCosmoMMF-0.0.7 → PyCosmoMMF-0.0.8}/PyCosmoMMF/__init__.py

@@ -10,7 +10,14 @@ import pkg_resources
 
 test_field = np.load(pkg_resources.resource_filename(__name__, 'test_density_cube.npy'))
 
-__version__ = '0.0.5'
+def version():
+    """
+    Print the version of the package.
+    """
+    __version__ = '0.0.8'
+
+    print("PyCosmoMMF version: ", __version__)
+
 
 
 

{PyCosmoMMF-0.0.7 → PyCosmoMMF-0.0.8}/PyCosmoMMF/signatures.py

@@ -59,20 +59,25 @@ def signatures_from_hessian(hessian):
                 ])
                 
                 e1, e2, e3 = np.sort(np.real(np.linalg.eigvals(hes_slice)))
-            
-                sigs[i,j,k,0] = (np.abs(e3 / e1) *
-                    np.abs(e3)) * (θ(-e1) * θ(-e2) * θ(-e3))
-                
-                sigs[i,j,k,1] = (np.abs(e2 / e1) *
-                    xθ(1-np.abs(e3 / e1))) * (np.abs(e2) * θ(-e1) * θ(-e2) )
-                
-                sigs[i,j,k,2] = (xθ(1-np.abs(e2 / e1)) *
-                    xθ(1-np.abs(e3 / e1))) * (np.abs(e1) * θ(-e1))
+
+                if e1 == 0:
+                    sigs[i,j,k,0] = 0
+                    sigs[i,j,k,1] = 0
+                    sigs[i,j,k,2] = 0
+                else:
+                    sigs[i,j,k,0] = (np.abs(e3 / e1) *
+                        np.abs(e3)) * (θ(-e1) * θ(-e2) * θ(-e3))
+                    
+                    sigs[i,j,k,1] = (np.abs(e2 / e1) *
+                        xθ(1-np.abs(e3 / e1))) * (np.abs(e2) * θ(-e1) * θ(-e2) )
+                    
+                    sigs[i,j,k,2] = (xθ(1-np.abs(e2 / e1)) *
+                        xθ(1-np.abs(e3 / e1))) * (np.abs(e1) * θ(-e1))
     
     return sigs
 
 
-def maximum_signature(Rs, field, alg='NEXUSPLUS'):
+def maximum_signature(Rs, field, alg='NEXUSPLUS', eps=1e-8):
     """
     Compute the maximum signatures across all scales Rs.
 
@@ -85,6 +90,8 @@ def maximum_signature(Rs, field, alg='NEXUSPLUS'):
     alg : string, optional
         The algorithm to use. Can be either 'NEXUS' or 'NEXUSPLUS'.
         Default is 'NEXUSPLUS'.
+    eps : float, optional
+        A small number to avoid division by zero errors and whatnot. Default is 1e-8.
 
     Returns
     -------
@@ -95,7 +102,6 @@ def maximum_signature(Rs, field, alg='NEXUSPLUS'):
         raise ValueError("alg must be either 'NEXUS' or 'NEXUSPLUS'")
     
     nx, ny, nz = field.shape
-    eps = 1e-5
     
     # Make sure the field has no 0 values
     field = field + eps

{PyCosmoMMF-0.0.7 → PyCosmoMMF-0.0.8}/PyCosmoMMF.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PyCosmoMMF
-Version: 0.0.7
+Version: 0.0.8
 Summary: A package for detecting structures in the Cosmic Web.
 Home-page: https://github.com/James11222/PyCosmoMMF/
 Author: James Sunseri

{PyCosmoMMF-0.0.7 → PyCosmoMMF-0.0.8}/setup.py

@@ -12,7 +12,7 @@ here = pathlib.Path(__file__).parent.resolve()
 
 setup(
     name='PyCosmoMMF',  # Required
-    version='0.0.7',  # Required
+    version='0.0.8',  # Required
     description='A package for detecting structures in the Cosmic Web.',  # Optional
     long_description=(here / 'README.md').read_text(encoding='utf-8'),  # Optional
     long_description_content_type='text/markdown',  # Optional (see note above)