Merge branch 'main' into update_formatter

Author: efaulhaber

Project.toml

@@ -1,7 +1,7 @@
 name = "PointNeighbors"
 uuid = "1c4d5385-0a27-49de-8e2c-43b175c8985c"
 authors = ["Erik Faulhaber <erik.faulhaber@uni-koeln.de>"]
-version = "0.6.0"
+version = "0.6.1"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"

src/cell_lists/full_grid.jl

@@ -222,7 +222,9 @@ function max_points_per_cell(cells)
     return 100
 end
 
-@inline function check_cell_bounds(cell_list::FullGridCellList, cell::Tuple)
+@inline function check_cell_bounds(cell_list::FullGridCellList{<:DynamicVectorOfVectors{<:Any,
+                                                                                        <:Array}},
+                                   cell::Tuple)
     (; linear_indices) = cell_list
 
     # Make sure that points are not added to the outer padding layer, which is needed
@@ -235,10 +237,10 @@ end
     end
 end
 
-# On GPUs, we can't throw a proper error message because string interpolation is not allowed
-@inline function check_cell_bounds(cell_list::FullGridCellList{<:DynamicVectorOfVectors{<:Any,
-                                                                                        <:AbstractGPUArray}},
-                                   cell::Tuple)
+# On GPUs, we can't throw a proper error message because string interpolation is not
+# allowed. Note that we cannot dispatch on `AbstractGPUArray`, as we are inside a kernel,
+# so the array types are something like `CuDeviceArray`, which is not an `AbstractGPUArray`.
+@inline function check_cell_bounds(cell_list::FullGridCellList, cell::Tuple)
     (; linear_indices) = cell_list
 
     # Make sure that points are not added to the outer padding layer, which is needed