Comment #1 on issue 1861 by asmeurer: Simplify uses all available memory in
some cases
http://code.google.com/p/sympy/issues/detail?id=1861
What modification makes it work fine? And what version of sympy are you
using? It works fine for me:
[0, 1, 1/a1, -1/a1, 0]
[0, 0, -1/a1, 1/a1 + 1/a2 + 1/a3, -1/a2 - 1/a3]
[1, 0, 0, -1/a2 - 1/a3, 1/a2 + 1/a3]
[0, 0, 1, 0, 0]
[0, 0, 0, -e, 1]
[ 0]
[ 0]
[ 0]
[-d1]
[ 0]
det(M): -1/a1 - 1/a2 - 1/a3 + e/a2 + e/a3
[ -d1*((-1/a2 - 1/a3)/(a1*(1/a1 + 1/a2 + 1/a3)) + e*(1/a2 + 1/a3
- (-1/a2 - 1/a3)**2/(1/a1 + 1/a2 +
1/a3))/(a1*(1 + e*(-1/a2 - 1/a3)/(1/a1 + 1/a2 + 1/a3))*(1/a1 + 1/a2 +
1/a3)))]
[-d1*((1/a2 + 1/a3)/(a1*(1/a1 + 1/a2 + 1/a3)) + e*(-(1/a2 + 1/a3)*(-1/a2 -
1/a3)/(1/a1 + 1/a2 + 1/a3) - 1/a2
- 1/a3)/(a1*(1 + e*(-1/a2 - 1/a3)/(1/a1 + 1/a2 + 1/a3))*(1/a1 + 1/a2 +
1/a3)))]
[
d1]
[ -d1*(e*(-1/a2
- 1/a3)/(a1*(1 + e*(-1/a2 - 1/a3)/(1/a1 + 1/a2 + 1/a3))*(1/a1
+ 1/a2 + 1/a3)**2) - 1/(a1*(1/a1 + 1/a2 + 1/a3)))]
[
d1*e/(a1*(1
+ e*(-1/a2 - 1/a3)/(1/a1 + 1/a2 +
1/a3))*(1/a1 + 1/a2 + 1/a3))]
[(a2*d1 + a3*d1 - a2*d1*e - a3*d1*e)/(a1*a2 + a1*a3 + a2*a3 - a1*a2*e -
a1*a3*e), (-a2*d1 - a3*d1 + a2*d1*e
+ a3*d1*e)/(a1*a2 + a1*a3 + a2*a3 - a1*a2*e - a1*a3*e), d1, a2*a3*d1/(a1*a2
+ a1*a3 + a2*a3 - a1*a2*e -
a1*a3*e), a2*a3*d1*e/(a1*a2 + a1*a3 + a2*a3 - a1*a2*e - a1*a3*e)]
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--
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.