Status: Started
Owner: ronan.l...@gmail.com
Labels: Type-Defect Priority-Medium NeedsReview Milestone-Release0.7.0
WrongResult
New issue 1969 by ronan.l...@gmail.com: Matrix.LUdecompositionFF() fails
for some inputs
http://code.google.com/p/sympy/issues/detail?id=1969
Matrix.LUdecompositionFF() returns wrong results in some cases. For
instance:
In [1]: M = Matrix([[0,0,1], [2,3,0], [3,1,4]])
In [2]: P, L, d, U = M.LUdecompositionFF()
In [4]: P*M == L*d.inv()*U
Out[4]: False
I have a branch that fixes it:
[http://github.com/rlamy/sympy/commits/LUdec]. I have also used the
opportunity to cleanup the method by making the code more pythonic and
improving the docstring.
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to sympy-iss...@googlegroups.com.
To unsubscribe from this group, send email to
sympy-issues+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.