diff --git a/sim/pbc.py b/sim/pbc.py index dd23bca5212d74a41096f12ea8720021d07fa7c4..1f066be079ec89809b1ec3c8f4eb370823f6506c 100644 --- a/sim/pbc.py +++ b/sim/pbc.py @@ -89,14 +89,13 @@ class SetupPBC: npbc.set(0) for d in range(0, ndims): for i in For(sim, 0, nlocal + npbc): - #last_id = nlocal + npbc + last_id = nlocal + npbc # TODO: VecFilter? for _ in Filter(sim, positions[i][d] < grid.min(d) + cutneigh): for capacity_exceeded in Branch(sim, npbc >= pbc_capacity): if capacity_exceeded: resize.set(Select(sim, resize > npbc, resize + 1, npbc)) else: - last_id = nlocal + npbc pbc_map[npbc].set(i) pbc_mult[npbc][d].set(1) positions[last_id][d].set(positions[i][d] + grid.length(d)) @@ -112,7 +111,6 @@ class SetupPBC: if capacity_exceeded: resize.set(Select(sim, resize > npbc, resize + 1, npbc)) else: - last_id = nlocal + npbc pbc_map[npbc].set(i) pbc_mult[npbc][d].set(-1) positions[last_id][d].set(positions[i][d] - grid.length(d))