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.

Reply via email to