Generates the vertices of a unit hexacosichoron

hexacosichoron.py

import itertools
from sympy.combinatorics import Permutation

def get_hexacosichoron_vertices():
	phi = (1 + 5 ** .5) / 2
	U = [.5 * phi, .5, .5 / phi]

	permutation_list = [P for P in itertools.permutations(range(4)) if Permutation(P).is_even]

	for i in range(16):
		yield [.5 if i & 2 ** j == 0 else -.5 for j in range(4)]
			
	for i in range(4):
		for x in [-1, 1]:
			yield [x if i == j else 0 for j in range(4)]

	for i in range(8):
		A = [0] + [U[j] if i & 2 ** j == 0 else -U[j] for j in range(3)]
		for P in permutation_list:
			yield [A[j] for j in P]