Change pypol to linpy
[linpy.git] / examples / tesseract.py
1 #!/usr/bin/env python3
2
3 from pypol import *
4
5 x, y, z, t = symbols('x y z t')
6
7 tesseract = \
8 Le(0, x) & Le(x, 1) & \
9 Le(0, y) & Le(y, 1) & \
10 Le(0, z) & Le(z, 1) & \
11 Le(0, t) & Le(t, 1)
12
13 def faces(polyhedron):
14 for points in polyhedron.faces():
15 face = points[0].aspolyhedron()
16 face = face.union(*[point.aspolyhedron() for point in points[1:]])
17 face = face.aspolyhedron()
18 yield face
19
20 print('Faces of tesseract\n\n {}\n\nare:\n'.format(tesseract))
21 for face in faces(tesseract):
22 assert(len(face.vertices()) == 8)
23 print(' {}'.format(face))