From 2a87f9c46ac6248b5661ee3a7492cffb62a448c2 Mon Sep 17 00:00:00 2001
From: = <=>
Date: Tue, 6 Feb 2024 16:41:00 +0100
Subject: [PATCH] tentative first implementation

---
 python/polycube.py | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100644 python/polycube.py

diff --git a/python/polycube.py b/python/polycube.py
new file mode 100644
index 0000000..da23b0f
--- /dev/null
+++ b/python/polycube.py
@@ -0,0 +1,32 @@
+from python.tuple_tools import reorient_tuple, generate_neighbors
+from itertools import combinations
+
+class PolyCube:
+
+    #cubes              set{tuple}
+    #orientations       list{set{tuple{}}}
+    #last_additions     set{tuple}
+
+    def __init__(self, cubes_list, last_additions = []) -> None:
+        self.cubes = set(cubes_list)
+        self._compute_reorientations()
+        self.last_additions = last_additions
+    
+    def _compute_reorientations(self):
+        n_cubes = len(self.cubes)
+        n_orientations = 24
+        orientations = [set()]*n_orientations
+        for i in range(n_orientations):
+            for cube in self.cubes:
+                orientations[i].add(cube)
+        self.orientations = orientations
+    
+    def generate_children(self, max_length):
+        growth_candidates = set([])
+        for last_addition in self.last_additions:
+            growth_candidates = growth_candidates.union(generate_neighbors(*last_addition))
+        growth_candidates = growth_candidates.difference(self.cubes)
+        realizable_growth =  combinations(growth_candidates, max_length)
+        for addition in realizable_growth:
+            yield PolyCube(self.cubes.union(addition), addition)
+