LinearAlgebraError: SVD did not converge
Damian Menscher
menscher+python at uiuc.edu
Tue Aug 27 05:26:52 EDT 2002
I'm trying to get the SVD inverse of a matrix, so I'm calling
LinearAlgebra.generalized_inverse(M, rcond). But sometimes it
bombs out with the error "LinearAlgebraError: SVD did not
converge". Seems odd that it would die, considering Numerical
Recipes states "SVD... (theoretically) cannot fail...." (section
15.4). I'm guessing this is because lapack_lite uses an
iterative method rather than a deterministic one? (I'm using
version 2.1, in case it's relevant.)
Does anyone know of a package that has a *working* SVD inverse?
I'd try to debug the lapack_lite one, but all I've got is a .so
file....
I can provide the inverse I'm trying to invert if needed. (I
didn't include it here since it's kinda large -- 26x26.)
TIA for any help or suggestions,
Damian Menscher
--
-=#| Physics Grad Student & SysAdmin @ U Illinois Urbana-Champaign |#=-
-=#| 488 LLP, 1110 W. Green St, Urbana, IL 61801 Ofc:(217)333-0038 |#=-
-=#| 1412 DCL, Workstation Services Group, CITES Ofc:(217)244-3862 |#=-
-=#| <menscher at uiuc.edu> www.uiuc.edu/~menscher/ Fax:(217)333-9819 |#=-
More information about the Python-list
mailing list