Revision history [back]

Change function p to:

def p(*t):
   t = min( (t[1],t[0],t[3],t[2]), t )
   return SR.var(('p' + '_{}'*len(t)).format(*[str(i) for i in t]))

which will use the symmetry and use lexicographically smallest indices.