Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- c= Sum(u(i)*<u(i),y>/v(i) )
- k=np.array([[1,2,3],[2,3,4],[2,7,8]])
- y=np.array([1,4,5])
- w,u=np.linalg.eigh(k)
- w=array([ -2.02599523, 0.47346124, 13.552534 ])
- u=array([[-0.18897996, 0.95770742, 0.21698634],
- [ 0.82245177, 0.03363605, 0.5678395 ],
- [-0.53652554, -0.28577109, 0.79402471]])
- uDoty=np.dot(u,y)
- div=np.divide(y,w)
- div=np.divide(uDoty,w)
- r=np.tile(div,(len(u),1))
- a=u*r.T
- c=sum(a)
- U = np.dot(u, y)
- d = U/w
- a = u*d[:,None]
- c = a.sum()
Add Comment
Please, Sign In to add comment