[Numpy-discussion] SVD error in Numpy. Bug?
Wed Mar 19 08:28:03 CDT 2008
I tried sending this message yesterday, but it is
being held up because the MatrixMarket attachment is
too large. The moderator my release it to the group,
but I don't know so I am sending the original message
minus the attachment. If anyone wants the
MatrixMarket version of my problem matrix, just let me
know and I will send it them directly on email.
---- The original message:
The determinant of my matrix is
I expect it to be small near a solution to my problem
whose solution is the vector closest to the null space
of the original matrix. That's the reason I am using
The MatrixMarket file of the complex 36 x 36 matrix is
attached as requested.
FYI: I found a curious workaround. If I catch the
linalg.linalg.LinAlgError exception that svd throws
and then "square" the original matrix:
the SVD on newmat works fine and the square root of
the minimum singular value (which is what I am looking
for) appears correct. If condition number were the
problem in some way, I would expect newmat to be
worse. Maybe the newmat symmetric form is better
behaved. Why? Beats me.
Thanks for your help.
-- Lou Pecora, my views are my own.
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile. Try it now. http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ
More information about the Numpy-discussion