Hi again,
> - A more general solution is to do what Anton proposed. I include a > sample script below. > sym.add_site_family(lat.sublattices[0], other_vectors=[(-1, 2)]) > sym.add_site_family(lat.sublattices[1], other_vectors=[(-1, 2)]) Thanks for this. I had tried do do the following: sym.add_site_family(lat, other_vectors=[(-1, 2)]) instead of calling `add_site_family` for each sub-lattice, and it did not work. I therefore thought that my understanding was flawed, but it seems that for a polyatomic lattice one has to manually add each sub-lattice. Maybe a word in the documentation would be useful to prevent this potential pitfall Joe
signature.asc
Description: OpenPGP digital signature