Oct. 26, 2010
2:42 a.m.
On 25 October 2010 06:03, Pauli Virtanen <pav@iki.fi> wrote:
Sun, 24 Oct 2010 18:42:14 -0700, Hao Xiong wrote:
I am trying to compute the eigenvectors corresponding to the d+1 smallest eigenvalues of A=W.T*W. I started with W as a dense matrix and then W = sparse.csr_matrix(W) A = W.dot(W) # W.T * W
That is W*W and not (W.T)*W
Just quickly: is your matrix really symmetric? This could lead to all sorts of bizarre failures to converge. Anne