cube24.py

import itertools as it
import numpy as np

def cube24(points):
    points = np.asarray(points)
    pp = [(0, 1, 2), (1, 2, 0), (2, 0, 1)]
    sp = [(1, 1, 1), (-1, -1, 1), (-1, 1, -1), (1, -1, -1)]
    pn = [(0, 2, 1), (1, 0, 2), (2, 1, 0)]
    sn = [(-1, 1, 1), (1, -1, 1), (1, 1, -1), (-1, -1, -1)]
    for perm, sign in it.chain(it.product(pp, sp), it.product(pn,sn)):
        yield points[:,perm] * sign