Frederik Hennig
--- a/src/pystencils/gpu/periodicity.py 0 → 100644
+++ b/src/pystencils/gpu/periodicity.py 0 → 100644
+import numpy as np
+from itertools import product
+
+from pystencils import CreateKernelConfig, create_kernel
+from pystencils import Assignment, Field
+from pystencils.enums import Target
+from pystencils.slicing import get_periodic_boundary_src_dst_slices, normalize_slice
+
+
+def create_copy_kernel(domain_size, from_slice, to_slice, index_dimensions=0, index_dim_shape=1, dtype=np.float64):
-def create_copy_kernel(domain_size, from_slice, to_slice, index_dimensions=0, index_dim_shape=1, dtype=np.float64):
+def create_copy_kernel(domain_size, src_slice, dst_slice, index_dimensions=0, index_dim_shape=1, dtype=np.float64):
--- a/src/pystencils/gpu/periodicity.py
+++ b/src/pystencils/gpu/periodicity.py
+    src_slice,
+    dst_slice,
+    index_dimensions=0,
+    index_dim_shape=1,
+    dtype=np.float64,
+):
    """Copies a rectangular part of an array to another non-overlapping part"""

-    f = Field.create_generic("pdfs", len(domain_size), index_dimensions=index_dimensions, dtype=dtype)
-    normalized_from_slice = normalize_slice(from_slice, f.spatial_shape)
-    normalized_to_slice = normalize_slice(to_slice, f.spatial_shape)
+    field = Field.create_generic(
+        "field", len(domain_size), index_dimensions=index_dimensions, dtype=dtype
+    )
+    normalized_from_slice = normalize_slice(src_slice, field.spatial_shape)
+    normalized_to_slice = normalize_slice(dst_slice, field.spatial_shape)
-    normalized_from_slice = normalize_slice(src_slice, field.spatial_shape)
-    normalized_to_slice = normalize_slice(dst_slice, field.spatial_shape)
+    normalized_src_slice = normalize_slice(src_slice, field.spatial_shape)
+    normalized_dst_slice = normalize_slice(dst_slice, field.spatial_shape)
--- a/src/pystencils/gpu/periodicity.py
+++ b/src/pystencils/gpu/periodicity.py
+    field = Field.create_generic(
+        "field", len(domain_size), index_dimensions=index_dimensions, dtype=dtype
+    )
+    normalized_from_slice = normalize_slice(src_slice, field.spatial_shape)
+    normalized_to_slice = normalize_slice(dst_slice, field.spatial_shape)

-    offset = [s1.start - s2.start for s1, s2 in zip(normalized_from_slice, normalized_to_slice)]
-    assert offset == [s1.stop - s2.stop for s1, s2 in zip(normalized_from_slice, normalized_to_slice)], \
-        "Slices have to have same size"
+    offset = [
+        s1.start - s2.start
+        for s1, s2 in zip(normalized_from_slice, normalized_to_slice)
+    ]
+    assert offset == [
+        s1.stop - s2.stop for s1, s2 in zip(normalized_from_slice, normalized_to_slice)
+    ], "Slices have to have same size"
-    offset = [
-        s1.start - s2.start
-        for s1, s2 in zip(normalized_from_slice, normalized_to_slice)
-    ]
-    assert offset == [
-        s1.stop - s2.stop for s1, s2 in zip(normalized_from_slice, normalized_to_slice)
-    ], "Slices have to have same size"
+    offset = [
+        s1.start - s2.start
+        for s1, s2 in zip(normalized_src_slice, normalized_dst_slice)
+    ]
+    assert offset == [
+        s1.stop - s2.stop for s1, s2 in zip(normalized_src_slice, normalized_dst_slice)
+    ], "Slices have to have same size"