def assemble(L,n):
print("newloop : " + str(n))
sol=[]
for i in range(0,len(L)):
print("i="+str(i))
a=dimChecker(L[i])
print(a)
if a>n:
return []
if a==n:
return [L[i]]
else:
print("here")
temp=assemble(L,n-a)
print("temp=")
print(temp)
if len(temp)>=1:
return L[i]+temp
else :
return []
#In G2C3
print("In G2")
a=expand("[(2,0)/(0,2)/(1,3)/(0,2)^9/(0,2)^5/(1,1)^9/(1,1)^5/(0,0)^9/(0,0)^9/(0,0)^3]")
print(checker(a,basisA12,orderList2))
print(dimChecker(a))
basisA12.append(a)
#InC3
print("in C3")
b=expand("[(6,0)/(2,0)/(3,1)/(4,0)^4/(0,0)^4/(3,1)^4/(4,0)^3/(3,1)^3/(0,2)/(0,0)^3/(0,0)^8/(0,0)^6]")
c=expand("[(2,0)/(4,2)/(0,2)/(4,0)^4/(2,2)^4/(4,0)^3/(2,2)^3/(0,0)^3/(0,0)^8/(0,0)^3]")
d=expand("[(2,0)^3/(0,2)/(1,1)^2/(2,0)^4/(1,1)^8/(2,0)^3/(1,1)^6/(0,0)^3/(0,0)^9/(0,0)^9/(0,0)^6/(0,0)^9]")
e=expand("[(2,0)/(0,2)/(1,1)^4/(1,0)^8/(0,1)^8/(1,1)^4/(1,0)^8/(0,0)^8/(0,0)^8/(0,1)^8/(0,0)^3/(0,0)^9/(0,0)^3]")
if test[i][0]==1:
x=(test[i][1],)
for j in range(2,len(test[i])):
x=x+(test[i][j],)
soltest.append(x)
return(soltest)
###############################################################
###############################################################
###########A18
###############################################################
###############################################################
a=expand("[(2,0,0,0,0,0,0,0)/(0,2,0,0,0,0,0,0)/(0,0,2,0,0,0,0,0)/(0,0,0,2,0,0,0,0)/(0,0,0,0,2,0,0,0)/(0,0,0,0,0,2,0,0)/(0,0,0,0,0,0,2,0)/(0,0,0,0,0,0,0,2)/(1,1,1,1,0,0,0,0)/(1,1,0,0,1,1,0,0)/(1,1,0,0,0,0,1,1)/(1,0,1,0,1,0,0,1)/(1,0,1,0,0,1,1,0)/(1,0,0,1,1,0,1,0)/(1,0,0,1,0,1,0,1)/(0,0,1,1,1,1,0,0)/(0,0,1,1,0,0,1,1)/(0,0,0,0,1,1,1,1)/(0,1,1,0,1,0,0,1)/(0,1,1,0,0,1,1,0)/(0,1,0,1,1,0,1,0)/(0,1,0,1,0,1,0,1)]")
print("A18")
print(dimChecker(a))
print(latex(a))
basisA18=[[0]]
basisA18.append(a)
###############################################################
###############################################################
###########A17
###############################################################
###############################################################
a=[0,0,0,0,0,00,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,00,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,00,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,00,0,0,0,0,0,0,0,0,0,0]
orderList = order(7)
basisA17=[[0]]
#We make the basis from basis of A18
a[0]=diagA1(basisA18[1],1,2)
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:
orderList3=order(3)
basisA13=[[0]]
basisA13Lambda1=[[0]]
for i in range(1,38):
print("case i")
LE8[i]=[]
LE8[i]=concat([(1,1)],lambda1[i])
LE8[i]=LE8[i]+concat([(1,0)],lambda4[i])
LE8[i]=LE8[i]+[(2,0,0),(0,2,0)]
LE8[i]=LE8[i]+concat([(0,0)],lambda2[i])
LE8[i]=LE8[i]+concat([(0,1)],lambda5[i])
if dimChecker(LE8[i])!=248:
print("LE8 does not have dim 248 at i=")
print(i)
temp=checker(LE8[i],basisA13,orderList3)
print(temp)
if temp==-1:
print("New element, added to basis")
basisA13.append(LE8[i])
basisA13Lambda1.append(lambda1[i])
else:
print("conjugate")
t=diagA1(basisA18[1],2,7)
t2=removeA1(t,7)
#print(t2[1])
#print(checker(t2,basis2,orderList6))
basis2.append(t2)
#In A1B4 in A1F4
temp=expand("[(1,1,1,0,1,0)/(1,1,1,0,0,1)/(1,1,0,1,1,0)/(1,1,0,1,0,1)/(1,1,1,1,0,0)/(1,1,0,0,1,1)/(1,1,0,0,0,0)^2/(1,3,0,0,0,0)/(0,2,1,1,0,0)/(0,2,0,0,1,1)/(0,2,0,0,0,0) /(0,2,1,0,1,0)/(0,2,1,0,0,1)/(0,2,0,1,1,0)/(0,2,0,1,0,1)/(0,2,0,0,0,0)^2/(0,0,2,0,0,0)/(0,0,0,2,0,0)/(0,0,0,0,2,0)/(0,0,0,0,0,2)/(0,0,1,1,1,1)/(0,0,1,1,0,0)/(0,0,0,0,1,1)/ (0,0,1,0,1,0)/(0,0,1,0,0,1)/(0,0,0,1,1,0)/(0,0,0,1,0,1)/(2,0,0,0,0,0)]")
print("In A1B4 in A1F4")
print(dimChecker(temp))
print(checker(temp,basis2,orderList6))
# In G2C3
print("in G2C3")
temp2=expand("[ (1,0,2,1,0,0)/(1,1,1,1,0,0)/(1,0,2,0,1,0)/(1,1,1,0,1,0)/(1,0,2,0,0,1)/(1,1,1,0,0,1)/(1,0,0,1,1,1)/(1,0,0,1,0,0)/(1,0,0,0,1,0)/(1,0,0,0,0,1)/ (0,0,2,0,0,0)^2/(0,1,1,0,0,0)^2/(0,0,2,1,1,0)/(0,1,1,1,1,0)/(0,0,2,1,0,1)/(0,1,1,1,0,1)/(0,0,2,0,1,1)/(0,1,1,0,1,1)/(0,2,0,0,0,0)/(0,0,2,0,0,0)/ (0,1,3,0,0,0)/(0,0,0,2,0,0)/(0,0,0,0,2,0)/(0,0,0,0,0,2)/(0,0,0,1,1,0)/(0,0,0,1,0,1)/ (0,0,0,0,1,1)/(2,0,0,0,0,0)]")
print(dimChecker(temp2))
print(checker(temp2,basis2,orderList6))
print("A_1^3 in D5")
lam=[0]
lam.append(expand("[(1,1,0)/(0,1,1)/(0,0,0)^2]"))
lam.append(expand("[(1,1,0)/(0,0,1)^2/(0,0,0)^2]"))
lam.append(expand("[(1,1,0)/(0,0,2)/(0,0,0)^3]"))
lam.append(expand("[(1,1,0)/(0,0,2)^2]"))
lam.append(expand("[(1,1,0)/(2,0,0)/(0,0,2)]"))
lam.append(expand("[(1,1,0)/(0,0,4)/(0,0,0)]"))
lam.append(expand("[(1,0,0)^2/(0,2,0)/(0,0,2)]"))
lam.append(expand("[(2,0,0)/(0,2,0)/(0,0,2)/(0,0,0)]"))
print(len(lam))
for i in range(1,len(lam)):
print(dimChecker(lam[i]))
basisA13=[0]
basisA13.append(expand("[(1,4,1)/(1,2,1)/(1,0,3)/(2,0,0)/(0,4,2)/(0,4,0)/ (0,2,2)/(0,2,0)/(0,0,2)]"))
basisA13.append(expand("[(1,2,2)^2/(1,0,0)^2/(2,0,0)/(0,4,0)/(0,0,4)/(0,2,0)/(0,0,2)/(0,2,2)^2/ (0,0,0)]"))
basisA13.append(expand("[(1,2,1)^2/(1,2,0)^2/(1,0,0)^2/(2,0,0)/(0,4,0)/(0,2,1)^2/(0,2,0)^3/(0,0,2)/(0,0,1)^2/(0,0,0)^2]"))
basisA13.append(expand("[(1,1,1)^2/(1,3,0)/(1,1,0)/(1,1,2)/(2,0,0)/(0,2,2)/ (0,2,0)^2/(0,0,2)/(0,2,1)^2/(0,0,1)^2/(0,0,0)]"))
