Store neighbour indices in crystal_t
First step towards cluster splitting
And another step to splitting...
Implemented cluster splitting, seems to work