remvoed limitation that value must be > 0

themachinethatgoesping · Jan 3, 2025 · 16f54c0 · 16f54c0
1 parent ef24439
commit 16f54c0
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/meson.build b/meson.build
@@ -9,7 +9,7 @@ project(
   'cpp',
   license: 'MPL-2.0',
 
-  version: '0.2.0',
+  version: '0.2.1',
   default_options: ['warning_level=2', 'buildtype=release', 'cpp_std=c++20'],
   meson_version: '>=1.3.2' #first version with  clang-cl openmp support
 )

diff --git a/python/themachinethatgoesping/gridding/functions/gridfunctions.py b/python/themachinethatgoesping/gridding/functions/gridfunctions.py
@@ -210,7 +210,8 @@ def grd_weighted_mean(
                     continue
 
             # print(ix,iy,iz,v,w)
-            if v >= 0:
+            #if v >= 0:
+            if np.isfinite(v):
                 image_values[ix][iy][iz] += v * w
                 image_weights[ix][iy][iz] += w
 
@@ -276,7 +277,8 @@ def grd_block_mean(
             if abs(iz) >= nz:
                 continue
 
-        if v >= 0:
+        #if v >= 0:
+        if np.isfinite(v):
             image_values[ix][iy][iz] += v
             image_weights[ix][iy][iz] += 1