[NumPy-Tickets] [NumPy] #2192: Inverse of a singular matrix -Bug

[NumPy Trac]

#2192: Inverse of a singular matrix -Bug
--------------------------+-------------------------------------------------
 Reporter:  Amit001       |       Owner:  pv         
     Type:  defect        |      Status:  new        
 Priority:  normal        |   Milestone:  Unscheduled
Component:  numpy.linalg  |     Version:  1.6.1      
 Keywords:                |  
--------------------------+-------------------------------------------------
 Hi,

 Numpy is not handling singular matrices properly.

 >>> numpy.__version__
 '1.6.1'

 for eg. Numpy is giving the inverse of the below singular matrix.

 a = array([[ 1,  2,  3],
        [ 5,  7,  9],
        [ 6,  9, 12]])

 >>> det(a)
 -5.3290705182007443e-015
 >>> rank(a)
 2
 >>> inv(a)
 array([[ -5.62949953e+14,  -5.62949953e+14,   5.62949953e+14],
        [  1.12589991e+15,   1.12589991e+15,  -1.12589991e+15],
        [ -5.62949953e+14,  -5.62949953e+14,   5.62949953e+14]])

-- 
Ticket URL: <http://projects.scipy.org/numpy/ticket/2192>
NumPy <http://projects.scipy.org/numpy>
My example project