Hi Matt, Well, I got the typos fixed and some spacing errors and was able to run the test_matrices.py. It crapped out on the matrix
A = Matrix(4,4,(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16)) This matrix has zero determinant by the other methods (I checked) but fails with my LU decomp method because LUdecomposition() itself fails. I gather that LU decomp has a problem for any square singular matrix. So, this method is not going to work for a general square matrix of numbers. However, I have not seen an example of a symbolic matrix for which it has failed. Hence, perhaps with the caveat that it is not ok for matrices which have det = 0 and probably works for most symbolic matrices, I think I could proceed with a pull request. What do you think? If the other developers are going to recommend killing the idea I will not bother to do a pull request. I do have the is_anti_symmetric method and the dual of a square matrix which can perhaps still go. What do you think? Comer -- You received this message because you are subscribed to the Google Groups "sympy" group. To post to this group, send email to [email protected]. To unsubscribe from this group, send email to [email protected]. For more options, visit this group at http://groups.google.com/group/sympy?hl=en.
