Comment #19 on issue 2956 by pr...@goodok.ru: Simetimes the
test_hyperexpand.py (test_meijerg_shift_operators) fails in master.
http://code.google.com/p/sympy/issues/detail?id=2956
To be sure, we must create reproducible failure without usage hyper(), e.g.
replace it with complex value (float complex).
I have tried to do this, but failed in it:
>>> c1 = hyper((3, 4, 5/2), (5, 9/2), z).subs(z, z1).n()
>>> c2 = hyper((2, 3, 3/2), (4, 7/2), z).subs(z, z1).n()
>>> c3 = hyper((1, 2, 1/2), (3, 5/2), z).subs(z, z1).n()
>>> c1
-0.894828455904212 + 1.86631565359559*I
>>> c2
-1.00803197294486 - 1.63879132343505*I
>>> c3
1.20830731062755 - 0.538594692525162*I
>>> h5 = S('z*(3*z*(-0.894828455904212 + 1.86631565359559*I)/35 +
2*(-1.00803197294486 - 1.63879132343505*I)/15)*((-65*z + 65)/(3*z) +
2/z*(-25*z + 25) + 2/z*(5*z - 5) + (20*z - 10)/(2*z) + 4/(3*z)*(50*z - 50)
+ (-70*z + 100)*(15*z - 15)/(18*z**2) + (3*z - 3)*(130*z - 250)/(9*z**2)) +
2*z*((-205*z + 395)/(6*z) + 2/z*(-25*z + 50) + (3*z - 3)*(205*z**2 - 720*z
+ 500)/(z*(9*z**2 - 9*z)) + 2/z*(5*z - 10) + (30*z - 65)/(2*z) +
4/(3*z)*(65*z - 100) + (-115*z + 200)*(15*z -
15)/(18*z**2))*(-1.00803197294486 - 1.63879132343505*I)/15 + (10 + (3*z -
3)*(70*z - 100)/(z*(9*z - 9)) - 20/(9*z)*(15*z - 15))*(1.20830731062755 -
0.538594692525162*I)')
>>> h6 = S(repr(h5)
>>> [h.subs(z,z1).n().expand() for h in [h5, h6]]
Out[32]: [-1.66600996337348 + 0.46405746708775⋅ⅈ, -1.66600996337348 +
0.464057467087749⋅ⅈ]
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to sympy-issues@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.
