Issue 1158: regression: roots of polys
http://code.google.com/p/sympy/issues/detail?id=1158
Comment #2 by goriccardo:
Thanks for the reply Mateusz.
The problem is that Matrix.eigenvals uses roots method and at the time it
doesn't
works in many cases with 3x3 or 4x4 matrices. 3x3 matrices works with the
cubics=True
keyword but often gives error calling eigenvects. It used to work in
0.5.15. BTW I'll
try to fix Matrices.
Here is another kind of problem, it seems easy to fix:
In [1]: p = x**4 - 6*x**3 - 21*x**2 + 118*x + 216
In [2]: roots(p,x,quartics=True)
---------------------------------------------------------------------------
NameError Traceback (most recent call last)
/home/riccardo/Sources/sympy/sympy/<ipython console> in <module>()
/home/riccardo/Sources/sympy/sympy/sympy/polys/rootfinding.pyc in roots(f,
*symbols,
**flags)
256
257 for i, h in enumerate(poly_sqf(g)):
--> 258 for zero in roots_trivial(h):
259 if zeros.has_key(zero):
260 zeros[zero] += i+1
/home/riccardo/Sources/sympy/sympy/sympy/polys/rootfinding.pyc in
roots_trivial(g)
226 elif n == 4:
227 if flags.get('quartics', False):
--> 228 zeros += roots_quartic(g)
229
230 return zeros
/home/riccardo/Sources/sympy/sympy/sympy/polys/rootfinding.pyc in
roots_quartic(f)
71 r = c**2 + a**2*d - a*b*c
72
---> 73 r = cubic(1, p, q, r)
74
75 u = r[1] + r[2] - r[0]
NameError: global name 'cubic' is not defined
--
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 [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
-~----------~----~----~----~------~----~------~--~---
