-
Notifications
You must be signed in to change notification settings - Fork 0
/
Frank_A13inD6.py
101 lines (82 loc) · 6.88 KB
/
Frank_A13inD6.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
import itertools
from operator import mul, add
from functools import reduce
from compare import checker
from compare import order
from compare import dimChecker
from compare import perm
def dim(x):
sol=1
for i in range(0,len(x)):
sol = sol*(x[i]+1)
return(sol)
def isZero(sol):
return reduce(mul, (reduce(add, x, 0) for x in zip(*sol)), 1) == 0
L=[]
d=12
#List of all possible weights
for i in range(0,d-3):
if i%2==0:
L.append([(i,0,0)])
else :
temp=[(i,0,0),(i,0,0)]
tempdim=dimChecker(temp)
if tempdim<=d-4:
L.append(temp)
o=order(3)
for i in range(0,len(L)):
for j in range(0,len(o)):
temp=perm(o[j],L[i])
if temp not in L:
L.append(temp)
Ladd = [[(1,1,0)],[(2,2,0)],[(1,3,0)],[(1,1,2)]]
for i in range(0,len(Ladd)):
for j in range(0,len(o)):
temp=perm(o[j],Ladd[i])
if temp not in Ladd:
Ladd.append(temp)
L= L+Ladd
print(L)
vs = [[(2,0,0)],[(1,1,0)],[(0,1,0),(0,1,0)],[(0,0,3),(0,0,3)]]
vs = [[(0, 0, 0, 0, 0)], [(1, 0, 0, 0, 0), (1, 0, 0, 0, 0)], [(2, 0, 0, 0, 0)], [(3, 0, 0, 0, 0), (3, 0, 0, 0, 0)], [(4, 0, 0, 0, 0)], [(5, 0, 0, 0, 0), (5, 0, 0, 0, 0)], [(6, 0, 0, 0, 0)], [(8, 0, 0, 0, 0)], [(10, 0, 0, 0, 0)], [(0, 1, 0, 0, 0), (0, 1, 0, 0, 0)], [(0, 0, 1, 0, 0), (0, 0, 1, 0, 0)], [(0, 0, 0, 1, 0), (0, 0, 0, 1, 0)], [(0, 0, 0, 0, 1), (0, 0, 0, 0, 1)], [(0, 2, 0, 0, 0)], [(0, 0, 2, 0, 0)], [(0, 0, 0, 2, 0)], [(0, 0, 0, 0, 2)], [(0, 3, 0, 0, 0), (0, 3, 0, 0, 0)], [(0, 0, 3, 0, 0), (0, 0, 3, 0, 0)], [(0, 0, 0, 3, 0), (0, 0, 0, 3, 0)], [(0, 0, 0, 0, 3), (0, 0, 0, 0, 3)], [(0, 4, 0, 0, 0)], [(0, 0, 4, 0, 0)], [(0, 0, 0, 4, 0)], [(0, 0, 0, 0, 4)], [(0, 5, 0, 0, 0), (0, 5, 0, 0, 0)], [(0, 0, 5, 0, 0), (0, 0, 5, 0, 0)], [(0, 0, 0, 5, 0), (0, 0, 0, 5, 0)], [(0, 0, 0, 0, 5), (0, 0, 0, 0, 5)], [(0, 6, 0, 0, 0)], [(0, 0, 6, 0, 0)], [(0, 0, 0, 6, 0)], [(0, 0, 0, 0, 6)], [(0, 8, 0, 0, 0)], [(0, 0, 8, 0, 0)], [(0, 0, 0, 8, 0)], [(0, 0, 0, 0, 8)], [(0, 10, 0, 0, 0)], [(0, 0, 10, 0, 0)], [(0, 0, 0, 10, 0)], [(0, 0, 0, 0, 10)], [(1, 1, 0, 0, 0)], [(2, 2, 0, 0, 0)], [(1, 3, 0, 0, 0)], [(1, 1, 2, 0, 0)], [(1, 0, 1, 0, 0)], [(1, 0, 0, 1, 0)], [(1, 0, 0, 0, 1)], [(0, 1, 1, 0, 0)], [(0, 1, 0, 1, 0)], [(0, 1, 0, 0, 1)], [(0, 0, 1, 1, 0)], [(0, 0, 1, 0, 1)], [(0, 0, 0, 1, 1)], [(2, 0, 2, 0, 0)], [(2, 0, 0, 2, 0)], [(2, 0, 0, 0, 2)], [(0, 2, 2, 0, 0)], [(0, 2, 0, 2, 0)], [(0, 2, 0, 0, 2)], [(0, 0, 2, 2, 0)], [(0, 0, 2, 0, 2)], [(0, 0, 0, 2, 2)], [(1, 0, 3, 0, 0)], [(1, 0, 0, 3, 0)], [(1, 0, 0, 0, 3)], [(3, 1, 0, 0, 0)], [(3, 0, 1, 0, 0)], [(3, 0, 0, 1, 0)], [(3, 0, 0, 0, 1)], [(0, 1, 3, 0, 0)], [(0, 1, 0, 3, 0)], [(0, 1, 0, 0, 3)], [(0, 3, 1, 0, 0)], [(0, 3, 0, 1, 0)], [(0, 3, 0, 0, 1)], [(0, 0, 1, 3, 0)], [(0, 0, 1, 0, 3)], [(0, 0, 3, 1, 0)], [(0, 0, 3, 0, 1)], [(0, 0, 0, 1, 3)], [(0, 0, 0, 3, 1)], [(1, 1, 0, 2, 0)], [(1, 1, 0, 0, 2)], [(1, 2, 1, 0, 0)], [(1, 2, 0, 1, 0)], [(1, 2, 0, 0, 1)], [(1, 0, 1, 2, 0)], [(1, 0, 1, 0, 2)], [(1, 0, 2, 1, 0)], [(1, 0, 2, 0, 1)], [(1, 0, 0, 1, 2)], [(1, 0, 0, 2, 1)], [(2, 1, 1, 0, 0)], [(2, 1, 0, 1, 0)], [(2, 1, 0, 0, 1)], [(2, 0, 1, 1, 0)], [(2, 0, 1, 0, 1)], [(2, 0, 0, 1, 1)], [(0, 1, 1, 2, 0)], [(0, 1, 1, 0, 2)], [(0, 1, 2, 1, 0)], [(0, 1, 2, 0, 1)], [(0, 1, 0, 1, 2)], [(0, 1, 0, 2, 1)], [(0, 2, 1, 1, 0)], [(0, 2, 1, 0, 1)], [(0, 2, 0, 1, 1)], [(0, 0, 1, 1, 2)], [(0, 0, 1, 2, 1)], [(0, 0, 2, 1, 1)]]
vs=[[(0, 0, 0, 0, 0)], [(1, 0, 0, 0, 0), (1, 0, 0, 0, 0)], [(2, 0, 0, 0, 0)], [(3, 0, 0, 0, 0), (3, 0, 0, 0, 0)], [(4, 0, 0, 0, 0)], [(5, 0, 0, 0, 0), (5, 0, 0, 0, 0)], [(6, 0, 0, 0, 0)], [(7, 0, 0, 0, 0), (7, 0, 0, 0, 0)], [(8, 0, 0, 0, 0)], [(10, 0, 0, 0, 0)], [(12, 0, 0, 0, 0)], [(14, 0, 0, 0, 0)], [(0, 1, 0, 0, 0), (0, 1, 0, 0, 0)], [(0, 0, 1, 0, 0), (0, 0, 1, 0, 0)], [(0, 0, 0, 1, 0), (0, 0, 0, 1, 0)], [(0, 0, 0, 0, 1), (0, 0, 0, 0, 1)], [(0, 2, 0, 0, 0)], [(0, 0, 2, 0, 0)], [(0, 0, 0, 2, 0)], [(0, 0, 0, 0, 2)], [(0, 3, 0, 0, 0), (0, 3, 0, 0, 0)], [(0, 0, 3, 0, 0), (0, 0, 3, 0, 0)], [(0, 0, 0, 3, 0), (0, 0, 0, 3, 0)], [(0, 0, 0, 0, 3), (0, 0, 0, 0, 3)], [(0, 4, 0, 0, 0)], [(0, 0, 4, 0, 0)], [(0, 0, 0, 4, 0)], [(0, 0, 0, 0, 4)], [(0, 5, 0, 0, 0), (0, 5, 0, 0, 0)], [(0, 0, 5, 0, 0), (0, 0, 5, 0, 0)], [(0, 0, 0, 5, 0), (0, 0, 0, 5, 0)], [(0, 0, 0, 0, 5), (0, 0, 0, 0, 5)], [(0, 6, 0, 0, 0)], [(0, 0, 6, 0, 0)], [(0, 0, 0, 6, 0)], [(0, 0, 0, 0, 6)], [(0, 7, 0, 0, 0), (0, 7, 0, 0, 0)], [(0, 0, 7, 0, 0), (0, 0, 7, 0, 0)], [(0, 0, 0, 7, 0), (0, 0, 0, 7, 0)], [(0, 0, 0, 0, 7), (0, 0, 0, 0, 7)], [(0, 8, 0, 0, 0)], [(0, 0, 8, 0, 0)], [(0, 0, 0, 8, 0)], [(0, 0, 0, 0, 8)], [(0, 10, 0, 0, 0)], [(0, 0, 10, 0, 0)], [(0, 0, 0, 10, 0)], [(0, 0, 0, 0, 10)], [(0, 12, 0, 0, 0)], [(0, 0, 12, 0, 0)], [(0, 0, 0, 12, 0)], [(0, 0, 0, 0, 12)], [(0, 14, 0, 0, 0)], [(0, 0, 14, 0, 0)], [(0, 0, 0, 14, 0)], [(0, 0, 0, 0, 14)], [(1, 1, 0, 0, 0)], [(2, 2, 0, 0, 0)], [(1, 3, 0, 0, 0)], [(1, 1, 2, 0, 0)], [(1, 0, 1, 0, 0)], [(1, 0, 0, 1, 0)], [(1, 0, 0, 0, 1)], [(0, 1, 1, 0, 0)], [(0, 1, 0, 1, 0)], [(0, 1, 0, 0, 1)], [(0, 0, 1, 1, 0)], [(0, 0, 1, 0, 1)], [(0, 0, 0, 1, 1)], [(2, 0, 2, 0, 0)], [(2, 0, 0, 2, 0)], [(2, 0, 0, 0, 2)], [(0, 2, 2, 0, 0)], [(0, 2, 0, 2, 0)], [(0, 2, 0, 0, 2)], [(0, 0, 2, 2, 0)], [(0, 0, 2, 0, 2)], [(0, 0, 0, 2, 2)], [(1, 0, 3, 0, 0)], [(1, 0, 0, 3, 0)], [(1, 0, 0, 0, 3)], [(3, 1, 0, 0, 0)], [(3, 0, 1, 0, 0)], [(3, 0, 0, 1, 0)], [(3, 0, 0, 0, 1)], [(0, 1, 3, 0, 0)], [(0, 1, 0, 3, 0)], [(0, 1, 0, 0, 3)], [(0, 3, 1, 0, 0)], [(0, 3, 0, 1, 0)], [(0, 3, 0, 0, 1)], [(0, 0, 1, 3, 0)], [(0, 0, 1, 0, 3)], [(0, 0, 3, 1, 0)], [(0, 0, 3, 0, 1)], [(0, 0, 0, 1, 3)], [(0, 0, 0, 3, 1)], [(1, 1, 0, 2, 0)], [(1, 1, 0, 0, 2)], [(1, 2, 1, 0, 0)], [(1, 2, 0, 1, 0)], [(1, 2, 0, 0, 1)], [(1, 0, 1, 2, 0)], [(1, 0, 1, 0, 2)], [(1, 0, 2, 1, 0)], [(1, 0, 2, 0, 1)], [(1, 0, 0, 1, 2)], [(1, 0, 0, 2, 1)], [(2, 1, 1, 0, 0)], [(2, 1, 0, 1, 0)], [(2, 1, 0, 0, 1)], [(2, 0, 1, 1, 0)], [(2, 0, 1, 0, 1)], [(2, 0, 0, 1, 1)], [(0, 1, 1, 2, 0)], [(0, 1, 1, 0, 2)], [(0, 1, 2, 1, 0)], [(0, 1, 2, 0, 1)], [(0, 1, 0, 1, 2)], [(0, 1, 0, 2, 1)], [(0, 2, 1, 1, 0)], [(0, 2, 1, 0, 1)], [(0, 2, 0, 1, 1)], [(0, 0, 1, 1, 2)], [(0, 0, 1, 2, 1)], [(0, 0, 2, 1, 1)]]
print(len(vs))
vs=L
print(len(vs))
solutions = set()
count = 0
def assemble(sol):
global count
count += 1
if count % 10000 == 0:
print(count, len(sol))
dim_sum = 0
for v in itertools.chain.from_iterable(sol):
dim_sum += dim(v)
if dim_sum == 12 and not isZero(list(itertools.chain.from_iterable(sol))):
return [sol]
if dim_sum >= 12:
return []
sols = []
for i in range(len(vs)):
sols += assemble(sol + [vs[i]])
return sols
solutions = assemble([])
print(list(solutions))
for sol in solutions:
print(sol)
print("Length of solution:")
print(len(solutions))
orderList3=order(3)
basis =[[0]]
for sol in solutions:
res=[]
for x in sol:
for y in x:
res.append(y)
if checker(res,basis,orderList3)==-1:
basis.append(res)
print("length of the basis:")
print(len(basis))
for b in basis:
print(b)