From 05903d266139c38e48deb24ac1e08837a7c9e0a8 Mon Sep 17 00:00:00 2001 From: Bryan Rumsey Date: Fri, 16 Dec 2022 13:06:33 -0500 Subject: [PATCH 1/2] Fixed issue with applying remove actions. --- spatialpy/core/domain.py | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/spatialpy/core/domain.py b/spatialpy/core/domain.py index c9062a48..b24267e7 100644 --- a/spatialpy/core/domain.py +++ b/spatialpy/core/domain.py @@ -388,8 +388,9 @@ def apply_remove_action(self, action): # remove the particles that fall within the defined region on_boundary = self.find_boundary_points(update=True) - for v_ndx in range(self.get_num_voxels()): - if action['geometry'].inside(self.coordinates()[v_ndx, :], on_boundary[v_ndx]): + for v_ndx, vertex in enumerate(self.vertices): + if action['geometry'].inside(vertex, on_boundary[v_ndx]): + self.vertices = numpy.delete(self.vertices, v_ndx) self.type_id = numpy.delete(self.type_id, v_ndx, 0) self.vol = numpy.delete(self.vol, v_ndx) self.mass = numpy.delete(self.mass, v_ndx) From d8bdc8a18581b218937c64429a1d87f137c6b99b Mon Sep 17 00:00:00 2001 From: Bryan Rumsey Date: Fri, 16 Dec 2022 13:27:15 -0500 Subject: [PATCH 2/2] Added axis to fix. --- spatialpy/core/domain.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/spatialpy/core/domain.py b/spatialpy/core/domain.py index b24267e7..81ea1f27 100644 --- a/spatialpy/core/domain.py +++ b/spatialpy/core/domain.py @@ -390,7 +390,7 @@ def apply_remove_action(self, action): on_boundary = self.find_boundary_points(update=True) for v_ndx, vertex in enumerate(self.vertices): if action['geometry'].inside(vertex, on_boundary[v_ndx]): - self.vertices = numpy.delete(self.vertices, v_ndx) + self.vertices = numpy.delete(self.vertices, v_ndx, 0) self.type_id = numpy.delete(self.type_id, v_ndx, 0) self.vol = numpy.delete(self.vol, v_ndx) self.mass = numpy.delete(self.mass, v_ndx)