fixed a bug with my usage of std set intersection. It is UNSAFE to use on pointer values and will fail in combination with periodic boundary conditions due to non sequential neighbors