>>> rcoord = beryl.get_scaled_positions()
>>> print(sg.unique_sites(rcoord))
[[ 0.33333333  0.66666667  0.25      ]
 [ 0.5         0.          0.25      ]
 [ 0.88        0.27        0.        ]
 [ 0.65        0.15        0.65      ]
 [ 0.24        0.31        0.5       ]]