reduce chunk size for mesh point calculations

mesoscope · Sep 14, 2023 · 4d6787a · 4d6787a
1 parent 91f4f17
commit 4d6787a
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/cellpack/autopack/MeshStore.py b/cellpack/autopack/MeshStore.py
@@ -4,7 +4,7 @@
 import trimesh
 from cellpack import autopack
 
-CHUNK_SIZE = 100000
+CHUNK_SIZE = 50000
 
 
 class MeshStore:
@@ -261,6 +261,7 @@ def contains_point(self, geomname, point):
 
     def contains_points_mesh(self, geomname, points):
         mesh = self.get_object(geomname)
+        inside = numpy.full(len(points), False)
         if mesh is not None:
             if len(points) <= CHUNK_SIZE:
                 return mesh.contains(points)  # TODO: check for memory leak
@@ -273,7 +274,7 @@ def contains_points_mesh(self, geomname, points):
                     else:
                         inside = numpy.append(inside, mesh.contains(chunk))
                 return inside
-        return numpy.full(len(points), False)
+        return inside
 
     def get_smallest_radius(self, geomname, center):
         mesh = self.get_object(geomname)