[sympy] Re: Matrix inversion bug

2008-11-06 Thread henrik johansson
Hi, Indeed. Could you please prepare a patch? We'll merge it. Right. That's why we should allow the user to pass his own simplifying functions. I think there will always be cases when simplify cannot do it, so that when the user can pass it's own zero finder, he can fix that easily. And

[sympy] Re: electric circuits solver

2008-10-25 Thread henrik johansson
Hi Ondrej and Travis, I have looked through Travis code and compared it to mine. And there are some differences in the way the matrices are constructed. And mine is more heavyweighted with more dependencies. Here are some features: * numeric and symbolic solutions * Pythonic interface for

[sympy] Re: Matrix inversion bug

2008-10-22 Thread henrik johansson
But there are probably cases where a simplify is not enough. I do not know if there is a more robust way of detecting if an expression is zero. Mathematically speaking, no. But in practise, I am sure it is. Try to research some literature, and also please create a new issue for each

[sympy] Re: Matrix inversion bug

2008-10-18 Thread henrik johansson
Hi Ondrej, Thanks for the help. I dug a little deeper and found that the problem was that both the LU and GaussElimination solver could not detect zeros in the pivot elements. The problems disappeared if the pivot elements were simplified first. For the GE solver I could just set the simplifed