so I just started using pymatgen and created a ‘Structure’-instance using the .from_spacegroup method. I put a species of atoms on a Wyckoff position and the atoms on symmetry-equivalent sites were created.
I also used the get_all_neighbors method to generate all the bonds up to a certain distance, which I find super convenient. However, is there a way to find all ‘symmetry-equivalent’ neighbors or bonds? As far as I can tell the code can only distinguish them by atom and distance. I guess I could search by brute-force, by iterating over all the bonds, applying all operations of the space group and then comparing and collecting them, but if this was already implemented somewhere, it would be very helpful to me.