X-Git-Url: https://svn.cri.ensmp.fr/git/linpy.git/blobdiff_plain/a8257bb17a1098f69625387a467170cac4b9f483..ce8dda550018a59520e2d632f10d3fe5650e9d0d:/examples/tesseract.py?ds=sidebyside diff --git a/examples/tesseract.py b/examples/tesseract.py index 383d7bd..d692e37 100755 --- a/examples/tesseract.py +++ b/examples/tesseract.py @@ -5,19 +5,22 @@ # the cube consists of 6 square faces, the hypersurface of the tesseract # consists of 8 cubical cells. -from linpy import * +from linpy import Le, symbols x, y, z, t = symbols('x y z t') tesseract = Le(0, x, 1) & Le(0, y, 1) & Le(0, z, 1) & Le(0, t, 1) + def faces(polyhedron): for points in polyhedron.faces(): face = points[0].aspolyhedron() - face = face.convex_union(*[point.aspolyhedron() for point in points[1:]]) + face = face.convex_union(*[point.aspolyhedron() + for point in points[1:]]) yield face + if __name__ == '__main__': print('Faces of tesseract\n\n {}\n\nare:\n'.format(tesseract)) for face in faces(tesseract):