def diagA1(L,index1,index2): n=len(L[0]) sol=[] ind1=index1-1 ind2=index2-1 for i in range(0,len(L)): #print("i= "+str(i)) x=tensor([(L[i][ind1],)],[(L[i][ind2],)]) #print("this is x:") #print(x) temp2=[] for t in range(0,len(x)): temp2.append((x[t],)) #for t in range(0,len(x)): # x=[(x[0],)] #print("This is temp2:") #print(temp2) s=() for j in range(0,n): if j!=ind1 and j!=ind2: s=s+(L[i][j],) #print("this is s:") #print(s) temp=[] temp.append(s) a=concat(temp2,temp) #print(a) for t in range(0,len(a)): sol.append(a[t]) #sol.append(a) return(sol)
lambda1=[[(0)]] lambda2=[[(0)]] l=[[(0)]] lambda1.append([(2,1)]) lambda1.append([(1,1),(1,0)]) lambda1.append([(1,1),(0,0),(0,0)]) lambda1.append([(3,0),(0,1)]) lambda1.append([(2,0),(0,2)]) lambda1.append([(2,0),(0,1),(0,0)]) lambda1.append([(1,0),(1,0),(0,1)]) lambda1.append([(1,0),(0,1),(0,0),(0,0)]) for i in range(1,len(lambda1)): lambda2.append(wedge2(lambda1[i])) l.append(tensor(lambda1[i], lambda1[i])) l[i].remove((0,0)) for i in range(1,len(lambda1)): Lie[i]=l[i] for j in range(0,6): Lie[i]=Lie[i]+lambda2[i] for j in range(0,8): Lie[i]=Lie[i]+[(0,0)] for i in range(1,len(lambda1)): conjugate = checker(Lie[i], basisA12,orderList2) # p=str(i)+" & "+latex(lambda1[i])+" & "+latex(Lie[i])+" & "+str(conjugate )+"\\\\ \\hline" # print(p)