[Scipy-tickets] [SciPy] #1458: lu_factor gives different results, depending on whether the input is C or Fortran contiguous (also affects lu_solve)

SciPy Trac scipy-tickets@scipy....
Thu Jun 9 11:06:34 CDT 2011


#1458: lu_factor gives different results, depending on whether the input is C or
Fortran contiguous (also affects lu_solve)
--------------------------+-------------------------------------------------
 Reporter:  mike.wimmer   |       Owner:  somebody
     Type:  defect        |      Status:  new     
 Priority:  normal        |   Milestone:          
Component:  scipy.linalg  |     Version:  0.9.0   
 Keywords:                |  
--------------------------+-------------------------------------------------

Comment(by warren.weckesser):

 This works for me, using numpy version 1.5.1, and for both scipy 0.9.0rc2
 (from EPD) and scipy trunk:
 {{{
 In [6]: A=np.array([[2,1],[0,1]])

 In [7]: la.lu_factor(A)
 Out[7]:
 (array([[ 2.,  1.],
        [ 0.,  1.]]), array([0, 1]))

 In [8]: la.lu_factor(np.asfortranarray(A))
 Out[8]:
 (array([[ 2.,  1.],
        [ 0.,  1.]]), array([0, 1]))
 }}}
 The ticket indicates version 0.9.0 of scipy--is that the version that you
 are actually using?

-- 
Ticket URL: <http://projects.scipy.org/scipy/ticket/1458#comment:2>
SciPy <http://www.scipy.org>
SciPy is open-source software for mathematics, science, and engineering.


More information about the Scipy-tickets mailing list