code: def _get_nsymop(self):
          
          if self.centrosymmetric:
              return 2 * len(self._rotations) * len(self._subtrans)
          else:
              return len(self._rotations) * len(self._subtrans)
docs: Returns total number of symmetry operations.