Python wedge2示例

示例#1

文件： A12inE7.py 项目： manonabiteboul/thesiscalculations

def checkerA12inE8(l,basis,orderList):
    print("Checking if conjugate to any of the complete ones:")
    print(checker(l,basis,orderList))

    #These are two two cases where we do not know lambda5: cases 28 and 29
    l28=[(4,), (4,), (0,), (0,)]
    l29=[(5,), (5,)]
    L28=[]
    L28=L28+concat([(0,)],wedge2(l28))
    L28=L28+[(2,0)]
    
    L29=[]
    L29=L29+concat([(0,)],wedge2(l29))
    L29=L29+[(2,0)]

    basisTemp=[[(0)]]
    basisTemp.append(l)
    print("Conjugate to case 28?")
    print(classifyIncomplete(L28,basisTemp,orderList2))
    if len(classifyIncomplete(L28,basisTemp,orderList2))>1:
        print("YES")
    print("Conjugate to case 29?")
    print(classifyIncomplete(L29,basisTemp,orderList2))
    if len(classifyIncomplete(L29,basisTemp,orderList2))>1:
        print("YES")

示例#2

文件： A12inE7.py 项目： manonabiteboul/thesiscalculations

#Check lambda5??
lambda1.append([(8,), (0,), (0,), (0,)])
lambda5.append(expand("[(10,)^2/(4,)^2]"))

#What is lambda5??
lambda1.append([(4,), (4,), (0,), (0,)])
lambda5.append(expand("[(6,)^2/(4,)^2/(2,)^2/(0,)^2]"))

#What is the spin here?
lambda1.append([(5,), (5,)])
lambda5.append(expand("[(8,)^2/(4,)^2/(0,)^4]"))


for i in range(1,len(lambda1)):
    lambda2.append(wedge2(lambda1[i]))
    #print(dimChecker(lambda2[i]))


orderList2=order(2)
L=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
for i in range(1,len(lambda1)):
    L[i]=[]
    L[i]=L[i]+concat([(0,)],lambda2[i])
    L[i]=L[i]+[(2,0)]
    L[i]=L[i]+concat([(1,)],lambda5[i])
    x=checker(L[i],basisA12,orderList2)
 #   print("///////////////////////////////")
 #   print(i)
 #   print(x)
    basisA12.append(L[i])

示例#3

文件： recreating_correctbasisA12inE7.py 项目： manonabiteboul/thesiscalculations

lambda1.append([(4,), (4,), (0,), (0,)])
lambda5.append(expand("[(6,)^2/(4,)^2/(2,)^2/(0,)^2]"))
lambda4.append(expand("[(6,)^2/(4,)^2/(2,)^2/(0,)^2]"))

#What is the spin here?
lambda1.append([(5,), (5,)])
lambda5.append(expand("[(8,)^2/(4,)^2/(0,)^4]"))
lambda4.append(expand("[(5,)^3/(9,)/(3,)]"))

#What is the spin here?
lambda1.append([(5,), (5,)])
lambda5.append(expand("[(5,)^3/(9,)/(3,)]"))
lambda4.append(expand("[(8,)^2/(4,)^2/(0,)^4]"))

for i in range(1,len(lambda1)):
    lambda2.append(wedge2(lambda1[i]))
    #print(dimChecker(lambda2[i]))



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])

示例#4

文件： E8simplified.py 项目： manonabiteboul/thesiscalculations

e=[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]


d[1]=expand("[(1,1,0,0,0,0)/(0,0,1,1,0,0)/(0,0,0,0,1,1)/(0,0,0,0,0,0)^4]")
d[2]=expand("[(1,1,0,0,0,0)/(0,0,1,1,0,0)/(0,0,0,0,1,1)^2]")
d[3]=expand("[(1,1,0,0,0,0)/(0,0,1,1,0,0)/(0,0,0,0,1,1)/(1,0,1,0,0,0)]")
d[4]=expand("[(1,1,0,0,0,0)/(0,0,1,1,0,0)/(0,0,0,0,1,1)/(0,0,0,0,0,1)^2]")
d[5]=expand("[(1,1,0,0,0,0)/(0,0,1,1,0,0)/(0,0,0,1,1,0)/(0,0,0,0,0,1)^2]")
d[6]=expand("[(1,1,0,0,0,0)/(0,0,1,1,0,0)/(0,0,0,0,1,0)^2/(0,0,0,0,0,1)^2]")
d[7]=expand("[(1,0,0,0,0,0)^2/(0,1,1,0,0,0)/(0,0,0,1,1,0)/(0,0,0,0,0,2)/(0,0,0,0,0,0)]")
d[8]=expand("[(1,1,0,0,0,0)/(0,1,1,0,0,0)/(0,0,0,1,1,0)/(0,0,0,0,0,2)/(0,0,0,0,0,0)]")
d[9]=expand("[(1,0,0,0,0,1)/(0,1,1,0,0,0)/(0,0,0,1,1,0)/(0,0,0,0,0,2)/(0,0,0,0,0,0)]")


for i in range(1,10):
        sol[i]= wedge2(d[i])
        
#Original : e[2]=expand("[(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,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)/()]")

e[1]=expand("[(1,0,1,0,1,0)^2/(1,0,1,0,0,1)^2/(1,0,0,1,1,0)^2/(1,0,0,1,0,1)^2/(0,1,1,0,1,0)^2/(0,1,1,0,0,1)^2/(0,1,0,1,1,0)^2/(0,1,0,1,0,1)^2]")
e[2]=expand("[(1,0,1,0,1,1)/(1,0,1,0,1,1)/(1,0,0,1,2,0)/(1,0,0,1,0,0)/(1,0,0,1,0,2)/(1,0,0,1,0,0)/(0,1,1,0,1,1)/(0,1,1,0,1,1)/(0,1,0,1,2,0)/(0,1,0,1,0,0)/(0,1,0,1,0,2)/(0,1,0,1,0,0)]")


e[3]=expand("[(1,0,2,0,1,0)/(1,0,0,0,1,0)/(2,0,1,0,0,1)/(0,0,1,0,0,1)/(2,0,0,1,1,0)/(0,0,0,1,1,0)/(1,0,1,1,0,1)/(0,1,2,0,1,0)/(0,1,0,0,1,0)/(1,1,1,0,0,1)/(1,1,0,1,1,0)/(0,1,1,1,0,1)]")



#e[4]=expand("[()/()/()/()/()/()/()/()]")

示例#5

文件： D5inE6.py 项目： manonabiteboul/thesiscalculations

from compare import classifyIncomplete
#from compare import classifyIncomplete2
from compare import classify
from compare import perm
from tensor_wedge import concat
from remove_tensor_A1 import removeFirstA1
from remove_tensor_A1 import removeA1
from remove_tensor_A1 import diagA1
from tensor_wedge import printLE8
from tensor_wedge import tensor


print("A_1^4 in D5")
a=[(1,1,0,0),(0,0,1,1),(0,0,0,0),(0,0,0,0)]
b=[(1,1,0,0),(0,0,2,0),(0,0,0,2),(0,0,0,0)]
print(latex(wedge2(a)))
print(latex(wedge2(b)))



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))