Previously it was never explicitly stated that the wavefunctions
returned by modes() are only defined over the unit cell of the
infinite system. This was confusing, as syst.sites contains
not only the unit cell sites, but also the connected sites in
neighboring cells.

Closes #178.