[SciPy-User] Eigenvectors of sparse symmetric matrix

Anne Archibald aarchiba@physics.mcgill...
Mon Oct 25 20:42:09 CDT 2010


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


More information about the SciPy-User mailing list