montimaj

Numpy_test

Jan 24th, 2018
126
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.83 KB | None | 0 0
  1. def get_terms(values, max_deg):
  2.     x=values[0]
  3.     y=values[1]
  4.     calc=[]
  5.     flag={}
  6.     d1=1
  7.     while d1<=max_deg:
  8.         d2=0
  9.         while d2<=max_deg:
  10.             if d1+d2<=max_deg:
  11.                 if (d2,d1) not in flag.keys():
  12.                     flag[(d1,d2)]=False
  13.                     print(d1,d2)
  14.                     calc.append((x**d1)*(y**d2))
  15.                     if d1!=d2:
  16.                         calc.append((y**d1)*(x**d2))
  17.             d2+=1
  18.         d1+=1
  19.     calc.append(1)
  20.     zero_val=[0]*len(calc)
  21.     calcx=calc+zero_val
  22.     calcy=zero_val+calc
  23.     return calcx, calcy
  24.  
  25. deg=5
  26. comp=[(1,2),(3,4)]
  27. users=[(2,3),(4,4)]
  28.  
  29. for point in comp:
  30.     terms=get_terms(point, deg)
  31.     print('Terms[Len: '+str(len(terms[0]))+']=',terms[0])
  32.     print('Terms[Len: '+str(len(terms[1]))+']=',terms[1],'\n')
Add Comment
Please, Sign In to add comment