[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