Advertisement
Guest User

Untitled

a guest
Nov 23rd, 2014
121
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.60 KB | None | 0 0
  1. Dict={}
  2.     start = time.time()
  3.     d=x.shape[1]
  4.     N=x.shape[0]
  5.     C=1
  6. #G[:N,d:d+N]=np.diag(np.diag(-1*np.ones((N,N))))
  7.     q=C*np.ones(N+d)
  8.     q[0:d]=0
  9.     P=np.zeros((N+d,N+d))
  10.     P[:d,:d]=np.diag(np.diag(np.ones((d,d))))
  11.     h=np.zeros(2*N)
  12.     h[:N]-=np.ones(N)
  13.     G=np.zeros((N+N, d+N))
  14.     G[:N,d:d+N]=np.diag(np.diag(-1*np.ones((N,N))))
  15.     G[N:2*N,d:d+N]=np.diag(np.diag(-1*np.ones((N,N))))
  16.     e=x.transpose()*y
  17.     G[:N,:d]=e.transpose()
  18.  
  19.     S=qp(matrix(P), matrix(q), matrix(G), matrix(h))
  20.    # print(S['primal objective'])
  21.    # print(S['dual objective'])
  22.     finish=time.time()
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement