Talk:LU decomposition

From Rosetta Code
Revision as of 16:18, 18 November 2014 by Waynes (talk | contribs)

The Python example has a divide by zero error for the matrix

b = [[1, 1, 1, 1], [1, 1, -1, -1], [1, -1, 0, 0], [0, 0, 1, -1]]

although a LUP decomposition exists: [1]

I guess the pivotization is at fault.