[SciPy-user] Bug in inverse matrix calculation for "scipy_core-0.8.4.win32-py2.3.exe"

Nils Wagner nwagner at mecha.uni-stuttgart.de
Wed Dec 14 09:01:13 CST 2005


Frank wrote:
>"On this note.  Is there a way of determining before hand in scipy if a
>matrix is singular, say by calculating its rank, e.g.:"
>
>if det(matrix) = 0, then the matrix is singular
>
>_______________________________________________
>SciPy-user mailing list
>SciPy-user at scipy.net
>http://www.scipy.net/mailman/listinfo/scipy-user
>  
This is not very reliable when using floating point operations.
See the book by Carl D. Meyer "Matrix analysis and applied linear
algebra" SIAM 2000
especially chapter 6 p. 466 Example 6.1.2

Nils



More information about the SciPy-user mailing list