docs: def get_rotations(self):
          """Return all rotations, including inversions for
          centrosymmetric crystals."""
          if self.centrosymmetric:
              return np.vstack((self.rotations, -self.rotations))
          else:
              return self.rotations