[sympy] Re: Question about LU decomposition for non-square and non-invertible matrices

2015-08-08 Thread AMiT Kumar
Hi Tom, Sorry for late reply, if you are interested in contributing code, then you can send us a Pull Request on Github: https://github.com/sympy/sympy/pulls/ BTW, Wolframalpha also doesn't calculates

[sympy] Re: Question about LU decomposition for non-square and non-invertible matrices

2015-08-08 Thread Tom H
Hi Amit, I submitted a pull request. Thanks for the guidance! Best, Tom On Saturday, August 8, 2015 at 9:07:57 AM UTC-7, AMiT Kumar wrote: Hi Tom, Sorry for late reply, if you are interested in contributing code, then you can send us a Pull Request on Github:

[sympy] Re: Question about LU decomposition for non-square and non-invertible matrices

2015-08-06 Thread Tom H
Hi Amit, Thanks for the speedy response! If I understand the code in gauss_jordan_solve correctly, solving A*x=b with multiple right hand sides requires a call to gauss_jordan_solve for each right hand side b, which means that the same Gaussian elimination steps are performed on all but the

[sympy] Re: Question about LU decomposition for non-square and non-invertible matrices

2015-08-04 Thread AMiT Kumar
Hi Tom, LUdecomposition() doesn't works for under determined systems, you can use , linsolve(), It supports all types of systems: In [1]: from sympy import * In [2]: from sympy.solvers.solveset import linsolve In [3]: A = Matrix([[1,2,3],[4,5,6],[7,8,9]]) In [4]: b = Matrix([4, 13, 22]) In