[SciPy-User] Eigen not working

[Vivek Lohani]

Hi,

I was trying to implement the code for diagonalizing sparse matrices through
scipy.sparse.linalg.eigen but it turns out that for N=9 (which leads to   
 2^N X 2^N matrix) and the subsequent input paramter= f, which sets J=-1,
I am not getting the correct output for -0.71<g<0.71. It also happens for
lower values of N(>=5) on my computer in a region -0.6<g<0.6. I am unable
to understand what is going wrong because i have a numpy routine too which
does the job correctly.


Vivek Lohani

-------------------------------------------------------------------------
M.Sc. Integrated Physics
Indian Institute of Technology
Kanpur 208016, Uttar Pradesh
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sparse_test.py
Type: text/x-python
Size: 3039 bytes
Desc: not available
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20110609/56c00e28/attachment.py>