Float rounding I don't like you.

smoia · Sep 8, 2023 · 8f541b4 · 8f541b4
1 parent fa58da2
commit 8f541b4
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/phys2cvr/phys2cvr.py b/phys2cvr/phys2cvr.py
@@ -759,6 +759,7 @@ def phys2cvr(
         )
 
     LGR.info("phys2cvr finished! Enjoy your outputs!")
+    LGR.warning("Due to float rounding, you might need to mask your output.")
 
 
 def _main(argv=None):

diff --git a/phys2cvr/stats.py b/phys2cvr/stats.py
@@ -388,9 +388,9 @@ def regression(
         pass
 
     # Assign betas, Rsquare and tstats to new volume
-    bout = np.zeros_like(mask, dtype=np.float64)
-    tout = np.zeros_like(mask, dtype=np.float64)
-    rout = np.zeros_like(mask, dtype=np.float64)
+    bout = np.zeros(mask.shape)
+    tout = np.zeros(mask.shape)
+    rout = np.zeros(mask.shape)
     bout[mask] = betas[-1, :]
     tout[mask] = tstats[-1, :]
     rout[mask] = r_square
-Original file line number
+Diff line change
@@ Expand Up / @@ -759,6 +759,7 @@ def phys2cvr( @@
             )
         LGR.info("phys2cvr finished! Enjoy your outputs!")
+        LGR.warning("Due to float rounding, you might need to mask your output.")
     def _main(argv=None):
@@ Expand Down @@