Status: New
Owner: ----
Labels: Type-Defect Priority-Medium
New issue 3216 by [email protected]: cse: substitution failure
http://code.google.com/p/sympy/issues/detail?id=3216
Dear All
Trying to cse a long expression (that I would like to input to fcode()
at the end), it fails at a substitution.
The script is attached and the output is:
eqs= [t11*(0.738558766382022*t15**(1/10)*(-0.104484691940934*exp(t33) +
0.104484691940934*exp(t38))/t2 + (0.738558766382022*t30 + t4)*Ei(t38) -
(0.738558766382022*t30*Ei(GRHO*t15**(3/10)*t2*(-0.0065302932463084*GRHO*a2*t14
- 2.75690433853138e-5*a5*t20 - 4.45572146470674e-6*a6*t17 -
0.00105542899627435*t19 - 0.000170578919531124*t21 -
7.20135751305192e-7*t22)/(0.0261211729852336*GRHO*b2*t14 +
0.161620459673995*b1*t15**(1/5)*t8**(1/3) + 0.00422171598509741*b3*t8*t9 +
0.000682315678124498*b4*t1*t12 + 0.000110276173541255*b5*t20 +
1.7822885858827e-5*b6*t17 + 2.88054300522077e-6*b7*t13*t6 +
4.65554684614493e-7*b8*t15**(8/5)*t8**(8/3) +
7.52431621307764e-8*b9*t15**(9/5)*t8**3 + 1.0)) + t4*Ei(t39) +
0.738558766382022*exp(t39)*((-40.7110354757367*t11*t2**2*t34 -
8.5073600005576*t18*t36 - 0.888888888888889*t36**(3/2))*(-1.2*A/t36**(3/2)
- 0.266666666666667*B/t36 - C*t32*(-0.000386980340521979*t16*t24 +
0.4*t16*t27*t3*(-0.000341157839062249*GRHO*a2*t14 -
1.44027150261038e-6*a5*t20 - 2.32777342307247e-7*a6*t17 -
5.51380867706276e-5*t19 - 8.91144292941349e-6*t21 -
3.76215810653882e-8*t22) + 0.4) +
t35*(-0.387889103217589*(t16*t27*(0.0261211729852336*a2*t16 +
0.000110276173541255*a5*t20 + 1.7822885858827e-5*a6*t17 +
0.00422171598509741*t19 + 0.000682315678124498*t21 +
2.88054300522077e-6*t22))**(1/2) + 2.4*(A + t30)**(1/2)) -
1.41796308072441*t35) - (0.444444444444444*B*t34 +
C*t36*(0.444444444444444*t18*t37 +
0.444444444444444)*(Sub(t16*t27*t3*(-0.000341157839062249*GRHO*a2*t14 -
1.44027150261038e-6*a5*t20 - 2.32777342307247e-7*a6*t17 -
5.51380867706276e-5*t19 - 8.91144292941349e-6*t21 -
3.76215810653882e-8*t22) + 1.0, 0.000967450851304948*t16*t24))))))]
<<<<<
so the "Sub(t16*t27*t3...." is my problem.
Any idea how to fix that?
Thanks
Valery
Attachments:
bug.py 27.0 KB
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" 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-issues?hl=en.
